Maritime Economic Modeling : A Survey , 2006-2016

We surveyed quite a number of papers (~60), where one third of them belonged to maritime economic dynamic systems, meaning multi variable models. The emphasis we placed was on works published the last 11 years or so (2006-2016). In addition, we provided a general picture of maritime economic research from 1973 to 2013. The survey had to be combined with advances in General Econometrics. As a result and despite long delays maritime economists (Glen et al. [1] applied the VAR—vector autoregressive (1987-1988) model due to Sims [2]. But this was not the only one as co-integration [3] [4] dominated the scene given that SEM-simultaneous equations model had defects in relation to exogenous variables. Moreover, VEC—the vector error correction model and a representation of VAR appeared based on Granger’s representation theorem. The models for 2006-2016 were 36 including 5 for maritime futures after 2003 and 7 on forecasting (2003-2016). A description of a shipping cycle has been given with the presence of time.


Introduction
A "survey" 1 is useful, we believe, because it traces the retrospective evolution of a science 2 .It also reveals the deficiencies remained after the research of econo-1 For maritime economy: Woo et al. [5]. 2 The papers appeared in the following journals ignored by Woo et al. [5]: "Maritime economics & logistics" (ME&L thereafter; published in 1999); "Shipping & transport logistics" (published in 2009); "Ocean Systems Management" (published in 2009).Maritime papers also appeared in transportation journals [7], these also ignored.

A. M. Goulielmos
the above general trend.
The IAME 2014 conference provided a list of the future research (conferences had to present what will be published tomorrow); 150 papers were accepted: "Shipping" 13 : 37%, with emphasis on modeling; "Ports" 14 : 30%; "Maritime logistics": 18%; "Environment":15 9%; "Maritime Geography": 4% and "Maritime Education/Training": 2% (=100%)."Maritime Geography" appeared there for the first time."Maritime logistics-ML emerged as a 3rd branch of maritime economics judged by the number of papers published; some [15] argue wrongly that ML is a discipline.ML deals with ports and port areas (hinterlands) with ships and other means of transport [18].
Modeling for the first time (2014) claimed a serious presence in an IAME conference.There is, however, a deficiency as much as "shipping forecasting" is concerned, which, despite its superior importance in business life, attracted a li-

Aim
This paper has 7 purposes: 1) to present briefly the (published) research production of maritime economists from 1973 (to 2012), on the economics of "dry cargo" and "tanker" shipping.2) To show selectively the ability of maritime economists to create seminal models from 1981 (to 2016), with special emphasis in last 11 years (2006-2016).3) To underline the religious devotion of maritime econometricians to GARCH model.4) To focus on most recent maritime econometric modeling since 2014 (to 2016).5) To record the efforts of maritime economists to forecast shipping markets and 6) to give a description of the mechanism of the "shipping economic cycle" taking into account time.

PAST Research Production of Maritime Economists (1973-2012)
This part presents the main points of a previous survey carried-out by Woo et al.
(2013) [6] over the last 40 years, restricted to the publications of one journal only (MP&M).
Worth noting is that whatever papers are published is not the unique responsibility of authors.We believe that the subjects of papers published in journals are firstly selected by their editors, and then by their reviewers.For example, the papers about "Marine Economics" 16 , when the late Prof. Goss R took over in 1981 as MP&M Editor, stopped.Moreover, journals seek a particular subject in selecting topics for "special" issues.Moreover, research is greatly helped by the degree of how much open (maritime) industry is to joint research with Academia 17 , as mentioned.
The progress of sciences in USA, we believe, is due to the joint research between Firms and Academia, especially at the level of doctoral theses.Shipping, unfortunately, is a "secretive" industry and its close cooperation with academia, though much improved, it is still desirable, especially for doctoral students 18 .

Frequency of Topics Appeared and Their Ranking
The papers published in the 2 journals MP&M and ME&L have been classified by the frequency they appeared, and according to their rank (Table 1).Table 1 shows the ranking of the first 10 only areas of research (2001-2012).
As shown, "Air emissions" appeared for the first time; MP&M shows that it is more dedicated to "seafarers' issues"; ME&L attracted more papers about "port performance", "supply chains" and "maritime security".1.1.1The research preference of maritime economists.
16 Goulielmos [11]; as a result MP&M focused thereafter to 3 industries instead of at least 10 previously.This may sound as a specialization.A doctoral student considers him-/her-self lucky if 50 international shipping firms respond to his/her questionnaire or to their interviews.19 Surveys: 9%; interviews: 3%; observations, reviews, contents' analysis and archives research: 2%; case-studies: 16%; total 30%.20 The share increased from 7% to 13% in 4 decades.One must take into account that Mathematics is used to capture an objective reality… A. M. Goulielmos    A. M. Goulielmos

Research Methodologies Used (till 2012)
The specific methodologies used among the 2 Journals MP&M and ME&L were in Table 2.
As shown, there are differences in % between the two main journals: 1) ME&L attracted more technical papers applying newer methods (41%) (=17% + 24%), against 20% in the other journal MP&M.Worth noting, is the high % of 53 in MP&M against 33% of ME&L over the papers characterized as "descriptive, legalistic and historical".This is important if this is taken as a quality assessment of papers.Notable is also the 24% in ME&L against 9% in MP&M over the "mathematical techniques" for optimization… Our comments here are based on the assumption that the more technical and mathematical is a paper, the more this is superior…

The Evolution of Maritime Econometric Models (1981-2016)
Modeling is the process of making a representation, (rather simple), of reality, by "observing" it.This definition, inspired by Small [23], focuses on the fact that the "reality we model, is the reality we observe".Models can be derived either: 1) from logic, and from known facts, or from assumptions about reality; e.g.Einstein's equation E = mc 2 is such a model; 2) by specifically mimic observations, especially when there is a lack of profound understanding of reality, or 3) from both.Shipping time series models fall in the second category… The start of modeling of "dry bulk" and "tanker" markets taking into account its topics, methods and data employed, can chronologically be placed by us in 1981.Worth noting is that shipping economics emerged in 1875 as a study of the 1st shipping cartel.

The Beenstock & Vergottis Model (1985-1993)
One of the 1st holistic works in maritime modeling is the book of "Beenstock & Vergottis" (BV thereafter) in 1993, based on 4 published papers between 1985 and 1993, and on Vergottis' doctoral thesis 21 (in 1988).The BV model is a "structural" and "SEM" model, involving "bulk 22 shipping" including prices of ships.They covered also "market's efficiency" (using the 3-stages least squares methodology).This "structural" approach, or "demand and supply" one, can be found also in Stopford [24] in relation to forecasting freight rates.Tsolakis et al. [25] were the ones to criticize BV.
However, the work of BV (1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993) put aside as the research had then to face the "stationarity" problem 23 , where shipping data were no exception.21 His PhD in London City University was on "(an aggregated) econometric model of the world shipping markets analyzing dry and wet, 2nd hand ship and freight rate markets, shipbuilding and scrap"… 22 Ships that transport in bulk one uniform cargo for their entire space, either liquid or dry.

23
A stochastic process is "strictly" stationary if its properties are unaffected by a change in its time origin [26] [27].Source: Compiled from Talley [14].
Stationarity means that a process is in a sense in a "statistical equilibrium" [28].
Based on the existence of trends, which also appear in shipping, research needs… differenced and the co-integrated time series to be stationary.
Moreover, the "rational expectations" hypothesis made by BV model was much more consistent with… VAR 24 .VAR was introduced in 1980s to develop certain operational tests, like "Granger's causality 25 " [1].Moreover, a natural extension of VAR was VARMA (= vector ARMA).Econometricians believe that VARs were good forecasting models, as argued by Qin [29], but not good structural ones [30].
"Maritime econometrics" committed also-as general econometrics did-the "spurious regression" error, because their explanatory variables were generated by a "unit root" AR(1) and/or by a "non-stationary" VAR.The above problem solved by "co-integration".In the late 1980s VAR found a new stimulus in the arrival of co-integration theory due to Granger in 1983 and Engle and Granger in 1987.A co-integration system can be represented by VAR [31], and thus it gave VAR a powerful link to long run equilibrium-based economic theories [32].
Notable is that the "single equation maritime economic models" pioneered by The presence of a "unit root" in an AR process creates a number of problems, one of which is the "spurious regression".Let variables Yt and Xt have common stochastic trends and that for simplicity Xt follows a RW.Let Yt = Bo + B 1 Xt + Ut, where Ut ~IID, (0, σ 2 ), meaning that the trend in Yt is a consequence of that in Xt or both; are I (1) (integrated of degree one).The linear combination (= co-integration) of Yt and Xt is I (0) and stationary.RW stands for random walk.A causal relationship between a Z sub-vector and an X one, in a VAR model, which possess autoregressive representations in a statement that X causes Z, and Z does not cause X, etc.
A. M. Goulielmos time ones, remain far away from reality if they exclude time from their models.
Moreover, VAR is superior as dynamic.Econometrics maintains the assumption, by now, that almost all variables in an economic system should be explicitly modeled by dynamic processes (e.g.VAR).Sims [1], in particular, was dissatisfied with SEM, because exogenous variables were explicitly left-out of the model [33]…

The Main Classes of Maritime Models, 1981-2016
The following 6 chronological classes of maritime models found during the last 35 years (Figure 2).
As shown, maritime economists started with structural models using structural equations, i.e. the original equations of a model describing a simultaneous system and where the endogenous variables are on RHS.Then the ARCH and CARCH models appeared in 1996 followed by VAR in 1981.In 2004 valuation models connected with assets like 2nd hand ships emerged.Indeed, the "Regression models" applying "ARCH", "ARIMA", and especially "GARCH", followed by "VAR" (Table 3), have dominated the scene.These models borrowed from "finance" and "commodity markets" and applied thereafter to maritime economy.
Table 3 shows analytically the trends shown in Figure 2, naming also the authors and their models involved.

Co-Integration
"Co-integration" analysis changed the focus of maritime research from "structural models" (BV) to the assumption that "the conditions of maritime markets (bulk)" are competitive... Maritime economists thereafter tried to model demand, expectations, conditional heteroscedasticity, seasonality and ship prices.Initials of: "Autoregressive-Integrated-Moving Average" process.A nonstationary stochastic process related to ARMA.If differenced d times (where d = an integer) it becomes stationary.ARMA is a stationary stochastic process mixed with AR and MA (moving average).AR is a stationary stochastic process where the current value of time series is related to past p values (p = an integer); if p = 1, then = AR(1) with ∞ variance; if p = 2, then AR(2) related to 2 previous values.Small [23] noted that his data of financial time series exhibited deterministic components assumed by ARIMA & GARCH to be stochastic.

27
ARIMA with an exogenous input.ARCH.Introduced by Engle [44]; stands for the 1st-order "Autoregressive Conditional Heteroskedastic" process.It requires the "conditional variance" (σ 2 ) always >0. "Autoregressive conditional" means that the changes in variability are controlled by data's past behavior.The variance is time-varying and conditional on past one.It exhibits frequency distributions with high peaks at the mean and fat-tails.29 This is an extension of ARCH, generalized, (i.e. made broader to accommodate more circumstances) due to Bollerslev 1986 [45] and in 1988.It is popular for modeling maritime models.GARCH stands for the "Generalized Auto-regressive Conditional Heteroscedasticity" model.It refers to a set of statistical tools to model data, whose variability changes with time (= heteroskedastic).It is used to model market volatility.Engle-sharing a Nobel in 2003 for this development based on the observation that "volatility clusters due to dependence".The model starts with a conventional Brownian model of price variation.When volatility jumps, new parameters are put in to make "bell" curve to grow, and the reverse also holds.The bell vibrates to fit to circumstances.It is used by options traders and finance directors in trying to model risk.But what makes bell to vibrate (?) [45] (Bell curve = normal distribution curve) remains a mystery… Maritime models were first in "reduced form 31 ", and then in "VAR", and later in "GARCH".
Maritime economists-apart from the models presented in Table 3-they also modeled "system dynamics" [20] [58] [59] [60] for a change.Engelen et al. [58] declared that "system dynamics" are totally well-suited for the detailed modeling of strategic and operational behavior in maritime businesses.Their model belonged to the "geometric Brownian motion" family of models.
Radical changes, however, take time to pass on to researchers they come after.
When first we applied chaos theory to model maritime economy in 2006 as a product of the doctoral thesis of my student Psifia M-E, no-one followed us, and the 1st massive quotation of our work appeared in 2016 [21].

Maritime Futures
Maritime economists showed also interest in models in the area of "maritime futures" (Table 4).

Maritime Models from 2006 to 2016
Adland et al. [65] tested the instantaneous equilibrium relationship between 2nd hand ship and contracting prices and freight rates, in a VEC model, and used the concept of the "time-varying delivery lag".Adland et al. [66] argued that there is

31
A general vector ARMA model is related to a dynamic SEM, where simultaneity can be removed.This is reduced form.It expresses each endogenous variable as a function of exogenous-lagged endogenous variables (the "predetermined").The directional forecasting accuracy of FFA and BIFFEX FFAs do not seem to be very accurate in revealing direction of future rates to a degree of 46% and 74% Alizadeh & Nomikos [50] Lead-lag relationships between spots and future rates Bidirectional causality in spots (all routes) No volatility spillovers found Kavussanos & Visvikis [62] Future rates are unbiased predictors of spots Kavussanos et al. [63] Futures trading They had a stability influence on spots' volatility Kavussanos et al. [64] no clear relationship between the size of the ship and her delivery lag.This is not realistic, we believe, as a larger ship needs more time to be constructed, other things being equal.Alizadeh and Nomikos [67] used co-integration between prices and earnings in the "2nd hand ship market"-in an effort to combine "technical" 32 (P/E ratios) and "fundamental" 33 analysis (use of historical data plus technical trading rules)-to find out the best investment or dis-investment timing, and thus the proper trading strategy.The effectiveness of this strategy tested against a benchmark-buy and hold-and the method found superior, despite "data snooping 34 bias" (a bootstrap 35 method used to counteract this bias).They used VEC.
Adland and Strandenes [68] applied "technical analysis" to freight markets considered to be appropriate by combining it with trading rules and profitability.The traditional EMH-efficient market hypothesis rejected on the grounds that freight rates are dynamic and shipping services cannot be stored or traded and replaced by a "semi-strong" efficiency.A trading strategy considered impossible to yield excess profits with no risk-premiums if the market is efficient… They used a "Kernel regression" for chartering decisions to be identified in peaks and troughs in the freight market cycles.32 This is the 2nd oldest form of analysis after "fundamental".This is a craft of recognizing patterns real or spurious of studying reams of price, volume, and indicator charts in search of clues to buy or sell.It expanded in the 1990s as thousands of neophytes took the Internet to trade stocks and insights, especially in currency markets.They use support points, trading ranges and other patterns in the tick-by-tick data.Chartists can at times be correct but they are financial astrologers.It is not a foundation on which to build a global risk-management system (Mandelbrot and Hudson, 2006 [45] pp.[8][9].This is the oldest and simplest method in financial markets, dealing with risk.If a stock is rising, seek the cause in a study of the company behind it, or the industry and economy.Predict stock's next move.The price of a stock, bond, derivative or currency moves because of a certain event or fact, mainly from outside the market.If one knows the cause, one can forecast the event, and manage risk; (Mandelbrot and Hudson, 2006 [45], pp.7-8).34 This is trying many variables in a regression without basing the selection of the candidate variables on an economic theory (Brooks, 2014) [69]).

A. M. Goulielmos
Syriopoulos and Roumpis [70] used an EGARCH and VAR model to identify the current relationships between returns and trade volume in the Handy and Aframax tankers.Price changes affect trading volume.Dikos et al. [60], based on Dikos' doctoral thesis, tried to model a structural relationship of supply and demand type.They examined tanker rates, 1980-2002, by using GARCH and EGARCH.Three functions, however, found nonlinear, while MSEs were low.
Alizadeh and Nomikos [71], examined whether certain chartering strategies based on "technical" analysis can generate excess returns by applying a "moving average-MA trading rule".This rule was to charter a vessel on a 12-month time, and charter it out on a 6-month, if the difference between the 2 rates is greater than the average difference over the n prior weeks.The MA approach showed the possibility of yielding excess profits using especially a rolling sample.
Goulielmos and Psifia [72] showed that the time series for the "trip and time charter indices" between 1968 and 2003 (428 months) and between 1971 and 2003 (391 months), were not "identically and independent distributed-IID".The above confirmed by Zheng and Lan [21].Moreover, they [21] showed also a nonlinear dependence on the basis of the BDS 36 test.A deviation of 5σ admittedly unexpected and high occurred not once in a century according to theory in the maritime industry.Time charters were less risky as their volatility was ≤3σ, except in 1974.Trip charters were more volatile with σ > 3 in more than one year case.

The Hurst Generalized Model of Randomness Due to Einstein (1905)
Einstein [74] studied extensively "Brownian motion 37 ", postulated for the 1st time in 1828 by Brown R, who observed it, but he could not prove it.This model became 38 the main model for "random walk" in the study of statistics.Einstein showed that the distance covered by a random particle undergoing random collisions from all sides, is directly related to the square root of time: R = the square root of time = k√T (1), where R is the distance covered, k is a constant and T time index.
Hurst [75] generalized [76] [77] Einstein's Equation (1) to include speeds adopted by time series below or above the square root of time.He used a dimensionless ratio: R/S = kT H (2) 39 , where R is the range (maximum value minus minimum value of time series), S is the (local) standard deviation, R/S is the 36 This is called BDS statistic and is based on the correlation integral that examines the probability that a purely random system could have the same scaling properties as the system under study.Named after its originators: Brock, Dechert and Scheikman [73].A computer code for BDS estimation is now widely available free of charge in the internet.37 Originally studied as the erratic movement of a small particle suspended in a fluid; this movement is due to water molecules colliding with the particle.38 Bachelier, five years before Einstein, in 1900, recognized the relationship between a "random walk" and the "Gaussian statistics".Weiner in 1976 modeled the "Brownian motion" as a "random walk" with underlying "Gaussian" statistical structure.Feder in 1988 explained this process (Peters, 1994, pp. 18-187 [76]).

A. M. Goulielmos
"rescaled range", k is a constant for time series, H is the Hurst exponent (a power law) and T is the index of time referring to the number of observations (sometimes denoted by n).
The advantage of R/S is that, by dividing R by S, the range is made independent of time.This means that data belonging to different times no matter how distant are comparable.We need no time adjustments like the ones made by Hampton [78].It also does not matter if time series has no characteristic time scale.
H can be equal to 0.50, which is Einstein's special case, or "random walk"; also, 0 ≤ H < 0.50, which is the case of "pink noise"/anti-persistence; the system in this case reverses more often; and 0.50 < H ≤ 1, which is the case of the "black noise", showing the "Joseph effect" (existence of trends and cycles in time series) and the "Noah effect", (the potential of sudden catastrophes in time series).Interesting is that almost all shipping time series, we have dealt with, showed "black noise".So these time series are probable candidates for sudden catastrophes, as happened in 1929-1933, in 1981-1987, in end 2008 and in early 2016 for the dry cargo market.

Garch Maritime Models
Jing et al. [17] examined the freight rate volatility in dry bulk shipping markets, recognizing the existence of an increased complexity, given the "incredible" Bessler et al. [80] tested the existence of autocorrelation and co-integration in the spot freight rates and the FFA prices.They found that spot freight rates are auto-correlated, while FFAs are not; and both are co-integrated.[81] were the first to apply the "value at risk" model-VaR 40 , (described below in forecasting models), to estimate risk in freight markets.They found that the simple non-parametric methods should be used to measure risk in the freight markets; e.g. the "Rescaled range analysis" is a non-parametric method.Alizadeh and Nomikos [82] analyzed the performance of the "regime-switching volatility models" in the freight market.They found that by incorporating regime changes into volatility models, they improved the precision of VaR estimates.[83] argued that ship prices are very volatile, given the observed "mean reversion" in freight earnings.They proposed a behavioral model where firms over-extrapolate exogenous d shocks and partially neglect the 40 VaR is written this way to avoid confusion with VAR.

A. M. Goulielmos
endogenous investment responses of their competitors.Their formal estimation of the model confirms that both types of expectation errors are required to account for empirical evidence.
Kalouptsidi [84] examined the impact of the "construction lag" and its extension in times of high investment activity by constructing a dynamic model of entry and exit from shipping markets.She found that moving from time-varying-to-constant-to-"no time" to build this reduces prices, while both increase the level and the volatility of investments.
Adland and Jia [85] tested whether "new building prices" are inelastic to "demand".They assumed that, for a shipowner, a "newbuilding contract, a resale case, (vessel sold during her construction), and an existing modern vessel" are "near perfect substitutes" (ignoring technological differences).They also assumed that the integration of the newbuilding, 2nd hand ship and freight market is feasible.This, we reckon, is a try to fit a model to theory, instead of fitting a model to reality… Nevertheless the model is interesting.
According to business economics [86], the one of the 3 conditions for goods to be close substitutes is to have the same, or similar, occasion for use [87].Given also the volatility of the freight market, the timing of entering into the freight market is of great importance (for us) e.g. a Cape earned $200,000/day in June 2008 and less than $5000 in Nov. 2008 (see graph 3 below for 2010).So, the delivery time should be a function of the level of freight rates… (for us).
Authors tried unsuccessfully, we believe to prove that despite subsidies to shipyards and/or to owners, which clearly affect newbuilding prices shipowners are indifferent between 3 cases: 1) ordering a "new building at her contract price 2) buying a resale, or 3) buying an existing vessel".They also assumed that the 2nd hand ship market for bulk carriers is closely co-integrated with those of the freight rate and new-buildings; moreover, they found no evidence of a "short term asset bubble" in the past boom years… Finally, authors assumed that newbuilding prices are not comparable across time due to the nature of the new contracts (which differ mainly on delivery time and on payments schedule).
They used a single equilibrium equation model, relating today's "prices" in all 4 ship markets (new-buildings, 2nd hand ships, freight rates and scrapping); it exists a "term structure of newbuilding prices" and a "time-varying delivery lag", which indicate that newbuilding prices and 2nd hand ship ones are different.
The model shows also how a "time varying delivery lag" positively correlated with freight rate will produce inelastic newbuilding prices and the market remaining efficient.
The model uses the classical "net present value" NPV method of a cash flow generated by a modern ship, between now (t) and future (τ): ( ) , where TC is the time-charter today and for τ-t years; r is the (risk-free) rate of interest 41 for the same period and OPEX 42 is the daily operating cost.Then the scrap value of the vessel is used as-

41
A risk premium?

A. M. Goulielmos
suming a linear depreciation method following Adland and Koekebakker [88] with life expectancy constant.They assumed also that investors do not expect the "scrap value" or the "expected useful economic value" to change...
The payments now to shipyard for a new building have the following present value: ( ) 5), where S τ,0 is the 2nd hand value of the vessel today.This equality may be not true as the amounts agreed to be paid to shipyard may not be equal to vessel's 2nd hand value on delivery.The 1st owner has to pay the price agreed.The 2nd owner faces a price on delivery, determined by market (supply and demand) then, which may be much higher or much lower than the one paid by 1st owner. Antiques , where S τ,τ is her 2nd hand price, Z t is her scrap value, and T τ is her expected life.In equilibrium, shipowners are assumed to be indifferent between the above 2 alternatives.Equation ( 7) is then inserted into Equation ( 6), and by using Equation ( 5), and assuming "myopic investors", i.e., Z τ = Z t and Τ τ = Τ t , one gets: Equation (8) gives the balance between "prices" of all 4 shipping markets today, where: S t,0 is the value of the new vessel today, TC is her time charter, Z t is her scrap value, and F t,τ is the value of her contract.If τ = 0, her delivery is immediate.Interesting is that this formula can be used to find out the newbuilding (contract) price F, given the price today of a new vessel S, her time charter TC and her scrap value Z.The delivery lag τ-t can also be determined as a function of time.They used also a numerical example using a Cape.
Worth noting is that the Cape market mentioned above in [85] suffered by extremely violent fluctuations (Figure 3).A finding useful for our "theory of depression" is that laid-up tonnage determines the short run fluctuations (ups and downs) of the freight rate from one month to next two say; the delivery lag affects supply after ~2 to ~4 years, depending on the phase of shipping cycle, while scrapping may need 5 to 9 years additional time to exert its full power on supply.

More Recent Models (2014-2016)
Zhang et al. [89] used a Brownian distance correlation method and a Granger causality analysis and applied them to the new ship market considering it relatively distant from 2nd hand ship and freight markets before and after end 2008.
Crisis made markets more closely related.
Zheng and Lan [21], investigated the dynamics of the spot rates for all main sizes of tankers, using a "multifractal de-trended fluctuation analysis" (MF-DFA; due to Kantelhardt J W et al. [90].Data: Jan. 27 th 1998 to Aug, 5 th 2013; from 2011 to 3884 days; using exponential ln freight rates and their 1st differences data found stationary; the IID hypothesis rejected by the BDS test.Normality is also rejected by the JB test.They concluded that the tanker market is more fractal than the crude oil commodity markets.What does this mean?
Original maritime research dealt with full scale econometric models, like those used by BV; later research, during the last 15 years or so, dealt with direct specifications of freight rate processes.Most studies made the assumption that tanker A. M. Goulielmos markets are "efficient" or "weakly efficient" in the sense of Fama in 1991.
Maritime economists almost exclusively used the traditional econometric models of ARIMA, VAR, VECM, and GARCH.They came to the conclusion that EMH does not apply to tanker freight markets, because… they are "untradeable services"… Finally, they [21] justified Goulielmos and Psifia [19] by saying that research reveals that the "tanker freight rates" are nonlinear, non-normal, non-stationary, and have memory effects!They [21] concluded also that the Hurst exponent of the daily rate returns showed that these markets are fractal but not multifractal 43 … The terms fractal-multifractal are used wrongly interchangeably by the authors.The paper, therefore, is not clear about fractality and multi-fractality.Fractal means that parts of time series are self-similar with the whole.Fractal is a pattern or object, whose parts echo the whole, only scaled down.But multi-fractality is different, because some parts of the object shrink quickly and others slowly.As argued by Mandelbrot & Hudson [45] "to improve almost any fractal model it is a good idea to replace it by a multifractal one".
The cycles found [21] were all short term, from 2.1 months to 6.6 months (on a 30 days calendar month).The Hurst exponent found between 0.4 and 1.It is not the same thing to have H < 0.50 and H > 0.50.Large tankers have relatively stronger multi-fractality and higher volatility, being commercially inflexible and having thin market conditions [21].
Most tanker markets have long term memory.Multi-fractality comes from both memory effects and fat tails, where the first is stronger.Finally, they found [21] that the crude oil markets are less multi-fractal than tanker markets; due mainly to fat tails; they exhibit longer non-periodic cycles.
Tsioumas and Papadimitriou [91] investigated the "excess returns" issue, between "time charter" and its corresponding "voyage" (using a TCE) in dry bulks (Cape and Panamax) for 2 t/c trips of BPI and BCI.They used weekly data: Jan.
2003 to Jan. 2014, to find-out a long term co-integration.TCE = time charter equivalent.
Their paper falls into "technical analysis" pioneered in 1980s by the "Norwegian School", i.e., by Norman [92], Adland [93], Adland and Koekebakker [56], Alizadeh and Nomikos [67] [71] and Adl and and Strandenes [68].Co-integration [3] [4] found in the long term.They then used the MMTM ("Modified Momentum Trading Model") due to Wilder in 1978 to formulate a chartering strategy, which is to get into voyage charters every week.This was compared with the benchmark strategy of always chartering in vessels on a t/c trip.The existence of excess returns confirmed.t/c = time charter.
They also used the "Relative Strength Index"-RSI at t + 1 using a 4 week moving average.The "voyage" is profitable when RSI moves above its MA (4), where MA(4) is a moving average of 4 periods.Authors could especially with Capes use 43 We believe that the way freight rate changes in shipping markets can cluster into zones and slow evolution.
A. M. Goulielmos monthly moving averages, which according to Figure 3 indicate the greatest fluctuations since June 2010.
Unfortunate for them is that their finding that all rates have right tales (+ skew) and kurtosis greater than 3 (leptokurtosis) and significant departures from normality… destroy the accuracy of their findings.Time series found non-stationary, if not in log-1st differences.Using MTMM one could improve excess returns from 15% to 31% for the Cape over Carter/Rotterdam; and from 6% to 10% for Panamax over Maracaibo/Rotterdam.The study confirmed that voyages are more volatile and that Cape had higher margins than Panamax [94].

Maritime Economic Models (2003-2016) in Forecasting
This area is one of the most important, but as Stopford [24] said: "maritime forecasting has a poor reputation".This opinion is shared also by the President of the shipping companies "Dynacom" and "Dynagas" Mr. George Procopiou G. [95].Maritime economists have indeed paid hard efforts to forecast freight markets and especially dry bulk freight rate markets, but in vain.
As we saw above (Table 4), Kavussanos and Nomikos [61] argued that the VECM generates the most accurate forecasts of spot prices, but it fails in forecasting future prices… It would be interesting to know why.
Batchelor et al. [96], compared ARIMA, VAR and VECM in predicting spot and forward freight rates, where ARIMA outperformed in FFAs; VAR and VECM outperformed, slightly, ARIMA in spots.[59] argued that in shipping, 2 cycles interplay, one 4 years and one 20 years in accordance with the "theory of waves" and on the capability to recognize patterns and construct indicator charts in search of clues for a move up or down; their forecasting for a turning point up in Aframax rates in 2010, proved, however, wrong.

Randers and Goluke
Chen et al. [97] applied ARIMA and VAR to predict freight rates of several dry bulk routes.VAR outperformed out of the sample ARIMA.Leonova and Nikolov [98] used a model based on "wavelet" and "neural network" in order to predict dry bulk freight rates.Bulut et al. [99] developed a "VAR fuzzy integrated logical" forecasting model for time charter rates.Duru et al. [100] used a "fuzzy-DELPHI adjustment method" to improve the accuracy of statistical forecasts in the dry bulk shipping index.
Chang et al. [101], applied the "Value-at-risk" (VaR 44 ) model to evaluate the risk of dry bulk freight rates (BDI, BPI, BCI, BHSI, BSI), showing an asymmetric long-memory volatility process, including expected shortfalls.They used the 44 In this method, you start off by deciding how "safe" you want to be, e.g. to be in 95% confidence level.So, in your investments you want a probability of 95% that any losses will stay below a danger point.In theory a 12% potential loss is advertised with an only 5% probability.The VaR model recognizes that actual loss could be greater than 12% due to the so called overhanging.After one financial crisis after the other, banks under the pressure of the Central banks resorted to more sophisticated risk models like "Extreme Value Theory" (*) found in the insurance market.It assumes that prices vary wildly with fat tails that scale, but bad news comes in flocks (Mandelbrot and Hudson, 2006, pp. 272-274) [45].(*) This theory has not as yet applied to shipping.

A. M. Goulielmos
FIGARCH (fractionally integrated), the Hyperbolic GARCH and FIAPARCH with 3 distributions.The precise VaR estimates are obtained from an asymmetric long-memory volatility system with skewed Student-t distribution.The asymmetric FIAPARCH model outperformed all other models in the out of the sample forecasting!The innovative element of this is the inclusion of a long-memory structure.As we said in the footnote, VaR has indeed ignored long memory/dependence in the time series originally.Santos et al. [102], attempted to forecast 1 month to 6 months ahead, the 1-year and 3-year "period charter rates" of the VLCC tankers, using "Artificial Neural Networks-ANNs" model, out-of-the sample.They covered: 1994 (Oct.) to 2010 (Sept.), or 192 observations.They found that the literature about both ANNs and period charter rates, vis-à-vis spots, is scarce.They used monthly time series, and their innovative element is the use of "Radial basis function"-RBF due to Haykin [103].They compared it with ARIMA.They found that ANNs plus RBF outperform ARIMA and NN-MLP (multi-layer perceptron).Dougherty [104] found only 2 models concerning shipping sector out of 52 (0.04%) using ANNs.
Looking at the results of Santos et al.Zeng et al. [109], forecast BDI, using the "empirical mode decomposition"-EMD (due to Premanode & Toumazou in 2013, and others before) and the "artificial neural networks-ANNs".They found that the "EMD-ANN" outperforms both ANN and VAR.The index decomposed into a number of "independent intrinsic mode functions"-IMF, consisting of 3 components: short term fluctuations, extreme events impact and long term trends.Authors recognized that in the dry bulk shipping market a great deal of non-linearity exists [110]; Goulielmos [111].As the authors argued, it is difficult for traditional econometric and statistical methods to capture the non-linear characteristics hidden in the dry bulk freight rates time series!ANN often, however, suffers from local minima.

A Model of the Shipping Cycle
Shipping cycle is the "outcome of the relationship of demand to supply", and vice versa.If only demand is cyclical and/or only supply is cyclical, then the freight rate will be cyclical.If we want to have no shipping cycles, then demand variations should be smoothly and timely addressed by a very elastic supply.But demand cannot be, and is not, adaptive to supply.The reverse is true.
Maritime economists are confused because, as other economists did, they kept time out of their analysis (except Shimojo as mentioned by McConville [112]. Let us use a graphical presentation in an effort to bring time into the picture (Figure 4).
As shown, as demand increased from D 1 to D 2 , supply increased by 0T 2 -0T 1 , due to ships that they were laid up and brought-in (T 1 T 2 ).This may take from 3 to 6 months for ships on average to be reactivated.The freight rate increased from F 1 to F 2 .As demand remained high, shipowners ordered ships, which when delivered increased supply by 0T 3 -0T 2 (T 2 T 3 ); this caused freight rate to fall to F 3 .This took 2 years on average.If demand does not increase to D 3 , then the new-buildings delivered will cause the fall in the freight rate to F 3 instead of going up to F 4 .If demand falls to D 1 then a surplus of Supply will be created equal to T 1 T 3 .This is a shipping cycle.This cycle took 27 -30 months to be accomplished, but this is quite symptomatic and it is dangerous to consider this as fixed and as valid for all cases.The "construction time 45 " can take from six months to up and near 4 years.The frequent and short cycles are due to the lay-up of ships depending on the tonnage accumulated there from last depression.New-buildings and the time these need to appear-given demand-is the decisive factor, as deliveries may be as much as 100m dwt for only one year (e.g. in 2012 for dry cargoes) and scrapping is unable to reduce it, but only by 33%.In shipping all are matters of time.So, time must return to maritime economy analysis to derive finally more sound conclusions.

Conclusions
Maritime economists since especially 2000 focused more to port economics.Moreover, a new branch formed titled "Maritime Logistics".Cruising and Shipbuilding have been neglected, though during the last 5 years or so we saw a certain number of papers dealing with them.
To renovate maritime research, and modeling, industry's journals have to direct research through planned "special issues".Moreover, non-listed shipping companies must boost joint research and create systematic links with Academia.
There is a serious progress made 1) in the number of maritime doctoral theses produced; 2) in the number of postgraduate degree courses including the e-learning ones 3) in the number of maritime conferences; 4) in the increasing number of published maritime papers and 5) in the number of maritime journals, which published the last 20 years or so.Quality has to surpass quantity though.This is the golden age of maritime research and maritime education, we believe.Journals' editors shape the various research directions.Maritime research produced (in only one journal), since 1973, almost 1000 papers by using 3 methods: "empirical" 1/3; "scientific" 1/3 and "conceptual" (34%).Econometric modeling and statistical analysis which drew our closer attention produced a slim number of 108 studies or ~3 papers per year.Moreover, market analysis, which we consider most important, had only 81 studies or 2 papers per year.Studies for shippers, transport chains.Maritime economists followed general econometrics but with long delays.
Modeling had a slim presence despite its rigor and power of mathematics but more serious is the lack of studies on maritime forecasting.We believe that if Academia wants to help the industry, it must improve forecasting.
Maritime economists are troubled, but at the same time helped, by general econometrics after long delays during this period.This period, stationarity emerged as a must and co-integration appeared as the proper solution not only econometrically, but also from the side of economic theory.We have been supplied with VAR (& VARMA) over the previous defected simultaneous equation models (SEMs)… Between 1981 and 2004, 19 models have been worked out by maritime economists using: co-integration; VAR; ARCH; GARCH; ARIMA; ARIMA-X; EGARCH and ECM-error correction model.There is devotion to linear models.
A pleasant exception was the "artificial neural networks" models with 5 studies.

A. M. Goulielmos
After 2006, the "Norwegian School" which was very active produced 5 -6 models using the "technical" and "fundamental" analyses.Moreover, Goulielmos and Psifia introduced for a change the nonlinear chaotic model of "Rescaled Range Analysis" in maritime forecasting in 2006 and 2007.Eight models that have appeared from 2008 till 2016 continued to use GARCH, EGARCH, VaR, technical analysis and co-integration.
In the area of forecasting, 8 models emerged since 2003, using VECM, ARIMA, VAR, Wavelet/ANN, VAR/fuzzy, fuzzy-DELPHI, VaR, AFIAPARCH and ANN.We are glad that maritime economists have discovered nonlinearity.They used "fractional integration" together with "empirical mode decomposition" and "ANN".
number of papers… The paper is organized as follows.Part I presents selectively the past research production of maritime economists since 1973 (to 2012).Part II presents briefly the evolution of maritime economic models since 1981 (to 2016).Part III presents separately the maritime economic models using GARCH as this is a milestone.Part IV focuses on the more recent maritime economic models since 2014 (to 2016).Part V focuses on models of forecasting maritime economic markets being of our special interest.Part VI provides a description of the shipping business cycle involving calendar time.Finally, we conclude.
Interesting is to see what a long boom can do (like the one or two in 2003-2014): 1) it prolongs vessel's scrapping age by ~145%, from near 21 years in 1993 to 30 years in 2004; and between 25 and 29 years from 2003 to 2014.Vessels therefore remain in service roughly additional 5 -9 years during a boom.This confirms our experience.This reinforces also our all along argument that scrapping is not a fast equilibrating market mechanism.2) It prolongs the delivery lag by ~200%; the lag from roughly 1.8 years in 1993, increased to 3.9 years in 2006, 3.8 years in 2009 and fell in 2014 to 1.9 years.So, "shipbuilding demand" cannot be satisfied immediately, but it needs a "construction plus a waiting time", (for us), which can be extended to almost double during boom years, i.e., from ~2 years to ~4 years(2003)(2004)(2005)(2006)(2007)(2008).The "delivery lag" shows a faster adaptation than scrapping.
, and it can be used to estimate the "term structure of newbuilding prices"; e.g. the newbuilding price of a Cape started at 75 m$ in 2008 and fell to $64 m in 2009, after a delivery lag of 2.4 years; finally, it fell to $45 m after a delivery lag of 5 years.
[102]: the 1 month forecasts usingNN-RBF (2004-2010), were quite close; at times, a $5000 maximum deviation occurred; similarly the 2 months had ~$30,000 deviations in 2008-2009 and in 2005 (3 year time charter) and 3 months forecasts with again ~$25,000 deviations (3 year time charter) in 2008-2009 and in 2005.Their 6 months forecasts were inaccurate in most years.ARIMA and the MLP model underperformed.It seems that the above model shows shortcomings in years of radical changes up, like 2005 and 2008-2009 (peak rates); also when forecasting period is longer than one month…

45 A
variable to be researched over a full shipping cycle.

Figure 4 .
Figure 4. Introducing time into the Freight market analysis.
' market or market of houses, are good examples of the above.This is why Greeks prefer buying used tonnage at rock bottom prices instead of ordering new ships, which during a boom are very expensive, and run the risk of oversupply, and under-demand, upon delivery...