Jonathan BEUVIER, Météo-France/ENSTA Florence SEVAULT, Météo-France Marine HERRMANN, Météo-France

1. Modelling the Mediterranean Seainterannual variability over the last 40 years:focus on the Eastern Mediterranean Transient (EMT) Jonathan BEUVIER, Météo-France/ENSTA Florence SEVAULT, Météo-France Marine HERRMANN, Météo-France Karine BÉRANGER, ENSTA Samuel SOMOT, Météo-France 2009 NEMO users meeting - Paris

2. Introduction • Eastern Mediterranean Transient: switch of the main source of dense waters in the Eastern Mediterranean in the early 1990’s, from the Adriatic Sea to the Aegean Sea. • Need of modelling to understand the EMT (Roether et al. 2007). • Interests: • variabilities at different time-scales: requires long and stable simulations, • good test for atmosphere and ocean models, • improves knowledge of possible past or future EMT, • need to be better simulated with realistic simulations (Samuel et al. 1999, Nittis et al. 2003, Bozec et al. 2006). • Questions: • are we able to reproduce the different phases of the EMT (winter deep convection, filling, overflow and spreading)? • what are the key processes that trigger the EMT? 2009 NEMO users meeting Paris

3. The NEMOMED8 configuration Adriatic Sea Gulf of Lions Otranto Strait Strait of Gibraltar Sicily Strait Aegean Sea Atlanticbuffer zone Ionian basin Levantine basin • Mediterranean version (Sevault et al. 2009), based on NEMO-v2 • Resolution of 1/8° x 1/8°cos(lat) (9 to 12 km with square meshes) • Grid tilted and stretched at Gibraltar (up to 6km resolution) • Z-coordinate partial steps (43 vertical Z-levels) • Atlantic buffer zone with 3D T-S damping (11°W to 7.5°W) • Explicit river forcing for 33 rivers + Black Sea (simulated as a river) 2009 NEMO users meeting Paris

4. Physics used in this study • Filtered free surface (with transfer of the evaporated water in the buffer zone). • TVD scheme for tracers. • Iso-neutral diffusion for tracers (laplacian operator). • Horizontal diffusion for momentum (bilaplacian operator). • Vertical diffusion based on TKE closure scheme. • EEN (energy and enstrophy conserving) scheme. • Feedback coefficient for SST damping: -40 W/m²/K. • No-slip condition for the lateral momentum boundary. • Non-linear bottom friction. 2009 NEMO users meeting Paris

5. Three 1961-2000 hindcast simulations • 15-year spin-up. • 40 years forced by ARPERA (dynamical downscaling of ERA40, 50 km of resolution over the Med). • SST relaxation (ERA40). • No SSS relaxation: monthly water flux correction. Horizontal grid and relief of ARPERA(Herrmann & Somot 2008) Climatological (---) and interannual(___) river and Black Sea runoffs Total Interannual Atlantic anomalies at 176m T S Black Sea Nile Rhone Po 2009 NEMO users meeting Paris

6. Total Med heat content Total Med salt content NM8-atl-rivNM8-rivNM8-clim Med salt content per layer Med heat content per layer 0-150m 150-600m 600m-bottom Global validation of the simulations • Comparison with the interannual climatology of Rixen et al., 2005 (mean and standard deviation): • T: good correlations (>0.7) despite global bias (+0.1°C), accurate surface layer, intermediate layer too warm (+0,2°C), trend in the bottom layer. • S: well simulated in average but not enough variability, surface layer too fresh, intermediate layer too salty. 2009 NEMO users meeting Paris

7. Agean winter convection in NM8-atl-riv Winter (NDJF) surface flux anomalies over the Aegean • Strong heat and water losses in winters 1991-92 and 1992-93, but not only. • Formation of dense and deep waters (> 29,2 kg/m3) during many successive winters in the 1970’s and the 1980’s. • Annual formation rate for σ>29.2kg/m3: 0,5 Sv in 1992 and 1,2 Sv in 1993. • In 1993, about 75% of the Aegean Sea filled by waters denser than 29,2 kg/m3. Net surface heat flux anomalies (W/m²) Net surface water flux anomalies (mm/day) Qtot R+P-E Monthly volume of Aegean dense waters (m3) Annual formation rate (Sv) ---- σ>29,2 kg/m3___ σ>29,3 kg/m3 ---- σ>29,2 kg/m3___ σ>29,3 kg/m3 2009 NEMO users meeting Paris

8. Potential density (kg/m3) on the Cretan Arc Straits sills Filling and overflowing of Aegean deep waters • Increase of potential density above the Cretan Arc Straits sills. • Overflow of warm, salty and dense waters toward the Ionian and Levantine seas. Location of the Cretan Arc Straits and paths of the outflowing waters Karpathos777m Antikithira534m Kassos542m 2009 NEMO users meeting Paris

9. May 1993 May 1994 Dispersion in the Eastern Mediterranean Salinity (colors, in psu) and depth (lines, in m) of the 29,165 kg/m3 isopycnal • Simulated EMT-waters warmer (0,3°C), saltier (0,05psu) and less dense (-0,03kg/m3) than the observed EMT-waters. • Simulated EMT-waters less dense than the bottom waters of the Eastern Mediterranean.=> they sink to a depth of 2200m (not to the bottom as observed). 2009 NEMO users meeting Paris

10. Impact of the interannual hydrological forcings • Main characteristics of the EMT not modified: same formation rates in 1992 and 1993, overflow, dispersion and sinking. • Affects mainly the chronology of the Aegean deep convection in the 1970’s and 1980’s.=> main motor of the EMT: atmospheric forcing. Monthly volume of Aegean dense waters (m3, left) and associated annual formation rate (Sv, right)---- σ>29,2 kg/m3 ___ σ>29,3 kg/m3 NM8-atl-rivNM8-rivNM8-clim NM8-atl-rivNM8-rivNM8-clim 2009 NEMO users meeting Paris

11. Conclusion • The EMT is well simulated: • realistic chronology (surface losses, winter deep convection), • good estimates of the dense water formation rates, • representation of an overflow and of a spreading. • The interannual hydrological forcings mostly impact the Aegean convection in the 1970’s and 1980’s. • Perspectives: • studies of the variability in other Mediterranean locations (e.g. the convection in the Gulf of Lions, the cascading at the Otranto Strait, …) • EMT modelling: • improvement of the initial conditions and of the spin-up, • test of new physical parametrisations (horizontal diffusion, vertical mixing, …), • use of higher resolution models (atmosphere and ocean). 2009 NEMO users meeting Paris