Abstract
Air-sea heat and freshwater water fluxes in the Mediterranean Sea play a crucial role in dense water formation. Here, we compare estimates of Mediterranean Sea heat and water budgets from a range of observational datasets and discuss the main differences between them. Taking into account the closure hypothesis at the Gibraltar Strait, we have built several observational estimates of water and heat budgets by combination of their different observational components. We provide then three estimates for water budget and one for heat budget that satisfy the closure hypothesis. We then use these observational estimates to assess the ability of an ensemble of ERA40-driven high resolution (25 km) Regional Climate Models (RCMs) from the FP6-EU ENSEMBLES database, to simulate the various components, and net values, of the water and heat budgets. Most of the RCM Mediterranean basin means are within the range spanned by the observational estimates of the different budget components, though in some cases the RCMs have a tendency to overestimate the latent heat flux (or evaporation) with respect to observations. The RCMs do not show significant improvements of the total water budget estimates comparing to ERA40. Moreover, given the large spread found in observational estimates of precipitation over the sea, it is difficult to draw conclusions on the performance of RCM for the freshwater budget and this underlines the need for better precipitation observations. The original ERA40 value for the basin mean net heat flux is −15 W/m2 which is 10 W/m2 less than the value of −5 W/m2 inferred from the transport measurements at Gibraltar Strait. The ensemble of heat budget values estimated from the models show that most of RCMs do not achieve heat budget closure. However, the ensemble mean value for the net heat flux is −7 ± 21 W/m2, which is close to the Gibraltar value, although the spread between the RCMs is large. Since the RCMs are forced by the same boundary conditions (ERA40 and sea surface temperatures) and have the same horizontal resolution and spatial domain, the reason for the large spread must reside in the physical parameterizations. To conclude, improvements are urgently required to physical parameterizations in state-of-the-art regional climate models, to reduce the large spread found in our analysis and to obtain better water and heat budget estimates over the Mediterranean Sea.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The Mediterranean Sea can be considered to be a thermodynamic machine that exchanges water and heat with the Atlantic Ocean through the Strait of Gibraltar and with the atmosphere through its surface. On average, the Mediterranean basin shows an excess of evaporation over freshwater inputs and a heat loss through air-sea interaction. It has an overall freshwater deficit, as the loss to the atmosphere by evaporation is larger than the gains by precipitation and runoff from the main rivers and input from the Black Sea. The total heat budget is negative, that is, the Mediterranean Sea loses more energy than it gains. These losses of freshwater and heat are compensated by the two-layer exchange at the Strait of Gibraltar comprising a relatively warm and fresh (15.4°C, 36.2 psu) upper water inflow and a relatively cooler and saltier (13°C, 38.4 psu) outflow to the Atlantic (Bryden et al. 1994; Tsimplis and Bryden 2000). On long timescales, the losses of freshwater and heat from the surface are compensated and the net salt flux is close to zero. Note that this closure hypothesis does not apply for short periods of time, as at monthly or yearly timescales (Castellari et al. 1998; Brankart and Pinardi 2000; Pettenuzzo et al. 2010). One main concern is to what extent the closure hypothesis is true under a changing climate. Some preliminary studies show a large drying of the Mediterranean area, leading to an increase of the freshwater loss at the end of the twenty first century (Somot et al. 2006; Mariotti et al. 2008; Sanchez-Gomez et al. 2009). The implications of these changes in the water cycle for the properties of the Mediterranean water masses and on the exchange at Gibraltar Strait need to be elucidated.
The circulation in the Mediterranean Sea is determined to a large extent by the air-sea exchanges of heat and freshwater, which depend on the meteorological conditions and ocean characteristics (Timplis et al. 2006). The water and heat fluxes play a crucial role in dense water formation, and hence in the Mediterranean Thermohaline Circulation (MTHC) (Béthoux et al. 1999). Consequently they affect the Mediterranean water mass characteristics (temperature, salinity, density) and then can potentially influence the Atlantic Ocean circulation by changing the properties of the Mediterranean Outflow water (MOW) (Béthoux et al. 1999; Potter and Lozier 2004; Artale et al. 2005; Millot et al. 2006). The Mediterranean water and heat budgets can also influence the atmospheric water content, the properties of the low level atmosphere and the occurrence of coastal intense precipitation events (Li 2006; Lebeaupin et al. 2006; Gimeno et al. 2010).
Improving our knowledge of the water and heat budgets and their variability is a challenging goal for the observational and modelling community of the Mediterranean region. Obtaining accurate estimates of every term in the water and heat budgets is crucial for understanding the Mediterranean ocean circulation and climate, and their evolution under climate change. Accurate modelling the Mediterranean air-sea fluxes will also provide the long-term atmospheric forcing for regional ocean models (Madec et al. 1991; Castellari et al. 2000; Somot et al. 2006; Beuvier et al. 2010) and will help to improve the design of fully-coupled Atmosphere–Ocean Regional Climate Models (Somot et al. 2008).
The Mediterranean basin mean water and heat flux estimates available in the literature over the last 30 years, vary depending on the processing techniques and datasets used. They can be estimated either from the atmospheric branch of the water cycle or from the oceanic branch, that is to say the Gibraltar net transport. Estimates obtained using these two independent methods are not fully consistent especially for the heat budget, and lead to a wide range of uncertainty. For the water budget, the atmospheric branch estimates range from a basin average of 520 to 950 mm/year (Castellari et al. 1998; Gilman and Garret 1994; Béthoux et al. 1999; Mariotti et al. 2002). This large uncertainty range shows the need to improve the dataset quality, since over the Mediterranean Sea, evaporation and precipitation observations are very sparse. Concerning the net water transport at the Strait of Gibraltar, large discrepancies also exist. Direct measurements of the inflow rate are extremely difficult to perform and currently there are very few. In the last 20 years, the estimations of the oceanic branch range from 0.04 to 0.09 Sv (1 Sv = 106 m3/s) that is to say equivalent to a basin average surface flux of between 515 and 1,150 mm/year (Bryden and Kinder 1991; Bryden et al. 1994; Tsimplis and Bryden 2000; Candela 2001; Baschek et al. 2001; Garcia de La Fuente et al. 2007). More recently estimates from the available reanalyses (NCEP/NCAR, Kalnay et al. 1996 or ERA40, Uppala et al. 2004) have also been proposed. They vary between 391 to 524 mm/year (Bouktir and Barnier 2000; Mariotti et al. 2002), showing an underestimation of the water loss by the Mediterranean sea surface. This underestimation may be due to the use of coarse spatial resolution models that are unable to resolve the complexity of the Mediterranean basin (orography, coast-line, islands). Thus, the use of higher resolution atmosphere models may allow resolving this problem (see Elguindi et al. 2009 for an impact of the resolution on the Mediterranean Sea water budget). The question of the potential added-value of high-resolution regional climate models (RCM) with respect to the ERA40 reanalysis is one of the scientific issues addressed in the current study. However dynamical downscaling techniques using RCMs is only one way of attempting to improve the low-resolution reanalysis. Another technique recently applied by Pettenuzzo et al. (2010) is to fit adjustment factors to the reanalysis variables in order to re-compute the various terms of the water (or heat) budget. The adjustment factors are applied to the wind, the downward radiation flux, the temperature or the humidity and are fitted with respect to observed estimates. Using this adjustment method, Pettenuzzo et al. (2010) provide an estimate of E–P freshwater loss of 640 mm/year over the 1958–2001 period with 1180 mm/year for the evaporation and 530 mm/year for the precipitation but do not provide estimate for the river and Black Sea freshwater inputs.
The net heat transport through the Strait of Gibraltar has been estimated by using mooring-based measurements. Results vary between 3 and 10 W/m2 (Béthoux 1979; Bunker et al. 1982; McDonald et al. 1994) with the most recent of these values being about 5 W/m2. The positive sign indicates a heat gain of the Mediterranean Sea. This gain should be compensated by an equivalent heat loss through the Mediterranean Sea surface, i.e. a basin mean net surface heat flux of around −5 W/m2, as the climate change signal remains weak up to now and so closure is expected to hold. However, earlier studies have shown that climatological estimates of the mean heat budget, determined from the surface heat fluxes, show discrepancies compared with the Strait of Gibraltar estimates of up to 20–30 W/m2 (Garrett et al. 1993; Gilman and Garret 1994; Artale et al. 2005; Tsimplis and Bryden 2000). These authors suggested that the bias is caused by a combination of overestimated shortwave gain and water vapour and underestimated longwave loss, latent and sensible heat fluxes. Subsequently a heat budget estimate over the Mediterranean Sea of +6 W/m2 has been obtained from a modified version of the ship-based NOC flux dataset (Josey et al. 1999), using radiative flux formulae appropriate for this basin as discussed later, which reduces the discrepancy to about 10 W/m2. In their analysis, Pettenuzzo et al. (2010) attempted to close the heat budget by making various plausible adjustments to the ERA40 meteorological fields (which have an imbalance of −15 W/m2). Using these adjustments with the transfer coefficient scheme of Kondo (1975) they are able to achieve closure with a mean net heat flux of −5 W/m2. However, this scheme is not supported by high quality direct measurements of the heat flux made in recent years by Fairall et al. (2003). Using the more up to date scheme of Kara et al. (2005) based on these measurements they instead obtain a value of 4 W/m2 which does not close the heat budget.
Mediterranean water and heat budgets have also been investigated recently using numerical simulations (Sotillo et al. 2005; Ruiz et al. 2008; Elguindi et al. 2009; Beuvier et al. 2010). The spatial resolution of typical general circulation models (GCMs) is unlikely to represent well the local thermal and dynamical processes occurring in the Mediterranean Sea. Such processes are very important in the Mediterranean area, since the air-sea exchanges (particularly the latent heat flux) strongly depend on a correct resolution for the temperature, humidity and wind effects. In a study using the ARPEGE atmosphere model, it has been shown (Elguindi et al. 2009) that increasing model resolution yields a more intense evaporation and hence water budget estimates that are closer to observed values. The surface heat fluxes have been analyzed (Ruiz et al. 2008) from a high resolution dataset issued from a dynamical downscaling performed with the regional climate model REMO, forced at its boundaries by the NCEP/NCAR reanalysis (HIPOCAS dataset, Sotillo et al. 2005). These authors show that the RCM experiment significantly improves the heat budget estimation when compared with the NCEP/NCAR reanalysis. This is a very promising result from the regional climate modeling community, and suggests that RCMs driven by reanalyses can constitute useful tools to study the air-sea fluxes over the Mediterranean region; the work reported here further develops this idea. We have used an ensemble of limited-area models to evaluate their performance in simulating the water and heat budgets in the Mediterranean Sea. Due to their reduced spatial domain, limited-area models provide an attractive approach allowing high spatial resolution climate simulations at an affordable computational cost. The database comes from the EU-FP6 ENSEMBLES project (Hewitt and Griggs 2004; Sanchez-Gomez et al. 2008) and consists of an ensemble of experiments performed by different RCMs with a spatial resolution of 25 km over the European-Mediterranean domain. This ensemble includes 12 RCMs driven by the ERA40 reanalysis for the period 1961–2000 at their lateral boundaries. A second ensemble is also available in the FP6 ENSEMBLES data base. This one includes the same RCMs but driven by different GCMs.
The main goals of this study are the following:
-
a.
To compare state-of-the-art observational datasets in terms of the heat and water budgets for the Mediterranean Sea;
-
b.
to evaluate the ability of an ensemble of ERA40-driven RCMs in simulating the components and net budgets of the heat and freshwater fluxes over the Mediterranean Sea. The RCM estimates will be compared to observed estimates available in the literature and those obtained in step (a). We will also determine whether the estimates of both heat and water budgets in the Mediterranean Sea provided by the RCMs are in accordance with the closure hypothesis at the Gibraltar Strait. The multi-model approach will allow us to assess the uncertainties associated with the water and heat budget estimates and
-
c.
to analyze for the several variables of water and heat budgets whether there is an added value of RCMs against the driving reanalysis (in our case ERA40).
While significant interannual variations have been observed in the water and heat budgets in the Mediterranean basin (Bouktir and Barnier 2000; Mariotti et al. 2002; Josey 2003; Struglia et al. 2004, Josey et al. 2011), and also long-term as a response to global warming (Mariotti et al. 2008; Sanchez-Gomez et al. 2009), in this work we focus only on the long-term annual estimates and on the seasonal cycle for the current climate.
The structure of the paper is as follows: In Sect. 2 we present the numerical RCM simulations and observational datasets, and a brief description of the methodology to compute the heat and water budgets. In Sect. 3 we present the water and heat budget estimates from different observational datasets and ERA40 and then the validation of RCMs according to these estimates. Finally, in Sect. 4 we summarize and draw our main conclusions.
2 Datasets and methodology
2.1 Regional climate models experiments
A summary of the main characteristics of the RCMs used in this work is presented in Table 1. More details about each individual model can be found in the ENSEMBLES project website: http://ensemblesrt3.dmi.dk. To produce this multi-model ensemble all RCM experiments have been performed for the time period 1961–2000 using six hourly lateral boundary conditions provided by the ERA40 reanalysis at 1.125° horizontal resolution. The sea surface temperature (SST) and sea-ice concentration are also from ERA40 dataset. All models are required to cover the entire Mediterranean Sea, though only a few cover the entire Black Sea basin. The RCMs used their own model setup as well as grid specifications like rotation and number of vertical levels, but a similar horizontal resolution of 25 km. The ENSEMBLES project has produced a second set of dynamical downscaling experiments (not used in this work) with the same horizontal resolution over the same geographical area, but in this case the RCMs are driven by diverse GCMs for the period 1950–2000.
2.2 Water Budget in the Mediterranean Sea
Following Mariotti et al. (2002), the long term mean of the Mediterranean water deficit is approximately equal to the net water flux at the Strait of Gibraltar over a long period of time:
where E is the evaporation, P the precipitation, R is the river discharge into the Mediterranean Sea, B the freshwater input from the Black Sea (B is actually equal to the E–P–R budget of the Black Sea over its own catchment area), and GW the Gibraltar Strait net water transport. In the following, all freshwater flux values are given in mm/year considering a Mediterranean Sea surface equal to 2.5 × 1012 m2. A summary of the observational datasets used to compute the estimates for the different water budget terms is presented in Table 2.
For the evaporation (E) we use estimates of the turbulent latent heat flux from three different products a.) the Objectively Analyzed air-sea Fluxes (OAFlux) (Yu and Weller 2007; b) the Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data set (HOAPS) (Andersson et al. 2007), and c.) a modified Mediterranean Sea version (referred to as NOC hereafter) of the National Oceanography Centre 1.1 (NOC1.1 - referred to as NOC hereafter) dataset (Josey et al. 1999). The OAFlux is a 50 year global dataset on a 1° × 1° grid for the period 1958–2008. This product is a result of merging satellite observations with surface moorings, ship reports, and atmospheric model reanalysed surface meteorology. The HOAPS data is derived from satellite measurements over the ice free global ocean and covers the period 1988–2005 with a resolution of 0.5° × 0.5°. The NOC dataset is based on Voluntary Observing Ships observations from the International Comprehensive Ocean–Atmosphere Data Set (ICOADS, Woodruff et al. 1998), and is presented on a 1° spatial grid for the period 1980–2004. The Mediterranean NOC dataset used in the current study has been developed by modifying the formulae employed to estimate the radiative flux components. In particular, the longwave flux is estimated using the formula of Bignami et al. (1995), which was developed using high quality radiometer measurements made in the Mediterranean Sea during several research ship cruises. The net longwave ocean heat loss obtained with the Bignami formula is typically stronger than that obtained with various other formulae developed over the North Atlantic (Josey et al. 1997, 2003). In addition, the shortwave flux has been corrected for aerosol loading following the method of Gilman and Garret (1994). The NOC dataset is only based on in situ data without satellite or model inputs, contrary to OAFlux or HOAPS. Note also that the higher resolution dataset HOAPS is also the shortest in time (19 years).
The observational datasets for precipitation (P) are the Global Precipitation Climatology Project (GPCP), the CMAP (CPC Merged Analysis of Precipitation) and the HOAPS datasets; each of which has its limitations. The GPCP (Adler et al. 2003) has been built by merging satellite and rain gauge data for the period 1979–2008. The horizontal resolution is 2.5° × 2.5° which is unlikely to be adequate for studying the Mediterranean Sea, especially because it merges land and sea points leading to a likely overestimation of the precipitation estimate. The CMAP precipitation consists of monthly values from January 1979 to July 2008. The spatial coverage is global with a 2.5° × 2.5° horizontal resolution. This dataset are obtained from gauge measurements, satellite data and reanalyses data (Xie and Arkin 1997 and http://daac.gsfc.nasa.gov/CAMPAIGN_DOCS/hydrology/hd_main.shtml). The HOAPS dataset has a higher spatial resolution and provides data only over the sea areas. However, it does not provide data over the Adriatic and Aegean seas because of issues with the land-sea mask in these regions. The HOAPS dataset is also known to underestimate the precipitation as it does not measure precipitation amounts less than 0.1 mm/h.
The river discharge R in the Mediterranean Sea has been provided by the dataset described in Ludwig et al. (2009). It consists of reconstructions of 40-year time series (1960–2000) of river discharges observed for the main rivers or reconstructed by using time series of temperature and precipitation. In our case we have computed the runoff for each of the 11 sub-regional catchments considered in Ludwig et al. (2009). The total runoff in the Mediterranean Sea results from the addition of the values from the 11 catchments. To evaluate the Black Sea discharge in the Mediterranean Sea simulated by the RCMs we have used hydrological data from Stanev et al. (2000). These data provide time series of precipitation, evaporation and runoff for the Black Sea basin over the 1923–1997 period. We recall that to obtain the term B in eq. (1) we compute the P + R–E budget of the Black Sea from Stanev et al. (2000) estimates. Table 2 summarizes all the datasets used for the water budget terms.
For the model runoff computations we have considered only the 7 RCMs containing the runoff field as an output (see Table 3). Most of models are non conservative, that is, the difference E–P over the river catchments is not equilibrated to the river discharge. In this case, it is preferable to use the runoff computed by the internal adjustments and hypotheses of the models. This is the case also for ERA40, whose runoff is not conservative and results in negative precipitation minus evaporation values over several catchments (Hagemman et al. 2005). This is a well known behaviour in the reanalysis data, in which the soil moisture is corrected to decrease the bias of the 2 m temperature. ERA40 runoff presents also other deficiencies as an unrealistic separation of total runoff into surface runoff and drainage and the absence of river routing that gives rise to instantaneous runoff (Betts et al. 2003; Hagemann et al. 2005). Hence in order to avoid misleading interpretations, we have omitted the comparison between the ERA40 and RCMs runoff estimates from our analysis and we only consider the comparisons for the precipitation and evaporation variables.
We have taken into account the missing drainage area in the models’ spatial domain for two of the river catchments: the Nile and the Black Sea. For the former we have replaced the runoff from the model by the observed Nile discharge (≈8 mm/year) in all the models, since the missing drainage area is too large to make a good estimate. In addition since the Aswan dam building at the beginning of the 60 s and the intensive use of irrigation in Egypt, the Nile river can be considered as mainly anthropogenically-driven (see Skliris and Lascaratos 2004; Ludwig et al. 2009, for a discussion of the problem). For the Black Sea we have estimated the missing contribution from the fraction of missing area as in Struglia et al. (2004). To compute the river discharge over each of the Mediterranean catchments we have interpolated the drainage area mask (available at a resolution of 0.5° × 0.5°) to the RCM native grids.
2.3 Heat budget in the Mediterranean Sea
The total heat budget is determined by the radiative and turbulent heat flux components. The radiative terms are the net solar (shortwave) flux at the surface (Qsw) and the net infrared radiation (longwave, QLW). The turbulent components are the latent (QLH) and the sensible (QSH) heat fluxes. For long time scales (several decades), we assume that the net heat transport through the Strait of Gibraltar (GH) is balanced by the Mediterranean Sea surface heat fluxes (Matsoukas et al. 2005; Ruiz et al. 2008). This hypothesis implies that the water masses do not show a significant warming. Rixen et al. (2005) found a warming equal to [1.3–1.5] × 1021 J for the 1950–2000 period (that is to say a trend of [0.0018–0.0020] °C/year) that corresponds to an imbalance less than 0.4 W/m2, so our hypothesis is reasonable. The total heat budget is given by the expression:
positive sign denotes a heat gain by the sea and negative sign indicates a heat loss from the sea. We estimate the different terms in eq. (2) from the RCMs and ERA40 data. The Black Sea heat budget is not considered in our computations since its contribution can be neglected (Garrett et al. 1993). In the first stage of our analysis, the simulated Mediterranean heat budget is evaluated for the ERA40 driven RCM runs and compared with the observational datasets described in the next paragraph and summarized in Table 2.
We consider various combinations of the flux components from different datasets. The turbulent terms are taken from the OAFlux, HOAPS and NOC datasets. The observations for the radiative fluxes come from the NOC data and also from the International Satellite Cloud Climatology project (ISCCP 2) dataset (Zhang et al. 2004). The ISCCP data are based on satellite data measurements and gridded with horizontal resolution of 2.5° × 2.5°. It starts at 1984 and currently it is planned to continue through 2010.
3 Results
In this study an evaluation of different water and heat budget components obtained from the observational datasets described above has been carried out. By combining different datasets and taking into account the Gibraltar constraint, we provide three observational estimates for water budget and only one for heat budget. We also analyze the E and P estimates of eq. (1) and all the terms in eq. (2) from the ERA40 data. For model data an assessment of whether the RCMs are in closer agreement with the observations than the low resolution reanalysis (i.e. whether the RCMs have ‘added value’) is performed for multi-annual means and seasonal cycle.
Note that the time periods covered by the various observational datasets and RCM experiments are different. To evaluate the uncertainties associated with the interannual variability in the budget estimates, we provide an error bar estimated as std × tα(n−1)/√(n), where std is the standard deviation of the interannual time series, and n is the number of years. Here we consider the value of the t-student test tα(n−1) for an error corresponding to 5%. For RCMs the multi-annual ensemble mean is also computed and the spread among the model obtained as std × tα(n−1)/√(n), where std is the inter-model standard deviation and n the number of models.
3.1 The Mediterranean Sea water budget
3.1.1 Multi-annual averages
First, we study the spatial structure of the two main water budget components simulated by the RCMs. Figure 1 presents the multi-annual average of E–P for the RCM ensemble mean (Fig. 1a) and the inter-model spread (Fig. 1b) calculated as the standard deviation of the climatological annual means computed for each individual model. In Fig. 1c the multi-annual mean of E–P for ERA40 is also represented. The multi-model mean (Fig. 1a) shows the well known west-east gradient of freshwater deficit in the Mediterranean Sea (Bouktir and Barnier 2000). The largest freshwater losses occur in the north coast of Libya, the Levantine basin and in the Aegean Sea. In the Levantine basin intense evaporation is produced by high sea surface temperatures. In the Aegean Sea, evaporation is mainly driven by the local wind system (Etesian winds). There is a secondary maximum located in the Gulf of Lyon, a deep convection zone in the Mediterranean Sea, also as a result of local wind forcing (Mistral winds). Precipitation larger than evaporation (negative values) only occurs over the coastal regions, where inflowing maritime air interacts with orography to yield strong precipitation. The west-east gradient is also observed in ERA40 (Fig. 1c), whereas the finer spatial structures described above are not always present because of the coarser horizontal resolution and of the land-sea mask. It can be observed that areas of strong E–P values are missing in ERA40 (Gulf of Lyon, Aegean Sea) as well as a major part of the coastal effects. Note that for ERA40 the Adriatic and Aegean basins are not fully included in the computations of water and heat budgets.
The inter-model spread (Fig. 1b) shows a maximum over the coast, where the orography influence is strongest, and in the Aegean Sea, a land-enclosed basin where it is more difficult to accurately represent processes related to the local wind fields and orography. This suggests that future improvements could be achieved by increasing the model resolution in the Aegean basin. A further cause of the larger spread in the Aegean basin may be internal variability of the RCMs, a process recently studied by Lucas-Picher et al. (2008). The internal variability is the inter-member spread of an ensemble generated using an RCM with the same lateral boundary forcing but different initial conditions. In a recent paper, Sanchez-Gomez et al. (2008), show that the internal variability and the inter-model spread are related. They found that the inter-model spread and also internal variability are larger in the eastern part of the ENSEMBLES domain, with a maximum located over the Balkan Peninsula, since the control exerted by the lateral boundary forcing (westerly flow) decreases with distance moved eastwards. The large spread observed over the Aegean Sea in Fig. 1b may also reflect the signature of internal variability.
Now we compare estimates of the water budget and its components obtained from the observational datasets, from the RCM ensemble and from ERA40. The basin averaged annual means of evaporation, precipitation, rivers and Black Sea discharges are shown in Table 3.
The observational datasets give a consistent picture of the Mediterranean basin evaporation term with a range going from 1,095 ± 80 (OAFlux) to 1,137 ± 90 mm/year (HOAPS), the highest spatial resolution dataset giving the highest water loss by evaporation. The ERA40 mean E value is 1,167 ± 50 mm/year which is within the error of the others estimates. Note that the ERA40 correction made by Pettenuzzo et al. (2010) leads to a consistent value (1,180 mm/year for the period 1958–2001) which however strongly depends on the bulk formula used to re-compute the evaporation terms from the atmosphere and ocean variables.
The RCM simulated annual mean estimates for the evaporation range between 1,066 mm/year (MPI model) and 1,618 mm/year (ICTP model). The interannual variability for the RCMs is typically 50–60 mm/year which is close to ERA40 but smaller than the observations which have values of about 80 mm/year. The ensemble mean value is 1,254 mm/year which is 7% higher than ERA40 and stands above all the observational estimates of evaporation. 7 out of the 12 RCMs present evaporation values higher than the largest observational estimates (HOAPS). The high evaporation values found for several RCMs can be explained by stronger local winds channelled by a better representation of the Mediterranean surrounding orography. This explanation is sustained by recent studies focusing on the Mediterranean winds and air-sea fluxes (Sotillo et al. 2005; Ruti et al. 2007; Herrmann and Somot 2008; Schroeder et al. 2010). Note that the inter-model spread of 164 mm/year constitutes 13% of the ensemble mean, indicating large differences among the regional models, forced by the same lateral and surface boundary conditions. Note also that in the ERA40 driven run all the RCMs used the same SST dataset.
The annual mean precipitation amount over the Mediterranean Sea is not coherent when comparing the GPCP (594 ± 56 mm/year), CMAP (467 ± 44 mm/year) and HOAPS (256 ± 44 mm/year) datasets, the difference exceeding 100% of the lower value. GPCP is known to overestimate the precipitation over the sea as it mixes land and sea points. Moreover its low resolution is not appropriated to the Mediterranean precipitation study. CMAP lies between GPCP and HOAPS values. On the contrary HOAPS is known to underestimate the precipitation as explained in the previous section. The “true” value is perhaps around 400 mm/year which is in the range 331–447 mm/year given by Mariotti et al. (2002). ERA40 gives a value of 386 ± 80 mm/year in agreement with the observation range, and the RCM ensemble mean a value of 442 mm/year which indicates stronger precipitation on average. We can thus characterise the difference of the 25-km dynamical downscaling mean with respect to the driver as an increase in precipitation by 15%. Contrary to evaporation, the interannual variability of precipitation is generally higher in ERA40, than the RCMs and most of the RCMs provide values that are similar to observations. None of the models has a value as low as HOAPS, and only one model (OURANOS) has a value even higher than GPCP. Note that the ERA40 correction applied by Pettenuzzo et al. (2010) with respect to the CMAP dataset gives a value of 530 mm/year; lying in the observation range.
The observed annual mean runoff is 142 ± 22 mm/year (Ludwig et al. 2009), which is larger than previous river discharge estimates given by Struglia et al. (2004) and Mariotti et al. (2002) who obtained a value of 102 mm/year from other observational datasets. Very recently a new estimate (116 - 135 mm/year) has been provided by Bouraoui et al. (2010) for the 1980–2000 period using a reconstruction method. The KNMI and METNO models have similar values of 146 and 140 mm/year respectively and the ensemble mean is 124 mm/year with a spread of 46 mm/year. The interannual variability of R is greater than the observed range for most of the RCMs with the closest being CNRM.
For the Black Sea discharge, the RCM ensemble mean is 87 mm/year showing a good agreement with the observational value of 80 mm/year by Stanev et al. (2000). In this computation, we have omitted the METNO model because its value (negative) is not realistic. However this variable presents the largest discrepancies among the RCMs, with an inter-model spread of 60 mm/year. Indeed the Black Sea input is mainly estimated from the river runoff term (P + R − E) of the Black Sea catchment, which is very large and covers a relatively flat area. As mentioned in a previous paragraph, the inter-model spread and also internal variability of RCMs are larger in the eastern part of the domain, since the control exerted by ERA40 decreases with distance as we moved eastwards. This can contribute to the large inter-model spread observed for the B term.
The use of different observation sources for evaporation and precipitation leads to five different estimates for the Mediterranean Sea water budget determined from eq. (1): WB1 (OAFlux for E, GPCP for P), WB2 (OAFlux for E and HOAPS for P), WB3 (HOAPS for both E and P), WB4 (NOC for E, HOAPS for P) and WP5 (NOC for E and CMAP for P). In all cases R is given by Ludwig et al. (2009) and B by Stanev et al. (2000). The choice for P drives the difference obtained between WB1 (280 mm/year) and the four other estimates (617, 660, 612 and 426 mm/year). We decided to discard the very low WB1 because it largely underestimates the lowest observation estimate of the Gibraltar Strait net transport (GW, which is between 515 and 1,150 mm/year as discussed earlier). WB2, WB3 and WB4 are considered to be more realistic on this basis as they lie within the GW range. WP5 is not coherent with Gibraltar Strait measurements. Note that the most recently cited value in the literature is 0.05 Sv for the Gibraltar Strait net transport (Baschek et al. 2001; Garcia de la Fuente et al. 2007 for example), that corresponds to about 630 mm/year (with a Mediterranean Sea surface equal to 2.5 × 1012 m2).
Because of the physics consistency in a climate model, it is interesting to look at the E–P budget too. Indeed a model having a high evaporation over the sea also tends to have a high precipitation. Despite this effect, the ERA40 E–P mean of 781 mm/year, lies slightly below the RCM ensemble mean of 812 mm/year. The E–P values from the models are closer to those provided by the observational WB2–WB3–WB4, estimates than ERA40 (Table 3). The WB1 and WP5 estimates of E–P seems to be unrealistic mainly due to the high GPCP and CMAP precipitation values.
The ERA40 water budget estimate has not been analyzed in this paper, since ERA40 based runoff values are not realistic (Hagemann et al. 2005). For the RCMs, the annual averages of the total Mediterranean water budget can be computed for only 7 models (CNRM, DMI, ETHZ, KNMI, MPI, OURANOS, UCLM, excluding METNO due to its unphysical Black Sea value). The ensemble mean is equal to 540 mm/year, which is below the observational estimates by approximately 100 mm/year. The RCMs range from a freshwater deficit of 420 mm/year for OURANOS to 720 mm/year for DMI. These are also the models with the highest and lowest E–P balance respectively. DMI has the highest E term and OURANOS the highest P term of the remaining 7 models. The CNRM model gives the best representation of the water budget with 635 mm/year (close to the Gibraltar 0.05 Sv value and the WB2 estimate). The ETHZ, MPI, KNMI and UCLM models stand surprisingly very close to each other between 470 and 506 mm/year.
Overall, our results demonstrate that regional climate models show some deficiencies to provide realistic water budget estimates, consequently some improvements may be required for the dynamical downscaling estimates in order to decrease the large discrepancies observed between different RCMs. The Mediterranean water budget strongly depends on the choices of physics model parameterizations and this study suggests that more effort needs to be made in this regard. As regards the observational datasets, our results also reinforce the need for improving the quality of data over the sea, especially the spatial resolution. Increasing resolution is crucial to a better representation of local phenomena that in the case of the Mediterranean Sea, play an important role in the water budget such as orography-driven local winds, orography-driven precipitation over land, land-sea contrast, wind-driven intense evaporation events and Mediterranean cyclones induced precipitation.
3.1.2 Seasonal cycle
We now consider the climatological annual cycles of evaporation, precipitation, river discharge and the Black Sea input, for both the RCMs and observations, which are shown in Fig. 2. In this figure only observations used in the net water budget estimates coherent with the Gibraltar Strait constraint (WB2–WB3–WB4, see Table 3) have been considered. These values have been averaged over the entire Mediterranean basin. We observe that in general, despite some spread among the models, the seasonal cycles for the various water budget terms are well represented by the RCMs.
Climatological seasonal cycle for the components of the water budget averaged over the whole Mediterranean basin. a Evaporation from RCMs, OAFlux and HOAPS data. b Precipitation from RCMs, GPCP and HOAPS data. c River discharge from RCMs and Ludwig et al. 2009. d Black Sea input from RCMs and Stanev et al. 2000
The Mediterranean evaporation (Fig. 2a) presents a minimum in late spring (May) and is most intense in autumn (November). The amplitude of the seasonal cycle for evaporation is very large. Consequently evaporation is the leading term of the water budget seasonal cycle as will be seen in Fig. 3 (discussed below). The seasonal cycle amplitude for E (Fig. 2a) is 1,052 mm/year for OAFlux, 950 mm/year for HOAPS and 980 mm/year for NOC. These values are comparable to the ERA40 value (1,004 mm/year) and to the RCM ensemble mean (1,012 mm/year), showing a good agreement between all the data sources. The HOAPS data presents larger evaporation values during summer and autumn than OAFLUX and NOC. The observed seasonal cycle is well represented by the RCMs although the ICTP model is a noticeable outlier with extremely strong evaporation.
Climatological seasonal cycle for the area-averaged Mediterranean freshwater deficit estimates for the period 1979–2000. The RCMs estimates are represented by lines in colors. The estimate WB1 has been obtained by combining the OAFlux with the GPCP datasets. The estimate WB2 has been obtained from HOAPS data, and the WB3 from NOC for E and HOAPS for P. R is given by Ludwig et al. (2009) and B by Stanev et al. (2000) for all estimates
The annual cycle of the area-averaged precipitation over the Mediterranean Sea (Fig. 2b) is also well simulated by the RCMs. ERA40 lies above the HOAPS observed precipitation with closer values from May to September. The RCMs mainly stay also above the observational dataset. The difference between the maximum and minimum value for the precipitation is 570 mm/year for HOAPS, 716 mm/year for ERA40, and 664 mm/year for the RCMs ensemble mean. These estimates are in agreement with the values obtained by Mariotti et al. (2002). In this case, we observe fewer discrepancies among the estimates of the amplitude of the seasonal cycle than in the annual values.
For the observed river discharge term minimum values are reached in summer and maximum values in winter (Fig. 2c), following the precipitation seasonal cycle. A second maximum is reached in March–April during the snow-melting season. For the models, the river discharge term has been analyzed only for the 7 RCMs having the runoff as an output. The seasonal cycle for the basin-integrated Mediterranean river discharge is reasonably well simulated by all the models, as they capture the main timing characteristics of the observed seasonal cycle with the minimum in summer and the maximum in winter with the snow-melting peak. However, the minimum in February between the two peaks is rarely represented by the models. Note also that throughout the year the ETHZ model overestimates the observational value, yielding to a large annual mean (see Table 3).
Although the timing of the river discharge seasonal cycle extremes are in reasonable agreement across observations and models, we observe large discrepancies concerning its amplitude. We obtain for the observations an amplitude of 82 mm/year. The ensemble mean underestimates the amplitude with 50 mm/year. Our results here support the suggestion of Elguindi et al. (2009), that higher resolution than 50 km is needed to adequately simulate the river runoff for an area of complex orography such as the surrounding of the Mediterranean Sea. Up to know, only some preliminary tests have been performed with 10 km-resolution RCMs for ERA40 dynamical downscaling since longer experiments at such high resolution are still far to be implemented. This new configuration often covers sub-regions of Europe as in the European project CECILIA or in Déqué and Somot (2008a, b).
The Black Sea input to the Mediterranean Sea has been also calculated only for the 7 models having the runoff as a diagnostic variable. This term is computed as P + R − E for the Black Sea. Indeed this is equivalent to the freshwater amount that flows into the Mediterranean Sea through the Dardanelles Strait assuming that the Black Sea level is constant over a long period of time. Positive values indicate a freshwater gain for the Black Sea and therefore for the Mediterranean Sea, whereas negative values correspond to a freshwater deficit (excess of evaporation against precipitation and runoff). The annual cycle of the Black Sea input (Fig. 2d) reaches the maximum value in spring (March–April), when the river discharge over the Black Sea drainage area is most intense due to the snowmelt. For the Black Sea input in the Mediterranean Sea, a lag is expected between the maximum in freshwater and the discharge in the Aegean Sea (Stanev et al. 2000). Indeed a hydraulic control is imposed by the Bosphorus Strait—Marmara Sea—Dardanelles Strait system, which leads to a maximum flow between the Black Sea and the Aegean Sea. During spring, a sea level increase is observed in the Black Sea due to a freshwater excess. Then the freshwater is transported through the Bosphorus Strait during late spring and summer, when the water balance is negative in the Black Sea. This delay is not taken into account here.
In April, the Black Sea input from the observations is 255 mm/year. The simulated values range between 91 mm/year for METNO model, which clearly underestimates the Black Sea water budget, and 500 mm/year for OURANOS. The freshwater deficit is largest in August, when the evaporation rates are most intense, with observed values of –91 mm/year. The RCM estimates for this minimum value vary between −168 mm/year (DMI model) and −20 mm/year (ETHZ model). The amplitude of the annual cycle for the Black Sea input is 350 mm/year for observations. This amplitude is about 20% greater for the RCM ensemble mean, which has a value of 423 mm/year.
The seasonal cycle for the total water budget in the Mediterranean Sea is displayed in Fig. 3 for models and the three WB estimates that verify the closure hypothesis. Note that the time lag between the the Black Sea E–P–R budget and the Black Sea input to the Aegean Sea is neglected (see Stanev et al. 2000 for a discussion). As mentioned above, we have combined the OAFlux, the HOAPS precipitation and river discharge data to obtain WB2 estimate. WB3 has been built from HOAPS and river runoff datasets, and WB4 is based on NOC evaporation, HOAPS precipitation and the river runoff datasets. For all three observational based estimates, the freshwater loss presents a maximum in the late summer, between August and September. During these months, the evaporation rate over the basin largely exceeds the freshwater input from precipitation and river discharge. In August the freshwater deficit is approximately 1,040 for WB2, 1,400 mm/year for WB3 and 1,060 mm/year for WB4. The water balance reaches its minimum in spring, from April to May. For WB2,WB3 and WB4 the water budget is close to zero because of the minimum of evaporation and a maximum for the river and Black Sea inputs. For some models, the water budget becomes negative, indicating that the freshwater input from precipitation and especially from river discharge is more important than evaporation. During the spring season (April–May), freshwater excess in the Mediterranean Sea is on average approximately +30 mm/year for WB2, WB3 and WB4.
3.2 The Mediterranean Sea heat budget
3.2.1 Multi-annual averages
As for the water budget, we first study the spatial structure of the total heat budget simulated by the RCMs. Figure 4a shows the multi-model annual mean over the whole period together with the spread among the models (Fig. 4b). In Fig. 4c the total heat budget for ERA40 is also represented. The spatial distribution of the multi-model total heat budget (Fig. 4a) is characterized by maxima of heat loss over the deep convection zones which are the Gulf of Lyon, and the Adriatic and Aegean Seas. Heat gain occurs over the Alboran Sea primarily as a result of reduced latent heat loss in this region (due to a smaller sea-air humidity difference and weaker winds). The spatial distribution of ERA40 total heat budget (Fig. 4c) differs from the RCM ensemble mean. Differences can be found in particular over the Gulf of Lyon, where the heat loss of the models is stronger than for ERA40. This is due to the higher resolution of the RCMs that results in stronger local winds and then enhanced evaporation over the region (Sotillo et al. 2005; Ruti et al. 2007; Herrmann and Somot 2008). Some differences are also observed in the Levantine basin, where the small scale coastal phenomena are not present in ERA40 and in the southern Lybian coast.
The inter-model spread (Fig. 4b) shows maxima over the coastal regions and also over the Alboran Sea, the Gulf of Lyon, Adriatic and Aegean Seas, where the influence of complex orography-related processes and local wind systems (e.g. Samuel et al. 1999) are dominant. Note that in the case of the total heat budget the signature of internal variability is not as evident as for the E–P budget. One reason for this is that the precipitation field may be more affected by internal variability than the radiative and turbulent fluxes.
The annual long term means for the heat budget terms from the different RCMs are displayed in Table 4. In the rows of Table 4 are indicated the estimates obtained with the RCMs for the ERA40 driven experiments and also the ensemble means values. The corresponding observational estimates and ERA40 are in Table 5. The net shortwave absorbed by the sea is 187 ± 3 W/m2 for ISCCP2 and 185 ± 2 W/m2 for NOC. The NOC estimate is close to the 183 W/m2 value provided in Gilman and Garret (1994). ERA40 has a lower value of 165 ± 3 W/m2 which is likely due to an overestimation of the cloud cover. The RCMs shortwave estimates show a wide range from 154 W/m2 for DMI to 214 W/m2 for METOHC, with an RCM ensemble mean of 181 ± 18 W/m2 which is slightly lower than, but in reasonable agreement with, the ISCCP2 and NOC estimates.
The annual mean estimate of the longwave radiative flux from ISCCP2 is 76 ± 4 W/m2, while the NOC longwave is considerably larger, 84 ± 1 W/m2. Contrary to the shortwave behaviour, the ERA40 longwave estimate of 78 ± 3 W/m2 lies between the two observation based values. Except for MPI and METNO that have very low values, all the other RCMs give values close to or within the range spanned by the ISCCP2 and NOC datasets. The RCM ensemble mean is 75 ± 6 W/m2, which is close to ERA40. However individually, RCMs provide values that are quite different from the driving field, showing that the downscaling has a major impact on the longwave flux.
The latent heat flux annual mean is given either by OAFlux (HB1), NOC (HB2) or HOAPS (HB3) with an observation range from 88 to 90 W/m2 (Table 5). The ERA40 value of 93 ± 4 W/m2 is just above the upper limit of this range, while the RCMs tend to lie well above the observation range (as noted in the previous section with the evaporation term of the water budget). The RCM ensemble mean estimate for the latent heat flux is 100 ± 13 W/m2 with values going from 85 to 128 W/m2. The RCMs at the lower end of this range (CNRM, KNMI, MPI, SMHI, UCLM) are clearly in better agreement with the observations than those at the upper values, although recall that closure of the water budget implies the latent heat flux may be higher than the observation range noted here.
For the sensible heat fluxes, the OAFlux and HOAPS datasets provides an estimate of 14 ± 2 W/m2 for the averaged basin and NOC a smaller estimate of 7 ± 1 W/m2. ERA40 (9 ± 1 W/m2) lies within the observed range as the model ensemble mean (13 ± 5 W/m2). However 5 models over 12 lie out of the observed range.
Overall, if we summarize the behaviour of the models compared to observations for the four different components of the heat budget, ERA40 performs well for QLW, QLH and QSH but strongly underestimates QSW. The RCM ensemble mean generally performs well but with a potential overestimation of QLH and a significant model spread. Very few models perform well for all the components and this underlies the need for model physics improvements as regards the air-sea exchange. To conclude, the dynamical downscaling of ERA40 seems to have a major impact on all the components of the Mediterranean heat budget, since the RCMs provide values quite different from the driving reanalysis. In view of the large inter-model spread observed, the reasons of these discrepancies are not likely related with the increase of horizontal resolution. Since the boundary conditions are identical for all models, the differences are probably related to the choice physical parameterisations used for each model.
Concerning the total heat balance, we have combined the different available datasets to obtain three values (see Table 5): HB1 uses ISCCP for QSW and QLW and OAFlux for QLH and QSH; HB2 uses NOC estimates for all the fluxes and HB3 uses NOC for QSW and QLW, HOAPS for QLH and QSH. The HB1 and HB2 estimates are positive (+9, +5 W/m2), indicating a non realistic heating of the Mediterranean Sea. Recall that the Gibraltar Strait net heat transport ranges from +3 to +10 W/m2 (see the Introduction) indicating an equivalent cooling of the Mediterranean Sea through its surface. We have also combined in HB3 the lowest estimate of QSW (NOC, 185 W/m2) and the highest for QLW (NOC, 84 W/m2), with the high HOAPS estimates for the turbulent fluxes (90 and 14 W/m2). This leads to a more realistic negative heat budget equal to −3 ± 8 W/m2, for which the closure hypothesis is satisfied within the stated error bounds, although the mean value of −3 W/m2 is greater than the −5 W/m2 found by McDonald et al. (1994) on the basis of Gibraltar Strait heat transport measurements.
As a result of its low QSW value, the ERA40 net heat flux estimate of −15 W/m2 is outside the range implied by the net heat transport at the Gibraltar Strait. Regarding the RCMs, we find a large uncertainty around the heat budget estimates. The values vary between −40 W/m2 for DMI and +21 W/m2 for METOHC, with 5 models simulating a heat gain by the Mediterranean Sea and 7 a heat loss. This result shows that individual RCMs have difficulties in simulating the net Mediterranean Sea heat budget. Moreover, the RCM ensemble mean of −9 W/m2 underestimated by 4 W/m2 the McDonald et al. (1994) value. In comparison, Pettenuzzo et al. (2010) obtain basin mean values between −5 and 4 W/m2 from their adjusted ERA40 fields depending on the choice of transfer coefficient scheme employed (see Sect. 1).
The model results show that, although closure is obtained with the ensemble mean, the individual RCMs are often not balanced in terms of the surface heat fluxes over the Mediterranean Sea. This is an important issue for the regional modelling community and suggests that more efforts are needed to improve the heat budget estimates. This is also an important message for the regional ocean modelling community, which is starting to evaluate the skill of dynamical downscaling of reanalysis over the Mediterranean Sea (Sotillo et al. 2005; Ruti et al. 2007), and to use them to force Mediterranean regional ocean models for climate-scale simulations (Herrmann and Somot 2008; Somot and Colin 2008; Sevault et al. 2009; Tsimplis et al. 2008; Beuvier et al. 2010; Herrmann et al. 2010). Long-term well balanced water and heat budgets are required skills of atmosphere forcing.
3.2.2 Seasonal cycle
Figure 5 shows the climatological annual cycle for the radiative fluxes (short and longwave) and for the turbulent fluxes (latent and sensible) averaged over the whole Mediterranean basin. For these figures only observations used in HB3 are considered (see Table 5). The factor with the largest contribution to the heat budget is the net QSW radiation (Fig. 5.a), which has a pronounced seasonal cycle. The net shortwave flux is in general well simulated by the models, though a larger inter-model spread is observed during the summer months, when the short wave absorbed by the sea is maximum, with values ranging from 240 to 320 W/m2 for July–August. The ERA40 estimates lie below the NOC observational dataset during the whole year, with a maximum underestimation during the summer months.
Climatological seasonal cycle for the components of the heat budget averaged over the whole Mediterranean basin. a Shortwave from the RCMs, ISCCP and NOC data. b Longwave from the RCMs, ISCCP and NOC data. c Latent heat flux from the RCMs, OAFlux, HOAPS and NOC data. d Sensible heat flux from the RCMs, OAFLUX, HOAPS and NOC data
The net longwave radiation (Fig. 5b) emitted by the sea has a seasonality much less pronounced than QSW (note the major difference in scales between Figs. 5a and 5b). These values indicate a heat loss from the sea to the atmosphere. According to this, the heat loss for NOC is especially large in autumn, winter and spring. There is a minimum of QLW during late spring (May). Most of the RCMs also show this minimum in May, though the inter-model spread is large for QLW. In May, QLW is 78 W/m2 for NOC. ERA40 is very close to NOC with 79 W/m2 and the RCM estimates vary between 60 W/m2 and 76 W/m2. QLW is maximum in summer, when cloud cover is at a minimum. For August the value is 80 W/m2 for NOC. ERA40 with 82 W/m2 is slightly higher than NOC and the RCM values range from 64 to 83 W/m2. In general, differences between the observational dataset and models for the longwave are less important than those for shortwave given the smaller magnitude of the longwave flux. Note the discrepancy between the NOC and ERA40 seasonal cycles, ERA40 showing a very flat seasonal cycle.
The annual cycle of latent heat flux will be only briefly discussed in this section, since its behaviour is the same as for evaporation in Sect. 3.1. QLH always shows a heat loss from the sea by evaporation in Fig. 5c. The annual cycle has a minimum in May and reaches a maximum in late autumn. In November during the maximum we obtain estimates 117 W/m2 for HOAPS and 131 W/m2 for ERA40. A range of 110–160 W/m2 is observed for the models at this time. In late spring, QLH drops to approximately 40 W/m2 for HOAPS and ERA40 and 55–77 W/m2 for the RCMs. The RCMs have a tendency to overestimate the latent heat flux compared to observations in all seasons. Note also that a subset of the models reproduces the high June–September anomaly values obtained from HOAPS. However OAFLUX and NOCS do not show this June–September anomaly values (not shown).
The sensible heat flux has the smallest contribution to the heat budget (again note the major difference in scales between Figs. 5c and 5d). Nevertheless it shows a pronounced seasonal cycle with a minimum in late spring and summer and a maximum from January to December (Fig. 5d). The values of QSH are in general an ocean heat loss, however during the minimum, the values of QSH for some models indicate a heat gain for the ocean, which is colder than the atmosphere. Again HOAPS shows an anomaly of the seasonal cycle for the June–September period compared to the other observational datasets (not shown). ERA40 lye below the observed estimate and the annual variations in QSH are reasonably well captured by the RCMs with ICTP an obvious outlier. None of the models shows the HOAPS anomaly. The amplitude of the seasonal cycle is on average 26 W/m2 for the RCM ensemble mean, compared with approximately 18 W/m2 for HOAPS and 25 W/m2 for ERA40.
The total heat budget seasonal cycle, as computed following eq. (2) is represented in Fig. 6. The annual mean values for HB3 have been discussed in the previous section; we recall here that only HB3 has a Mediterranean Sea yearly-mean negative heat loss, and hence is more consistent with the Gibraltar net transport. Considering the seasonal cycle, the heat loss from the Mediterranean Sea is maximum in winter; with values of 120 W/m2 for HB3. This is a result of weak shortwave and strong turbulent fluxes due to intense wind events and the cold and dry atmosphere. This loss is compensated by a heat gain of 120 W/m2 during summer, with a peak reached in June for both observations and models. In general, the models reproduce the seasonal cycle seen in the observations reasonably well with no major anomalies.
Climatological seasonal cycle for the area-averaged Mediterranean heat budget for the period 1983–2000. The estimate HB1 has been obtained by combining the OAFlux for QLH and QSH with the ISCCP for QSW and QLW. The estimate HB2 has been obtained entirely from NOC data. And HB3 is obtained by merging HOAPS for QLH and QSH, and NOC for QSW and QLW
4 Summary and discussion
This main goal of this work has been to assess the performance of an ensemble of regional models in terms of heat and water budgets in the Mediterranean Sea. The RCMs have been driven by the ERA40 reanalysis for the 1960–2000 period and with a spatial resolution of 25 km. These high resolution model datasets have been compared to different lower resolution observational datasets and also to the ERA40 reanalysis. We have used a range of recent observational datasets over the Mediterranean Sea to evaluate the RCMs. The focus has been on the area-averaged climatological seasonal cycle and the long-term annual mean. The main conclusions can be summarized as follows:
There are still large uncertainties concerning the observations of the various terms of the heat and water budget of the Mediterranean Sea. Improvements of the observational estimates over the sea, especially for the latent heat flux, precipitation and radiative fluxes, are required to provide tighter constraints on the models. Precipitation is the most problematic term because of the large difference between the GPCP, CMAP and HOAPS datasets.
Difficulties in obtaining agreement between the estimates of the water and heat budgets and transport based values obtained from the Gibraltar Strait suggest that the latent heat flux may have been underestimated in all of the observational datasets. For example, the use of the HOAPS low estimate of the precipitation over sea is required to obtain a realistic total freshwater loss (WB2, WB3, WB4). For the heat budget, a combination of NOC for shortwave and longwave and HOAPS for sensible heat and latent heat (HB3) is required to obtain a net heat loss of –3 W/m2 which is close to the value of –5 W/m2 implied by measurements of the Gibraltar heat transport. Other combinations of the data sources give unrealistically positive values for the net heat flux. The estimates obtained in the current study are: 1,095–1,115 mm/year for E, 256 mm/year for P, 142 mm/year for R, 80 mm/year for B, 612–660 mm/year for the E–P–R–B total water budget and 185 W/m2 for SW, −84 W/m2 for LW, 88–90 W/m2 for LH, 14 W/m2 for SH and −3 W/m2 for the total heat budget.
The ERA40 reanalysis latent heat flux or evaporation is slightly higher than the observational values and the net short wave radiation is significantly lower than the observations.
The ERA40-driven RCM ensemble is typically in good agreement with the observations of the shortwave, longwave, sensible heat, evaporation-precipitation, river runoff and Black Sea inputs. However, the RCMs overestimate the latent heat flux (or evaporation) and precipitation with respect to the observations (although as noted above the observation based latent heat flux may be an underestimate). For the water budget, good river runoff and Black Sea discharge estimates are obtained with the ensemble of RCMs. This change is likely to reflect the high spatial resolution which enables a more realistic representation of the complex orography (Elguindi et al. 2009). The long-term seasonal cycle of the water budget components (evaporation, precipitation and river discharge) is well simulated by the RCMs.
For the heat budget terms, and in view of our results, giving an overall conclusion about the added value of the dynamical downscaling is difficult. We note that the RCM ensemble mean value of −9 W/m2 for the net heat flux is closer to the basin-wide average of −5 W/m2 expected from the Gibraltar transport measurements than the original ERA40 value of −15 W/m2. Though the ensemble mean estimate is correct, the ensemble spread is large (+21 W/m2), indicating the discrepancies of the state-of-the-art RCMs. We also note that previous studies have shown that dynamical downscaling has an added value with respect the coarse horizontal scale driving field in local and extreme phenomena (winds, air-sea exchanges) as shown in Sotillo et al. (2005), Ruti et al. (2007) and Herrmann and Somot (2008), but this is not the case for this work
There is a tendency for a large spread between RCMs in the terms of both the water and heat budgets. For example, the mean spread among the models is 164 mm/year for evaporation. The simulated precipitation also has a large inter-model spread of approximately 84 mm/year. Note that this spread is even stronger in the observational precipitation datasets between GPCP and HOAPS which demonstrates the difficulty in using them to constrain the freshwater budget. The inter-model spread is also large for the river discharge and the Black Sea input. As mentioned above, the uncertainties associated with the heat budget estimates for individual RCMs remain large with basin mean values from –40 W/m2 to +21 W/m2. The extreme high and low values are likely due to the use of inappropriate physical parameterization of turbulent and radiative heat fluxes in some cases in the individual RCMs, or perhaps to the use of non coupled models.
To conclude, our results show that only a few of the RCMs considered provide balanced heat and water budgets when compared with observations, and also not significants improvements have been obtained with respect to the original ERA40 forcing. Previous studies have shown that the dynamical downscaling in the Mediterranean region improve the representation of air-sea exchanges (Sotillo et al. 2005; Ruti et al. 2007; Ruiz et al. 2008; Herrmann and Somot 2008) in several zones (Gulf of Lyon, Adriatic Sea). However we show here that for a basin wide average, the dynamical downscaling does not always yield to significant improvements of the water and heat budget estimates respect to ERA40.
From the information provided in this paper, it is possible to select an RCM that is suitable on budget grounds for use as a high resolution atmospheric forcing for an ocean model of the Mediterranean Sea. High resolution for the atmospheric forcing is vital to correctly represent the deep water formation and the Mediterranean thermohaline circulation and our new results indicate that more progress can be made in this area in the near future.
Acknowledgements: The financial support of this work has been provided by the European Union Sixth Framework programme project ENSEMBLES, under contract GOCE-CT-2006-037005. The authors wish to thank Alain Braun (CNRM), A. Farda and P. Stepanek (CHMI), O·B. Christensen (DMI), C. Schär (ETHZ), B. Rockel (GKSS), E. Buonomo (METO-HC), F. Giorgi (ICTP), E. van Meijgaard (KNMI), J.E. Haugen (METNO), D. Jacob (MPI), M. Rummukainen (SMHI), E. Sanchez (UCLM) and D. Paquin (OURANOS) for performing the RCMs experiments and providing the model outputs. We would like to acknowledge J.H. Christensen and M. Rummukainen for coordinating the RT3 activities and O.B. Christensen (DMI) for preparing and maintaining the website from the RT3 database. Special thanks are given to Romain Bourdalle-Badie (MERCATOR) for the discussions about the Mediterranean Sea circulation.
References
Adler RF, Huffman GJ, Chang A, Ferraro R, Xie P, Janowiak J, Rudolf B, Schneider U, Curtis S, Bolvin D, Gruber A, Susskind J, Arkin P (2003) The version 2 global precipitation climatology project (GPCP) monthly precipitation analysis (1979-present). J Hydrometeor 4:1147–1167
Andersson A, Bakan S, Fennig K, Grassl H, Klepp CP, Schulz J (2007) Hamburg ocean atmosphere parameters and fluxes from satellite data—HOAPS-3—monthly mean. World Data Center for Climate. doi:10.1594/WDCC/HOAPS3_MONTHLY
Artale V, Calmanti S, Malanotte-Rizzoli P, Pisacane G, Rupolo W, Tsimplis M (2005) Mediterranean climate variability. In: Lionello P, Malanotte-Rizzoli P, Boscolo R (eds) Elsevier, pp 282–323
Baschek B, Send U, Garcia de la Fuente J, Candela J (2001) Transport estimates in the Strait of Gibraltar with a tidal inverse model. J Geophys Res 112:31033–31044
Béthoux J (1979) Budgets of the Mediterranean Sea. Their dependence on the local climate and on the characteristics of the Atlantic waters. Oceanol Acta 2:157–163
Béthoux J, Gentili B, Taillez D (1999) Warming and freshwater budget change in the Mediterranean since the 1940s, their possible relation to the greenhouse effect. Geophys Res Lett 25:1023–1026
Betts AK, Ball JH, Viterbo P (2003) Evaluation of the ERA40 surface water budget and surface temperature for the Mackenzie river basin. ECMWF ERA40 Project Report Series, http://www.ecmwf.int/publications/
Beuvier J, Sevault F, Herrmann M, Kontoyiannis H, Ludwig W, Rixen M, Stanev E, Béranger K, Somot S (2010) Modelling the Mediterranean Sea interannual variability over the last 40 years: focus on the EMT. J Geophys Res (Ocean) 115. doi:10.1029/2009JC005950
Bhöm U, Küchen M, Ahrens W, Block A, Hauffe D, Keuler K, Rockel B, Will A (2006) CLM-the climate version of LM: brief description and long-term applications. COSMO Newslett n6
Bignami F, Marullo S, Santoleri R, Schiano ME (1995) Longwave radiation budget in the Mediterranean Sea. J Geophys Res 100:2501–2514
Bouktir M, Barnier B (2000) Seasonal and Interannual variations in the surface freshwater flux in the Mediterranean Sea from ECMWF reanalysis project. J Mar Sys 24:343–354
Bouraoui F, Grizzetti B, Aloe A. (2010) Estimation of water fluxes into the Mediterranean Sea. J Geophys Res. doi:10.1029/2009JD013451
Brankart JM, Pinardi N (2000) Abrupt cooling of the Mediterranean levantine intermediate water at the beginning of the 1980s: observational evidence and model simulation. J Phys Oceanogr 31:2307–2320
Bryden HL, Kinder TH (1991) Steady two-layer exchange through the strait of Gibraltar. Deep Sea Res 38:445–463
Bryden HL, Candela J, Kinder TH (1994) Exchange through the strait of Gibraltar. Prog Oceanogr 33:201–248
Bunker AF, Charnock H, Goldsmith RA (1982) A note of the heat balance of the Mediterranean and Red Seas. J Mar Res 40:73–84
Candela J (2001) The Mediterranean water and the Global circulation, observing and modelling the global ocean. In: Sielder G, Church J, Gould G (eds) Academic, San Diego, pp 419–429
Castellari S, Pinardi N, Leaman K (1998) A model study of air-sea interactions in the Mediterranean Sea. J Mar Sys 18:89–114
Castellari S, Pinardi N, Leaman K (2000) Simulation of water mass formation processes in the Mediterranean Sea: influence of the time frequency of the atmospheric forcing. J Geophys Res 105:24157–24181
Christensen JH, Christensen OB, Lopez P, van Meijgaard E, Botzet M (1996) The HIRHAM4 regional atmospheric climate model. Scientific Report DMI, Copenhagen 96–4, 51
Collins M, Booth BBB, Harris GR, Murphy JM, Sexton DMH, Webb MJ (2006) Towards quantifying uncertainty in transient climate change. Clim Dyn 27:127–147. doi:10.1007/s00382-006-0121-0
Déqué M, Somot S (2008a) Added value of high resolution for ALADIN Regional Climate Model. Research activities in atmospheric and oceanic modelling. CAS/JSC Working group on numerical experimentation. Report No. 38. http://collaboration.cmc.ec.gc.ca/science/wgne/index.html
Déqué M, Somot S (2008b) Extreme precipitation and high resolution with Aladin. Idöjaras Quaterly Journal of the Hungarian Meteorological Service 112:179–190
Elguindi N, Somot S, Déqué M, Ludwig W (2009) Climate change evolution of the hydrological balance of the Mediterranean, Black and Caspian Seas: impact of climate model resolution. Clim Dyn 36:205–228. doi:10.1007/s00382-009-0715-4
Fairall CW, Bradley EF, Hare JE, Grachev AA, Edson JB (2003) Bulk parameterization of air‐sea fluxes: updates and verification for the COARE algorithm. J Clim 16:571–591. doi:10.1175/1520-0442
Garcia de la Fuente J, Sanchez-Roman A, Diaz del Rio G, Sannino G, Sanchez Garrido JC (2007) Recent observations of seasonal variability in the Mediterranean outflow in the strait of Gibraltar. J Geophys Res 112:1–11
Garrett C, Outerbridge R, Thompson K (1993) Interannual variability in Mediterranean heat and buoyancy fluxes. J Climate 6:900–910
Gilman C, Garret C (1994) Heat flux parameterizations for the Mediterranean Sea: the role of atmospheric aerosols and constraints from the water budget. J Geophys Res 99:5119–5134
Gimeno L, Drumond A, Nieto R, Trigo R, Stohl A (2010) On the origin of continental precipitation. Geophys Res Lett. doi:10.1029/2010GL043712
Giorgi F, Mearns LO (1999) Introduction to special section: regional climate modelling revisited. J Geophys Res (Atmos) 104:6335–6352. doi:10.1029/98JD02072
Hagemann S, Arpe K, Bengtsson L (2005) Validation of the hydrological cycle of ERA40, ECMWF ERA40 Project Report Series, http://www.ecmwf.int/publications/
Haugen JE, Haakensatd H (2006) Validation of HIRHAM version 2 with 50 km and 25 km resolution. RegClim General Technical Report, No. 9, pp 159–173
Herrmann M, Somot S (2008) Relevance of ERA40 dynamical downscaling for modelling deep convection in the North-Western Mediterranean Sea. Geophys Res Let. doi:10.1029/2007GL032442
Herrmann M, Sevault F, Beuvier J, Somot S (2010) What induced the exceptional 2005 convection event in the North-western Mediterranean basin ? Answers from a modelling study. J Geophys Res (Oceans). doi:10.1029/2010JC006162
Hewitt CD and Griggs DJ (2004) Ensembles-based predictions of climate changes and their impacts. EOS, 85, 566 pp
Jacob D (2001) A note to the simulation of the annual and inter-annual variability of the water budget over the Baltic Sea drainage basin. Meteorology Atmospheric Physics 77:61–73
Josey SA (2003) Changes in the heat and freshwater forcing on the Eastern Mediterranean and their influence on deep water formation J Geophys Res. doi:10.1029/2003JC001778
Josey SA, Oakley D, Pascal RW (1997) On estimating the atmospheric longwave flux at the ocean surface from ship meteorological reports. J Geophys Res 102:27961–27972
Josey SA, Kent E, Taylor P (1999) New insights into the ocean heat budget closure problem from analysis of the SOC air-sea flux climatology. J Clim 12:2856–2880
Josey SA, Pascal RW, Taylor PK, Yelland MJ (2003) A new formula for determining the atmospheric longwave flux at the ocean surface at mid-high latitudes. J Geophys Res. doi:10.1029/2002JC001418
Josey SA, Somot S, Tsimplis M (2011) Impacts of atmospheric modes of variability on Mediterranean Sea surface heat exchange. J Geophys Res (Oceans), in press
Kalnay E, Kanamitsu M, Kistler R, Collins W, Deaven D, Gandin L, Iredell M, Saha S, White G, Woolen J, Zhu Y, Leetmaa A, Reynolds R (1996) The NCEP/NCAR 40 years reanalysis project. Bull Amer Meteor Soc 77:437–471
Kara AB, Hurlburt HE, Wallcraft AJ (2005) Stability‐dependent exchange coefficients for air-sea fluxes. J Atmos Oceanic Technol 22:1080–1094. doi:10.1175/JTECH1747.1
Kjellström E, Bärring L, Gollvik S, Hansson U, Jones C, Samuelsson P, Rummukainen M, Ullersig A, Willen U, Wyser K (2005) A 140-year simulation of European climate with the new version of the Rossby Centre regional atmospheric climate model (RCA3). Reports Meteorology and Climatology, 108, SMHI, SE-60176 Norrköping, Sweden, 54 pp
Kondo J (1975) Air‐sea bulk transfer coefficients in diabatic condition. Boundary Layer Meteorol 9:91–112. doi:10.1007/BF00232256
Lebeaupin C, Ducrocq V, Giordani H (2006) Sensitivity of Mediterranean torrential rain events to the sea surface temperature based on high-resolution numerical forecasts. J Geophys Res. doi:10.1029/2005JD006541
Lenderink G, van der Hurk B, van Meijgaard E, van Ulden A, Cujipers H (2003) Simulation of present day climate in RACMO2: first results and model developments. Technical report no 252, KNMI, 24 pp
Li L (2006) Atmospheric GCM response to an idealized anomaly of the Mediterranean Sea surface temperature. Clim Dyn. doi:10.1007/s00382-006-0152-6
Lucas-Picher P, Caya D, de Elia R, Laprise R (2008) Investigation of regional climate models’ internal variability with a ten-member ensemble of ten-year simulations over a large domain. Clim Dyn. doi:10.1007/s00382-008-0384-8
Ludwig W, Dumont E, Meybeck M, Heussner S (2009) River discharges of water and nutrients to the Mediterranean Sea: major drivers for ecosystem changes during past and future decades? Prog Oceanogr 80:199–217
Madec G, Chartier M, Crépon M (1991) The effect of the thermohaline forcing variability on deep water formation in the western Mediterranean Sea: a high-resolution three-dimensional numerical study. Dyn Atm Oceans 15:301–332
Mariotti A, Struglia MV, Zeng N, Lau KM (2002) The hydrological cycle in the Mediterranean Region and implications for the water budget of the Mediterranean Sea. J Climate 15:1674–1690
Mariotti A, Zeng N, Yoon J, Artale V, Navarra A, Alpert P, Li L (2008) Mediterranean water cycle changes: transition to drier 21st century conditions in observations and CMIP3 simulations. Env Res Lett. doi:10.1088/1748-9326/3/044001
Matsoukas C, Banks AC, Hatzianastassiou N, Pavlakis KG, Hatzidimitriu D, Drakakis E, Stackhouse PW, Vardavas I (2005) Seasonal heat budget of the Mediterranean Sea. J Geophys Res. doi:10.1029/2004JC002566
McDonald AM, Candela J, Bryden HL (1994) An estimate of the net heat transport flux trough the strait of Gibraltar, seasonal and interannual variability of the Western Mediterranean Sea, Coastal Estuarine Study. In: La Violette (ed) 46, AGU, Washigton DC, pp 12–32
Millot C, Candela J, Fuda JL, Tber Y (2006) Large warming and salinification of the Mediterranean outflow due to changes in its composition. Deep-Sea Res 53:656–666
Pettenuzzo D, Large WG, Pinardi N (2010) On the corrections of ERA40 surface flux products consistent with they Mediterranean heat and water budgets and the connection between basin surface total heat flux and NAO. J Geophys Res. doi:10.1029/2009JC005631
Plummer D, Caya D, Coté H, Frigon A, Biner S, Giguère M, Paquin D, Harvey R, de Elia R (2006) Climate and climate change over North America as simulated by the Canadian Regional Climate Model. J Climate 19:3112–3132
Potter R, Lozier S (2004) On the warming and salinification of the Mediterranean outflow waters in the North Atlantic. Geophys Res Lett. doi:10.1029/2003GL018161
Radu R, Déqué M, Somot S (2008) Spectral nudging in a spectral regional climate model. Tellus A 60:898–910
Rixen M, et al (2005) The Western Mediterranean deep water: a proxy for climate change. Geophys Res Lett. doi:10.1029/2005GL022702
Ruiz S, Gomis D, Sotillo MG, Josey SA (2008) Characterization of surface heat fluxes in the Mediterranean Sea from a 44-year high-resolution atmospheric dataset. Glob Plan Chan 63:258–274
Ruti P, Marullo S, D’Ortenzio F, Tremant M (2007) Comparison of analyzed and measured winds speeds in the perspective of oceanic simulations over the Mediterranean basin: Analyses of QuickScat and buoy data. J Mar Sys. doi:10.1016/2007.02.026
Samuel SL, Haines K, Josey SA, Myers PG (1999) Response of the Mediterranean Sea thermohaline circulation to observed changes in the winter wind stress field in the period 1980–1993. J Geophys Res 104:7771–7784
Sanchez E, Gallardo C, Gaertner MA, Arribas A, Castro M (2004) Future climate extreme events in the Mediterranean simulated by a regional climate model: a first approach. Glob Plan Chan 44:163–180
Sanchez-Gomez E, Somot S, Deque M (2008) Ability of an ensemble of regional climate models to reproduce weather regimes over Europe-Atlantic during the period 1961–2000. Clim Dyn. doi:10.1007/s00382-008-0502-7
Sanchez-Gomez E, Somot S, Mariotti A (2009) Future changes in the Mediterranean water budget projected by an ensemble of Regional Climate Models. Geophys Res Lett. doi:10.1029/2009GL040120
Schroeder K, Josey SA, Herrmann M, Grignon L, Gasparini GP, Bryden HL (2010) Abrupt warming and salting of the Western Mediterranean Deep Water: atmospheric forcings and lateral advection. J Geophys Res (Oceans), in press
Sevault F, Somot S, Beuvier J (2009) A regional version of the NEMO ocean engine on the Mediterranean Sea: NEMOMED8 user’s guide. Note de centre no 107, Groupe de Météorologie de Grande Echelle et Climat, CNRM
Skliris N, Lascaratos A (2004) Impacts of the Nile River damming on the thermohaline circulation and water mass characteristics of the Mediterranean Sea. J Mar Sys. doi:10.1016/j.jmarsys.2004.02.005
Somot S, Colin J (2008) First step towards a multi-decadal high resolution Mediterranean Sea reanalysis using dynamical downscaling of ERA40. Research activities in atmospheric and oceanic modelling. CAS/JSC Working Group on numerical experimentation, Report No. 38, http://collaboration.cmc.ec.gm.ca/science/wgne/index.html
Somot S, Sevault F, Déqué M (2006) Transient climate change scenario simulation of the Mediterranean Sea for the 21st century using a high-resolution ocean circulation model. Clim Dyn. 27:851–879. doi:10.1007/s00382-006-0167-z
Somot S, Sevault F, Déqué M, Crépon M (2008) 21st century climate change scenario for the Mediterranean using a coupled atmosphere-ocean regional climate model. Glob Plan Chan 63:112–126. doi:10.1016/j.gloplacha.2007.10.003
Sotillo MG, Ratsimandresy AW, Carretero JC, Bentamy A, Valero F, González-Rouco F (2005) A high-resolution 44-year atmospheric hindcast for the Mediterranean Basin: contribution to the regional improvement of global reanalysis. Clim Dyn. doi:10.1007/s00382-005-0030-7
Stanev EV, Le Traon PY, Peneve EL (2000) Sea level variations and their dependency on meteorological and hydrological forcing: analysis of altimeter and surface data for the Black Sea. J Geophys Res 105:17203–17216
Struglia MV, Mariotti A, Filograsso A (2004) River discharge into the Mediterranean Sea: climatology and aspects of the observed variability. J Clim 17:4740–4751
Timplis M, Zervakis V, Josey SA, Peneva EL, Struglia MV, Stanev E, Teocharis A, Lionello P, Malanotte-Rizzoli P, Artale V, Tragou E, Oguz T (2006) Changes in the oceanography of the Mediterranean Sea and their link to climate variability. In: Lionello P, Malanotte-Rizzoli P, Boscolo R (eds) Elsevier, pp 226–281
Tsimplis MN, Bryden HL (2000) Estimation of the transport through the strait of Gibraltar. Deep Sea Res Part I 47:2219–2242
Tsimplis M, Marcos M, Somot S, Barnier B (2008) The forcing of the sea level in the Mediterranean Sea between 1960–2000. Glob Plan Chan 65:325–332. doi:10.1016/j.gloplacha.2008.07.004
Uppala S, Kallberg P, Hernandez A, Saarinen S, Fiorino M, Li X, Onogi K, Sokka N, Andrae U, da Costa Bechtold V (2004) ERA-40: ECMWF 45-year reanalysis of the global atmosphere and surface conditions 1957–2002. ECMWF Newslett 101:2–21
Woodruff SD, Diaz HF, Elms JD, Worley SJ (1998) COADS release 2 and metadata enhancements of marine surface flux fields. Phys Chem Earth 23:517–526
Xie P, Arkin PA (1997) Global precipitation: A 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs. Bull Amer Meteor Soc 78:2539–2558
Yu L, Weller RA (2007) Objectively analyzed air-sea heat fluxes (OAFlux) for the global oceans. Bull Ameri Meteor Soc 88:527–539
Zhang Y, Rossow WB, Lacis AA, Oinas V, Mischenki MI (2004) Calculation of radiative fluxes from the surface to top of atmosphere based on ISCCP and other global data sets: refinements of the radiative transfer model and the input data. J Geophys Res. doi:10.1029/2003JD004457
Open Access
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This is an open access article distributed under the terms of the Creative Commons Attribution Noncommercial License (https://creativecommons.org/licenses/by-nc/2.0), which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
About this article
Cite this article
Sanchez-Gomez, E., Somot, S., Josey, S.A. et al. Evaluation of Mediterranean Sea water and heat budgets simulated by an ensemble of high resolution regional climate models. Clim Dyn 37, 2067–2086 (2011). https://doi.org/10.1007/s00382-011-1012-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00382-011-1012-6