The Italian Group of Operational Oceanography (GNOO) has developed and maintains the Adriatic Forecasting System (AFS) operational since April 2003 (Oddo et al. 2005, 2006, Guarnieri et al. 2008).

The ocean model used is POM (Princeton Ocean Model). It is a three-dimensional finite difference, free surface, primitive equation numerical model, utilizing the Boussinesque and the hydrostatic approximation and a split mode time step. The model contains a second-order turbulence closure submodel providing the vertical mixing coefficients (Mellor and Yamada, 1982). Horizontal diffusion is parameterized following the scheme of Smagorinsky (1993), as coded into POM by Mellor and Blumberg (1985). Density is calculated by an adaptation of the UNESCO equation of state devised by Mellor.

Compared to the standard version of POM an important change in the model structure has been implemented into this version of POM. The default POM centred difference scheme for the advection of tracers has been substituted with the MUSCL scheme (Estubier and Lévy, 2000).

The POM model has been implemented in the Adriatic Sea (this implementation is called AREG: Adriatic REGional model), and its domain encompasses the whole Adriatic basin and extends south of the Otranto channel into the northern Ionian Sea, where the only open boundary is located. Grid points located over the Thyrrenian Sea have been masked out. The AREG grid has a horizontal resolution of about 1/45° (about 2.2 km), on 31 σ-layers. The bottom topography was obtained from the U.S. Navy 1/60° bathymetric database DBDB1, by bilinear interpolation of the depth data into the model grid. All the depths shallower than 10 metres have been flattened to 10. The model domain and the bathymetry are shown in below.


The surface fluxes are interactively computed using model predicted sea surface temperature and realistic atmospheric data provided by the European Centre for Medium Range Weather Forecast (ECMWF) with a frequency of six hours and a resolution of 0.125°. Realistic fresh water has been implemented in the surface boundary condition for the vertical velocity. The river input into the basin has been implemented through river climatology (Raicich, 1994) for all the rivers except for the Po, which is a very important forcing for the Adriatic Sea. The Po runoff mplemented is on a daily basis. For what concerns the simulations the runoff values are calculated from the observed level of the river at the closure cross section of Pontelagoscuro, while for what concerns the forecast, the last observed value is persisted for the whole period. The runoff is partitioned over several grid points approximately representing the proportion of the fresh water discharged through the mouth of the delta (Provini et al. 1992). The precipitation comes from the global climatological monthly means of Legates and Wilmott, 1990.

The initial and lateral boundary conditions for temperature, salinity and velocity come from the Mediterranean Forecasting System, MFS (Pinardi et al, 2003, Tonani et al, 2008). The lateral boundary conditions are taken on a daily basis. Since December 2008 the tidal signal has been introduced in the model through the southern lateral open boundary conditions on barotropic velocities as proposed by Flather in 1976.

A nine day ocean forecast for all the Adriatic Basin is done and released via ftp and on a web bulletin every day.

In order to have the best estimate of the state of the sea for the production of a new forecast, the daily cycle of the system is combined once a week with a weekly cycle, necessary in order not to have a degeneration of the system. In fact, since the Adriatic system is nested into the Mediterranean System (, and this latter is based on a weekly cycle of assimilation, also the Adriatic System needs to be rewinded once a week to have the best estimate of the sea state for the forecast runs. This rewind of the simulation cycle takes place every Wednesday morning, and it provides the initial conditions for the 9-day forecast. During the other days of the week the forecast is based on a daily simulation which provides the initial condition for the next 9-day forecast.

A scheme of the system is shown in the panel below.



Rita Lecci

  • Blumberg, A.F. & G.L. Mellor. 1987. A description of a three-dimensional coastal ocean circulation model. In: N.S. Heaps (Editor). Three-dimensional coastal ocean models. American Geophysical Union, Washington D.C., pp. 1-16.
  • Estubier A., and Lévy M. Quel schéma numérique pour le transport d'organismes biologiques par la circulation océanique. Note Techniques du Pôle de modélisation, Institut Pierre-Simon Laplace; 2000 : 81pp.
  • Flather, R.A., 1976. A tidal model of the northwest European continental shelf. Memories de la Societe Royale des Sciences de Liege 6 (10), 141–164.
  • Guarnieri, A. P. Oddo, M. Pastore, N. Pinardi. 2008. “The Adriatic Basin Forecasting System new model and system development”. Coastal to Global Operational Oceanography: Achievements and Challenges. Eds. H. Dahlin, N.C Fleming, and S.E. Petersson. Proceeding of 5th EuroGOOS Conference, Exeter (accepted)
  • Legates, D. R. & C. J. Wilmott. 1990. Mean seasonal and spatial variability in a gauge corrected global precipitation. Int. J. Climatol., 10: 121-127.
  • Mellor, G. L. & T. Yamada. 1982. Development of a turbulence closure submodel for geophysical fluid problems. Rev. Geophys. Space Phys., 20: 851-875.
  • Oddo, P., N. Pinardi & M. Zavatarelli. 2005. A numerical study of the interannual variability of the Adriatic Sea (2000-2002). Sci. Total Environ., 353 (2005); pp 39-56.
  • Oddo, P., N. Pinardi, M. Zavatarelli and A. Colucelli. The Adriatic Basin forecasting system, 2006, Acta Adriatica, 47(Suppl):169-184.
  • Oddo P. and N. Pinardi 2007. Lateral open boundary conditions for nested limited area models: A scale selective approach. Ocean Modelling 20 (2008) 134–156.
  • Pinardi, N., I. Allen, E. Demirov, P. De Mey, G. Korres, A. Lascaratos, P.Y. Traon, C. Maillard, G. Manzella & C. Tziavos. 2003. The Mediterranean ocean forecasting system: first phase of implementation (1998-2001). Ann. Geophys., 21: 3-20.
  • Provini, A., G. Crosa & R. Marchetti. 1992. Nutrient export from the Po and Adige river basin over the last 20 years. Sci. Total Environ., Suppl.: 291-313.
  • Raicich, F. 1994. Note on flow rates of the Adriatic rivers. Technical Report. CNR Istituto Talassografico Sperimentale, Trieste. RF 02/94, 8 pp.
  • Smagorinsky, J. 1993. Some historical remarks on the use of nonlinear viscosities. In: Galperin, B. & S.A. Orszag (Editors). Large eddy simulations of complex engineering and geophysical flows. Cambridge Univ. Press, New York, pp. 3-36.
  • Tonani M., N.Pinardi, S. Dobricic, I. Pujol and C. Fratianni, (2008). A High Resolution Free Surface Model on the Mediterranean Sea. Ocean Science, 4, 1-14.
  • Tonani M., N.Pinardi, M.Adani, A.Bonazzi, G.Coppini, M.DeDominicis, S.Dobricic, M.Drudi, N.Fabbroni, C.Fratianni, A.Grandi, S.Lyubartsev, P.Oddo, D.Pettenuzzo, J.Pistoia, and I.Pujol: “The Mediterranean Ocean Forecasting System”, Coastal to Global Operational Oceanography: Achievements and Challenges. Eds. H. Dahlin, N.C Fleming, and S.E. Petersson. Proceeding of 5th EuroGOOS Conference, Exeter, 2008 (Accepted).