Regional and global effects of southern ocean constraints in a global model

Global ocean circulation models do not usually take high-latitude processes into account in an adequate form due to a limited model domain or insufficient resolution. Without the processes in key areas contributing to the lower part of the global thermohaline circulation, the characteristics and flow of deep and bottom waters often remain unrealistic in these models. In this study, various sections of the Bremerhaven Regional Ice Ocean Simulation model results are combined with a global inverse model by using temperature, salinity, and velocity constraints for the Hamburg Large Scale Geostrophic ocean general circulation model. The differences between the global model with and without additional constraints from the regional model demonstrate that the Weddell Sea circulation exerts a significant influence on the course of the Antarctic Circumpolar Current with consequences for Southern Ocean water mass characteristics and the spreading of deep and bottom waters in the South Atlantic. The influence of the Ross Sea is found to be less important in terms of global influences. However, regional changes in the Pacific sector of the Southern Ocean are found to be of Ross Sea origin. The additional constraints change the hydrographic conditions of the global model in the vicinity of the Antarctic Circumpolar Current in such a way that transport values, e.g., in Drake Passage no longer need to be prescribed to obtain observed transports. These changes not only improve the path and transport of the Antarctic Circumpolar Current but affect the meso- and large-scale circulation. With a higher (lower) mean Drake Passage transport, the mean Weddell Gyre transport is lower (higher). Furthermore, an increase (decrease) in the Antarctic Circumpolar Current leads to a decrease (increase) of the circum-Australian flow, i.e., a decrease (increase) of the Indonesian Throughflow.     

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.     

Regional and interannual variability of ecosystem dynamics in the Southern Ocean

To investigate regional and interannual variability of the ecosystem in the Southern Ocean, a coupled circumpolar ice–ocean–plankton model has been developed. The ice–ocean component (known as BRIOS-2) is based on a modified version of the s-coordinate primitive equation model (SPEM) coupled to a dynamic–thermodynamic sea-ice model. The biological model (BIMAP) comprises two biogeochemical cycles – silica and nitrogen – and a prognostic iron compartment to include possible effects of micronutrient limitation. Simulations with the coupled ice–ocean–plankton model indicate that the physical–biological interaction is not limited to the effect of a varying surface mixed-layer depth. In the Pacific sector, large anomalies in winter mixed-layer depth cause an increased iron supply and enhance primary production and plankton biomass in the following summer, whereas in the Atlantic sector variability in primary production is caused mainly by fluctuations of oceanic upwelling. Thus, the Antarctic Circumpolar Wave (ACW) induces regional oscillations of phytoplankton biomass in both sectors, but not a propagating signal. Furthermore, interannual variability in plankton biomass and primary production is strong in the Coastal and Continental Shelf Zone and the Seasonal Ice Zone around the Antarctic continent. Interannual variability induced by the ACW has large effects on the regional scale, but the associated variability in biogenic carbon fluxes is small compared to the long-term carbon sequestration of the Southern Ocean.     

Passive tracers in a general circulation model of the Southern Ocean

Passive tracers are used in an off-line version of the United Kingdom Fine Resolution Antarctic Model (FRAM) to highlight features of the circulation and provide information on the inter-ocean exchange of water masses. The use of passive tracers allows a picture to be built up of the deep circulation which is not readily apparent from examination of the veloCity or density fields. Comparison of observations with FRAM results gives good agreement for many features of the Southern Ocean circulation. Tracer distributions are consistent with the concept of a global “conveyor belt” with a return path via the Agulhas retroflection region for the replenishment of North Atlantic Deep Water.     

Effect of wave-induced mixing on Antarctic sea ice in a high-resolution ocean model

Ocean Dynamics - In this study, we investigate how sea ice extent and thickness are affected by additional turbulent mixing from unbroken surface waves. This ‘surface wave mixing’ (SWM)...     

Interannual variability and characteristics of the East India Coastal Current associated with Indian Ocean Dipole events using a high resolution regional ocean model

The characteristics and variability of the East India Coastal Current (EICC), the western boundary current in the Bay of Bengal (BoB) during the Indian Ocean Dipole (IOD) years between 2006 and 2012 have been investigated using the high-resolution Regional Ocean Modeling System (ROMS). The evolution of temperature, mixed layer depth (MLD), and seasonal basin scale circulation in the upper ocean simulated by the model agrees well with the observations. The EICC in BoB is characterized by a seasonal reversal flow: the poleward EICC during February?May and the equatorward EICC during August?December. A long-term simulation from 2006 to 2012 suggest that the circulation pattern, boundary current structure, and transport in the western BoB are completely different in positive and negative IOD years. As IOD is mainly phase-locked to the seasonal cycle with most significant influence in the Borel autumn, the equatorward EICC is affected during the IOD years. It is found that the strength of this EICC is ~?5 Sv in October 2010 and a weaker EICC dominated by the presence of eddies is observed in October 2006. We also quantified the local and remote forcing effects on the variability of EICC and found that the seasonal coastal Kelvin waves (KWs) play a dominant role in the development of the EICC. During positive IOD year 2006, due the absence of second downwelling KW, the EICC is completely disorganized and dominated by the eddies, whereas in the negative IOD year 2010, the strong second downwelling KW plays a key role in developing organized and stable EICC in the western BoB.     

Intraseasonal variability of subsurface ocean temperature anomalies in the Indian Ocean during the summer monsoon season

Ocean Dynamics - The present study focuses on the variability of subsurface ocean temperature and associated planetary waves (oceanic Kelvin and Rossby waves) in the Indian Ocean during the boreal...     

Microplastics in the Southern Ocean

A field survey to collect microplastics with sizes < 5 mm was conducted in the Southern Ocean in 2016. We performed five net-tows and collected 44 pieces of plastic. Total particle counts of the entire water column, which is free of vertical mixing, were computed using the surface concentration (particle count per unit seawater volume) of microplastics, wind speed, and significant wave height during the observation period. Total particle counts at two stations near Antarctica were estimated to be in the order of 100,000 pieces km^? 2.     

On the influence of adequate Weddell Sea characteristics in a large-scale global ocean circulation model

Global ocean circulation models usually lack an adequate consideration of high-latitude processes due to a limited model domain or insufficient resolution. Without the processes in key areas of the global thermohaline circulation, the characteristics and flow of deep and bottom waters cannot be modeled realistically. In this study, a high-resolution (~20 km) ocean model focused on the Weddell Sea sector of the Southern Ocean is combined with a low-resolution (2^° × 2^°) global ocean model applying the state estimation technique. Temperature, salinity, and velocity data on two Weddell Sea sections from the regional model are used as constraints for the large-scale model in addition to satellite altimetry and sea-surface temperatures. The differences between the model with additional constraints and without document that the Weddell Sea circulation exerts significant influence on the course of the Antarctic Circumpolar Current with consequences for Southern Ocean water mass characteristics and the spreading of deep and bottom waters in the South Atlantic. Furthermore, a warming trend in the period 1993–2001 was found in the Weddell Sea and adjacent basins in agreement with float measurements in the upper Southern Ocean. Teleconnections to the North Atlantic are suggested but need further studies to demonstrate their statistical significance.     

A non-breaking-wave-generated turbulence mixing scheme for a global ocean general circulation model

Ocean Dynamics - A novel vertical mixing scheme to describe the influence of the non-breaking surface waves in ocean general circulation models is proposed based on the second-order turbulence...     

Ocean in peril: Reforming the management of global ocean living resources in areas beyond national jurisdiction

This article presents the outcome of research aimed at assisting governments in meeting their commitments and legal obligations for sustainable fisheries, based on increasing evidence that global fisheries are in crisis. The article assesses the effectiveness of the existing legal and institutional framework for high seas living resources. It focuses on: (1) the role of regional fisheries management organizations (RFMOs); (2) tools for compliance and enforcement to stem illegal fishing; and (3) mechanisms for habitat protection. The article further highlights a variety of options for addressing key weaknesses and gaps in current ocean governance, including United Nations General Assembly (UNGA) resolutions, reforms at the regional level, as well as a possible new legal instrument, with a view to informing international discussions on ways to ensure the sustainable use of high seas resources without compromising the health of the marine environment.     

Spatiotemporal variability of the ocean since 1900: testing a new analysis approach using global sea level reconstruction

Ocean Dynamics - A new approach for analysis of spatiotemporal variability across ocean basins was tested using global 1° × 1° monthly sea level reconstruction (RecSL)...     

Argo data assimilation in ocean general circulation model of Northwest Pacific Ocean

The Argo temperature and salinity profiles in 2005–2009 are assimilated into a coastal ocean general circulation model of the Northwest Pacific Ocean using the ensemble adjustment Kalman filter (EAKF). Three numerical tests, including the control run (CTL) (without data assimilation, which serves as the reference experiment), ensemble free run (EnFR) (without data assimilation), and EAKF experiment (with Argo data assimilation using EAKF), are carried out to examine the performance of this system. Using the restarts of different years as the initial conditions of the ensemble integrations, the ensemble spreads from EnFR and EAKF are all kept at a finite value after a sharp decreasing in the first few months because of the sensitive of the model to the initial conditions, and the reducing of the ensemble spread due to Argo data assimilation is not much. The ensemble samples obtained in this way can well represent the probabilities of the real ocean states, and no ensemble inflation is necessary for this EAKF experiment. Different experiment results are compared with satellite sea surface temperature (SST) data and the Global Temperature-Salinity Profile Program (GTSPP) data. The comparison of SST shows that modeled SST errors are reduced after data assimilation; the error reduction percentage after assimilating the Argo profiles is about 10?% on average. The comparison against the GTSPP profiles, which are independent of the Argo profiles, shows improvements in both temperature and salinity. The comparison results indicated a great error reduction in all vertical layers relative to CTL and the ensemble mean of EnFR; the maximum value for temperature and salinity reaches to 85?% and 80?%, respectively. The standard deviations of sea surface height are employed to examine the simulation ability, and it is shown that the mesoscale variability is improved after Argo data assimilation, especially in the Kuroshio extension area and along the section of 10°N. All these results suggest that this system is potentially useful for improving the simulation ability of oceanic numerical models.     

The roles of surface heat flux and ocean heat transport convergence in determining Atlantic Ocean temperature variability

The temperature variability of the Atlantic Ocean is investigated using an eddy-permitting (1/4°) global ocean model (ORCA-025) forced with historical surface meteorological fields from 1958 to 2001. The simulation of volume-averaged temperature and the vertical structure of the zonally averaged temperature trends are compared with those from observations. In regions with a high number of observations, in particular above a depth of 500 m and between 22° N and 65° N, the model simulation and the dataset are in good agreement. The relative contribution of variability in ocean heat transport (OHT) convergence and net surface heat flux to changes in ocean heat content is investigated with a focus on three regions: the subpolar and subtropical gyres and the tropics. The surface heat flux plays a relatively minor role in year-to-year changes in the subpolar and subtropical regions, but in the tropical North Atlantic, its role is of similar significance to the ocean heat transport convergence. The strongest signal during the study period is a cooling of the subpolar gyre between 1970 and 1990, which subsequently reversed as the mid-latitude OHT convergence transitioned from an anomalously weak to an anomalously strong state. We also explore whether model OHT anomalies can be linked to surface flux anomalies through a Hovmöller analysis of the Atlantic sector. At low latitudes, increased ocean heat gain coincides with anomalously strong northward transport, whereas at mid-high latitudes, reduced ocean heat loss is associated with anomalously weak heat transport.     

Low-frequency sea level variability and impact of recent sea ice decline on the sea level trend in the Arctic Ocean from a high-resolution simulation

Ocean Dynamics - The Arctic Ocean is undergoing significant changes, with rapid sea ice decline, unprecedented freshwater accumulation, and pronounced regional sea level rise. In this paper, we...     

Wave effects in global ocean modeling: parametrizations vs. forcing from a wave model

One of the main challenges of the Copernicus Marine Service is the implementation of coupled ocean/waves systems that accurately estimate the momentum and energy fluxes provided by the atmosphere to the ocean. This study aims to investigate the impact of forcing the Nucleus for European Modelling of the Ocean (NEMO) ocean model with forecasts from the wave model of Météo-France (MFWAM) to improve classical air-sea flux parametrizations, these latter being mostly driven by the 10-m wind. Three wave-related processes, namely, wave-state-dependent stress, Stokes drift-related effects (Stokes-Coriolis force, Stokes drift advection on tracers and on mass), and wave-state-dependent surface turbulence, are examined at a global scale with a horizontal resolution of 0.25°. Three years of sensitivity simulations (2014–2016) show positive feedback on sea surface temperature (SST) and currents when the wave model is used. A significant reduction in SST bias is observed in the tropical Atlantic Ocean. This is mainly due to the more realistic momentum flux provided by the wave model. In mid-latitudes, the most interesting impact occurs during the summer stratification, when the wind is low and the wave model produces a reduction in the turbulence linked with wave breaking. Magnitudes of the large-scale currents in the equatorial region are also improved by 10% compared to observations. In general, it is shown that using the wave model reduces on average the momentum and energy fluxes to the ocean in tropical regions, but increases them in mid-latitudes. These differences are in the order of 10 to 20% compared with the classical parametrizations found in stand-alone ocean models.     

The computation of the free barotropic oscillations of a global ocean model including friction and loading effects

A set of 61 normal modes with periods between 7.8 and 133.1 h has been calculated, using a 1° model of the global ocean, including the Arctic Ocean. The model explicitly considers frictional forces and ocean self-attraction and loading effects. The latter effects have generally been taken into account by parameterization, but for some modes the effects have also been considered fully. Due to friction, the computed eigenfrequencies are complex, exhibiting also the varying dissipative properties of the modes and their dependence on the distribution of potential and kinetic energies over the oceanic regions. In detail, gravity modes having periods less then 80 h and dominating the semi-diurnal and the diurnal tides, topographically controlled vorticity modes with periods longer than diurnal, and two planetary vorticity modes with periods of 96.8 and 119.4 h have been identified. These planetary vorticity modes have their energies distributed over Pacific, Atlantic, and Indian Oceans, while the other modes with periods longer than 80 h, as vorticity modes, have their energies concentrated on topographic structures of restricted extension. The modes are discussed with respect to their wave properties, e.g., concerning quasi-standing-wave resonances and to the appearance of Kelvin waves of different orders and trapped by different coastlines. In particular, the relevance of specific modes for the development of the fields of the most important semi-diurnal and diurnal tidal constituents is investigated.     

Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy-permitting resolution

Series of sensitivity tests were performed with a z-coordinate, global eddy-permitting (1/4°) ocean/sea-ice model (the ORCA-R025 model configuration developed for the DRAKKAR project) to carefully evaluate the impact of recent state-of-the-art numerical schemes on model solutions. The combination of an energy–enstrophy conserving (EEN) scheme for momentum advection with a partial step (PS) representation of the bottom topography yields significant improvements in the mean circulation. Well known biases in the representation of western boundary currents, such as in the Atlantic the detachment of the Gulf Stream, the path of the North Atlantic Current, the location of the Confluence, and the strength of the Zapiola Eddy in the south Atlantic, are partly corrected. Similar improvements are found in the Pacific, Indian, and Southern Oceans, and characteristics of the mean flow are generally much closer to observations. Comparisons with other state-of-the-art models show that the ORCA-R025 configuration generally performs better at similar resolution. In addition, the model solution is often comparable to solutions obtained at 1/6 or 1/10° resolution in some aspects concerning mean flow patterns and distribution of eddy kinetic energy. Although the reasons for these improvements are not analyzed in detail in this paper, evidence is shown that the combination of EEN with PS reduces numerical noise near the bottom, which is likely to affect current–topography interactions in a systematic way. We conclude that significant corrections of the mean biases presently seen in general circulation model solutions at eddy-permitting resolution can still be expected from the development of numerical methods, which represent an alternative to increasing resolution.