Global Ocean Simulation Research Articles

Modulation of Cross-Isothermal Velocities with ENSO in the Tropical Pacific Cold Tongue

AbstractThe tropical Pacific Ocean cold tongue (CT) plays a major role in the global climate system. The strength of the CT sets the zonal temperature gradient in the Pacific that couples with the atmospheric Walker circulation. This coupling is an essential component of the El Niño–Southern Oscillation (ENSO). The CT is supplied with cold water by the Equatorial Undercurrent that follows the thermocline as it shoals toward the east, adiabatically transporting cold water toward the surface. As the thermocline shoals, its water is transformed through diabatic processes, producing water mass transformation (WMT) that allows water to cross mean isotherms. Here, we examine WMT in the cold-tongue region from a global high-resolution ocean simulation with saved budget terms that close its heat budget exactly. Using the terms of the heat budget, we quantify each individual component of WMT (vertical mixing, horizontal mixing, eddy fluxes, and solar penetration) and find that vertical mixing is the single most important contribution in the thermocline and solar heating dominates close to the surface. Horizontal diffusion is much smaller. During El Niño events, vertical mixing, and hence cross-isothermal flow as a whole, are much reduced, whereas, during La Niña periods, strong vertical mixing leads to strong WMT, thereby cooling the surface. This analysis demonstrates the enhancement of diabatic processes during cold events, which in turn enhances cooling of the CT from below the surface.

Winter surface salinity in the northeastern Gulf of Maine from five years of SMAP satellite data

Wintertime freshwater transport into the Gulf of Maine (GoM) is typically controlled by a seasonal velocity increase in the fresh upstream Nova Scotia Current (NSC). Repeat satellite observations from the Soil Moisture Active Passive mission have mapped significant GoM surface salinity anomalies in four of five recent winters. These satellite data are used in combination with Jordan basin buoy and model datasets to investigate the likelihood that variable wind-forcing of the NSC contributed to these anomalies. This stems from regional ocean circulation studies suggesting that strengthening of alongshore southwesterly winds on the Scotian shelf weakens NSC transport into the GoM, while cross-shore southeasterly winds may also contribute to NSC weakening, and vice versa. A 17-year time series analysis of GoM buoy and satellite data shows that near-surface salinity in the eastern GoM can indeed be modulated by both alongshore and cross-shore winds through their impact on the NSC. The NSC geostrophic current modulation correlates with buoy-observed surface salinity anomalies when using a one month advective lead time. For a shorter 5-year period, a SMAP-derived salinity anomaly index in the eastern GoM indicates a similar relation to NSC variation and also correlates with Scotian Shelf wind anomalies. The relationships between winter wind and NSC transport variability are confirmed using output from a 19-year high resolution global ocean simulation. Attribution of these local wind anomalies to basin-scale atmospheric patterns shows that most years with strong winter GoM freshening coincide with weakening of a North Atlantic Oscillation-like wind pattern. But its reversals do not always correspond with a saline GoM. This contrast suggests that wind forcing more directly controls fresh winter GoM anomalies but not salty. Instead, they are partly due to density-driven advection from neighboring warm and saline Atlantic slope water including episodic Gulf Stream instabilities. • SMAP detects significant eastern GoM SSSA (up to 1 psu) in four of five winters. • Jordan basin SSSA correlation with Nova Scotia Current geostrophic anomaly is ~0.7. • Only half of the above correlation is explained by wind-driven current variations. • Eastern GoM satellite SSSA index is better indicator of fresh Scotian Shelf inflows.

A Kinematic Kinetic Energy Backscatter Parametrization: From Implementation to Global Ocean Simulations

AbstractOcean models at eddy‐permitting resolution are generally overdissipative, damping the intensity of the mesoscale eddy field. To reduce overdissipation, we propose a simplified, kinematic energy backscatter parametrization built into the viscosity operator in conjunction with a new flow‐dependent coefficient of viscosity based on nearest neighbor velocity differences. The new scheme mitigates excessive dissipation of energy and improves global ocean simulations at eddy‐permitting resolution. We find that kinematic backscatter substantially raises simulated eddy kinetic energy, similar to an alternative, previously proposed dynamic backscatter parametrization. While dynamic backscatter is scale‐aware and energetically more consistent, its implementation is more complex. Furthermore, it turns out to be computationally more expensive, as it applies, among other things, an additional prognostic subgrid energy equation. The kinematic backscatter proposed here, by contrast, comes at no additional computational cost, following the principle of simplicity. Our primary focus is the discretization on triangular unstructured meshes with cell placement of velocities (an analog of B‐grids), as employed by the Finite‐volumE Sea ice‐Ocean Model (FESOM2). The kinematic backscatter scheme with the new viscosity coefficient is implemented in FESOM2 and tested in the simplified geometry of a zonally reentrant channel as well as in a global ocean simulation on a 1/4° mesh. This first version of the new kinematic backscatter needs to be tuned to the specific resolution regime of the simulation. However, the tuning relies on a single parameter, emphasizing the overall practicality of the approach.

Antarctic Circumpolar Current Transport Through Drake Passage: What Can We Learn From Comparing High‐Resolution Model Results to Observations?

AbstractUncertainty exists in the time‐mean total transport of the Antarctic Circumpolar Current (ACC), the world's strongest ocean current. The two most recent observational programs in Drake Passage, DRAKE and cDrake, yielded transports of 141 and 173.3 Sv, respectively. In this paper, we use a realistic 1/12° global ocean simulation to interpret these observational estimates and reconcile their differences. We first show that the modeled ACC transport in the upper 1,000 m is in excellent agreement with repeat shipboard acoustic Doppler current profiler (SADCP) transects and that the exponentially decaying transport profile in the model is consistent with the profile derived from repeat hydrographic data. By further comparing the model results to the cDrake and DRAKE observations, we argue that the modeled 157.3 Sv transport, that is, approximately the average of the cDrake and DRAKE estimates, is actually representative of the time‐mean ACC transport through the Drake Passage. The cDrake experiment overestimated the barotropic contribution in part because the array undersampled the deep recirculation southwest of the Shackleton Fracture Zone, whereas the surface geostrophic currents used in the DRAKE estimate yielded a weaker near‐surface transport than implied by the SADCP data. We also find that the modeled baroclinic and barotropic transports are not correlated; thus, monitoring either baroclinic or barotropic transport alone may be insufficient to assess the temporal variability of the total ACC transport.

On the interplay between horizontal resolution and wave drag and their effect on tidal baroclinic mode waves in realistic global ocean simulations

The effects of horizontal resolution and wave drag damping on the semidiurnal M2 tidal energetics are studied for two realistically-forced global HYbrid Coordinate Ocean Model (HYCOM) simulations with 41 layers and horizontal resolutions of 8 km (1∕12.5∘; H12) and 4 km (1∕25∘; H25). In both simulations, the surface tidal error is minimized by tuning the strength of the linear wave drag, which is a parameterization of the surface-tide energy conversion to the unresolved baroclinic wave modes. In both simulations the M2 surface tide error with TPXO8-atlas, an altimetry constrained model, is 2.6 cm. Compared to H12, the surface tide energy conversion to the resolved vertical modes is increased by 50% in H25. This coincides with an equivalent reduction in the tuned loss of energy from the surface tide to the wave drag. For the configurations studied here, the horizontal and not the vertical resolution is the factor limiting the number of vertical modes that are resolved in most of the global ocean: modes 1–2 in H12 and modes 1–5 in H25. The wave drag also dampens the resolved internal tides. The 40% reduction in wave-drag strength does not result in a proportional increase in the mode-1 energy density in H25. In the higher-resolution simulations, topographic mode-scattering and wave–wave interactions are better resolved. This allows for an energy flux out of mode 1 to the higher modes, mitigating the need for an internal tide damping term. The HYCOM simulations are validated with analytical conversion models and altimetry-inferred sea-surface height, fluxes, and surface tide dissipation. H25 agrees best with these data sets to within ∼10%. To facilitate the comparison of stationary tide signals extracted from time series with different durations, we successfully apply a spatially-varying correction factor.

Impact of observation-based snow albedo parameterization on global ocean simulation results

Albedo parameterization is of fundamental importance for accurate representation of high-latitude climate variability by modeling studies. Field observations show that near-infrared snow albedo decreases dramatically when surface air temperature exceeds −2 °C. This can influence reproduction of sea ice simulations taking into consideration the importance of the drastic change in albedo in early melt season for the seasonal change of sea ice extent. Therefore, we conducted global ocean data-assimilative simulation experiments using a modified snow albedo parameterization. The modified parameterization reduced the albedo directly and achieved a comparable indirect reduction via changes in the modeled snow and sea ice distributions (ice–albedo feedback). As a result, sea ice thickness was reduced by more than 0.4–1 cm over most of the central Arctic Ocean. Sea ice velocities were also reduced by enhanced ocean drag with weakened surface ocean circulation in the Beaufort Gyre. In the Southern Ocean, the modified parameterization caused snow thicknesses to be decreased by up to 2 cm in the Weddell Sea. These impacts, which were generally larger than the spread of ensemble experiment results and therefore robust, at least in our model, provide useful information for quantifying the results of albedo modification in climate modeling studies.

Ocean Kinetic Energy Backscatter Parametrization on Unstructured Grids: Impact on Global Eddy‐Permitting Simulations

AbstractIn this study we demonstrate the potential of a kinetic energy backscatter scheme for use in global ocean simulations. Ocean models commonly employ (bi)harmonic eddy viscosities causing excessive dissipation of kinetic energy in eddy‐permitting simulations. Overdissipation not only affects the smallest resolved scales but also the generation of eddies through baroclinic instabilities, impacting the entire wave number spectrum. The backscatter scheme returns part of this overdissipated energy back into the resolved flow. We employ backscatter in the FESOM2 multiresolution ocean model with a quasi‐uniform 1/4° mesh. In multidecadal ocean simulations, backscatter increases eddy activity by a factor 2 or more, moving the simulation closer to observational estimates of sea surface height variability. Moreover, mean sea surface height, temperature, and salinity biases are reduced. This amounts to a globally averaged bias reduction of around 10% for each field, which is even larger in the Antarctic Circumpolar Current. However, in some regions such as the coastal Kuroshio, backscatter leads to a slight overenergizing of the flow and, in the Antarctic, to an unrealistic reduction of sea ice. Some of the bias increases can be reduced by a retuning of the model, and we suggest related adjustments to the backscatter scheme. The backscatter simulation is about 2.5 times as expensive as a simulation without backscatter. Most of the increased cost is due to a halving of the time step to accommodate higher simulated velocities.

Trends of Coastal Sea Level Between 1993 and 2015: Imprints of Atmospheric Forcing and Oceanic Chaos

The observation and simulation of the variability of coastal sea level are impacted by various uncertainties, such as measurement errors and sampling biases, unresolved processes, and model and forcing biases. Ocean model simulations suggest that another uncertainty should be taken into account for the attribution of sea-level changes. Global ocean simulations indeed show that resolving mesoscale turbulence (even partly) promotes the emergence of low-frequency (LF) chaotic intrinsic variability (CIV) which causes substantial random fluctuations of sea level up to multiple decades in eddy-active regions of the world ocean. This random LFCIV is superimposed on the atmospherically forced (or simply “forced”) fluctuations, which are directly controlled by the atmospheric variability. We show from a large ensemble of global oceanic hindcasts that this multi-decadal LFCIV leaves a substantial imprint on the long-term trends (1993–2015) of coastal sea level: over 17–20% of the global ocean coastal area, in particular along the coastlines of the northwestern Pacific and Indian Oceans, and around the Gulf of Mexico, random sea-level trends may blur their atmospherically forced counterpart, such that simulated (and potentially observed) coastal sea-level trends cannot be unambiguously attributed to atmospheric or anthropic causes. The steric and manometric sea-level change contributions of these uncertainties are discussed, suggesting that they mostly come from the manometric sea-level trends near the coasts.

Surface Kinetic Energy Distributions in the Global Oceans From a High‐Resolution Numerical Model and Surface Drifter Observations

AbstractThe surface kinetic energy of a 1/48° global ocean simulation and its distribution as a function of frequency and location are compared with the one estimated from 15,329 globally distributed surface drifter observations at hourly resolution. These distributions follow similar patterns with a dominant low‐frequency component and well‐defined tidal and near‐inertial peaks globally. Quantitative differences are identified with deficits of low‐frequency energy near the equator (factor 2) and at near‐inertial frequencies (factor 3) and an excess of energy at semidiurnal frequencies (factor 4) for the model. Owing to its hourly resolution and its near‐global spatial coverage, the array of surface drifters is an invaluable tool to evaluate the realism of tide‐resolving high‐resolution ocean simulations used in observing system simulation experiments. Sources of bias between model and drifter data are discussed, and associated leads for future work highlighted.

On the Long-Wavelength Validation of the SWOT KaRIn Measurement

AbstractThe Surface Water and Ocean Topography (SWOT) mission will measure the sea surface height (SSH) using a Ka-band radar interferometer (KaRIn) over a swath off the nadir of the satellite tracks. The mission requires calibration and validation (CalVal) of the SSH wavenumber spectrum at wavelengths between 15 and 1000 km. The CalVal in the short-wavelength range (15–150 km) requires in situ observations. In the long-wavelength range (150–1000 km), the CalVal will use the onboard Jason-class nadir altimeter. Using a high-resolution global ocean simulation, this study identifies the spatial scales beyond which the nadir and off-nadir observations can be considered comparable. Our results suggest that the ocean signals at nadir can represent off-nadir ocean signals at wavelengths longer than 120 and 70 km along the midswath and the inner edge of the KaRIn grid, respectively, indicating that the nadir altimeter is able to fulfill its goal to validate the long-wavelength KaRIn measurement. The wavelength along the inner edge is limited around 70 km because the onboard nadir altimeter cannot resolve spatial scales longer than ~70 km. These wavelengths provide a reference point for the required spatial coverage of the SWOT SSH in situ CalVal.

Improving surface tidal accuracy through two-way nesting in a global ocean model

In global ocean simulations, forward (non-data-assimilative) tide models generally feature large sea-surface-height errors near Hudson Strait in the North Atlantic Ocean with respect to altimetry-constrained tidal solutions. These errors may be associated with tidal resonances that are not well resolved by the complex coastal-shelf bathymetry in low-resolution simulations. An online two-way nesting framework has been implemented to improve global surface tides in the HYbrid Coordinate Ocean Model (HYCOM). In this framework, a high-resolution child domain, covering Hudson Strait, is coupled with a relatively low-resolution parent domain for computational efficiency. Data such as barotropic pressure and velocity are exchanged between the child and parent domains with the external coupler OASIS3-MCT. The developed nesting framework is validated with semi-idealized basin-scale model simulations. The M2 sea-surface heights show very good accuracy in the one-way and two-way nesting simulations in Hudson Strait, where large tidal elevations are observed. In addition, the mass and tidal energy flux are not adversely impacted at the nesting boundaries in the semi-idealized simulations. In a next step, the nesting framework is applied to a realistic global tide simulation. In this simulation, the resolution of the child domain (1/75°) is three times as fine as that of the parent domain (1/25°). The M2 sea-surface-height root-mean-square errors with tide gauge data and the altimetry-constrained global FES2014 and TPXO9-atlas tidal solutions are evaluated for the nesting and no-nesting solutions. The better resolved coastal bathymetry and the finer grid in the child domain improve the local tides in Hudson Strait and Bay, and the back-effect of the coastal tides induces an improvement of the barotropic tides in the open ocean of the Atlantic.

On the Spatial Scales to be Resolved by the Surface Water and Ocean Topography Ka-Band Radar Interferometer

AbstractThe Surface Water and Ocean Topography (SWOT) mission aims to measure the sea surface height (SSH) at a high spatial resolution using a Ka-band radar interferometer (KaRIn). The primary oceanographic objective is to characterize the ocean eddies at a spatial resolution of 15 km for 68% of the ocean surface. This resolution is derived from the ratio between the wavenumber spectrum of the conventional altimeter (projected to submesoscale) and the SWOT SSH errors. While the 15-km threshold is useful as a global approximation of the spatial scales resolved by SWOT (SWOT scale), it can be misleading for regional studies. Here we revisit the problem using a high-resolution (~2-km horizontal grid spacing) tide-resolving global ocean simulation and map the SWOT scale as a function of location and season. The results show that the SWOT scale increases, in general, from about 15 km at low latitudes to ~30–45 km at mid- and high latitudes but with a large geographical dependence. A SWOT scale smaller than 30 km is expected in the high-latitude energetic regions. The SWOT scale varies seasonally as a result of the seasonality in both the noise and ocean signals. The seasonality also has a geographical dependence. Both eddies and internal gravity waves/tides contribute significantly to the SWOT scale variation. Our analysis provides model predictions for interpreting the anticipated observations from SWOT and guidance for the development of analysis methodologies.

Impact of Tropical Cyclones on the Global Ocean: Results from Multidecadal Global Ocean Simulations Isolating Tropical Cyclone Forcing

Tropical cyclone (TC)-induced ocean vertical mixing can alter the upper-ocean temperature structure, influencing ocean heat content variability and meridional ocean heat transport. TC–ocean interactions can influence tropical variability on seasonal to interannual time scales. Here the impacts of TCs on the global ocean and the associated feedbacks are investigated using a hierarchy of high-resolution global ocean model simulations featuring the Community Earth System Model (CESM). The aim is to understand the potential impact of the model’s self-generated transient TC events on the modeled global ocean. Two ocean-only simulations are performed using the atmosphere boundary conditions from a fully coupled preindustrial CESM simulation configured with 0.25° atmosphere resolution and the nominal 1° ocean resolution (with ~0.25° meridional resolution in the tropics). The high-resolution coupled model is capable of directly simulating TC events with wind structure and climatology generally consistent with observations. TC effects at the ocean–atmosphere boundary are filtered out in one of the ocean simulations (OCN_FILT) while fully retained in the other (OCN_TC) in order to isolate the effect of the TCs on regional and global ocean variability across multiple time scales (from intraseasonal to interdecadal). Results show that the model-simulated TCs can 1) alter surface and subsurface ocean temperature patterns and variability; 2) affect ocean energetics, including increasing ocean mixed layer depth and strengthening subtropical gyre and meridional overturning circulations; and 3) influence ocean meridional heat transport and ocean heat content from seasonal to interannual time scales. Results help provide insights into the model behavior and the physical nature of the effect of TCs within the Earth system.

Mesoscale Effects on Carbon Export: A Global Perspective

Carbon export from the surface to the deep ocean is a primary control on global carbon budgets and is mediated by plankton that are sensitive to physical forcing. Earth system models generally do not resolve ocean mesoscale circulation ( (10–100) km), scales that strongly affect transport of nutrients and plankton. The role of mesoscale circulation in modulating export is evaluated by comparing global ocean simulations conducted at 1° and 0.1° horizontal resolution. Mesoscale resolution produces a small reduction in globally integrated export production (<2%); however, the impact on local export production can be large (±50%), with compensating effects in different ocean basins. With mesoscale resolution, improved representation of coastal jets block off‐shelf transport, leading to lower export in regions where shelf‐derived nutrients fuel production. Export is further reduced in these regions by resolution of mesoscale turbulence, which restricts the spatial area of production. Maximum mixed layer depths are narrower and deeper across the Subantarctic at higher resolution, driving locally stronger nutrient entrainment and enhanced summer export production. In energetic regions with seasonal blooms, such as the Subantarctic and North Pacific, internally generated mesoscale variability drives substantial interannual variation in local export production. These results suggest that biogeochemical tracer dynamics show different sensitivities to transport biases than temperature and salinity, which should be considered in the formulation and validation of physical parameterizations. Efforts to compare estimates of export production from observations and models should account for large variability in space and time expected for regions strongly affected by mesoscale circulation.

An Observing System Simulation Experiment for the Calibration and Validation of the Surface Water Ocean Topography Sea Surface Height Measurement Using In Situ Platforms

AbstractThe wavenumber spectrum of sea surface height (SSH) is an important indicator of the dynamics of the ocean interior. While the SSH wavenumber spectrum has been well studied at mesoscale wavelengths and longer, using both in situ oceanographic measurements and satellite altimetry, it remains largely unknown for wavelengths less than ~70 km. The Surface Water Ocean Topography (SWOT) satellite mission aims to resolve the SSH wavenumber spectrum at 15–150-km wavelengths, which is specified as one of the mission requirements. The mission calibration and validation (CalVal) requires the ground truth of a synoptic SSH field to resolve the targeted wavelengths, but no existing observational network is able to fulfill the task. A high-resolution global ocean simulation is used to conduct an observing system simulation experiment (OSSE) to identify the suitable oceanographic in situ measurements for SWOT SSH CalVal. After fixing 20 measuring locations (the minimum number for resolving 15–150-km wavelengths) along the SWOT swath, four instrument platforms were tested: pressure-sensor-equipped inverted echo sounders (PIES), underway conductivity–temperature–depth (UCTD) sensors, instrumented moorings, and underwater gliders. In the context of the OSSE, PIES was found to be an unsuitable tool for the target region and for SSH scales 15–70 km; the slowness of a single UCTD leads to significant aliasing by high-frequency motions at short wavelengths below ~30 km; an array of station-keeping gliders may meet the requirement; and an array of moorings is the most effective system among the four tested instruments for meeting the mission’s requirement. The results shown here warrant a prelaunch field campaign to further test the performance of station-keeping gliders.

Benchmarking the mesoscale variability in global ocean eddy-permitting numerical systems

The role of data assimilation procedures on representing ocean mesoscale variability is assessed by applying eddy statistics to a state-of-the-art global ocean reanalysis (C-GLORS), a free global ocean simulation (performed with the NEMO system) and an observation-based dataset (ARMOR3D) used as an independent benchmark. Numerical results are computed on a 1/4 ∘ horizontal grid (ORCA025) and share the same resolution with ARMOR3D dataset. This “eddy-permitting” resolution is sufficient to allow ocean eddies to form. Further to assessing the eddy statistics from three different datasets, a global three-dimensional eddy detection system is implemented in order to bypass the need of regional-dependent definition of thresholds, typical of commonly adopted eddy detection algorithms. It thus provides full three-dimensional eddy statistics segmenting vertical profiles from local rotational velocities. This criterion is crucial for discerning real eddies from transient surface noise that inevitably affects any two-dimensional algorithm. Data assimilation enhances and corrects mesoscale variability on a wide range of features that cannot be well reproduced otherwise. The free simulation fairly reproduces eddies emerging from western boundary currents and deep baroclinic instabilities, while underestimates shallower vortexes that populate the full basin. The ocean reanalysis recovers most of the missing turbulence, shown by satellite products , that is not generated by the model itself and consistently projects surface variability deep into the water column. The comparison with the statistically reconstructed vertical profiles from ARMOR3D show that ocean data assimilation is able to embed variability into the model dynamics, constraining eddies with in situ and altimetry observation and generating them consistently with local environment.

The relationship between Arabian Sea upwelling and Indian Monsoon revisited in a high resolution ocean simulation

Studies based on sediment records, sea-surface temperature and wind suggest that upwelling along the western coast of Arabian Sea is strongly affected by the Indian summer Monsoon. We examine this relationship directly in an eddy-resolving global ocean simulation STORM driven by atmospheric reanalysis over the last 61 years. With its very high spatial resolution (10 km), STORM allows us to identify characteristics of the upwelling system. We analyse the co-variability between upwelling and meteorological and oceanic variables from 1950 to 2010. The analysis reveals high interannual correlations between coastal upwelling and along-shore wind-stress (r = 0.73) as well as with sea-surface temperature (r = −0.83). However, the correlation between the upwelling and the Monsoon is small. We find an atmospheric circulation pattern different from the one that drives the Monsoon as the main modulator of the upwelling variability. In spite of this, the patterns of temperature anomalies that are either linked to Arabian Sea upwelling or to the Monsoon are spatially quite similar, although the physical mechanisms of these links are different. In addition, no long-term trend is detected in our modelled upwelling in the Arabian Sea.

Objectively estimating the temporal evolution of accuracy and skill in a global ocean reanalysis

ABSTRACTThe accuracy of ocean reanalyses is acknowledged to have increased dramatically during the last few decades, as a consequence of the evolution of the oceanic and atmospheric observing networks. Because of evolution of the atmospheric observing system, a similar accuracy increase also occurs in assimilation‐free global ocean simulations. However, the temporal evolution of traditional accuracy metrics based on comparing reanalyses with observations is affected by the change in the number and spatial sampling of the verifying observations themselves. To disentangle the accuracy evolution and the changes in the verifying dataset, procedures insensitive to the latter are recommended. A simple and easy‐to‐implement methodology inspired by bootstrapping is proposed, specifically conceived for objectively estimating the accuracy and skill evolution of reanalyses over long (e.g. inter‐decadal) periods. This allows the artificial impact of the observational network evolution on the verification metrics, which generally led to under‐estimated accuracy at the beginning of the reanalysis period and much too sharp an increase later, to be overcome. Applied to the assessment of ocean temperature accuracy and skill in a state‐of‐the‐art eddy‐permitting global ocean reanalysis covering the period 1982–2012, the methodology illuminates the crucial impact on the upper ocean (top 100 m) reanalysis skill of both the atmospheric and oceanic observing network evolution, the former leading to a quicker skill increase. Conversely, the fundamental effect of the oceanic observing network below 100 m emerges from the study, as the assimilation‐free simulation shows no skill at any time in the comparison against observations.

Designing variable ocean model resolution based on the observed ocean variability

AbstractIf unstructured meshes are refined to locally represent eddy dynamics in ocean circulation models, a practical question arises on how to vary the resolution and where to deploy the refinement. We propose to use the observed sea surface height variability as the refinement criterion. We explore the utility of this method (i) in a suite of idealized experiments simulating a wind‐driven double gyre flow in a stratified circular basin and (ii) in simulations of global ocean circulation performed with FESOM. Two practical approaches of mesh refinement are compared. In the first approach the uniform refinement is confined within the areas where the observed variability exceeds a given threshold. In the second one the refinement varies linearly following the observed variability. The resolution is fixed in time. For the double gyre case it is shown that the variability obtained in a high‐resolution reference run can be well captured on variable‐resolution meshes if they are refined where the variability is high and additionally upstream the jet separation point. The second approach of mesh refinement proves to be more beneficial in terms of improvement downstream the midlatitude jet. Similarly, in global ocean simulations the mesh refinement based on the observed variability helps the model to simulate high variability at correct locations. The refinement also leads to a reduced bias in the upper‐ocean temperature.

Energy budget-based backscatter in an eddy permitting primitive equation model

Increasing computational resources are starting to allow global ocean simulations at so-called “eddy-permitting” resolutions, at which the largest mesoscale eddies can be resolved explicitly. However, an adequate parameterization of the interactions with the unresolved part of the eddy energy spectrum remains crucial. Hyperviscous closures, which are commonly applied in eddy-permitting ocean models, cause spurious energy dissipation at these resolutions, leading to low levels of eddy kinetic energy (EKE) and weak eddy induced transports. It has recently been proposed to counteract the spurious energy dissipation of hyperviscous closures by an additional forcing term, which represents “backscatter” of energy from the un-resolved scales to the resolved scales. This study proposes a parameterization of energy backscatter based on an explicit sub-grid EKE budget. Energy dissipated by hyperviscosity acting on the resolved flow is added to the sub-grid EKE, while a backscatter term transfers energy back from the sub-grid EKE to the resolved flow. The backscatter term is formulated deterministically via a negative viscosity, which returns energy at somewhat larger scales than the hyperviscous dissipation, thus ensuring dissipation of enstrophy. The parameterization is tested in an idealized configuration of a primitive equation ocean model, and is shown to significantly improve the solutions of simulations at typical eddy-permitting resolutions.