A Mean Dynamic Topography of the Mediterranean Sea computed from altimetric data, in-situ measurements and a general circulation model

In the Mediterranean Sea, where the mean circulation is largely unknown and characterized by smaller scales and less intensity than in the open ocean, the interpretation of altimetric Sea Level Anomalies (SLA) is rather difficult. In the context of operational systems such as MFS (Mediterranean Forecasting System) or MERCATOR, that assimilate the altimetric information, the estimation of a realistic Mean Dynamic Topography (MDT) consistent with altimetric SLA to be used to reconstruct absolute sea level is a crucial issue. A method is developed here to estimate the required MDT combining oceanic observations as altimetric and in-situ measurements and outputs from an ocean general circulation model (OGCM).In a first step, the average over the 1993–1999 period of dynamic topography outputs from MFS OGCM provides a first guess for the computation of the MDT. Then, in a second step, drifting buoy velocities and altimetric data are combined using a synthetic method to obtain local estimates of the mean geostrophic circulation which are then used to improve the first guess through an inverse technique and map the MDT field (hereafter the Synthetic Mean Dynamic Topography or SMDT) on a 1/8° resolution grid.Many interesting current patterns and cyclonic/anticyclonic structures are visible on the SMDT obtained. The main Mediterranean coastal currents are well marked (as the Algerian Current or the Liguro–Provenço–Catalan Current). East of the Sicily channel, the Atlantic Ionian Stream divides into several main branches crossing the Ionian Sea at various latitudes before joining at 19°E into a unique Mid-Mediterranean Jet. Also, strong signatures of the main Mediterranean eddies are obtained (as for instance the Alboran gyre, the Pelops, Ierapetra, Mersa-Matruh or Shikmona anticyclones and the Cretan, Rhodes or West Cyprius cyclones). Independent in-situ measurements from Sea Campaigns NORBAL in the North Balearic Sea and the North Tyrrhenian Sea and SYMPLEX in the Sicily channel are used to validate locally the SMDT: deduced absolute altimetric dynamic topography compares well with in-situ observations. Finally, the SMDT is used to compute absolute altimetric maps in the Alboran Sea and the Algerian Current. The use of absolute altimetric signal allows to accurately follow the formation and propagation of cyclonic and anticyclonic eddies in both areas.     

Variability in hydrographic conditions to the east and northwest of South Georgia, 1996–2001

Six years of high-resolution hydrographic data from the eastern and northwestern sides of South Georgia (southwest Atlantic) are used to study the changing circulation and water mass properties of the region. One year of data from these locations was used previously to describe the oceanographic conditions at those times; using the much greater volume of data now available, we identify which features appear temporally robust and which are transient, and begin addressing topics relating to the forcing of the inter-annual variability and the potential consequences for the local ecosystem. Waters on the shelf and those over the adjacent deep ocean invariably have different hydrographic properties, though the transition between them can be abrupt or gradual. The onshelf/offshelf differences vary greatly from year to year, due to the combined influences of local and remote processes. There are several instances of strong physical coupling between the eastern and northwestern sides of South Georgia; this offers potential for distinguishing physically-induced ecosystems changes separately from biologically-induced ecosystems changes. On the northeast side of the shelf, close to Cumberland Bay, there is evidence of an often intense, but variable, cyclonic circulation that is the result of interaction with the local bathymetry. This may act as a retention mechanism, and enhance local productivity. Two examples of extreme cold anomalies are present in the series of measurements. One of these (in 2000/2001) affected a limited area at the eastern side of the region surveyed and was due to an intrusion of the Southern Antarctic Circumpolar Current Front. The other (in early 1998) was due to the combined effects of the passage of a large-scale ocean anomaly that had its origins upstream in the Pacific Ocean, and strong air/sea interaction. Both of these were associated with the strong 1997/1998 El Niño event. Whilst previous studies have observed a link between El Niño forcing and ocean response around South Georgia with a temporal lag of around 3 years, we observe a much more rapid response to the extremely strong 1997/1998 El Niño event. This indicates that the ocean and ecosystem around South Georgia are more immediately susceptible to extreme instances of remote climatic forcing than had been supposed.     

Identification of potential sites for deep-ocean waste isolation with a geographic site-selection model

Identification of optimal sites for the isolation of waste on the abyssal seafloor was performed with two approaches: by the traditional method of map overlays of relevant attributes, and by a specially developed, automated Site-Selection Model (SSM). Five initial, Surrogate Sites, identified with the map-overlay approach, were then compared with the more rigorously produced scores from the SSM. The SSM, a process for optimization of site locations, accepts subjective, expert-based judgments and transforms them into a quantitative, reproducible, and documented product. The SSM is adaptable to any siting scenario. Forty-one factors relevant to the isolation scenario, including 21 weightable factors having a total of 123 scorable categories, have been entered into the SSM. Factors are grouped under project definition, unique environments, anthropogenic, geologic, biologic, weather, oceanographic and distance criteria. The factor scores are linked to a georeferenced database array of all factors, corresponding to 1°×1° latitude–longitude squares. The SSM includes a total of 2241 one-degree squares within 1000 n.m. of the U.S. coasts, including the western North Atlantic, the Gulf of Mexico, and the eastern North Pacific. Under a carefully weighted and scored scenario of isolation, the most favorable sites identified with the SSM are on the Hatteras and Nares Abyssal Plains in the Atlantic. High-scoring sites are also located in the Pacific abyssal hills province between the Murray and Molokai Fracture Zones. Acceptable 1° squares in the Gulf of Mexico are few and of lower quality, with the optimum location on the northern Sigsbee Abyssal Plain. Two of the five Surrogate Site locations, on the Hatteras and Sigsbee Abyssal Plains, correspond to the best SSM sites in each ocean area. Two Pacific and a second Atlantic Surrogate Site are located in low-scoring regions or excluded by the SSM. Site-selection results from the SSM, although robust, are an initial attempt to quantify the site-selection process. The SSM database exposes a significant lack of high-quality information for many areally mappable attributes on the abyssal seafloor, particularly bottom-current speed and measures of biologic productivity and flux. Terminologies and classifications of some measures, such as sediment types, suffer from parochialism and vary by ocean. Considerable research is needed even for a broad understanding of the environmental measures required to make sound societal decisions about use of the abyssal seafloor for disposal or other purposes.     

The West African coastal upwelling filaments and cross-frontal water exchange conditioned by them

One of the important problems in the oceanography of the wind-driven upwelling regions of the Ocean is the investigation of water exchange processes in the coastal zone. Satellite data (thermal and colour imagery) have changed our view on these processes after the relatively recent discovery of cold, chlorophyll-rich, narrow (< 50 km wide) offshore flowing filaments off the west coasts of North America, North and South Africa. On the basis of satellite IR images and oceanographic original and archive data we investigated systems of filaments in the northwest and southwest African upwelling regions. The spatial distribution of filaments was analyzed. It was found that seasonal variability of the filaments' location depends upon the general intensity of upwelling motion along the coast during the year. The main statistical characteristics of filaments were obtained. An example of the three-dimensional velocity structure of the filament was presented. In order to estimate the intensity of cross-frontal exchange process due to filaments, a special procedure was proposed. The values of the velocity of cross-frontal water exchange produced by filaments in these upwelling zones are obtained. It was shown that upwelling filaments play an important role in water exchange between the coastal zone and the open ocean. A non-dimensional parameter which characterizes the permeability of oceanic fronts for water exchange due to mesoscale dynamical structures was introduced and estimated.     

High-resolution synthetic monitoring by a 4-dimensional variational data assimilation system in the northwestern North Pacific

A data assimilation system is applied to the integrated monitoring of oceanic states in the northwestern North Pacific by combining a high resolution ocean general circulation model with an adjoint method. A comparison of assimilation results with observations shows that the system is better able to represent synoptic features of ocean circulation than do models or data alone. Furthermore, meso-scale features associated with frontal structures and eddies, which are often seen in the Kuroshio and Oyashio extension regions and the Sea of Japan, are better defined in the assimilation results. These features suggest that our 4D-VAR high-resolution data assimilation system is capable of providing time series data which satisfy the model physics and fit the observations, and hence the ocean state derived from our system has greater information content than that obtained from earlier methods.     

Containerization, inter-port competition, and port selection

This paper explores the criteria shippers employ in the port selection process. By focusing on the containerized traffic between the North American Mid-West and Western Europe, the factors considered by exporters and freight forwarders are examined. The findings suggest that decision-makers are influenced more by price and service considerations of land and ocean carriers than by perceived differences in the ports of entry and exit. Port infrastructures do not appear to play an important role in the routeing decisions made by an important group of independent businesses involved in the North Atlantic container trade.     

Mesoscale subduction at the Almeria–Oran front: Part 1: Ageostrophic flow

This paper presents a detailed diagnostic analysis of hydrographic and current meter data from three, rapidly repeated, fine-scale surveys of the Almeria–Oran front. Instability of the frontal boundary, between surface waters of Atlantic and Mediterranean origin, is shown to provide a mechanism for significant heat transfer from the surface layers to the deep ocean in winter. The data were collected during the second observational phase of the EU funded OMEGA project on RRS Discovery cruise 224 during December 1996. High resolution hydrographic measurements using the towed undulating CTD vehicle, SeaSoar, traced the subduction of Mediterranean Surface Water across the Almeria–Oran front. This subduction is shown to result from a significant baroclinic component to the instability of the frontal jet. The Q-vector formulation of the omega equation is combined with a scale analysis to quantitatively diagnose vertical transport resulting from mesoscale ageostrophic circulation. The analyses are presented and discussed in the presence of satellite and airborne remotely sensed data; which provide the basis for a thorough and novel approach to the determination of observational error.     

Interactions of wind and density driven currents in North Sea ROFIs—a model study

Three different versions of a baroclinic three-dimensional circulation model of the North Sea are used to obtain information on the wind and density interactions in the North Sea ROFIs (Regions Of Freshwater Influence): the standard version with fully prognostic treatment of salinity and temperature is compared to a barotropic model run on the same grid on the one hand and to an also fully prognostic model run on a four times coarser grid on the other hand. In order to gain knowledge on the wind and density interactions, two opposing wind directions are chosen for investigation, namely a time of strong north wind, 21st–28th April 1982, and a time of strong southwest wind, 22nd–24th May 1982. In the April case the effect of the salinity gradients on the border of the ROFIs of Rhine, Weser, Ems and Elbe, i.e. along the continental shore, is shown to lead to a clear enhancement of the mean surface currents. In May this result is partly disguised by the additional effect of the thermocline in the deeper parts of the North Sea, i.e. in the classical shelf sea regime region. Nevertheless, the same pattern of enhanced mean surface currents along the coast is detected and is of the same order of magnitude as in the April case. It is thus concluded that although the circulation in the North Sea is reversed by the wind, the density induced component of the general circulation is modified only slightly.     

Validation of a hybrid optimal interpolation and Kalman filter scheme for sea surface temperature assimilation

A hybrid data assimilation scheme designed for operational assimilation of satellite sea surface temperatures (SST) into an ocean model has been developed and validated against in-situ observations. The scheme consists of an optimal interpolation (OI) part and a greatly simplified Kalman filter (KF) part.The OI is performed only in the longitudinal and latitudinal directions. A climatological field is used as a background field for the interpolation. It is constructed by fitting daily averages of satellite SST to the annual mean, annual, and semiannual harmonics in a 20 km by 20 km grid. The background error covariance is approximated by a spatially varying two-dimensional exponential covariance model. The parameters of the covariance model are fitted to the deviations of the satellite data from the background field using data from a full year.The simplified KF uses ocean model forecasts as a background field. It is based on the assumption that it is possible to neglect horizontal SST covariances in the filter and that the typical time scale for vertical mixing in the mixed layer is much shorter than the average time between observations. We therefore assume that the error variance in a column of water is evenly spread out throughout the mixed layer. The result of these simplifications is a computationally very efficient KF.A one year validation of the scheme is performed for year 2001 using an operational eddy resolving ocean model covering the North Sea and the Baltic Sea. It is found that assimilation of sea surface temperature data reduces the model root mean square error from 1.13 °C to 0.70 °C. The hybrid scheme is found to reduce the root mean square error slightly more than the simplified KF without OI to 0.66 °C. The inclusion of spatially varying satellite error variances does not improve the performance of the scheme significantly.     

CO2 partial pressure in Northeast Atlantic and adjacent shelf waters: Processes and seasonal variability

CO₂ partial pressure in surface water was measured in the Northeast Atlantic and in the Hebride Shelf/North Sea area during a cruise with R.V. Poseidon in June 1991. A mean p_CO2 of 303 μatm was found in the Atlantic between 50°N and 60°N. For an atmospheric CO₂ content of 357.5 ppm(v) this corresponds to a partial pressure difference of −55 μatm. This supports the view that the subarctic Atlantic is a significant sink within the CO₂ cycle between the ocean and the atmosphere. A comparison of our measurements with other data reveals that the p_CO2 distribution changes significantly during May/June. This explained by seasonal warming, CO₂ exchange with the atmosphere and biomass production. The contribution by each of these processes to the seasonal variations is calculated. It was found that during a plankton bloom the production of biomass is the dominating factor and may lower seawater p_CO2 by almost 100 μatm. The shelf areas are charactrized by strong p_CO2 gradients which are explained by water exchange with the Atlantic, temperature effects and biomass production.     

Rates and mechanisms of climate change: Implications for ports and harbours

The paper explores some of the key ocean-atmosphere processes relevant to the themes of ‘global warming“ and ‘sea level rise“ addressing some of issues that form the basis of the science. It is shown that both processes as well as the data used to characterise such processes demonstrate that a) there is no such thing as a single global sea level and b) changes in extreme storminess are not uniform worldwide. Instead, trends in storminess (including Asian typoon and Atlantic hurricane activity) exhibit marked regional variability with strong regional signature linked to episodes of major El Niño and La Niña activity. Similarly, satellite data demonstrates that whereas some ocean areas have recently experienced sea level rise caused principally by thermal expansion effects, other ocean areas have been subject to cooling and sea level lowering. This account provides an introductory account that discusses the rationale for understanding the regional variability in the nature of coastal flood risk to the world’s ports and harbours.     

Oceanic CO2 uptake and biogeochemical variability during the formation of the Eastern North Atlantic Central water under two contrasting NAO scenarios

Long-term variability of the biogeochemical properties during the formation of central waters in the Eastern North Atlantic were analyzed between 42–47°N and 10–20°W from the dataset gathered during the Galicia VII (GVII) and C. Darwin 58/59 (CD58/59) cruises. These cruises that showed important changes in the thermohaline properties and the nutrient abundance of the upper layers were carried out under contrasting conditions of the North Atlantic Oscillation (NAO) index. The different climate forcing led a meridional shift of the transition zone between the formation regions of subpolar and subtropical Eastern North Atlantic Central Water (ENACWp and ENACWt, respectively). This displacement conditioned the presence of each ENACW in the study region and so the thermohaline and biogeochemical properties. The effect of the observed variability at decadal scale on the air–sea CO₂ gradient (ΔfCO₂) and exchange (F_CO2) was analyzed using 1D model approach throughout 11 weekly-steps that simulated the development of a spring bloom during the shoaling of the mixed layer. The outputs of the model showed an intensification of the ocean CO₂ uptake due to higher biological CO₂ drawdown, during positive NAO conditions and its weakening under negative NAO influence.     

Evolving the subspace of the three-dimensional multiscale ocean variability: Massachusetts Bay

A data and dynamics driven approach to estimate, decompose, organize and analyze the evolving three-dimensional variability of ocean fields is outlined. Variability refers here to the statistics of the differences between ocean states and a reference state. In general, these statistics evolve in time and space. For a first endeavor, the variability subspace defined by the dominant eigendecomposition of a normalized form of the variability covariance is evolved. A multiscale methodology for its initialization and forecast is outlined. It combines data and primitive equation dynamics within a Monte-Carlo approach.The methodology is applied to part of a multidisciplinary experiment that occurred in Massachusetts Bay in late summer and early fall of 1998. For a 4-day time period, the three-dimensional and multivariate properties of the variability standard deviations and dominant eigenvectors are studied. Two variability patterns are discussed in detail. One relates to a displacement of the Gulf of Maine coastal current offshore from Cape Ann, with the creation of adjacent mesoscale recirculation cells. The other relates to a Bay-wide coastal upwelling mode from Barnstable Harbor to Gloucester in response to strong southerly winds. Snapshots and tendencies of physical fields and trajectories of simulated Lagrangian drifters are employed to diagnose and illustrate the use of the dominant variability covariance. The variability subspace is shown to guide the dynamical analysis of the physical fields. For the stratified conditions, it is found that strong wind events can alter the structures of the buoyancy flow and that circulation features are more variable than previously described, on multiple scales. In several locations, the factors estimated to be important include some or all of the atmospheric and surface pressure forcings, and associated Ekman transports and downwelling/upwelling processes, the Coriolis force, the pressure force, inertia and mixing.     

Multivariate assimilation in MERCATOR project: New statistical parameters from forecast error estimation

The new operational prototype of Mercator (french Global Ocean Data Assimilation Experiment contribution) is composed of a North Atlantic primitive equation ocean model OPA (Ocean Parallel Algorithm between 20°S and 70°N, [Madec, G., P. Delecluse, M. Imbard and C. Lévy (1998). OPA8.1 ocean general circulation model reference manuel. Notes du pôle de modélisation IPSL. n°11: 91p]) and of a multivariate and multidata assimilation scheme [De Mey, P. and M. Benkiran (2002). “A multivariate reduced-order optimal interpolation method and its application in Mediterranean basin-scale circulation.” Ocean Forecasting : Conceptual basis and application, Pinardi, N., Springer Verlag.] This system has already given some significant improvements from previous Mercator configurations (M. Benkiran, personal communication). However some biases on ocean state still remain in the tropics where the reduced-order optimal interpolation scheme is suspected to be ill-parameted in the model forecast error. Indeed the guess error covariance matrix is decomposed into an error variance value and a spatio-temporal correlation function which are assumed to have some “good” properties (spatial homogeneity of the correlation function, constant ratio between signal and error variance). This study shows how we can use ensemble methods to validate these assumptions. We can see that the correlation function can reach negative values locally, mostly in regions of high variability contradictory with the homogeneous hypothesis. The reduced space used in the operational configuration is based on the signal seasonal Empirical Orthogonal Functions (EOFs). An empirical relationship between signal and error variance has been set and the correlation function is the same on every dimension of the reduced space. By projection of the estimated guess error variance onto the reduced space, we find a repartition of this quantity quite different to what was set in the system. The error statistics is found to be inhomogeneous compared to hypothesis made in the assimilation scheme. These two new parameters tested separately in the assimilation scheme gives significant improvements of the forecast and analysis results. This is particularly obvious in the tropics. But relationship between signal and error statistics (as assumed in the optimal interpolation) is found to be complex.     

Container shipping in 1990

By 1990 express shuttle container ships will dominate the USA/North European and USA/Pacific trades. The pattern is already apparent in the Pacific. These will touch only a few ports such as Halifax, New York and Virginia on the Atlantic and Los Angeles/Long Beach, Oakland and Seattle/Tacoma in the West. Their low costs and good services will capture over 80% of available general cargo. Round the world systems will be left with lucrative fringe business in the Mediterranean, US Gulf/South Atlantic, Singapore, India and Arab countries.

Progress will be sensitive to twelve trends. Most important of these will be labour resistance and management concerns about rapid obsolescence of ships and terminals. Rail and ship costs have been rationalized and are not as flexible as terminal costs. The ultimate terminal will consume about one third of the time and money presently needed.

This paper is chiefly intended to suggest a forecasting methodology.     

On the key role of nutrient data to constrain a coupled physical–biogeochemical assimilative model of the North Atlantic Ocean

A sequential assimilative system has been implemented into a coupled physical–biogeochemical model (CPBM) of the North Atlantic basin at eddy-permitting resolution (1/4°), with the long-term goal of estimating the basin scale patterns of the oceanic primary production and their seasonal variability. The assimilation system, which is based on the SEEK filter [Brasseur, P., Verron, J., 2006. The SEEK filter method for data assimilation in oceanography: a synthesis. Ocean Dynamics. doi: 10.1007/s10236-006-0080-3], has been adapted to this CPBM in order to control the physical and biogeochemical components of the coupled model separately or in combination. The assimilated data are the satellite Topex/Poseidon and ERS altimetric data, the AVHRR Sea Surface Temperature observations, and the Levitus climatology for salinity, temperature and nitrate.In the present study, different assimilation experiments are conducted to assess the relative usefulness of the assimilated data to improve the representation of the primary production by the CPBM. Consistently with the results obtained by Berline et al. [Berline, L., Brankart, J-M., Brasseur, P., Ourmières, Y., Verron, J., 2007. Improving the physics of a coupled physical–biogeochemical model of the North Atlantic through data assimilation: impact on the ecosystem. J. Mar. Syst. 64 (1–4), 153–172] with a comparable assimilative model, it is shown that the assimilation of physical data alone can improve the representation of the mixed layer depth, but the impact on the ecosystem is rather weak. In some situations, the physical data assimilation can even worsen the ecosystem response for areas where the prior nutrient distribution is significantly incorrect. However, these experiments also show that the combined assimilation of physical and nutrient data has a positive impact on the phytoplankton patterns by comparison with SeaWiFS ocean colour data, demonstrating the good complementarity between SST, altimetry and in situ nutrient data. These results suggest that more intensive in situ measurements of biogeochemical nutrients are urgently needed at basin scale to initiate a permanent monitoring of oceanic ecosystems.     

Using the SKAGEX dataset for evaluation of ocean model skills

Numerical ocean models are now being applied in numerous oceanographic studies. However, the qualities of the model results are often uncertain and there is a great need for standards and procedures for evaluation of the skills of numerical general circulation models. In this paper measurements from repeated hydrographical sections across Skagerrak taken in 1990, the SKAGEX dataset, are used to evaluate the skills of two σ-coordinate ocean models and to study the sensitivity of these models to model parameters. A methodology for quantification of model skills based on observations from repeated hydrographical sections in general is suggested. Area averages of absolute differences are for Skagerrak completely dominated by the discrepancies in the upper few meters of the ocean and may not be used to assess models' abilities to reproduce the fields in the larger and deeper part of the ocean. Therefore, discrepancies between average values in time from the observed fields and time averaged values from model outputs are related to the natural variability of the fields. The numbers produced with the suggested measure are relative numbers that will be specific for each section and for each series of observation. Ideally we would therefore like to see the measures computed for a number of sections for various models and choices of model parameters in order to assess model skills. The value of the SKAGEX dataset as a tool for model improvements is demonstrated. Evidence to support the importance of applying non-oscillatory, gradient preserving advection schemes in areas with sharp density fronts is given. The method is used to identify that the forcing/initial values/boundary values for the temperature field are inferior to the corresponding values for the salinity field. With the present coarse resolution, 11 layers in the vertical, it is shown that it is far from obvious that the quality of the model results improve when replacing simple Richardson number formulations for vertical mixing processes with higher order turbulence closure in the Skagerrak area.     

Comparison between a reanalyzed product by 3-dimensional variational assimilation technique and observations in the Ulleung Basin of the East/Japan Sea

Reanalyzed products from a MOM3-based East Sea Regional Ocean Model with a 3-dimentional variational data assimilation module (DA-ESROM), have been compared with the observed hydrographic and current datasets in the Ulleung Basin (UB) of the East/Japan Sea (EJS). Satellite-borne sea surface temperature and sea surface height data, and in-situ temperature profiles have been assimilated into the DA-ESROM. The performance of the DA-ESROM appears to be efficient enough to be used in an operational ocean forecast system.Comparing with the results from Mitchell et al. [Mitchell, D. A., Watts, D. R., Wimbush, M., Teague, W.J., Tracey, K. L., Book, J. W., Chang, K.-I., Suk, M.-S., Yoon, J.-H., 2005a. Upper circulation patterns in the Ulleung Basin. Deep-Sea Res. II, 52, 1617-1638.], the DA-ESROM fairly well simulates the high variability of the Ulleung Warm Eddy and Dok Cold Eddy as well as the branching of the Tsushima Warm Current in the UB. The overall root-mean-square error between 100 m temperature field reproduced by the DA-ESROM and the observed 100-dbar temperature field is 2.1 °C, and the spatially averaged grid-to-grid correlation between the two temperature fields is high with a mean value of 0.79 for the inter-comparison period.The DA-ESROM reproduces the development of strong southward North Korean Cold Current (NKCC) in summer consistent with the observational results, which is thought to be an improvement of the previous numerical models in the EJS. The reanalyzed products show that the NKCC is about 35 km wide, and flows southward along the Korean coast from spring to summer with maximum monthly mean volume transport of about 0.8 Sv in August–September.     

Evaluating seaport operation and capacity analysis—preliminary methodology

The marine terminal in general and the container terminal in particular is the physical link between ocean and land modes of transportation. Port capacity is commonly defined as the amount of cargo that can be handled by a port per time period, usually a year. For containers it is the number of processed containers per year, for bulk cargo-tons or pallets per year and for Ro-Ro cargoes-autos per year. As the number of ships and the amount of cargo passing through the port increase, a point is reached at which some elements of the port system are functioning near or beyond their maximum productivity rate. The concept that should direct the terminal capacity analysis is controlled by the terminal activity that determines the lowest capacity rate per unit of time, during one cycle of operation. Four performance measures are introduced in order to establish a recommended methodology that quantifies the port's quality of operation and also functions as a tool for decision makers to justify a required investment. The methodology should become a standard that might turn into a requirement for examining every port's level of service.