1. INTRODUCTION
An accurate and reliable description of the weather that a ship will experience in its life is essential, not only for design purposes, but also for fleet management, to schedule mantainance and to provide recomendations to the master and the route planner.
Recently several global wave climate studies, using wave hindcasts or reanalysis, have been published (Uppala et al., Reference Uppala, Kållberg, Simmons, Andrae, da Costa Bechtold, Fiorino, Gibson, Haseler, Hernandez, Kelly, Li, Onogi, Saarinen, Sokka, Allan, Andersson, Arpe, Balmaseda, Beljaars, van de Berg, Bidlot, Bormann, Caires, Chevallier, Dethof, Dragosavac, Fisher, Fuentes, Hagemann, Hólm, Hoskins, Isaksen, Janssen, Jenne, McNally, Mahfouf, Morcrette, Rayner, Saunders, Simon, Sterl, Trenberth, Untch, Vasiljevic, Viterbo and Woollen2005; Chawla et al., Reference Chawla, Spindler and Tolman2012; Stopa et al., Reference Stopa, Cheung, Tolman and Chawla2013; Fan et al., Reference Fan, Lin, Held, Yu and Tolman2012). Also some long-term regional wave climate studies using downscaled wave hindcasts were pursued (e.g., Weisse and Günther, Reference Weisse and Günther2007; Pilar et al., Reference Pilar, Guedes Soares and Carretero2008; Lionello and Galati, Reference Lionello and Galati2008). This attempt to investigate sea waves in more detail is not only impelled by scientific interest. Wind waves have a significant impact on offshore and coastal infrastructures, ship design and routing, beach erosion and sediment transport, for example, and are an important element in extreme events, like storm surges and flooding in coastal areas.
In ship design the North Atlantic weather is widely used as a reference for the establishment of the reference values of extreme wave induced loads in ships and it is the one adopted by the International Association of Classification Societies (IACS) to establish the design conditions for ships. The attempts to use different wave databases to determine the most probable wave loads a ship will be subjected to, show a significant scatter in the predictions (Guedes Soares and Viana, Reference Guedes Soares, Viana, Schreffler and Zienkiewicz1988; Guedes Soares and Trovão, Reference Guedes Soares, Trovão, Price, Temarel and Keane1991; Reference Guedes Soares and Trovão1992; Bitner-Gregersen and Guedes Soares, Reference Bitner-Gregersen, Guedes Soares, Guedes Soares and Das2007), and also shows differences from the values along different routes (Guedes Soares and Moan, Reference Guedes Soares and Moan1991; Queirós et al., Reference Queirós, Bernardino, Guedes Soares, Guedes Soares, Garbatov, Sutulo and Santos2012). However the use of the various databases normally leads to higher values than prescribed in the design rules and the explanation is that ships tend to avoid the roughest weather and thus are exposed to somewhat lower sea states. However, the effect of storm avoidance has not yet been fully quantified and this paper makes a step in that direction.
Besides the importance of the IACS recommendations, especially to take into account the extreme events that a ship can be asked to withstand maybe a few times in its life, it is useful to represent the actual environment that ships usually experience in their voyages across the Atlantic. While in the design phase a conservative approach is still recommended, when the schedule of a ship during the year or of an entire fleet has to be planned by assessing operability, knowledge of what conditions can be expected along the actual routes would significantly improve the results. This becomes particularly relevant in the present conditions in which energy efficiency, fuel consumption and emissions are being more considered as a priority (Papanikolaou et al., Reference Papanikolaou, Zaraphonitis, Bitner-Gregersen, Shigunov, El Moctar, Guedes Soares, Reddy and Sprenger2014; Prpic-Oršic et al., Reference Prpic-Oršic, Vettor, Guedes Soares, Faltinsen, Guedes Soares and Santos2015).
As was shown in a previous work (Vettor and Guedes Soares, Reference Vettor and Guedes Soares2015), for geographic and economic reasons, most of the trades in the North Atlantic are carried out along a few specific tracks, so it is there that more effort has been made to provide voluntary observations of weather conditions.
Sophisticated systems are nowadays used on board to choose the route that fulfils the targets, and onshore to provide recommendations to avoid dangerous conditions. This allows ships to avoid some rough weather, an aspect that is inherently contained in the Visual Observing Ships (VOS) wave data. Ships, in fact, provide reports only relative to the location and moment in which they are passing. Consequently, if ships tend to avoid storms, which means avoiding high Significant Wave Heights (SWHs), the VOS database will not contain those events, or at least they will be in a lower percentage than they actually occur in climatological databases. This paper attempts to quantify this effect along the above-mentioned North Atlantic routes.
Section 2 describes the main North Atlantic routes. In Section 3 an overview of the three wave databases used in this study is given. In Section 4 the different climate, dependent on the covered route is described with an attempt to provide all useful elements to adopt as an input for ship performance and response analysis. In this perspective, besides the average significant wave height, also the scatter diagram and the relative wave heading probability are presented. In Section 5 a discussion on the bias encountered in the VOS data and a simple correction is applied, highlighting the effect of storm avoidance on statistics from VOS reports. The concept of storm avoidance correction is introduced in Section 6 and an equation for its estimation is proposed. Conclusions and final remarks are made in Section 7.
2. MAIN NORTH ATLANTIC ROUTES
In Vettor and Guedes Soares (Reference Vettor and Guedes Soares2015) the tracks of the main North Atlantic routes were identified from the VOS database following the area with a higher density of reports and with some consideration of global economy and geography.
Six main trans-oceanic passages were detected as follows and depicted in Figure 1:
• Route 1: Channel - Puerto Rico (North), [Ch_PR-N]
• Route 2: Channel - Puerto Rico (South), [Ch_PR-S]
• Route 3: Channel – Virginia, [Ch_VA]
• Route 4: Strait - Virginia (North), [St_VA-N]
• Route 5: Strait - Virginia (South), [St_VA-S]
• Route 6: Strait – Miami, [St_MIA]
It was shown that the most used stream is the one from Northern Europe (the Channel) to the Caribbean Sea (Ch_PR-N and Ch_PR-S) with 36% of the traffic in the considered areas (higher percentages in winter). More than one-third of these ships prefer the northern and shorter orthodrome.
The routes between the Strait of Gibraltar and the North-East Coast of USA (St_VA-N and St_VA-S) denote a similar concentration with about 32% of the trades, while the routes from the Channel to the Virginia area (Ch_VA) and from Gibraltar to Miami (St_MIA) contribute 21% and 10% respectively.
Rough weather avoidance was found to influence the choice of the routes and it is mainly highlighted on the northern passages, which are always less frequented in winter. Overall, May is the month with the highest presence of ships, whereas in December it reduces by almost 20%.
3. WAVE DATABASES
A wave database contains climatological data in a spatial grid covering a specific area, which can be obtained from simulation models, measurements or observations. Modern phase averaging third generation wave models (Komen et al., Reference Komen, Cavaleri, Donelan, Hasselmann, Hasselmann and Janssen1994; WISE Group, Reference Cavaleri, Alves, Ardhuin, Babanin, Banner, Belibassakis, Benoit, Donelan, Groeneweg, Herbers, Hwang, Janssen, Janssen, Lavrenov, Magne, Monbaliu, Onorato, Polnikov, Resio, Rogers, Sheremet, McKee Smith, Tolman, van Vledder, Wolf and Young2007) predict the spatial and temporal evolution of the directional spectrum solving the spectral energy equation.
Numerical models as well as data assimilation schemes have gone through various improvements in recent decades causing variability in the long-term data not related to natural variability. The reanalysis aim is to overcome the inhomogeneities related to the model (Sterl and Caires, Reference Sterl and Caires2005): the best available model is used to repeat the analysis procedure throughout the period in consideration. Inhomogeneities related to the changes in data coverage remain.
Two numerical databases, namely HIPOCAS (Hindcast of Dynamic Processes of the Ocean and Coastal Areas of Europe) and ERA40 (ECMWF Re-analysis-40), and the International Comprehensive Ocean-Atmosphere Data Set (ICOADS) archive containing the voluntary observations have been used and compared in this work.
The first hindcast database derives from the EU project HIPOCAS (Guedes Soares, Reference Guedes Soares2008) aimed at producing a 44 year database of wave, wind and sea level for all seas around Europe with a 2° × 2° resolution. In addition a data set from the whole North Atlantic was also produced in order to determine the swell conditions that would be arriving on the Atlantic coasts of Europe (Pilar et al., Reference Pilar, Guedes Soares and Carretero2008). The initial wind field was from the NCEP reanalysis (Kalnay et al., Reference Kalnay, Kanamitsu, Kistler, Collins, Deaven, Gandin, Iredell, Saha, White, Woollen, Zhu, Chelliah, Ebisuzaki, Higgins, Janowiak, Mo, Ropelewski, Wang, Leetmaa, Reynolds, Jenne and Joseph1996), which was used to force a regional atmospheric model REMO (Jacob and Podzun, Reference Jacob and Podzun1997; von Storch et al., Reference von Storch, Langenberg and Feser2000), which finally forced the WAM wave model (WAMDI, Reference Hasselmann, Hasselmann, Bauer, Janssen, Komen, Bertotti, Lionello, Guillaume, Cardone, Greenwood, Reistad, Zambresky and Ewing1988). The data was validated with buoy measurements establishing its general adequacy (Pilar et al., Reference Pilar, Guedes Soares and Carretero2008). A comparison with the database from the National Oceanic and Atmospheric Administration (NOAA) showed the general consistency between these two databases (Campos and Guedes Soares, Reference Campos, Guedes Soares, Guedes Soares, Garbatov, Sutulo and Santos2012).
The WAM model is also the base of the European Centre for Medium-Range Weather Forecast (ECMWF) atmospheric ERA-40 reanalysis (Uppala et al., Reference Uppala, Kållberg, Simmons, Andrae, da Costa Bechtold, Fiorino, Gibson, Haseler, Hernandez, Kelly, Li, Onogi, Saarinen, Sokka, Allan, Andersson, Arpe, Balmaseda, Beljaars, van de Berg, Bidlot, Bormann, Caires, Chevallier, Dethof, Dragosavac, Fisher, Fuentes, Hagemann, Hólm, Hoskins, Isaksen, Janssen, Jenne, McNally, Mahfouf, Morcrette, Rayner, Saunders, Simon, Sterl, Trenberth, Untch, Vasiljevic, Viterbo and Woollen2005), which runs from September 1957 to August 2002. The dataset consists of six-hourly global fields with a 1·5° × 1·5° grid resolution. ERA-40 was produced using the ECMWF Integrated Forecasting System, a two-way coupled atmosphere–wave model system and a three-dimensional Variational Data Assimilation (3DVAR) scheme. Observations of ocean wind speeds from VOS, buoys, satellite scatterometer and (when available) satellite altimeter SWHs were assimilated in the analysis process.
The ICOADS is a database resulting from the Voluntary Observing Ships’ (VOS) scheme (Kent et al., Reference Kent, Taylor, Truscott and Hopkins1993), which is an international program by which ships are recruited by National Meteorological Services for taking and transmitting meteorological observations. The ICOADS database has records of visually observed height, period and direction for wind-waves and swell allowing the evaluation of ocean wave climate (Gulev et al., Reference Gulev, Grigorieva, Sterl and Woolf2003). However, only 40–60% of the total number of reports include wave information (Gulev and Hasse, Reference Gulev and Hasse1998). This type of data has some uncertainties; Guedes Soares (Reference Guedes Soares1986a; Reference Guedes Soares1986b) and Gulev et al. (Reference Gulev, Cotton and Sterl1998) among others have concluded that wave height is always overestimated in VOS observations in comparison with remote sensed data and numerical models in the areas with small SWH. In this work, thanks to the fact that the analysis is specifically conducted along the ship routes and the behaviour of ships and of masters in facing heavy seas is taken into account, a new interpretation of this bias will be given.
VOS data have the irreplaceable value to be the only source of meteorological data along the ship paths and the ones that have the rough weather avoidance built into them. Understanding how this influences the statistics from the VOS database is fundamental.
4. WAVE CLIMATE ALONG THE ROUTES
In order to analyse the wave climate along the routes, it is firstly necessary to select the wave data associated with each route from the databases. The three databases considered are not homogenous in resolution. In particular the VOS data by their nature do not have a resolution, but are unevenly distributed depending on the ships’ paths. To bypass this problem, the data are gridded in a squared area of 2° × 2° around a number of previously defined points (see Vettor and Guedes Soares, Reference Vettor and Guedes Soares2015).
Figure 2 shows the path and all the coordinates selected from the different databases for the northern route from the Channel to Puerto Rico (Ch_PR-N). Even if the data are not exacly located at the same points, they are close to each other, reflecting the climate of the same area. Thus the derived statistics can be considered consistent.
4.1. Significant Wave Height
A straightforward indication of the severity of the seastate along the routes is given in Figure 3 where the evolution of the average significant wave height by dependency on the longitude are provided for the extreme meteorological seasons: summer (June-July-August) and winter (December-January-February).
ERA40 and HIPOCAS results are comparable both in tendency and amplitude. In the westernmost part of the North Atlantic Ocean (latitudes smaller than 50°N) HIPOCAS tends to heavily underestimate the SWH in winter, especially in the northern part (routes ending in the Virginia area: Ch_VA, St_VA-N and St_VA-S). This can be explained considering the limitations to 64°N latitude in the domain used to produce the hindcast database. This might have caused the loss of some information relative to the storms generated closer to the US coasts. VOS averages show a strong bias all along the routes in summer and more in general when the SWH is lower. This issue will be discussed in Section 5.
Besides the quantitative differences between the databases, some considerations are evident. In winter ships are asked to cross a critical area between 35°W and 25°W longitude with an average SWH that can reach 5 m, which is in fact near to the North Atlantic SWH climatological maximum, located south-west of Iceland. In contrast, in the area close to the US coast up to 50°W longitude, the sea conditions are calmer, around 3·5 m. In summer the weather is more homogeneous and no noteworthy differences were noted in the evolution with a SWH varying between 1·2 m and 2·2 m depending on the source data.
4.2. Scatter diagrams
While the evolution of the SWH along the route gives a clear qualitative view to estimate how heavy the navigation can be, when ship motion has to be analysed more in detail, the knowledge of the wave periods is essential as it greatly influences the amplitude of ship responses. Moreover, when one is dealing with routes and the life-time ship loads have to be estimated, a long-term study is required. These types of studies are perfomed as an extension of the short term formulations and take into account the occurrences of the various seastates. As an input a scatter diagram has then to be considered.
Traditionally two sources are adopted depending on the scope. When one is verifying the operability of the ship, the global wave statistics database (Hogben, et al., Reference Hogben, Da Cunha and Olliver1986) is the most common reference; it divides the world into large zones allowing consideration of the long term statistics of the areas where the ship will operate. In the design phase a conservative approach is recommended and a standard scatter diagram is adopted (IACS, 2001; Figure 4a), computed considering only the data (modified according to Bitner-Gregersen et al., Reference Bitner-Gregersen, Cramer and Korbijin1995) of the zones relative to the North Atlantic, which is widely considered as a reference for highly demanding weather conditions.
They are both descriptions of wide areas, which can include zones in which the ship will never operate. If in a design phase it can be accepted both for safety reasons and to allow a possible conversion of the ship, when the operability analysis and schedule management are the goals, more specific input data would improve the quality and reliability of the results. The scatter diagrams proposed here are computed from the ICOADS database only taking into account the data in a narrow belt near the path of the main routes; they are thus more suitable for long-term operational analysis on those specific routes. As a drawback, they need to be re-computed if another route needs to be investigated.
Figure 4 compares the standard scatter diagram with the ones relative to the Ch_VA route, the northernmost and closer to the climatologial maximum, and the St_MIA route, which is the southermost one and shows generally lower waves and higher probability of swell, due to its distance from the area of generation of the extra-tropical storms (Vettor et al., Reference Vettor, Semedo, Guedes Soares, Breivik, Sterl, Reistad, Guedes Soares and Lopez Pena2014).
4.3. Relative wave heading
Another parameter that has great influence on the ship motions and loads is the relative wave heading (Guedes Soares, Reference Guedes Soares1995). The influence of this parameter is often either neglected considering only the head wave condition (usually, but not always, the most demanding one), or assuming a uniform probability for all the different headings.
In the North Atlantic ocean, waves tend to assume a west-east direction, generally more north-wise above the 40th parallel and south-wise below. On the other hand, the routes can be either eastward when starting from US or westward when starting from Europe. This combination allows assessment of the relative wave heading probability: at each waypoint the main ship course (towards the destination) is compared with the hindcast wave direction (from HIPOCAS database) for the whole period generating an array of relative wave headings from which the histogram is eventually computed.
The relative wave heading probabilities for westward passages for all the routes is shown in Figure 5 (eastward probabilities are obviously symmetrical around 90°). In general it can be noted that due to the North Atlantic wave climate and the route direction, ships directed to the Virginia area mostly encounter quartering-head waves (relative heading angle θ ⩾ 135°).
5. BIAS IN VISUAL OBSERVATIONS
Figure 3 shows a significant bias in the average SWH computed using the VOS database when compared with the ones derived from numerical models. This bias is clearly evident in the summer and generally when the average SWH is low. Figure 6 shows the relationship between the seasonal mean SWH (from ERA40) and the trajectory of the considered routes.
The presence of a seasonal bias in VOS data has been already noted in Vettor and Guedes Soares (Reference Vettor and Guedes Soares2012) where just one point of coordinates (50°N, 23°W) along the route between the Channel and Virginia (Ch_VA), was analysed. The present study provides evidence that the bias has two main causes: the visual error and the difference due to storm avoidance.
The reliability of the ICOADS database has been discussed in the literature (e.g. Gulev et al., Reference Gulev, Grigorieva, Selemenov and Zolina2003) and studies of the errors in the statistics derived from these data can be found in Guedes Soares (Reference Guedes Soares1986a; Reference Guedes Soares1986b) and Gulev et al. (Reference Gulev and Hasse1998) among others. Generally an overestimation in areas with smaller SWH was identified and some correction has been proposed. The data used in this work have been preprocessed accordingly to correct for those differences (see Vettor and Guedes Soares, Reference Vettor and Guedes Soares2012). Nevertheless the correction only affects the very small waves (lower than 0·5 m) and a bias still persists as can clearly be seen in Figure 3, especially in the summer cases, and requires further investigation.
The solution of the problem would be to understand how this bias works, keeping in mind how the source of the observations can influence the data. A possibility could be to adapt a nonparametric regression method, already applied to wave data (Caires and Sterl, Reference Caires and Sterl2005), in order to correct the ICOADS database, using only the more reliable reports. However, since a correction of the database is out of the scope of this paper, a simple method based on logical assumptions will be adopted here.
The main characteristic of the VOS dataset is to be the only source of meteorological data along the ships’ paths, which reflects the weather condition in which the ships operate. For this reason, even if the single observation were precise and perfectly reliable, how the masters act when a storm has to be faced would affect the statistics derived from those observations.
In summer and in low SWH situations there is no storm avoidance as there is no storm. This means that statistics can be more significant and reliable, depending on the number of observations, and other factors, but they are not affected by the weather. Instead the differences between the visual results and the numerical databases are due to the visual observation process, which shows as an almost constant value across all longitudes. Determining the average value of this systematic difference allows correction of the data both in summer and in winter, when a separation of visual data and the effect of rough weather would otherwise be hardly possible.
During winter, in addition to the observational bias the storm avoidance bias is also included. After correcting the winter data with the average observational bias determined from the summer, the only bias left is the one related to the rough weather avoidance which can then be determined.
The observational bias computed considering only the summer data is shown in Table 1. It can be seen that it does not significantly differ among different routes, especially when the percentage value is considered. Thus the average percentage value of 20·83% can be considered representative for all the routes and used to correct the observational component of the bias.
After the correction, the evolution of the mean SWHs along the routes can be seen in Figure 7. Once the observational bias due to errors in observation is corrected, the resulting statistics are assumed to be ruled by the behaviour of ships in waves and the decisions of the masters regarding whether to face or not the forecasted weather conditions. Thus, after this correction, the strong winter bias present in all the routes in the longitude range [50°W, 10°W] represents the effect of storm avoidance on wave statistics from VOS reports, while during the summer there are no significant differences among all three databases. The reduction of the mean SWH from 0·8 m for the southernmost routes (St-MIA and St-VA-S) up to 1·3 m for the northernmost ones (Ch-PR-N and Ch-VA) is shown.
6. STORM AVOIDANCE CORRECTION
Looking at Figure 7, some results are clear. In general the bias increases going from west to east and in some cases this trend seems more regular and almost linear. As can be seen in Figure 6b and as the scatters of the bias versus the longitude (not shown here) demonstrated, this behaviour is accidentally due to the fact that the climatological maximum of the SWH in winter is located in the north east part of the North Atlantic basin, thus the increase of the bias with the longitude actually corresponds to the increase in the SWH.
More interesting is to determine how storm avoidance is correlated with the intensity of the weather. In Figure 8 the scattering between the winter storm avoidance correction (mean of numerical minus corrected observed SWH) is shown. In the first four routes the correlation is clear and the bias increases linearly with the SWH. In areas characterised by sea states with an average SWH of 4 m, ships experience conditions with SWH reduced of about 1 m.
In the routes between the Strait of Gibraltar and the north-east side of the USA (St_VA-N and St_VA-S) the scatter is wider, especially in the second/southern one. This is due to the westernmost part of the routes close to the Virginia area (highlighted in the red circles in Figure 8 and also evident in the scatter of the SWH with the longitude, not shown here) where the effect of storm avoidance is low (between 0 and 20 cm) even in quite demanding sea-states (between 3·0 to 3·5 m SWH). It can be noticed that even the route between the Channel and the Virginia (Ch_VA) area is affected by this circumstance which contributes to make the slope of the regression higher.
To find an explanation it is necessary to take into account the mean peak periods. In winter, due to the short fetch, the wave climate in the northwest part of the North Atlantic basin presents a local maximum of the wind-sea waves (lower peak periods), whereas the swell waves' amplitude is in line with the global climate of the basin (see Vettor et al., Reference Vettor, Semedo, Guedes Soares, Breivik, Sterl, Reistad, Guedes Soares and Lopez Pena2014). This means that this area is characterised by local storms, which are more difficult to predict for the master than the global climate and the swell waves, as demonstrated by the high rate of weather related accidents offshore the US Atlantic coast (Guedes Soares et al., Reference Guedes Soares, Bitner-Gregersen and Antão2001), especially for ships sailing westward, which leave the European coast five to seven days before reaching the destination and have to plan the voyage using medium-long term weather forecast.
These considerations suggest also including the wave period in the analysis. The mean SWH, although being an important factor (with a global correlation coefficient equal to 0·88), is not fully sufficient to understand the effect of weather routing and storm avoidance on the weather conditions experienced by the ships. The confirmation is given in Figure 9 where the correlation between the storm avoidance correction and the mean peak period is shown. The red-circled outliers of Figure 8 do not appear in this case, suggesting that they were indeed caused by the wave steepness (wave height wave period ratio) and how the ships face those conditions. All in all the correlation coefficient between the storm avoidance correction and the peak period is equal to 0·82.
The correlation is strong and adequately represented by a linear regression in all cases. Nevertheless, as shown in Figure 10, the gradient as well as the magnitude of the regression differs from one route to the other. The former is similar for the three routes including respectively the Channel and the Strait of Gibraltar, while the latter is step by step higher as the route moves further north.
From the previous considerations it appears fundamental to take into account both parameters as in Figure 11. In this case a bilinear regression (Equation (1)) seems to be suitable and to keep in consideration the two most important variables ruling storm avoidance. As a result one can extend the meaning of the equation considering that, if in a region or along a route the mean wave climate is known in terms of SWH and peak period, thanks to the storm avoidance effect the average SWH that a ship will encounter can be estimated from:
The validity of the equation is limited by the boundaries of the input values [2 m, 5 m] for the SWH and [8 s, 13 s] for the peak period. This is fundamental mainly for the peak period since changing the encounter frequencies the ship dynamic dramatically changes, while a certain extrapolation might be reasonable for higher SWHs. Table 2 shows the mean SWH that a ship, which follows the common safety recommendations, can expect to encounter while sailing a route characterised by the sea-state described by the given SWHs and peak periods. Of course, being a general approach it does not take into account the characteristics of the vessel, but it can probably be a good term of reference for trans-oceanic ships, especially if operating along the same route for a considerable period, as might often happen to a container-ship.
7. CONCLUSIONS
An analysis of the wave climate along the main North Atlantic routes has been presented considering three different sources: two databases generated by numerical wave models (HIPOCAS and ERA40) and weather reports from the voluntary observations. A description of the mean seasonal significant wave height as well as the joint probability of significant wave height and wave period was provided. It showed a mean SWH between 1 m and 2 m in summer (lower limit for the southwest portions of the routes and higher limit for the northeast ones) and between 2 m (approaching to Puerto Rico) and 5 m (close to the Channel) in winter with a general tendency of increasing from west to east, except for the routes operating in the northwest-southeast direction. Considering its influence in ship response, the relative wave heading was also computed. For the case of westward passages, in the routes directed to the Virginia area ships mostly experience head-quartering seas, whereas in the route from the Strait of Gibraltar to Miami/Gulf of Mexico following-quartering sea is predominant.
The evolution of the SWH along the routes highlighted a strong bias on the visual observation in summer and more in general in calmer areas. A calibration of the data taking into account the tendency of ships to avoid heavy weather and its consequences on statistics from the VOS dataset led to a new interpretation of the voluntary meteorological reports, which have been considered to overestimate the SWH by about 20% of the actual value. This correction allowed the identification of the effect of storm avoidance in the waves that ships experience along a route. The impact of crossing the area of generation of the extra-tropical storms, where the mean SWH is around 4·5 m in winter, is clear in the routes starting from/directed to the Channel where the mean SWH is up to 1·3 m less than the one computed by numerical codes.
The concept of storm avoidance correction, indicating the difference between the actual mean SWH and the mean SWH that the ships experience has been introduced. Analysing the correlations between this bias and the commonly used sea-states parameters (SWH and peak period), a general equation to estimate the storm avoidance correction has been proposed.
ACKNOWLEDGEMENTS
This work was performed within the project SHOPERA - Energy Efficient Safe SHip OPERAtion, which was partially funded by the EU under contract 605221. The first author was supported by the Portuguese Foundation for Science and Technology (FCT - Fundação para a Ciência e Tecnologia, Portugal) under the contract no. SFRH/BD/89476/2012.