Nucleus For European Modelling Of The Ocean Research Articles

Mixed layer depth calculation in deep convection regions in ocean numerical models

Mixed Layer Depths (MLDs) diagnosed by conventional numerical models are generally based on a density difference with the surface (e.g., 0.01 kg.m−3). However, the temperature-salinity compensation and the lack of vertical resolution contribute to over-estimated MLD, especially in regions of deep convection. In the present work, we examined the diagnostic MLD, associated with the deep convection of the Labrador Sea Water (LSW), calculated with a simple density difference criterion. The over-estimated MLD led us to develop a new tool, based on an observational approach, to recalculate MLD from model output. We used an eddy-permitting, 1/12° regional configuration of the Nucleus for European Modelling of the Ocean (NEMO) to test and discuss our newly defined MLD. We compared our new MLD with that from observations, and we showed a major improvement with our new algorithm. To show the new MLD is not dependent on a single model and its horizontal resolution, we extended our analysis to include 1/4° eddy-permitting simulations, and simulations using the Modular Ocean Model (MOM) model.

Explicit representation and parametrised impacts of under ice shelf seas in the <i>z</i><sup>∗</sup> coordinate ocean model NEMO 3.6

Abstract. Ice-shelf–ocean interactions are a major source of freshwater on the Antarctic continental shelf and have a strong impact on ocean properties, ocean circulation and sea ice. However, climate models based on the ocean–sea ice model NEMO (Nucleus for European Modelling of the Ocean) currently do not include these interactions in any detail. The capability of explicitly simulating the circulation beneath ice shelves is introduced in the non-linear free surface model NEMO. Its implementation into the NEMO framework and its assessment in an idealised and realistic circum-Antarctic configuration is described in this study. Compared with the current prescription of ice shelf melting (i.e. at the surface), inclusion of open sub-ice-shelf cavities leads to a decrease in sea ice thickness along the coast, a weakening of the ocean stratification on the shelf, a decrease in salinity of high-salinity shelf water on the Ross and Weddell sea shelves and an increase in the strength of the gyres that circulate within the over-deepened basins on the West Antarctic continental shelf. Mimicking the overturning circulation under the ice shelves by introducing a prescribed meltwater flux over the depth range of the ice shelf base, rather than at the surface, is also assessed. It yields similar improvements in the simulated ocean properties and circulation over the Antarctic continental shelf to those from the explicit ice shelf cavity representation. With the ice shelf cavities opened, the widely used three equation ice shelf melting formulation, which enables an interactive computation of melting, is tested. Comparison with observational estimates of ice shelf melting indicates realistic results for most ice shelves. However, melting rates for the Amery, Getz and George VI ice shelves are considerably overestimated.

Tuning without over-tuning: parametric uncertainty quantification for the NEMO ocean model

Abstract. In this paper we discuss climate model tuning and present an iterative automatic tuning method from the statistical science literature. The method, which we refer to here as iterative refocussing (though also known as history matching), avoids many of the common pitfalls of automatic tuning procedures that are based on optimisation of a cost function, principally the over-tuning of a climate model due to using only partial observations. This avoidance comes by seeking to rule out parameter choices that we are confident could not reproduce the observations, rather than seeking the model that is closest to them (a procedure that risks over-tuning). We comment on the state of climate model tuning and illustrate our approach through three waves of iterative refocussing of the NEMO (Nucleus for European Modelling of the Ocean) ORCA2 global ocean model run at 2° resolution. We show how at certain depths the anomalies of global mean temperature and salinity in a standard configuration of the model exceeds 10 standard deviations away from observations and show the extent to which this can be alleviated by iterative refocussing without compromising model performance spatially. We show how model improvements can be achieved by simultaneously perturbing multiple parameters, and illustrate the potential of using low-resolution ensembles to tune NEMO ORCA configurations at higher resolutions.

Comparing sea ice, hydrography and circulation between NEMO3.6 LIM3 and LIM2

Abstract. A set of hindcast simulations with the new version 3.6 of the Nucleus for European Modelling of the Ocean (NEMO) ocean–ice model in the ORCA1 configuration and forced by the DRAKKAR Forcing Set version 5.2 (DFS5.2) atmospheric data was performed from 1958 to 2012. Simulations differed in their sea-ice component: the old standard version Louvain-la-Neuve Sea Ice Model (LIM2) and its successor LIM3. Main differences between these sea-ice models are the parameterisations of sub-grid-scale sea-ice thickness distribution, ice deformation, thermodynamic processes, and sea-ice salinity. Our main objective was to analyse the response of the ocean–ice system sensitivity to the change in sea-ice physics. Additional sensitivity simulations were carried out for the attribution of observed differences between the two main simulations.In the Arctic, NEMO-LIM3 compares better with observations by realistically reproducing the sea-ice extent decline during the last few decades due to its multi-category sea-ice thickness. In the Antarctic, NEMO-LIM3 more realistically simulates the seasonal evolution of sea-ice extent than NEMO-LIM2. In terms of oceanic properties, improvements are not as evident, although NEMO-LIM3 reproduces a more realistic hydrography in the Labrador Sea and in the Arctic Ocean, including a reduced cold temperature bias of the Arctic Intermediate Water at 250 m. In the extra-polar regions, the oceanographic conditions of the two NEMO-LIM versions remain relatively similar, although they slowly drift apart over decades. This drift is probably due to a stronger deep water formation around Antarctica in LIM3.

Optimizing domain decomposition in an ocean model: the case of NEMO

Earth System Models are critical tools for the study of our climate and its future trends. These models are in constant evolution and their growing complexity entails an incrementing demand of the resources they require. Since the cost of using these state-of-the-art models is huge, looking closely at the factors that are able to impact their computational performance is mandatory. In the case of the state-of-the-art ocean model NEMO (Nucleus for European Modelling of the Ocean), used in many projects around the world, not enough attention has been given to the domain decomposition. In this work we show the impact that the selection of a particular domain decomposition can have on computational performance and how the proposed methodology substantially improves it.

Effects of wave-induced forcing on a circulation model of the North Sea

The effect of wind waves on water level and currents during two storms in the North Sea is investigated using a high-resolution Nucleus for European Modelling of the Ocean (NEMO) model forced with fluxes and fields from a high-resolution wave model. The additional terms accounting for wave-current interaction that are considered in this study are the Stokes-Coriolis force, the sea-state-dependent energy and momentum fluxes. The individual and collective role of these processes is quantified and the results are compared with a control run without wave effects as well as against current and water-level measurements from coastal stations. We find a better agreement with observations when the circulation model is forced by sea-state-dependent fluxes, especially in extreme events. The two extreme events, the storm Christian (25–27 October 2013), and about a month later, the storm Xaver (5–7 December 2013), induce different wave and surge conditions over the North Sea. Including the wave effects in the circulation model for the storm Xaver raises the modelled surge by more than 40 cm compared with the control run in the German Bight area. For the storm Christian, a difference of 20–30 cm in the surge level between the wave-forced and the stand-alone ocean model is found over the whole southern part of the North Sea. Moreover, the modelled vertical velocity profile fits the observations very well when the wave forcing is accounted for. The contribution of wave-induced forcing has been quantified indicating that this represents an important mechanism for improving water-level and current predictions.

Effects of lateral processes on the seasonal water stratification of the Gulf of Finland: 3-D NEMO-based model study

Abstract. This paper aims to fill the gaps in knowledge of processes affecting the seasonal water stratification in the Gulf of Finland (GOF). We used a state-of-the-art modelling framework NEMO (Nucleus for European Modelling of the Ocean) designed for oceanographic research, operational oceanography, seasonal forecasting, and climate studies to build an eddy-resolving model of the GOF. To evaluate the model skill and performance, two different solutions were obtained on 0.5 km eddy-resolving and commonly used 2 km grids for a 1-year simulation. We also explore the efficacy of non-hydrostatic effect (convection) parameterizations available in NEMO for coastal application. It is found that the solutions resolving submesoscales have a more complex mixed layer structure in the regions of the GOF directly affected by the upwelling/downwelling and intrusions from the open Baltic Sea. Presented model estimations of the upper mixed layer depth are in good agreement with in situ CTD (BED) data. A number of model sensitivity tests to the vertical mixing parameterization confirm the model's robustness. Further progress in the submesoscale process simulation and understanding is apparently not connected mainly with the finer resolution of the grids, but with the use of non-hydrostatic models because of the failure of the hydrostatic approach at submesoscale.

Simulation of the 2014 Anomalous Warming in the Northeast Pacific

ABSTRACTThe Nucleus for European Modelling of the Ocean (NEMO) numerical model was used to simulate the North Pacific Ocean beginning in January 1960. The model had a horizontal resolution of 0.25o, 46 vertical levels, and employed a spectral nudging assimilation scheme that, unlike standard nudging, nudges only specific frequency and wavenumber bands. This simulation was nudged to the mean and monthly Levitus climatology of potential temperature and salinity. The model was forced with mean monthly winds, net heat flux, and precipitation from the National Centers for Environmental Prediction (NCEP). The simulation was used to study a recent intrusion of much warmer and less saline water than normal from the west into the Gulf of Alaska, beginning in December 2013 and lasting until at least the early summer of 2014. The observed surface temperature anomalies were more than 4 standard deviations above normal. The model reproduced these anomalies in both a qualitative and quantitative manner, reproducing the same scale of anomalies over the region. An anomalous increase in the North Pacific Current (NPC) was found in the model in 2012 and the beginning of 2013, in agreement with observations. This increase in the NPC is associated with the positive phase of the North Pacific Gyre Oscillation. The causes of the temperature anomalies in the Gulf of Alaska could be due to three key factors: (i) an anomalously high, positive, surface heat flux in 2013 over the greater North Pacific; (ii) a significant decrease in the eastward flow of the NPC starting in late 2013 (with an accompanying decrease in cold water advection) after a period of historically strong eastward flow; and (iii) weaker winds throughout most of 2013 accompanied, however, by a shift to stronger northward winds (with an accompanying increase in warm water advection) in late 2013.

Results of the first Wave Glider experiment in the southern Tyrrhenian Sea

<p>A wave-propelled autonomous vehicle (Wave Glider) instrumented with a variety of oceanographic and meteorological sensors was launched from Gulf of Naples on the 12<sup>th </sup>of September 2012 for a two-week mission in the Southern Tyrrhenian Sea. The main objective of the mission was a preliminary evaluation of the potential of commercial autonomous platforms to provide reliable measurements of sea surface parameters which can complement existing satellite based products moving from the local to the synoptic scale. To this aim Wave Glider measurements were compared to equivalent, or near-equivalent, satellite products achieved from MODIS (Moderate Resolution Imaging Spectroradiometer) sensors onboard the EOS (Earth Observing System) satellite platforms and from AVISO (Archiving Validation and Interpretation of Satellite Oceanographic Data). Level-3 near real time and Level-4 reprocessed sea surface foundation temperature products provided by the CMEMS (Copernicus Marine Environment Monitoring Service) were also included in this study as well as high resolution model output supplied by NEMO (Nucleus for European Modelling of the Ocean). The Wave Glider was equipped with sensors to measure temperature, salinity, currents, as well as Colored Dissolved Organic Matter (CDOM), turbidity and refined fuels fluorescence. The achieved results confirmed the emerging value of Wave Gliders in the framework of multiplatform monitoring systems of the ocean surface parameters. In particular, they showed that Wave Glider measurements captured the southern Tyrrhenian Sea major surface oceanographic features, including the coast to open sea haline gradient and the presence of a cyclone-anticyclone system in the southeastern sub-region. The Wave Glider also had the capability to monitor upper ocean currents at finer spatial and temporal scales than satellite altimetric observations and model outputs. Nonetheless, results stressed the existence of several limits in the combined use of satellite and Wave Glider observations and the necessity of further analyses concerning the monitoring of the ocean optical properties. In fact, Wave Glider and satellite-based products agree in terms of sea surface temperature and currents patterns, while bio-optical properties turned out to be less well correlated. No significant traces of refined fuels have been detected along the WG track.</p>

The Louvain-La-Neuve sea ice model LIM3.6: global and regional capabilities

Abstract. The new 3.6 version of the Louvain-la-Neuve sea ice model (LIM) is presented, as integrated in the most recent stable release of Nucleus for European Modelling of the Ocean (NEMO) (3.6). The release will be used for the next Climate Model Inter-comparison Project (CMIP6). Developments focussed around three axes: improvements of robustness, versatility and sophistication of the code, which involved numerous changes. Robustness was improved by enforcing exact conservation through the inspection of the different processes driving the air–ice–ocean exchanges of heat, mass and salt. Versatility was enhanced by implementing lateral boundary conditions for sea ice and more flexible ice thickness categories. The latter includes a more practical computation of category boundaries, parameterizations to use LIM3.6 with a single ice category and a flux redistributor for coupling with atmospheric models that cannot handle multiple sub-grid fluxes. Sophistication was upgraded by including the effect of ice and snow weight on the sea surface. We illustrate some of the new capabilities of the code in two standard simulations. One is an ORCA2-LIM3 global simulation at a nominal 2° resolution, forced by reference atmospheric climatologies. The other one is a regional simulation at 2 km resolution around the Svalbard Archipelago in the Arctic Ocean, with open boundaries and tides. We show that the LIM3.6 forms a solid and flexible base for future scientific studies and model developments.

PISCES-v2: an ocean biogeochemical model for carbon and ecosystem studies

Abstract. PISCES-v2 (Pelagic Interactions Scheme for Carbon and Ecosystem Studies volume 2) is a biogeochemical model which simulates the lower trophic levels of marine ecosystems (phytoplankton, microzooplankton and mesozooplankton) and the biogeochemical cycles of carbon and of the main nutrients (P, N, Fe, and Si). The model is intended to be used for both regional and global configurations at high or low spatial resolutions as well as for short-term (seasonal, interannual) and long-term (climate change, paleoceanography) analyses. There are 24 prognostic variables (tracers) including two phytoplankton compartments (diatoms and nanophytoplankton), two zooplankton size classes (microzooplankton and mesozooplankton) and a description of the carbonate chemistry. Formulations in PISCES-v2 are based on a mixed Monod–quota formalism. On the one hand, stoichiometry of C / N / P is fixed and growth rate of phytoplankton is limited by the external availability in N, P and Si. On the other hand, the iron and silicon quotas are variable and the growth rate of phytoplankton is limited by the internal availability in Fe. Various parameterizations can be activated in PISCES-v2, setting, for instance, the complexity of iron chemistry or the description of particulate organic materials. So far, PISCES-v2 has been coupled to the Nucleus for European Modelling of the Ocean (NEMO) and Regional Ocean Modeling System (ROMS) systems. A full description of PISCES-v2 and of its optional functionalities is provided here. The results of a quasi-steady-state simulation are presented and evaluated against diverse observational and satellite-derived data. Finally, some of the new functionalities of PISCES-v2 are tested in a series of sensitivity experiments.

NEMO–ICB (v1.0): interactive icebergs in the NEMO ocean model globally configured at eddy-permitting resolution

Abstract. An established iceberg module, ICB, is used interactively with the Nucleus for European Modelling of the Ocean (NEMO) ocean model in a new implementation, NEMO–ICB (v1.0). A 30-year hindcast (1976–2005) simulation with an eddy-permitting (0.25°) global configuration of NEMO–ICB is undertaken to evaluate the influence of icebergs on sea ice, hydrography, mixed layer depths (MLDs), and ocean currents, through comparison with a control simulation in which the equivalent iceberg mass flux is applied as coastal runoff, a common forcing in ocean models. In the Southern Hemisphere (SH), drift and melting of icebergs are in balance after around 5 years, whereas the equilibration timescale for the Northern Hemisphere (NH) is 15–20 years. Iceberg drift patterns, and Southern Ocean iceberg mass, compare favourably with available observations. Freshwater forcing due to iceberg melting is most pronounced very locally, in the coastal zone around much of Antarctica, where it often exceeds in magnitude and opposes the negative freshwater fluxes associated with sea ice freezing. However, at most locations in the polar Southern Ocean, the annual-mean freshwater flux due to icebergs, if present, is typically an order of magnitude smaller than the contribution of sea ice melting and precipitation. A notable exception is the southwest Atlantic sector of the Southern Ocean, where iceberg melting reaches around 50% of net precipitation over a large area. Including icebergs in place of coastal runoff, sea ice concentration and thickness are notably decreased at most locations around Antarctica, by up to ~ 20% in the eastern Weddell Sea, with more limited increases, of up to ~ 10% in the Bellingshausen Sea. Antarctic sea ice mass decreases by 2.9%, overall. As a consequence of changes in net freshwater forcing and sea ice, salinity and temperature distributions are also substantially altered. Surface salinity increases by ~ 0.1 psu around much of Antarctica, due to suppressed coastal runoff, with extensive freshening at depth, extending to the greatest depths in the polar Southern Ocean where discernible effects on both salinity and temperature reach 2500 m in the Weddell Sea by the last pentad of the simulation. Substantial physical and dynamical responses to icebergs, throughout the global ocean, are explained by rapid propagation of density anomalies from high-to-low latitudes. Complementary to the baseline model used here, three prototype modifications to NEMO–ICB are also introduced and discussed.

Modelling of the anthropogenic tritium transient and its decay product helium-3 in the Mediterranean Sea using a high-resolution regional model

Abstract. 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.

CMCCCCESMMNEMO: Toward the New CMCC Earth System Model

This report describes the physical core of the new CMCC Earth System Model CMCC–CESM–NEMO. The model is largely based on the Community Earth System Model (CESM) project developed at the National Center for Atmospheric Research (NCAR). The ocean component of the CESM model has been replaced by the Nucleus for European Modelling of the Ocean (NEMO) general circulation model, developed by an European consortium of modelling centers of which the CMCC is an active member. An overview of the CMCC–CESM–NEMO model is provided, together with a detailed description of the integration of NEMO into CESM and of the design of the physical coupling interface.

Exploratory ensemble designs for environmental models using k-extended Latin Hypercubes.

In this paper we present a novel, flexible, and multi-purpose class of designs for initial exploration of the parameter spaces of computer models, such as those used to study many features of the environment. The idea applies existing technology aimed at expanding a Latin Hypercube (LHC) in order to generate initial LHC designs that are composed of many smaller LHCs. The resulting design and its component parts are designed so that each is approximately orthogonal and maximises a measure of coverage of the parameter space. Designs of the type advocated for in this paper are particularly useful when we want to simultaneously quantify parametric uncertainty and any uncertainty due to the initial conditions, boundary conditions, or forcing functions required to run the model. This makes the class of designs particularly suited to environmental models, such as climate models that contain all of these features. The proposed designs are particularly suited to initial exploratory ensembles whose goal is to guide the design of further ensembles aimed at, for example, calibrating the model. We introduce a new emulator diagnostic that exploits the structure of the advocated ensemble designs and allows for the assessment of structural weaknesses in the statistical modelling. We provide illustrations of the method through a simple example and describe a 400 member ensemble of the Nucleus for European Modelling of the Ocean (NEMO) ocean model designed using the method. We build an emulator for NEMO using the created design to illustrate the use of our emulator diagnostic test. © 2015 The Authors. Environmetrics published by John Wiley & Sons Ltd.

Major modes of short-term climate variability in the newly developed NUIST Earth System Model (NESM)

Abstract A coupled earth system model (ESM) has been developed at the Nanjing University of Information Science and Technology (NUIST) by using version 5.3 of the European Centre Hamburg Model (ECHAM), version 3.4 of the Nucleus for European Modelling of the Ocean (NEMO), and version 4.1 of the Los Alamos sea ice model (CICE). The model is referred to as NUIST ESM1 (NESM1). Comprehensive and quantitative metrics are used to assess the model’s major modes of climate variability most relevant to subseasonal-to-interannual climate prediction. The model’s assessment is placed in a multi-model framework. The model yields a realistic annual mean and annual cycle of equatorial SST, and a reasonably realistic precipitation climatology, but has difficulty in capturing the spring-fall asymmetry and monsoon precipitation domains. The ENSO mode is reproduced well with respect to its spatial structure, power spectrum, phase locking to the annual cycle, and spatial structures of the central Pacific (CP)-ENSO and eastern Pacific (EP)-ENSO; however, the equatorial SST variability, biennial component of ENSO, and the amplitude of CP-ENSO are overestimated. The model captures realistic intraseasonal variability patterns, the vertical-zonal structures of the first two leading predictable modes of Madden-Julian Oscillation (MJO), and its eastward propagation; but the simulated MJO speed is significantly slower than observed. Compared with the T42 version, the high resolution version (T159) demonstrates improved simulation with respect to the climatology, interannual variance, monsoon-ENSO lead-lag correlation, spatial structures of the leading mode of the Asian-Australian monsoon rainfall variability, and the eastward propagation of the MJO.

전지구 해양·해빙예측시스템 NEMO-CICE/NEMOVAR의 북극 영역 해빙초기조건 특성 분석

s: In this study, the characteristics of sea ice initial conditions generated from a global ocean and sea ice prediction system, the Nucleus for European Modeling of the Ocean (NEMO) - Los Alamos Sea Ice Model (CICE)/ NEMOVAR were analyzed for the period June 2013 to May 2014 over the Arctic region. For the purpose, the observed and reanalyzed data were used to compare with the sea ice initial conditions. Results indicated that the variability of the monthly sea ice extent and thickness in model initial conditions were well represented as compared to the observation, while it was found that the extent and thickness of Arctic sea ice in initial data were narrower and thinner than those in reanalysis and observation for the period. The reason for the narrower sea ice extent in model initial conditions seems to be due to the fact that the initial sea ice concentration at the boundary area of sea ice was about 20 percent less than the reanalysis data. Also, the reason for the thinner sea-ice thickness in the Arctic region is due to the underestimation of Arctic sea ice thickness (about 60 cm) of the model initial conditions in the Arctic Ocean area adjacent to Greenland and Arctic archipelago where thick sea ice appears all the year round.

IMarNet: an ocean biogeochemistry model intercomparison project within a common physical ocean modelling framework

Abstract. Ocean biogeochemistry (OBGC) models span a wide variety of complexities, including highly simplified nutrient-restoring schemes, nutrient–phytoplankton–zooplankton–detritus (NPZD) models that crudely represent the marine biota, models that represent a broader trophic structure by grouping organisms as plankton functional types (PFTs) based on their biogeochemical role (dynamic green ocean models) and ecosystem models that group organisms by ecological function and trait. OBGC models are now integral components of Earth system models (ESMs), but they compete for computing resources with higher resolution dynamical setups and with other components such as atmospheric chemistry and terrestrial vegetation schemes. As such, the choice of OBGC in ESMs needs to balance model complexity and realism alongside relative computing cost. Here we present an intercomparison of six OBGC models that were candidates for implementation within the next UK Earth system model (UKESM1). The models cover a large range of biological complexity (from 7 to 57 tracers) but all include representations of at least the nitrogen, carbon, alkalinity and oxygen cycles. Each OBGC model was coupled to the ocean general circulation model Nucleus for European Modelling of the Ocean (NEMO) and results from physically identical hindcast simulations were compared. Model skill was evaluated for biogeochemical metrics of global-scale bulk properties using conventional statistical techniques. The computing cost of each model was also measured in standardised tests run at two resource levels. No model is shown to consistently outperform all other models across all metrics. Nonetheless, the simpler models are broadly closer to observations across a number of fields and thus offer a high-efficiency option for ESMs that prioritise high-resolution climate dynamics. However, simpler models provide limited insight into more complex marine biogeochemical processes and ecosystem pathways, and a parallel approach of low-resolution climate dynamics and high-complexity biogeochemistry is desirable in order to provide additional insights into biogeochemistry–climate interactions.

Sensitivity of the Mediterranean sea level to atmospheric pressure and free surface elevation numerical formulation in NEMO

Abstract. The sensitivity of the dynamics of the Mediterranean Sea to atmospheric pressure and free surface elevation formulation using NEMO (Nucleus for European Modelling of the Ocean) was evaluated. Four different experiments were carried out in the Mediterranean Sea using filtered or explicit free surface numerical schemes and accounting for the effect of atmospheric pressure in addition to wind and buoyancy fluxes. Model results were evaluated by coherency and power spectrum analysis with tide gauge data. We found that atmospheric pressure plays an important role for periods shorter than 100 days. The free surface formulation is important to obtain the correct ocean response for periods shorter than 30 days. At frequencies higher than 15 days−1 the Mediterranean basin's response to atmospheric pressure was not coherent and the performance of the model strongly depended on the specific area considered. A large-amplitude seasonal oscillation observed in the experiments using a filtered free surface was not evident in the corresponding explicit free surface formulation case, which was due to a phase shift between mass fluxes in the Gibraltar Strait and at the surface. The configuration with time splitting and atmospheric pressure always performed best; the differences were enhanced at very high frequencies.

Recent development of the Met Office operational ocean forecasting system: an overview and assessment of the new Global FOAM forecasts

Abstract. The Forecast Ocean Assimilation Model (FOAM) is an operational ocean analysis and forecast system run daily at the Met Office. FOAM provides modelling capability in both deep ocean and coastal shelf sea regimes using the NEMO (Nucleus for European Modelling of the Ocean) ocean model as its dynamical core. The FOAM Deep Ocean suite produces analyses and 7-day forecasts of ocean tracers, currents and sea ice for the global ocean at 1/4° resolution. Satellite and in situ observations of temperature, salinity, sea level anomaly and sea ice concentration are assimilated by FOAM each day over a 48 h observation window. The FOAM Deep Ocean configurations have recently undergone a major upgrade which has involved the implementation of a new variational, first guess at appropriate time (FGAT) 3D-Var, assimilation scheme (NEMOVAR); coupling to a different, multi-thickness-category, sea ice model (CICE); the use of coordinated ocean-ice reference experiment (CORE) bulk formulae to specify the surface boundary condition; and an increased vertical resolution for the global model. In this paper the new FOAM Deep Ocean system is introduced and details of the recent changes are provided. Results are presented from 2-year reanalysis integrations of the Global FOAM configuration including an assessment of short-range ocean forecast accuracy. Comparisons are made with both the previous FOAM system and a non-assimilative FOAM system. Assessments reveal considerable improvements in the new system to the near-surface ocean and sea ice fields. However there is some degradation to sub-surface tracer fields and in equatorial regions which highlights specific areas upon which to focus future improvements.