Analysis and projection of global iron ore trade: a panel data gravity model approach

In this study, an empirical model of global trade in iron ore is developed and applied. The empirical specification is based on the trade gravity theory in which the trade is determined by the income of the trading countries, the distance between the countries, and other characteris-tics of the countries. The model is specified allowing for country-specific effects. The estimation is performed with panel data for global bilateral iron ore trade flows from 1980 to 2016 including 121 countries and almost 14,000 observations. The results indicate a strong support of the gravity model hypotheses. On average, the trade value is projected to increase by approximately 5% per year up until 2035. The trade potential of iron ore is estimated to 410 million USD per year. Applied to forecasting and policy analysis, the results represent another worthwhile source of information providing an alternative view of the global trade in iron ore that can be helpful for decision-makers.


Introduction
In 2015, almost 63% of the global iron production was exported to a value of more than 65.4 billion USD (UN 2017). This is an increase in real value by more than 500% compared to 1980. The production, consumption, and trade patterns of iron ore have also changed over the last decades. Since 1980, the focus of trade has shifted from Europe and North America markets towards Asian and South American markets. But even within these markets, there have been major changes. For instance, China has replaced Japan as the major importer of iron ore in Asia while Brazil has replaced Canada as a major exporter. In addition, most European countries have experienced reduced shares of the trade value. These insights give raise to several interesting and important research questions. What have been the main drivers behind the changes in trade patterns? How will the future trade pattern change? Which countries will be the relative winners of these changes? What is the trade potential in iron ore? The objective of this paper is to apply a trade gravity model on the global iron ore trade in order to project and assess future trade patterns and trade potentials for iron ore.
The trade gravity model is commonly used to evaluate and predict trade patterns. The model has been applied in numerous specifications and contexts and it has been argued that the gravity model is the leading empirical model in economics regarding international trade flow (Anderson 2011). The basic trade gravity equation states that the trade flow between countries depends primarily on the size of the economy of the two countries, the dist a n c e b e t w e e n t h e m . I t h a s b e e n a p p l i e d t o Electronic supplementary material The online version of this article (https://doi.org/10.1007/s13563-017-0125-8) contains supplementary material, which is available to authorized users. macroeconomic issues dealing with aggregate trade flows and, e.g., the effects of currency unions (e.g., Frankel and Rose 2002;Glick and Rose 2002;Berger and Nitsch 2008). It has also been applied to disaggregated trade flows, such as agricultural commodities trade (e.g., Zahaniser et al. 2002;Peterson et al. 2013) and forest products trade (e.g., Buongiorno 2015Buongiorno , 2016Olofsson et al. 2017). In addition, the gravity model has been discussed in the context of cultural factors affecting trade (e.g., Hofstede 1994;Kogut and Singh 1998) and political factors (e.g., Srivastava and Green 1986). However, despite the importance of trade for the iron ore industry and the explanatory power of the gravity model, only a few studies have been found that quantitatively analyze the effects of globally changing economic conditions using gravity equations. Robertson and Robitaille (2014) investigate the importance of geography for trade in iron ore by estimating a gravity model. They show that iron ore has high trade elasticity with respect to distance. Moreover, Babri et al. (2017) separate the trade flow of iron ore into a long-term contract bound part and into a spot market (choice theory) part, using a gravity model. Figure 1 depicts the development of trade value and quantity for iron ore between 1980 and 2016. In 2016, the global export trade value of iron ore was 54,319 million USD (constant 2005 value). This is an increase by almost 1000 million USD compared to the year before and more than a fivefold increase compared to 1980 (UN 2017). As seen in the figure, the trade value peaked 2011 after a sharp increase starting around 2004. Between 1980 and 2004 the trade value was fairly constant. The sharp increase and decline in trade value is contrasting with the steady increase in traded weight. During the whole period, the trade in iron ore increased from 215 to 1496 million ton per year.

Global iron ore trade
The decline in trade value after 2011 can be explained by a number of interacting factors. Firstly, until 2005 the price of iron ore was fairly constant with a system of annual contracts keeping the price volatility down. Now, the pricing system is based on more flexible monthly contracts and spot markets. As a consequence, price fluctuations that used to occur over a matter of months now occur over a matter of days. Secondly, the decline in trade value coincides with an equally large decline in iron ore price. In turn, the price decline can partly be explained by the sagging demand growth from the Chinese steel industry. Since its peak in 2011, the iron ore price has declined by more than two thirds.
In terms of export value, Australia is the by far largest exporter and accounted for almost 59% of the global export value. Brazil is on a remote second place and accounted for approximately 20%. South Africa, Canada, and Sweden follow with five, four, and 2%, respectively. In the import side, Chine accounts for almost 71% of the trade value followed by Japan, Korea, and Germany with nine, five, and 3%, respectively. Thus, both in terms of export and import there are only few countries that account for the majority of the trade.
Iron ore is usually considered a heterogeneous commodity. It is predominantly mined as magnetite (72.4% Fe) and hematite (69.9% Fe) and traded as sinter or pellets. Ores containing very high quantities of hematite or magnetite are known as Bnatural ore^or DSO Bdirect shipping ore^and can be fed directly into iron-making blast furnaces. Iron ore can also be pre-treated before traded. For instance, sinter feed is made from DSO and pellets are made by taking very fine iron ore and agglomerating them with moisture and a binder to form balls. Pellets are essentially artificial lump. The technical specification of the individual steel mills sets the restriction of the type of iron ore that can be used. But on an aggregate level, the demand of iron ore is still a derived demand suggesting that it is the demand for steel products that are the main determinants of trade.

Trade gravity model
The first application of the gravity model on trade flows was by Tinbergen (1962) but it was not until 1995 it entered mainstream of trade research. Head and Mayer (2013) describe the latter development as three distinct steps: First, the realization that previous models could not explain Bmissing^trade (Trefler 1995). Second, the research on multilateral resistance/fixed effects established theoretical relationships between gravity models and, e.g., monopolistic competition and Armington and thus connecting the gravity model to solid micro-foundations (Anderson and van Wincoop 2003). It became clear that fixed effects could be used to capture multilateral resistance terms that emerged from different theoretical models (Redding and Venables 2004). The final step was the convergence of the literature on gravity models and heterogeneous firms (Chaney 2008). This work suggested changes in how the gravity equations should be estimated and how the estimated coefficients should be interpreted. Gravity model, used in an economic context, were initially criticized for lacking theoretical foundations (Deardorff 1984). However, the theoretical foundation of the gravity model has been developed and can now be considered consistent with microeconomic foundations. This work started with Anderson (1979) and has continued since. The structure of the gravity model can be derived under different theoretical frameworks, e.g., Ricardian (Eaton and Kortum 2002); Heckscher-Ohlin (Evenett and Kneller 2002) and models based on increasing returns and product differentiation (Bergstrand 1989). Because of the flexibility of the gravity model, its empirical success does not per se provide support for any particular model of trade. As with most empirical models, the trade gravity model also has limitations. Most prominent among the limitations are that trade costs of third parties are not considered important and do not affect the trade between two trading partners. Another limitation is that trade creation and diversion cannot be captured by the gravity model alone.
In its most basic form the trade gravity model stipulates that trade flows between trading partners vary positively with the magnitude of the partners' economic activities and negatively with the distance between them. Thus, the basic relationship can be expressed as: where T is the trade flow in monetary terms for trading partner i and j, A is a constant (across cross sections), Y is the magnitude of the economic activity, D is the distance, and the β's are unknown parameters. In this specification, possible third-country effects, if any, must be captured indirectly via either the exporting or the importing country. Eq. 1 can be linearly estimated by transforming the equation into logarithmic form. 1 In general, and for previously evaluated traded commodities, the trade gravity model has made good predictions. Feenstra et al. (2001) conclude that with no entry barriers, export is more sensitive to changes in domestic income than on changes in trading partners' income. In a later study, Stay and Kulkarni (2016) test the validity of the gravity model by examine United Kingdom (UK) and its trading partners. More specifically, they investigate how the trade flow between UK and its trading partners is influenced by colonial history and membership in the European Union. They conclude that the gravity model is good at predicting the trade flows. McCallum (1995) used the gravity equation to refute the notion that national borders had lost their economic relevance and estimated the effects of trade integration policies and the understanding of Bborder effects^. There have also been a few attempts to apply the gravity model on forestry products. Among the first empirical studies are Kangas and Niskanen (2003) who analyzed trade in forest products between EU and central and eastern European countries (CEE) by using a gravity model. The data include trade between 15 EU and ten CEE countries and has 498 observations. The result showed that the trade and production of forestry products under the 1990s increased in CEE countries. According to estimates with the gravity model, the trade between EU countries and CEE countries did not reached its expected value. Later, Akyüz et al. (2010) examines the trade with forestry products between Turkey and the EU countries using a gravity model. They used panel data for the period 2000-2006.
The results indicate that Turkey is below its potential export with EU countries and they conclude that it would be beneficial for both Turkey and EU if Turkey joins EU. Buongiorno (2015) analyzed if the introduction of the European monetary union had any effect on the trade flow of forestry products. He applies a gravity model on three forest products using a panel data on bilateral trade between 12 Euro-countries for the time period 1988-2013. By applying the gravity model in differential form, the time-invariant effects, e.g., distance between countries and common borders, was eliminated. The model was estimated using fixed effects. He found that trade increased for all products and countries after the introduction of the monetary union. In a later study, Buongiorno (2016) studied forestry product trade flow and made forecasts of the value of trade between member countries in Trans-Pacific Partnership (TPP) for the same forest products as in Buongiorno (2015). The data used is a panel of TPP participant for the period 2005-2014. The model was estimated using both fixed and random effects, resulting in similar results. The trade in all three forest products were positively affect by the TPP.

Data
The Observations with a trade value less than 1000 USD has been removed from the sample. The GDP and distance variables are obtained from CEPII (2017). Observations and projections for the GDP variables are summarized from the EconMap 2.4 Reference (Fouré et al. 2012(Fouré et al. , 2013. Three GDP variables are available: GPD in million constant 2005 USD; GDP in 2005 PPP per capita and GDP in million current USD. The distances between any two countries are summarized in the bilateral GeoDist database (Mayer and Zignago 2011). The database contains two simple and two weighted distance variables. The two simple measures are based on (1) latitudes and longitudes of the most important cities (in terms of population) and (2) the distance between the geographical coordinates of the capital cities. From the same database information on common borders, language and colonial history is obtained, as dummy variables. The common border variable is set to unity if two countries are contiguous. Two common language measures are defined: (1) if the two countries have the same official language and (2) if at least 9% of the populations speak the same language. The final dummy variable is if two countries ever had a colonial link. The descriptive statistics and the correlation matrix for the continuous variables are presented in Table 1.

Econometric issues
In order to comply with the theoretical conditions of the gravity model, the multilateral resistance term need to be accounted for empirically. Initially, a remoteness variable was used as a proxy, but as the theory developed it became clear that it was a too weak proxy for the resistance term (e.g., Wei 1996;Helliwell 1998). Instead, it is now common practice to used fixed effects estimation. This does not impose strong structural assumptions on the model. Another, more practical advantage of using fixed effects is that intermediate trade destinations can be accounted for. That is, countries such as the Netherlands and Singapore have large imports but they re-export. Ideally, the exporter should also be the producer and the importer the consumer. But the world is more complicated than that. However, fixed effects can control for this since they account for any unobservable that contributes to shift of the overall level of exports or imports of a country. Equation 1 is linearized by taking the logarithm and estimated with panel data from multiple bilateral trade flows between countries and years using fixed effects. Each observation is represented by the export value of iron ore between a county pair in a given year. By using the fixed effects approach it is assumed that the disturbance term can be expressed as u i, j, t = a i, j + ε i, j, t where a i, j is a constant for each bilateral trade flow. That is, a i, j measures the difference in exports for a particular country pair due to other variables than those explicitly in the model. The residual ε i, j, t is assumed i.d.d. with zero mean and constant variance across trade flows and years.
To test the validity of the assumptions, the fixed effects estimation is compared to random effects and OLS. The OLS model, under the assumption that the fixed effects are nil, can be rejected using an F test. A Hausman test is conducted to evaluate random vs. fixed effects. The Hausman test rejects the null hypothesis that random effects are appropriate. Thus, both test support the use of the fixed effect approach. In addition, the countries being observed is not a random sample from a larger population,  (Baier and Bergstrand 2007). The estimation is applied with robust estimation of the standard errors to account for possible heteroscedasticity or within-group correlation. Moreover, it is carried out with the Stata software using the xtreg, fixed effects option.

Results
The results are presented in three parts. First, the basic econometric results based on Eq. 1 using fixed effects are presented and discussed. Second, the estimated coefficients are used to simulate trade flow projections of iron ore up until the year 2035. Third, the trade potential is calculated. Table 2 presents the regression results based on five model specifications (I-V). Models I-III is based on the basic specification, which is the logarithmic linearization of Eq. 1. Model II uses a squared distance variable compared to Model I that is using a normal distance variable. Model III is using GDP per capita instead of GDP as in Model I. The specification is expanded in Model IV and V by introducing multilateral resistance terms. In Model IV, a dummy variable for common language, common border and colonial links are added. Two definitions of common language are tested (individually) in the specification of Model IV: (1) a common official language and (2) if a common language is spoken by at least 9% of the population in both countries. The results of Model IV are based on the official common language definition (the other definition produced statistically insignificant results as well and is not presented). Finally, in Model V, common language and colonial links are dropped. Of the five specifications, we chose to continue the analysis using Model V. It might be pointed out that although most of the iron ore trade is seaborne, suggesting that a common land border might not be an important determinant for trade, the data consist of all bilateral trade flows of iron ore. In fact, the descriptive statistics of the common border dummy variable suggest that almost 17% of the observations (bilateral trade flows) occur between countries that share a common border. The inclusion of the common border dummy variable is also supported by findings in the previous literature.
The results based on Model V support the basic structure of the trade gravity model. The trade increases with the size of the trading partners' economies and decreases with the distance between them. More specifically, if the GDP of the exporting country of iron ore increases by 1%, the trade flow value will increase by 0.64%. Similarly, of the GDP of the importing country increases by 1%, the trade flow value will increase by 0.48%. The further apart the countries are, the less trade will they engage in. For each percentage increase in distance, the trade flow will be reduced by 0.32%. Thus, in relative terms, the results suggest that the size of the economy has a larger effect on trade flow value than distance. Finally, the coefficient for the common border dummy suggests that if the trading countries share a common border (or more specifically, they constitute a continuous land mass) the trade flow value of iron ore will be 1.95 million USD (e 0.67 ) higher per year than otherwise. All continuous variables in logarithmic form * , ** , *** represent statistical significance at 10, 5, and 1%, respectively a Analyzed model In general, these results can be compared with previous studies using the trade gravity model, albeit not for iron ore specifically since previous studies are lacking for iron ore. Head and Mayer (2013) make a metaanalysis of gravity models estimates based on 159 studies and more than 2500 estimates. Their results can be contrasted with the results for iron ore trade presented herein. For the structural gravity models they conclude that the average effect of origin GDP (exporter GDP) is 0.74 and of destination GDP (importer GDP) is 0.58. These are similar results as for iron ore, both in magnitude and in relative ranking (i.e., exporter GDP has a larger effect than importer GDP). These results are in line with Feenstra et al. (2001) argument of home market effects. The average distance elasticity reported by Head and Mayer (2013) is − 1.1, which is lower than the results for iron ore trade. A plausible explanation could be that iron ore is mostly traded by sea transport making it less sensitive to distance compared to trade aggregates (or other specific goods and commodities) using land transportation modes. The contiguity effect (common border) is similar to that reported by Head and Mayer (2013). They estimate an average metaresult on 0.66, while the result for iron ore is 0.67. Overall, the results for iron ore seem to be in line with previous literature.

Iron ore trade projections
The iron ore trade projections are based on the estimated coefficients from Model V. Table 3 shows the projected 5-year incremental growth rates from 2015 to 2035 for the global iron ore trade. The table lists the total export from specific countries. The GDP projections are obtained from CEPII (2017). In addition, the observed export value of iron ore (million USD) for 2015 is presented in the table. The projections are based on bilateral trade observations higher than 30 million USD and are under the assumption that countries keep the same trading partners as it had during 2015. This delimitation is made in order to simply the calculations. The detailed trade projection calculations are presented in Table 5 in the Appendix.
As suggested by Table 3, we can expect the trade value of iron ore to increase in the future. In total, the annual average growth rate is expected to 5% initially, only to slowly decline over time and end up around 4.4% in 2035. However, even with a declining growth rate, the trade value is expected to increase over time. From 2020, the annual trade value change is 3.19 billion USD, which increase to 3.52 billion in 2025, 3.89 billion in 2030 and 4.25 billion in 2035. This represents an accumulated increase in trade value on 85 billion USD for the time period. Labson (1997) made the projection that the iron ore demand would increase by approximately 122 million metric tons per year until 2000 of which China alone was predicted to account for 46% of the increase. He also concluded that Australia, Brazil and India were expected to further expand their production and export of iron ore. Australian production was projected to increase by 27.7 million tons, Brazil by 22.3 million tons and India by 15.1 million tons.
Mongolia is projected to have the highest rate of increase in trade value (albeit from low levels). Of the six countries that had an iron ore trade value higher than 1000 million USD in 2015 (Australia, Brazil, Canada, South Africa, Ukraine, and Sweden), Table 3 indicate that all but Ukraine will experience an increasing but diminishing growth rate. However, if their growth rate is expressed in relative terms to Mongolia's growth rate another picture emerges (see Table 6). Australia, Canada, South Africa, and Ukraine will have an increase growth rate in iron ore trade value relative Mongolia, while Brazil and Sweden will have a steady growth rate. This disparity can partly be explained by the GDP projections. Firstly, the economies are expected to increase suggesting that more iron ore will be trade to a higher total value. Second, the increase in GDP is not uniform; some countries are expected to have a fast GDP growth compared to other countries. This is reflected in the relative growth rate in iron ore trade value.

Iron ore trade potentials
Having estimated the gravity model for global bilateral trade flows of iron ore, we proceed to estimate the global trade potential for iron ore. The regression estimates give the prediction of the average trade level. This suggests that some trade flows trade more than average while other flows trade less. By applying the estimated coefficients, it is possible to calculate the potential trade by analyzing the difference between the observed trade flows and the estimated (i.e., the size of the error term), with the caveat that no new trade routes are established (c.f. Batra 2004). The interpretation of a large negative trade potential suggests that a country could be trading more based on their economic and geographical fundamentals. Something the model cannot explain is holding trade back. However, the trade potential needs to be interpreted carefully since the error term also include statistical noise and possible measurement errors. That is, it is important not to overemphasize the trade potential but rather use it as a first indication on particular trade relationships.
Tables 4, 5, 6 present the calculated trade potential based on 2015 statistics. Countries with a trade potential less than 1 million USD has been excluded from the table due to space limitation. Overall, the trade in iron ore has the potential to increase by 410 million USD per year. For individual countries, Germany, France, the USA, and the UK have the largest potential to increase their trade in iron ore.

Conclusions
The objective of this paper is to apply a trade gravity model on the global iron ore trade in order to project and assess future trade patterns and trade potentials for iron ore. In general, the gravity model for iron ore performed well and is in line with the theoretical prediction of the model. That is, countries in proximity to each other will trade more with each other than with counties further away and large economies tend to have a larger trade than small ones. The empirical results conform to the metaresults for gravity models, suggesting that the approach can be used to analyze global iron ore trade.
The estimated coefficients are expected values over a very large number of trade flows and years of observation. There was strong statistical evidence supporting the gravity hypothesis, but as indicated by the standard errors, there were substantial variations for individual trade flows and years. The development of the home market has a larger effect on the trade flow value compared to the development of foreign markets. This suggests that policies towards developing home markets will also have a positive effect on the trade value of iron ore.
The trade projections for iron ore suggest a continuous increase in the trade value. This is explained by the expected economic development in both the exporting and importing countries. On average, the trade value is projected to increase by approximately 5% per year up until 2035. However, like most econometric projections, such forecasts must be viewed cautiously since they are subject to large potential errors. The trade potential of iron ore is estimated to 410 million USD per year. Conclusions made from this result must also be made carefully. Rather, it might give an indication of the possible trade potential.