Search

We present the very-high-resolution (VHR) version of the EC-Earth global climate model, EC-Earth3P-VHR, developed for HighResMIP. The model features an atmospheric resolution of ∼16 km and an oceanic resolution of 1/12° (∼8 km), which makes it one of the finest combined resolutions ever used to complete historical and scenario-like CMIP6 simulations. To evaluate the influence of numerical resolution on the simulated climate, EC-Earth3P-VHR is compared with two configurations of the same model at lower resolution: the ∼100 km grid EC-Earth3P-LR (LR) and the ∼25 km grid EC-Earth3P-HR (HR). Of the three configurations, VHR shows the smallest drift in the global mean ocean temperature and salinity at the end of a 100-year 1950s control simulation, which points to a faster equilibrating phase than in LR and HR. In terms of model biases, we compare the historical simulations against observations over the period 1980–2014. In contrast to LR and HR, VHR shows a reduced equatorial Pacific cold tongue bias, an improved Gulf Stream representation with a reduced coastal warm bias and a reduced subpolar North Atlantic cold bias, and more realistic orographic precipitation over mountain ranges. By contrast, VHR shows a larger warm bias and overly low sea ice extent over the Southern Ocean. Such biases in surface temperature have an impact on the atmospheric circulation aloft, connected with a more realistic storm track over the North Atlantic yet a less realistic storm track over the Southern Ocean compared to the lower-resolution model versions. Other biases persist or worsen with increased resolution from LR to VHR, such as the warm bias over the tropical upwelling region and the associated cloud cover underestimation, a precipitation excess over the tropical South Atlantic and North Pacific, and overly thick sea ice and an excess in oceanic mixing in the Arctic. VHR shows improved air–sea coupling over the tropical region, although it tends to overestimate the oceanic influence on the atmospheric variability at midlatitudes compared to observations and LR and HR. Together, these results highlight the potential for improved simulated climate in key regions, such as the Gulf Stream and the Equator, when the atmospheric and oceanic resolutions are finer than 25 km in both the ocean and atmosphere. Thanks to its unprecedented resolution, EC-Earth3P-VHR offers a new opportunity to study climate variability and change of such areas on regional and local spatial scales, in line with regional climate models.

The relative importance of river solid discharge, deposited sediment remobilisation, and atmospheric dust as sources of neodymium (Nd) to the ocean is the subject of ongoing debate, the magnitudes of these fluxes being associated with a significant uncertainty. The Mediterranean basin is a specific basin; it receives a vast amount of emissions from different sources and is surrounded by continental margins, with a significant input of dust as compared to the global ocean. Furthermore, it is largely impacted by the Atlantic water inflow via the Strait of Gibraltar. Here, we present the first simulation of dissolved Nd concentration ([Nd]) and Nd isotopic composition (εNd) using a high-resolution regional model (NEMO/MED12/PISCES) with an explicit representation of all Nd inputs, and the internal cycle, i.e. the interactions between the particulate and dissolved phases. The high resolution of the oceanic model (at 1/12∘), essential to the simulation of a realistic Mediterranean circulation in present-day conditions, gives a unique opportunity to better apprehend the processes governing the Nd distribution in the marine environment. The model succeeds in simulating the main features of εNd and produces a realistic distribution of [Nd] in the Mediterranean Sea. We estimated the boundary exchange (BE, which represents the transfer of elements from the margin to the sea and their removal by scavenging) flux at 89.43 × 106 g(Nd) yr−1, representing ∼84.4 % of the total external Nd source to the Mediterranean basin. The river discharge provided 3.66 × 106 g(Nd) yr−1, or 3.5 % of the total Nd flow into the Mediterranean. The flux of Nd from partially dissolved atmospheric dusts was estimated at 5.2 × 106 g(Nd) yr−1, representing 5 % of the total Nd input, and 7.62 × 106 g(Nd) yr−1 comes from the Atlantic across the Strait of Gibraltar, i.e. 7.1 % of the total Nd input. The total quantity of Nd in the Mediterranean Sea was estimated to 7.28 × 109 g(Nd); this leads to a new calculated Nd residence time of ∼68 year. This work highlights that the impact of river discharge on [Nd] is localised near the catchments of the main rivers. In contrast, the atmospheric dust input has a basin-wide influence, correcting for a too-radiogenic εNd when only the BE input is considered and improving the agreement of simulated dissolved Nd concentration with field data. This work also suggests that εNd is sensitive to the spatial distribution of Nd in the atmospheric dust, and that the parameterisation of the vertical cycling (scavenging/remineralisation) considerably constrains the ability of the model to simulate the vertical profile of εNd.

We present the very-high-resolution (VHR) version of the EC-Earth global climate model, EC-Earth3P-VHR, developed for HighResMIP. The model features an atmospheric resolution of ∼16 km and an oceanic resolution of 1/12 ° (∼8 km), which makes it one of the finest combined resolutions ever used to complete historical and scenario-like CMIP6 simulations. To evaluate the influence of numerical resolution on the simulated climate, EC-Earth3P-VHR is compared with two configurations of the same model at lower resolution: the ∼100 km grid EC-Earth3P-LR (LR) and the ∼25 km grid EC-Earth3P-HR (HR). Of the three configurations, VHR shows the smallest drift in the global mean ocean temperature and salinity at the end of a 100-year 1950s control simulation, which points to a faster equilibrating phase than in LR and HR. In terms of model biases, we compare the historical simulations against observations over the period 1980–2014. In contrast to LR and HR, VHR shows a reduced equatorial Pacific cold tongue bias, an improved Gulf Stream representation with a reduced coastal warm bias and a reduced subpolar North Atlantic cold bias, and more realistic orographic precipitation over mountain ranges. By contrast, VHR shows a larger warm bias and overly low sea ice extent over the Southern Ocean. Such biases in surface temperature have an impact on the atmospheric circulation aloft, connected with a more realistic storm track over the North Atlantic yet a less realistic storm track over the Southern Ocean compared to the lower-resolution model versions. Other biases persist or worsen with increased resolution from LR to VHR, such as the warm bias over the tropical upwelling region and the associated cloud cover underestimation, a precipitation excess over the tropical South Atlantic and North Pacific, and overly thick sea ice and an excess in oceanic mixing in the Arctic. VHR shows improved air–sea coupling over the tropical region, although it tends to overestimate the oceanic influence on the atmospheric variability at midlatitudes compared to observations and LR and HR. Together, these results highlight the potential for improved simulated climate in key regions, such as the Gulf Stream and the Equator, when the atmospheric and oceanic resolutions are finer than 25 km in both the ocean and atmosphere. Thanks to its unprecedented resolution, EC-Earth3P-VHR offers a new opportunity to study climate variability and change of such areas on regional and local spatial scales, in line with regional climate models. [ABSTRACT FROM AUTHOR]

The Earth system model EC-Earth3 for contributions to CMIP6 is documented here, with its flexible coupling framework, major model configurations, a methodology for ensuring the simulations are comparable across different high-performance computing (HPC) systems, and with the physical performance of base configurations over the historical period. The variety of possible configurations and sub-models reflects the broad interests in the EC-Earth community. EC-Earth3 key performance metrics demonstrate physical behavior and biases well within the frame known from recent CMIP models. With improved physical and dynamic features, new Earth system model (ESM) components, community tools, and largely improved physical performance compared to the CMIP5 version, EC-Earth3 represents a clear step forward for the only European community ESM. We demonstrate here that EC-Earth3 is suited for a range of tasks in CMIP6 and beyond.

In this paper, we present and evaluate the skill of an EC-Earth3.3 decadal prediction system contributing to the Decadal Climate Prediction Project – Component A (DCPP-A). This prediction system is capable of skilfully simulating past global mean surface temperature variations at interannual and decadal forecast times as well as the local surface temperature in regions such as the tropical Atlantic, the Indian Ocean and most of the continental areas, although most of the skill comes from the representation of the external radiative forcings. A benefit of initialization in the predictive skill is evident in some areas of the tropical Pacific and North Atlantic oceans in the first forecast years, an added value that is mostly confined to the south-east tropical Pacific and the eastern subpolar North Atlantic at the longest forecast times (6–10 years). The central subpolar North Atlantic shows poor predictive skill and a detrimental effect of initialization that leads to a quick collapse in Labrador Sea convection, followed by a weakening of the Atlantic Meridional Overturning Circulation (AMOC) and excessive local sea ice growth. The shutdown in Labrador Sea convection responds to a gradual increase in the local density stratification in the first years of the forecast, ultimately related to the different paces at which surface and subsurface temperature and salinity drift towards their preferred mean state. This transition happens rapidly at the surface and more slowly in the subsurface, where, by the 10th forecast year, the model is still far from the typical mean states in the corresponding ensemble of historical simulations with EC-Earth3. Thus, our study highlights the Labrador Sea as a region that can be sensitive to full-field initialization and hamper the final prediction skill, a problem that can be alleviated by improving the regional model biases through model development and by identifying more optimal initialization strategies.

We present a new framework for global ocean–sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean–sea-ice models (JRA55-do). We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean–ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean–sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80 % of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP-2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP-2. For example, the sea surface temperatures of the OMIP-2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating process-level responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean–sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.

We here present the very-high resolution version of the EC-Earth global climate model, EC-Earth3P-VHR, developed for HighResMIP. The model features an atmospheric resolution of ~16 km and an oceanic resolution of 1/12° (~8 km), which makes it one of the finest combined resolutions ever used to complete historical and scenario-like CMIP6 simulations. To evaluate the influence of numerical resolution on the simulated climate, EC-Earth3P-VHR is compared with two configurations of the same model at lower resolution: the ~100-km-grid EC-Earth3P-LR, and the ~25-km-grid EC-Earth3P-HR. The models' biases are evaluated against observations over the period 1980–2014. Compared to LR and HR, VHR shows a reduced equatorial Pacific cold tongue bias, an improved Gulf Stream representation with a reduced coastal warm bias and a reduced subpolar North Atlantic cold bias, and more realistic orographic precipitation over mountain ranges. By contrast, VHR shows a larger warm bias and overly low sea ice extent over the Southern Ocean. Such biases in surface temperature have an impact on the atmospheric circulation aloft, with improved stormtrack over the North Atlantic, yet worsened stormtrack over the Southern Ocean compared to the lower resolution model versions. Other biases persist with increased resolution from LR to VHR, such as the warm bias over the tropical upwelling region and the associated cloud cover underestimation, and the precipitation excess over the tropical South Atlantic and North Pacific. VHR shows improved air–sea coupling over the tropical region, although it tends to overestimate the oceanic influence on the atmospheric variability at mid-latitudes compared to observations and LR and HR. Together, these results highlight the potential for improved simulated climate in key regions, such as the Gulf Stream and the Equator, when the atmospheric and oceanic resolutions are finer than 25 km in both the ocean and atmosphere. Thanks to its unprecedented resolution, EC-Earth3P-VHR offers a new opportunity to study climate variability and change of such areas on regional/local spatial scales, in line with regional climate models. [ABSTRACT FROM AUTHOR]

The river routing scheme (RRS) in the Organising Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) land surface model is a valuable tool for closing the water cycle in a coupled environment and for validating the model performance. This study presents a revision of the RRS of the ORCHIDEE model that aims to benefit from the high-resolution topography provided by the Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales (HydroSHEDS), which is processed to a resolution of approximately 1 km. Adapting a new algorithm to construct river networks, the new RRS in ORCHIDEE allows for the preservation of as much of the hydrological information from HydroSHEDS as the user requires. The evaluation focuses on 12 rivers of contrasting size and climate which contribute freshwater to the Mediterranean Sea. First, the numerical aspect of the new RRS is investigated, in order to identify the practical configuration offering the best trade-off between computational cost and simulation quality for ensuing validations. Second, the performance of the new scheme is evaluated against observations at both monthly and daily timescales. The new RRS satisfactorily captures the seasonal variability of river discharge, although important biases stem from the water budget simulated by the ORCHIDEE model. The results highlight that realistic streamflow simulations require accurate precipitation forcing data and a precise river catchment description over a wide range of scales, as permitted by the new RRS. Detailed analyses at the daily timescale show the promising performance of this high-resolution RRS with respect to replicating river flow variation at various frequencies. Furthermore, this RRS may also eventually be well adapted for further developments in the ORCHIDEE land surface model to assess anthropogenic impacts on river processes (e.g. damming for irrigation operation).

Abstract The current feedback to the atmosphere (CFB) contributes to the oceanic circulation by damping eddies. In an ocean‐atmosphere coupled model, CFB can be correctly accounted for by using the wind relative to the oceanic current. However, its implementation in a forced oceanic model is less straightforward as CFB also enhances the 10‐m wind. Wind products based on observations have seen real currents that will not necessarily correspond to model currents, whereas meteorological reanalyses often neglect surface currents or use surface currents that, again, will differ from the surface currents of the forced oceanic simulation. In this study, we use a set of quasi‐global oceanic simulations, coupled or not with the atmosphere, to (i) quantify the error associated with the different existing strategies of forcing an oceanic model, (ii) test different parameterizations of the CFB, and (iii) propose the best strategy to account for CFB in forced oceanic simulation. We show that scatterometer wind or stress are not suitable to properly represent the CFB in forced oceanic simulation. We furthermore demonstrate that a parameterization of CFB based on a wind‐predicted coupling coefficient between the surface current and the stress allows us to reproduce the main characteristics of a coupled simulation. Such a parameterization can be used with any forcing set, including future coupled reanalyses, assuming that the associated oceanic surface currents are known. A further assessment of the thermal feedback of the surface wind in response to oceanic surface temperature gradients shows a weak forcing effect on oceanic currents.

An extensive compilation of published neodymium (Nd) concentrations and isotopic compositions (Nd IC) was realized in order to establish a new database and a map (using a high-resolution geological map of the area) of the distribution of these parameters for all the Mediterranean margins. Data were extracted from different kinds of samples: river solid discharge deposited on the shelf, sedimentary material collected on the margin or geological material outcropping above or close to a margin. Additional analyses of surface sediments were done in order to improve this data set in key areas (e.g. Sicilian strait). The Mediterranean margin Nd isotopic signatures vary from non-radiogenic values around the Gulf of Lion, (εNd values ∼ −11) to radiogenic values around the Aegean and the Levantine sub-basins up to +6. Using a high-resolution regional oceanic model (1/12° of horizontal-resolution), εNd distribution was simulated for the first time in the Mediterranean Sea. The high resolution of the model provides a unique opportunity to represent a realistic thermohaline circulation in the basin and thus apprehend the processes governing the Nd isotope distribution in the marine environment. Results are consistent with the preceding conclusions on boundary exchange (BE) as an important process in the Nd oceanic cycle. Nevertheless this approach simulates a too-radiogenic value in the Mediterranean Sea; this bias will likely be corrected once the dust and river inputs will be included in the model. This work highlights that a significant interannual variability of εNd distribution in seawater could occur. In particular, important hydrological events such as the Eastern Mediterranean Transient (EMT), associated with deep water formed in the Aegean sub-basin, could induce a shift in εNd at deep/intermediate depths that could be noticeable in the eastern part of the basin. This underlines that the temporal and geographical variations of εNd could represent an interesting insight of Nd as tracer of the Mediterranean Sea circulation, in particular in the context of palaeo-oceanographic applications.

The French research community in the Mediterranean Sea modeling and the French operational ocean forecasting center Mercator Océan have gathered their skill and expertise in physical oceanography, ocean modeling, atmospheric forcings and data assimilation to carry out a MEDiterranean sea ReanalYsiS (MEDRYS) at high resolution for the period 1992–2013. The ocean model used is NEMOMED12, a Mediterranean configuration of NEMO with a 1∕12° ( ∼ 7 km) horizontal resolution and 75 vertical z levels with partial steps. At the surface, it is forced by a new atmospheric-forcing data set (ALDERA), coming from a dynamical downscaling of the ERA-Interim atmospheric reanalysis by the regional climate model ALADIN-Climate with a 12 km horizontal and 3 h temporal resolutions. This configuration is used to carry a 34-year hindcast simulation over the period 1979–2013 (NM12-FREE), which is the initial state of the reanalysis in October 1992. MEDRYS uses the existing Mercator Océan data assimilation system SAM2 that is based on a reduced-order Kalman filter with a three-dimensional (3-D) multivariate modal decomposition of the forecast error. Altimeter data, satellite sea surface temperature (SST) and temperature and salinity vertical profiles are jointly assimilated. This paper describes the configuration we used to perform MEDRYS. We then validate the skills of the data assimilation system. It is shown that the data assimilation restores a good average temperature and salinity at intermediate layers compared to the hindcast. No particular biases are identified in the bottom layers. However, the reanalysis shows slight positive biases of 0.02 psu and 0.15 °C above 150 m depth. In the validation stage, it is also shown that the assimilation allows one to better reproduce water, heat and salt transports through the Strait of Gibraltar. Finally, the ability of the reanalysis to represent the sea surface high-frequency variability is shown.

The Mediterranean Sea is one of the most oligotrophic regions of the oceans, and nutrients have been shown to limit both phytoplankton and bacterial activities, resulting in a potential major role of dissolved organic carbon (DOC) export in the biological pump. Strong DOC accumulation in surface waters is already well documented, though measurements of DOC stocks and export flux are still sparse and associated with major uncertainties. This study provides the first basin-scale overview and analysis of organic carbon stocks and export fluxes in the Mediterranean Sea through a modeling approach based on a coupled model combining a mechanistic biogeochemical model (Eco3M-MED) and a high-resolution (eddy-resolving) hydrodynamic simulation (NEMO-MED12). The model is shown to reproduce the main spatial and seasonal biogeochemical characteristics of the Mediterranean Sea. Model estimations of carbon export are also of the same order of magnitude as estimations from in situ observations, and their respective spatial patterns are mutually consistent. Strong differences between the western and eastern basins are evidenced by the model for organic carbon export. Though less oligotrophic than the eastern basin, the western basin only supports 39 % of organic carbon (particulate and dissolved) export. Another major result is that except for the Alboran Sea, the DOC contribution to organic carbon export is higher than that of particulate organic carbon (POC) throughout the Mediterranean Sea, especially in the eastern basin. This paper also investigates the seasonality of DOC and POC exports as well as the differences in the processes involved in DOC and POC exports in light of intracellular quotas. Finally, according to the model, strong phosphate limitation of both bacteria and phytoplankton growth is one of the main drivers of DOC accumulation and therefore of export.

This numerical study provides the first simulation of the anthropogenic tritium invasion and its decay product helium-3 (3He) in the Mediterranean Sea. The simulation covers the entire tritium (3H) transient generated by the atmospheric nuclear weapons tests performed in the 1950s and early 1960s and is run till 2011. Tritium, helium-3 and their derived age estimates are particularly suitable for studying intermediate and deep-water ventilation and spreading of water masses at intermediate/deep levels. The simulation is made using a high-resolution regional model NEMO (Nucleus for European Modelling of the Ocean), in a regional configuration for the Mediterranean Sea called MED12, forced at the surface with prescribed tritium evolution derived from observations. The simulation is compared to measurements of tritium and helium-3 performed along large-scale transects in the Mediterranean Sea during the last few decades on cruises of R/V Meteor: M5/6, M31/1, M44/4, M51/2, M84/3, and R/V Poseidon: 234. The results show that the input function used for the tritium generates a realistic distribution of the main hydrographic features of the Mediterranean Sea circulation. In the eastern basin, the results highlight the weak formation of Adriatic Deep Water in the model, which explains its weak contribution to the Eastern Mediterranean Deep Water (EMDW) in the Ionian sub-basin. It produces a realistic representation of the Eastern Mediterranean Transient (EMT) signal, simulating a deep-water formation in the Aegean sub-basin at the beginning of 1993, with a realistic timing of deep-water renewal in the eastern basin.

Heavy precipitation events (HPEs) are frequent in southern France in autumn. An HPE results from landward transport of low-level moisture from the Western Mediterranean: large potential instability is then released by local convergence and/or orography. In the upstream zone, the sea surface temperature (SST) undergoes significant variations at the submonthly time scale primarily driven by episodic highly energetic events of relatively cold outflows from the neighbouring mountain ranges (the Mistral and Tramontane winds). Here, we study the HPE of 22–23 September 1994 which is preceded by a strong SST cooling due to the Mistral and Tramontane winds. This case confirms that the location of the precipitation is modulated by the SST in the upstream zone. In fact, changes in latent and sensible heat fluxes due to SST changes induce pressure and stratification changes which affect the low-level dynamics. Using three companion regional climate simulations running from 1989 to 2009, this article statistically shows that anomalies in the HPEs significantly correlate with the SST anomalies in the Western Mediterranean, and hence with the prior history of Mistral and Tramontane winds. In such cases, the role of the ocean as an integrator of the effect of past wind events over one or several weeks does indeed have an impact on HPEs in southern France.