ALBERT

Description: Phytoplankton blooms in the northwestern Mediterranean Sea are seasonal events that mainly occur in a specific area comprising the Gulf of Lion and the Provençal basin, where they are promoted by a general cyclonic circulation, strong wind-driven mixing and subsequent re-stratification of the water column. At the southern boundary of this area, a persistent density front known as the north Balearic front can be found. The front is presumed to cause an early phytoplankton bloom in its vicinity because (a) it enhances the transport of nutrients into the euphotic layer and (b) it promotes the speedy re-stratification of the water column (through frontal instabilities). In February and March 2013, a glider, equipped with a CTD (conductivity, temperature, and depth device) and a fluorometer, was deployed on a mission that took it from the Balearic Islands to Sardinia and back. The frontal zone was crossed twice, once during the outbound leg and the once on the return leg. The data provided by the glider clearly showed the onset of a bloom soon after a decrease in wind-driven turbulent convection and mixing. The in situ observations were supported and confirmed by satellite imagery. It is shown that frontal dynamics play a key role in the promotion and acceleration of re-stratification, which is a necessary pre-conditioning factor for the onset of blooms much like other relevant processes such as an enhanced biological pump. Swift re-stratification stimulates new production by inhibiting mixing. Finally, viewing the blooming phenomenon from a regional perspective, it seems that Sverdrup's critical depth model applies in the northern well-mixed area whereas, in the south, front-related re-stratification seems to be the principal cause.

Description: The impact of the assimilation of MyOcean sea level anomalies along-track data on the analyses of the Sicily Channel Regional Model was studied. The numerical model has a resolution of 1/32° degrees and is capable to reproduce mesoscale and sub-mesoscale features. The impact of the SLA assimilation is studied by comparing a simulation (SIM, which does not assimilate data) with an analysis (AN) assimilating SLA along-track multi-mission data produced in the framework of MyOcean project. The quality of the analysis was evaluated by computing RMSE of the misfits between analysis background and observations (sea level) before assimilation. A qualitative evaluation of the ability of the analyses to reproduce mesoscale structures is accomplished by comparing model results with ocean colour and SST satellite data, able to detect such features on the ocean surface. CTD profiles allowed to evaluate the impact of the SLA assimilation along the water column. We found a significant improvement for AN solution in terms of SLA RMSE with respect to SIM (the averaged RMSE of AN SLA misfits over 2 years is about 0.5 cm smaller than SIM). Comparison with CTD data shows a questionable improvement produced by the assimilation process in terms of vertical features: AN is better in temperature while for salinity it gets worse than SIM at the surface. This suggests that a better a-priori description of the vertical error covariances would be desirable. The qualitative comparison of simulation and analyses with synoptic satellite independent data proves that SLA assimilation allows to correctly reproduce some dynamical features (above all the circulation in the Ionian portion of the domain) and mesoscale structures otherwise misplaced or neglected by SIM. Such mesoscale changes also infer that the eddy momentum fluxes (i.e. Reynolds stresses) show major changes in the Ionian area. Changes in Reynolds stresses reflect a different pumping of eastward momentum from the eddy to the mean flow, in turn influencing transports through the channel.

Description: The effects of the 2003 European heatwave on the sea surface layer of the Central Mediterranean were studied using a regional 3-D ocean model. The model was used to simulate the period 2000 to 2004 and its performance was validated using remotely-sensed and in situ data. Analysis of the results focused on changes in the Sea Surface Temperature (SST) and on changes to the surface and sub-surface current field. This permitted us to identify and quantify the anomalies of atmospheric and sea surface parameters that accompanied the heatwave. The dominant annual cycle in each variable was first removed and a wavelet analysis then used to locate anomalies in the time-frequency domain. We found that the excess heating affecting the sea surface in the summer of 2003 was related to a significant increase in air temperature, a decrease in wind stress and reduction of all components of the upward heat flux. The monthly averages of the model SST were found to be in good agreement with remotely-sensed data during the period studied, although the ocean model tended to underestimate extreme events. The spatial distribution of SST anomalies as well as their time-frequency location was similar for both the remotely-sensed and model temperatures. We also found, on the basis of the period of the observed anomaly, that the event was not limited to the few summer months of 2003 but was part of a longer phenomenon. Both the model results and experimental data suggest the anomalous heating mainly affected the top 15 m of ocean and was associated with strong surface stratification and low mixing. The skill of the model to reproduce the sub-surface hydrographic features during the heatwave was checked by comparison with temperature and salinity measurements. This showed that the model was generally in good agreement with observations. The model and observations showed that the anomalous warming also modified the currents in the region, most noticeably the Atlantic Ionian Stream (AIS) and the Atlantic Tunisian Current (ATC). The AIS was reduced in intensity and showed less meandering, mainly due to the reduced density gradient and low winds, while the ATC was enhanced in strength, the two currents appearing to modulate each other in order to conserve the total transport of Modified Atlantic Water.

Description: The spatial and temporal variability of eddy and mean kinetic energy of the Central Mediterranean region has been investigated, from January 2008 to December 2010, by mean of a numerical simulation mainly to quantify the mesoscale dynamics and their relationships with physical forcing. In order to understand the energy redistribution processes, the baroclinic energy conversion has been analysed, suggesting hypotheses about the drivers of the mesoscale activity in this area. The ocean model used is based on the Princeton Ocean Model implemented at 1/32° horizontal resolution. Surface momentum and buoyancy fluxes are interactively computed by mean of standard bulk formulae using predicted model Sea Surface Temperature and atmospheric variables provided by the European Centre for Medium Range Weather Forecast operational analyses. At its lateral boundaries the model is one-way nested within the Mediterranean Forecasting System operational products. The model domain has been subdivided in four sub-regions: Sardinia channel and southern Tyrrhenian Sea, Sicily channel, eastern Tunisian shelf and Libyan Sea. Temporal evolution of eddy and mean kinetic energy has been analysed, on each of the four sub-regions, showing different behaviours. On annual scales and within the first 5 m depth, the eddy kinetic energy represents approximately the 60 % of the total kinetic energy over the whole domain, confirming the strong mesoscale nature of the surface current flows in this area. The analyses show that the model well reproduces the path and the temporal behaviour of the main known sub-basin circulation features. New mesoscale structures have been also identified, from numerical results and direct observations, for the first time as the Pantelleria Vortex and the Medina Gyre. The classical kinetic energy decomposition (eddy and mean) allowed to depict and to quantify the permanent and fluctuating parts of the circulation in the region, and to differentiate the four sub-regions as function of relative and absolute strength of the mesoscale activity. Furthermore the Baroclinic Energy Conversion term shows that in the Sardinia Channel the mesoscale activity, due to baroclinic instabilities, is significantly larger than in the other sub-regions, while a negative sign of the energy conversion, meaning a transfer of energy from the Eddy Kinetic Energy to the Eddy Available Potential Energy, has been recorded only for the surface layers of the Sicily Channel during summer.

Description: As a part of the project Mediterranean Network to Assess and Upgrade Monitoring and Forecasting Activity in the Region (MAMA) we implemented a high resolution nested hydrodynamic model (1/40° horizontal grid, 16 sigma levels) for the coastal, shelf and open sea areas off the Lebanese coast, East Levantine Basin of the Eastern Mediterranean Sea. The Lebanese Shelf Model (LSM) is a version of the Princeton Ocean Model (POM). It is nested in a coarse resolution model the Aegean Levantine Eddy Resolving Model (1/20° horizontal grid, 25 sigma levels), ALERMO, that covers the Eastern Mediterranean. The nesting is one way so that velocity, temperature, and salinity along the open boundaries are interpolated from the relevant coarse model variables. Numerical simulations have been carried out under climatological surface and lateral forcing. Due to the relatively small domain, the results closely follow the simulation of the intermediate model with more details especially over the narrow shelf region. Simulations reproduce main circulation features and coastal circulation characteristics over the eastern Levantine shelf. This paper describes the modeling system setup, compares the simulations with the corresponding results of the coarse model ALERMO, and with the observed climatological circulation characteristics in the Levantine Basin off the Lebanese coast.

Description: The effects of anomalous weather conditions on the sea surface layer over the Central Mediterranean were studied with an eddy resolving regional ocean model by performing a 5-year long simulation from 2000 to 2004. The focus was on surface heat fluxes, temperature and dynamics. The analysis of the time series of the selected variables permitted us to identify and quantify the anomalies of the analysed parameters. In order to separate the part of variability not related to the annual cycle and to locate the anomalies in the time-frequency domain, we performed a wavelet analysis of anomalies time series. We found the strongest anomalous event was the overheating affecting the sea surface in the summer of 2003. This anomaly was strictly related to a strong increase of air temperature, a decrease of both wind stress and upward heat fluxes in all their components. The simulated monthly averages of the sea surface temperature were in a good agreement with the remotely-sensed data, although the ocean regional model tended to underestimate the extreme events. We also found, on the basis of the long-wave period of the observed anomaly, this event was not limited to the few summer months, but it was probably part of a longer signal, which also includes negative perturbations of the involved variables. The atmospheric parameters responsible for the overheating of the sea surface also influenced the regional surface and sub-surface dynamics, especially in the Atlantic Ionian Stream and the African Modified Atlantic Water current, in which flows seem to be deeply modified in that period.

Description: The Sicily Channel Regional Model forecasting system was tested using an optimization package for the initial and lateral boundary conditions. Spurious high frequency oscillations during the spin-up time were successfully reduced both in duration and magnitude by optimizing the time tendency of the free surface elevation using the Variational Initialization and Forcing Platform method developed in the framework of the Mediterranean Forecasting System Towards the Environmental Prediction project. The effect of optimization was most profound for the free surface elevation, where all oscillations with periods shorter than 4 h were suppressed. The overall forecast skill was assessed on a 5 day case study starting on 6 April 2005, characterized by a fast passage of a deepening atmospheric low-pressure field with strong winds and marked wind direction change. We compared the predicted variables with in-situ and remotely sensed data. The forecasts of temperature, including the sea surface temperature, and salinity were quite successful, while the forecasted currents, especially within the surface layer, were not in good agreement with the measurements.

Description: A twin numerical experiment was conducted in the seas of Sardinia (Western Mediterranean) to assess the impact, at coastal scales, of the use of relative winds (i.e. taking into account ocean surface currents) in the computation of heat and momentum fluxes through bulk formulas. The model, the Regional Ocean Modeling System (ROMS), was implemented at 2 km of resolution in order to well resolve (sub-)mesoscale dynamics. Small changes (1–2%) in terms of spatially-averaged fluxes correspond to quite large spatial differences of such quantities (up to 15–20%) and to comparably significant differences in terms of mean velocities of the surface currents. Wind power input of the wind stress to the ocean surface P results also reduced by a 15%, especially where surface currents are stronger. Quantitative validation with satellite SST suggests that such a modification on the fluxes improves the model solution especially in areas of cyclonic circulation, where the heat fluxes correction is predominant in respect to the dynamical correction. Surface currents changes above all in their fluctuating part, while the stable part of the flow show changes mainly in magnitude and less in its path. Both total and eddy kinetic energies of the surface current field results reduced in the experiment where fluxes took into account for surface currents. Dynamically, the largest correction is observed in the SW area where anticyclonic eddies approach the continental slope. This reduction also impacts the vertical dynamics and specifically the local upwelling that results diminished both in spatial extension as well in magnitude. Simulations suggest that, even at local scales and in temperate regions, it is preferable to take into account for such a component in fluxes computation. Results also confirm the tight relationship between local coastal upwelling and eddy-slope interactions in the area.