Assessment of freight traffic flows and harmful emissions in euro-mediterranean context: scenario analyses based on a gravity model

Maritime transport contributes significantly to environmental pollution. The Mediterranean Sea Area is particularly affected by marine emissions from particulate matter, black carbon, nitrogen oxides and sulphur oxides. In addition to international and crossing traffic, the Mediterranean Sea is affected by the freight flows moved between its shores (in the North–South direction and vice versa), but also by transverse freight flows (East–West); these are freights transported through Short Sea Shipping. The Mediterranean Sea is the area where there is the highest concentration of short sea shipping in the EU-27. There are different types of Short Sea transport: container and bulk handling, general cargo ships and Ro-Ro transport. In such scenario, the analysis of trade flows as well as their spatial and geographical distribution, becomes fundamental. The paper proposes a gravity model for estimating trade flows, considering 18 countries boarding the Euro-Mediterranean Sea, in 2019. The proposed gravity model assumes that one of the main factors affecting trade is the economic dimension of a country which is directly related to the volume of imports and exports. In the paper, after a literature review on the gravity models, the illustration of the different phases investigation and construction of the database (specification, calibration and validation) to the definition and implementation of the proposed model is proposed. Scenario analyses are therefore proposed for assessing the environmental impacts generated by maritime transport in the Mediterranean basin as the freight flows vary. The analyses are carried out using a simulation approach which, starting from hypotheses on the economic and social development trends of the countries of the southern shore, made it possible to evaluate the variations in terms of freight flows and environmental impacts in terms of pollution.


Introduction
The Mediterranean Sea area, place of the most of European shipping, is strongly affected by particulate matter emissions, black carbon, nitrogen oxides and sulphur oxides. More than 600 freight yards are in the Euro-Mediterranean area (more than 20 located in the European Union) and 500,000 employees are in the maritime transport sector, generating an added value (GVA) equal to 27 billion of Euro. The Mediterranean fleet is made up of approximately 8,000 ships, representing 210 million gross tonnages (Randone et al. 2019).
In addition to international and crossing traffic (Suez Canal-Strait of Gibraltar), the Mediterranean Sea is affected by the freight flows moving among its shores (in the North-South direction and vice versa), but also by cross freight flows (direction East-West and vice versa); these are typically Short Sea Shipping flows.
There are different types of Short Sea transport as containers and bulk handling, general cargo, Ro-Ro ships, etc. An important role is played by the traffic of liquid and solid bulk, with a share equal to approximately 58% of the total goods moved in the Short Sea segment.
In the last 20 years, the Ro-Ro traffic in the Mediterranean has exceeded the threshold of 4 million units, with a growth of 255%, so that, transport by ferry has reached the containers transport. In 2017, 70,000 Ro-Ro vehicles were detected in the Mediterranean, an increase of 7.4% compared to 2012 (SRM 2018). Eurostat data show that in the EuroMed area Ro-Ro traffic in 2018 handled about 81 million tons of goods (about 14% of the total). 22 countries bordering the Mediterranean Sea have begun negotiations for the establishment of a low-emission zone (MED ECA), such as the one in force in the Baltic Sea, the North Sea and the English Channel for the protection of the whole eco-system.
In this scenario, the analysis of trade flows and their spatial and geographical distribution becomes fundamental in order to evaluate the current and predict future traffic, and to estimate the emissions generated by maritime transport in the Euro-Mediterranean area.
The paper proposes a gravity model for estimating the trade flows within the Mediterranean basin with specific reference to maritime transport. The gravity model assumes that one of the main factors affecting trade is the economic dimension of a country which is directly related to the volume of imports and exports: economically "stronger" countries on the one hand produce more goods and services and therefore can sell more on the foreign market, on the other hand they generate more income and therefore their residents are able to import more.
The gravity model is widely used in the literature for the estimation of freight flows with reference to limited geographical areas. The work proposes a calibration of the model with reference to a wide area such as the Euro-Mediterranean one, which has never been considered until now.
The model wants to be more agile, faster and more friendly than the literature models that consider a set (often numerous) of proxy variables with the aim of explicitly taking into account environmental conditions (political situation, conflict, social status, trade agreements) that could affect the freight flow among states.
After a brief state of the art about the gravity models applied in the literature, the different phases of the model construction will be introduced: survey and building of the database, specification, calibration and validation.
Therefore, some scenario analysis will be proposed for the assessment of the environmental impacts generated by maritime transport in the Mediterranean basin according to the variation of the freight flows. The analysis has been carried out through a simulation approach; some hypotheses have been made on the economic and social development trends of the countries of the southern shore, to evaluate the variations in terms of freight flows and environmental pollution.
Future scenarios include an option based on the use of "green" ships with clean technologies (zero-carbon fuels) through a mix of solutions: LNG, electrification, liquid hydrogen/ammonia, hybrid engines.

Literature review
The availability of a good, in a given area, for the final consumption is made possible thanks to its production in another area: the resulting movement between the production and the consumption area can therefore be represented by a matrix of demand for goods. In the specialized literature, the simulation of demand exchanges is made according to two different approaches. In the first case, the export and import demand are assessed separately using a generation-distribution model and an attraction-acquisition model respectively. In the second case, the joint definition of the import/export demand is carried out using a gravity model.
In world trade theories, the universal gravity model gives a positive connection among the international trade flows and the country size, measured by its GDP and the distance between the import country and the export country. The economic dimension is important in the model because larger countries generate more income from the sale of goods and services, therefore their residents are able to import more. The distance, instead, affects the ability to create contacts and facilitate communications, which affect the trade.
The model was firstly formulated by the economist Tinbergen in 1962 (after reviewed and consolidated by Pulliainen in 1963 andBergstrand in 1985), and it is based on the idea that the geographical proximity facilitates trade for several reasons, from low transport costs to institutional and language similarities among states.
Other authors have used the gravity model to evaluate the impact of some economic policy choices, both in its basic and extended version to include variables of different nature. Such models are sometimes complex as they consider a set (often numerous) of proxy variables with the aim of explicitly take account surrounding conditions (political situation, conflict, social state, specific commercial agreements) that could influence the freight flow among states. Table 1 shows some literature studies related to the estimation of trade flows by a gravity model. The reference areas and temporal horizon of data are specified.  1975,1980,1985,1990,1995,1999 Kandogan (

Gravity models for the freight transport. The methodology
The proposed methodology follows the Timbergen's Gravity Model: where F ij is the export flow from country i to country j; G is a constant; M i and M j represent the economic mass of the exporting and importing countries respectively (generally represented by GDP and, sometimes, by the population or by both variables); D ij is the distance between the countries; β k are the model parameters.
The logarithmic transformation of the previous equation translates in the following linear function: where α is the intercept and the sum of β n is the slop of the straight line. Several hypotheses have been formulated regarding the determination of the optimal distance measure, including the geographical distance among the countries "centers of gravity" or among the country's capitals. Other proposals consist in measuring the distances directly through transport costs. Other variables can be combined with the distance in determining the impedance factors, such as the technological distance, understood as the difference in the technological level between the two countries.
According to Head (2000), β 1 and β 2 should be equal to one but, usually, in the several empirical works the values are between 0.7 and 1.1. The comparison between β 1 and β 2 can be used for some economic interpretations: according to some authors (Feenstra et al. 2001) if β 1 ≥ β 2 , there is the so-called home market effect which occurs when the output supply of a country grows more than the domestic demand, thus stimulating export flows. β 3 should be negative (β 3 < 0) because the distance is used as a "proxy" variable that is as an approximation of the transport costs and other factors such as the perishable nature of the goods, the communication and transaction costs and the "cultural distance", that constitute an obstacle to international trade. In fact, in Head's study, the average value of β 3 is equal to -0.94: this can be interpreted in the sense that a doubling of the distance would have the effect of halving trade between the two countries (Head 2000).

Gravity model for freight flows in the euro-mediterranean context
A gravity model for the estimation of the freight flows in the Mediterranean Sea has been specified and calibrated. It is simple to be easily applied to a large study area but, at the same time, it retains some relevant modelling characteristics such as to allow a not excessively simplified representation of reality. The gravity model appears suitable for modelling demand flows for the Euro-Mediterranean area. The specific reference is to maritime transport and, in particular, dry bulk, liquid bulk, containers, Ro-Ro and other cargo transport.
In the phase of model calibration, the import/export flows in the Euro-Mediterranean area among 18 countries (8 on the northern shore and 10 on the southern shore) have been considered (Fig. 1).

Fig. 1 Context of model calibration
The information and data useful for specifying and calibrating the model have been obtained through a complex procedure that required the acquisition, processing and aggregation of data from several online databases. In particular, Word Bank Group, World Trade Organization and EUROSTAT database have been consulted. On the basis of the information obtained, the commercial profiles of the States bordering the Mediterranean Sea have been carried out, indicating, for each of them, the GDP, the population, the trade flows in terms of exports and imports moved by 5 different types of maritime transport: liquid bulk, dry bulk, containers, Ro-Ro and general cargo. The data have been systematized and organized in a synthetic, modular, potentially and periodically updatable database.
The freight exchanges in the area have been evaluated through the composition of appropriate Origin/Destination matrices. Table 2 shows the O/D matrix of export flows of container transport in 2019.
The reference model is a gravity model in its logarithmic shape in which the dependent variable (F ij ) is the logarithm of the export trade flow between country i and country j; the three independent variables are represented by the logarithms of the economic mass of the origin country i (M i ) and of the destination country j (M j ) and by the logarithm of the distance between them (D ij ): The economic mass has been assessed by referring to the GDP of the single country; the distance has been evaluated as the intermodal temporal distance, that is, the time necessary to ensure the connection between the exporting country and the importing one by intermodal road-sea transport (Fig. 2).
The distance has been assessed by taking as reference the main container port of each country; the origin and destination of the commercial flow have been in barycentric position of each country area, considering that, in the transport sector the barycentre is chosen in relation to the activities and the population density of the area. In African countries, the centroid is located towards the Mediterranean coast.
The temporal distance T ij has been defined as: where • T io : time to reach the port o of the exporting country by road starting from the barycentre i of the same country; • T od : time for the maritime connection between the port o of the exporting country and that one d of the importing country, including manoeuvre and hoteling times; • T dj : time to reach the final destination (that is the centre of activity j) starting from the port d by road mode.
The calculation of the travel times on the road has been carried out by an aggregate model:  where L i is the distance by road travelled in the country i and v i is the commercial speed varying in relation to the specific country.
The travel time by sea (T od ) has been assessed considering the navigation time and the time lost at the port for manoeuvring, loading/unloading and retooling (hoteling time). Table 3 shows the average speeds and the average times at node for different types of ships.
Different time distances matrices have been obtained in relation to the type of the maritime service. By way of example, Table 4 shows the matrix of the time distances for container transport.
The individuation of the model variables has been carried out on the basis of statistical considerations through the analysis of regression matrices, scatter plots describing all the possible correlations between pairs of variables.
The matrices construction has highlighted the type and degree of relationship existing between the response variable, that is the flow in export from i to j, and the explanatory variables of the problem; in addition, the matrices made it possible to assess the degree of relationship between the independent variables examined, allowing for checks on their independence.  The model calibration, for the different types of maritime service, has been carried out with reference to a multiple linear regression model and using the least squares method for estimating the unknown parameters. The calibration results are shown in Table 5.
For goods moved by liquid bulk, container and general cargo transports, β 1 is greater than β 2 for which there is monopolistic competition that produces the socalled home-market effect. Monopolistic competition is a market structure that combines elements of monopoly and competitive markets. Country i (origin) is configured as a net exporter, i.e., the commercial supply produced generally exceeds domestic demand with an accentuation of export flows, therefore the exports of country i are more sensitive to its own income than to that of the destination country.
For goods moved by dry bulk and Ro-Ro transport, this relation is reversed (β 1 < β 2 ) and the exports of country i (origin) are more influenced by the income of the destination country j.    The intermodal distance weighs in a decisive way in the development of commercial relations between the countries considered; in fact, β 3 assumes values higher than or close to one (except in container transport) demonstrating that the doubling of the distance corresponds to a halving of export trade flows between the countries.

GR ES FR HR IT CY MT SI AL TR DZ LY MR TN EG IL LB SY
For the validation of the models, in the first instance, the evaluation of the determination coefficient R 2 and of the multiple coefficient of determination R have been carried out. Table 6 shows the validation results.
The linear regression models proposed for the different maritime transport provide multiple R values close to 1, demonstrating the strong linear relation between the observed values and the corresponding estimated.
Even the values of R 2 are close to one, it is therefore possible to state that the variability of the export flow is satisfactorily explained by the explanatory variables GDP exporting country, GDP importing country and intermodal temporal distance through the regression model with a good linear correlation.
There are also low standard errors ranging from 1.2 to 1.5. After evaluating the "goodness of fit" of the calibrated models, through the F test (Table 7), the presence of a significant relation between the dependent variable and the set of explanatory variables considered has been verified.
The significance value of F is well below 5% therefore the results obtained are statistically significant (i.e., reliable), so the model used is good.
Despite the goodness of the model, the analysis and comparison between the estimated and real values led to consider a correction factor G in the transition from logarithmic to exponential notation (Table 8).

Emissions estimation from maritime shipping
Maritime transport generates numerous impacts on the marine ecosystem due, for example, to the engines, the spillage of oils and other contaminating liquids into the water, the loss or discharge of solid materials, the emission of air pollutants. In the paper, attention is paid to emissions into the air (NO x ; SO 2 ; CO 2 ; HC). CO 2 contributes to global warming by trapping heat in the atmosphere and negatively affects marine ecosystems by increasing the acidity of sea water. Sulphur dioxide (SO x ) emitted by the maritime sector contributes to acid rain, which has a negative impact on health. Nitrogen oxides (NO x ) are gases that can cause acidification and eutrophication of water and soil; in fact, NO x emissions on the one hand increase the presence of nutrients in sea water leading to abnormal growth of algae, on the other hand they favour the formation of particulate matter and ground-level ozone.
The evaluation of the emissions generated by a pollutant p (E p ) can be evaluated using the following expression: where EF p is the emission factor of the pollutant p in g/KWh; T od is the time of navigation by sea from the port o to the port d; P s is the engine power for the s ship type; N sy is the number of s ships type offering the maritime transport service in year y.
Tables 9 and 10 give respectively the values of the emission factor pollutant and the average power values of the main engines by type of ship.
The time T od can be evaluated as described previously.
The number of s ships type for the transport service s in year y can be evaluated as: E p = EF p · T od · P s · N sy N sy = F sy /C s  where F sy is the flow of goods moved by the transport service s in year y; C s is the average capacity of the ships operating the maritime service s. This average capacity (Table 11) has been assessed with reference to a large database built by a detailed survey which included the monitoring of ships in the Mediterranean Sea and the detection of their characteristics.
Using the model, the polluting emissions generated by intra-Mediterranean maritime transport have been evaluated, in 2019 ( Figs. 3 and 4). Intra-Mediterranean maritime traffic flows generated the emission into the air of about 1,800,000 tons of pollutants of different types. About 96% of emissions concerns the production of CO 2 , 2.35% the NO x , 1.63% SO 2 and 0.08% HC.

Case study and scenario analysis
The specified, calibrated and validated model has been used to analyze full-scale scenario assessments. The analysis has been carried out using a simulation approach (what if ) and making some hypotheses on the economic and social development trends of the southern Mediterranean countries. The variations of the freight flows and environmental impacts in terms of pollution have evaluated. Considering a temporal horizon of 40 years, three different scenarios have been envisaged: • Scenario 0: current state; • Scenario 1A: GDP increase of 20% for the countries on the Mediterranean southern shore; • Scenario 1B: evolution of Scenario 1A with renewal of 20% of the maritime fleet using "green" ships with clean technologies (zero-carbon fuels) and a mix of solutions as: LNG, electrification, liquid hydrogen/ammonia, hybrid engines. The use of "green" ships would allow an average reduction of 40% of the emission factor of each kind of ships.
The time horizon of 40 years has been chosen in relation to the evolutionary dynamics of the market; some changes in the structure are expressed in a medium-long time horizon. The scenario hypotheses are from the direct experience of the authors and from the bibliometric analysis; in particular, many analysts point to a higher GDP growth trend for countries on the south eastern shore than for European countries. Cautiously, an annual growth of 0.5% has been assumed. Starting from the concept that an average life span of a fleet is around 50-60 years, a share of the vessels could be newly built, a 20% renovation is reasonable.
Finally, 40% reduction of ships emission has been defined in relation to some bibliographic studies that foresee a drastic reduction of emissions in relation to the technological innovation of ships engines. The 20% increase in GDP of the countries on the southern shore of the Mediterranean area generates an increase in goods flows different in relation to the maritime service considered (Fig. 5). For Liquid Bulk and Dry Bulk transport the increase in export flows is 4%, for container traffic it is about 2%, for General Cargo it is about 13%. Ro-Ro traffics show the lowest percentage increase, which is around 1%.
The increase in traffic flows determines a consequent increase in the number of ships moving in the Mediterranean basin in a year. In scenario 1A, a significant increase in polluting emissions emerges on average equal to 5% for Liquid Bulk Dry Bulk, Ro-Ro and Container services around 60% for General Cargo services.
In scenario 1B, the renewal of 20% of the fleet by using innovative and green ships leads to a drastic reduction in emissions compared to scenario 0 on average equal to 8.5% for all types of transport except for the general cargo service (average reductions of 48%). Compared to scenario 1A, the reduction of emissions is on average equal to 8% for all types of pollutants and services (Figs. 6 and 7).

Conclusions
The paper has focused attention on trade flows in the Euro-Mediterranean area. The specific objective of the research has been to quantitatively analyze the extent of commercial traffic within the Mediterranean basin between the countries bordering it. By the use of data from various sources, the commercial profiles of each country have been constructed, focusing with particular attention on maritime import/export flows (liquid bulk, dry bulk, container, Ro-Ro, general cargo). In detail, the export exchange within the area has been evaluated through the composition of appropriate Origin/ Destination matrices. The collected and systematized data have been used for the specification, calibration and validation of a log-linear gravity model useful for estimating commercial traffic in the Euro-Mediterranean area.
The gravity model had been reliable in giving the actual trade flows; the results have showed a good adaptation of the model to the data of commercial use of exchanges between the countries of the Euro-Mediterranean area, confirming the importance of the variables used in the phase. specification of the same (GDP and intermodal distance). Assessments have been carried out in terms of polluting emissions into the air generated by maritime transport which made it possible to highlight the importance of projecting the maritime sector towards a green turning point that can allow a significant reduction in emissions while still guaranteeing the possibility of conveying flows by sea important goods.