Seasonal evolution of both surface signature and subsurface structure of a Mediterranean mesoscale anticyclones is assessed using the Coastal and Regional Ocean Community high‐resolution numerical ...model with realistic background stratification and fluxes. In good agreement with remote‐sensing and in‐situ observations, our numerical simulations capture the seasonal cycle of the anomalies induced by the anticyclone, both in the sea surface temperature (SST) and in the mixed layer depth (MLD). The eddy signature on the SST shifts from warm‐core in winter to cold‐core in summer, while the MLD deepens significantly in the core of the anticyclone in late winter. Our sensitivity analysis shows that the eddy SST anomaly can be accurately reproduced only if the vertical resolution is high enough (∼4 m in near surface) and if the atmospheric forcing contains high‐frequency. In summer with this configuration, the vertical mixing parameterized by the k − ϵ closure scheme is three times higher inside the eddy than outside the eddy, and leads to an anticyclonic cold core SST anomaly. This differential mixing is explained by near‐inertial waves, triggered by the high‐frequency atmospheric forcing. Near‐inertial waves propagate more energy inside the eddy because of the lower effective Coriolis parameter in the anticyclone core. On the other hand, eddy MLD anomaly appears more sensitive to horizontal resolution, and requires SST retroaction on air‐sea fluxes. These results detail the need of high frequency forcing, high vertical and horizontal resolutions to accurately reproduce the evolution of a mesoscale eddy.
Plain Language Summary
Mesoscale eddies are turbulent structures present in every regions of the world ocean, and accounting for a significant part of its kinetic energy budget. These structures can be tracked in time and recently revealed a seasonal cycle from in situ data. An anticyclone (clockwise rotating eddy in the northern hemisphere) is observed in the Mediterranean to be predominantly warm at the surface and to deepen the mixed layer in winter, but shifts to a cold‐core summer signature. This seasonal signal is not yet understood and studied in ocean models. In this study we assess the realism of an anticyclone seasonal evolution in high resolution numerical simulations. Eddy surface temperature seasonal shift is retrieved and is linked to an increased mixing at the eddy core spontaneously appearing at high vertical resolution (vertical grid size smaller than 4 m) in the presence of high frequency atmospheric forcing. This increased mixed is due to the preferred propagation of near‐inertial waves in the anticyclone due to its negative relative vorticity. Eddy‐induced mixed layer depth anomalies also appear to be triggered by sea surface temperature retroaction on air‐sea fluxes. These results suggest that present‐day operational ocean forecast models are too coarse to accurately retrieve mesoscale evolution.
Key Points
Enhanced mixing in anticyclones explain inverse eddy sea surface temperature (SST) signature
Vertical resolution is crucial to model eddy core mixing triggered by near‐inertial waves
Mixed layer anomaly is mainly driven by SST retroaction on air‐sea fluxes
Studying the larval dispersal of bottom-dwelling species is necessary to understand their population dynamics and optimize their management. The black-lip pearl oyster (Pinctada margaritifera) is ...cultured extensively to produce black pearls, especially in French Polynesia's atoll lagoons. This aquaculture relies on spat collection, a process that can be optimized by understanding which factors influence larval dispersal. Here, we investigate the sensitivity of P. margaritifera larval dispersal kernel to both physical and biological factors in the lagoon of Ahe atoll. Specifically, using a validated 3D larval dispersal model, the variability of lagoon-scale connectivity is investigated against wind forcing, depth and location of larval release, destination location, vertical swimming behavior and pelagic larval duration (PLD) factors. The potential connectivity was spatially weighted according to both the natural and cultivated broodstock densities to provide a realistic view of connectivity. We found that the mean pattern of potential connectivity was driven by the southwest and northeast main barotropic circulation structures, with high retention levels in both. Destination locations, spawning sites and PLD were the main drivers of potential connectivity, explaining respectively 26%, 59% and 5% of the variance. Differences between potential and realistic connectivity showed the significant contribution of the pearl oyster broodstock location to its own dynamics. Realistic connectivity showed larger larval supply in the western destination locations, which are preferentially used by farmers for spat collection. In addition, larval supply in the same sectors was enhanced during summer wind conditions. These results provide new cues to understanding the dynamics of bottom-dwelling populations in atoll lagoons, and show how to take advantage of numerical models for pearl oyster management.
Coastal waters are subject to great environmental and anthropogenic pressures. The diffusion and the transport of these waters are a key element for environmental, ecological and economic management. ...There are numerous indicators of hydrological characteristics based on theories of transport time scale. However, these indicators strongly depend on the geographical shape of the studied area and tend to give information after long integration time periods, generally on the order of weeks. Here, to qualify a coastal area’s dispersion more precisely, we combined two Lagrangian approaches and estimated a local diffusivity. This paper presents the numerical implementation and the results obtained over a tidally flushed, semi-enclosed water body located at mid-latitude. This new coefficient was estimated using the hydrodynamic model MARS 3D with a barycentric repositioning technique over a tidal period to ensure its reliability. We highlight the existing relationships between local diffusivity and both horizontal and vertical processes. Methodological aspects were analysed based on a reference case (number and distribution of particles, resolution, integration time period). The consistency and sensitivity of the coefficient were studied with different forcing conditions (hydrodynamical and meteorological regimes). In conclusion, our local diffusivity provides a new perspective for understanding the land–sea interface and coastal dispersion and holds potential for future studies of coastal marine ecosystems.
Tritium concentrations in oceans were compiled from the literature, online databases and original measurements in order to determine the global distribution of tritium concentrations according to ...latitude and depth in all oceans.
The total inventory of tritium decay corrected in 2016 has been estimated using evaluation of the natural and artificial contributions in 23 spatial subdivisions of the total ocean. It is determined equal to 26.8 ± 14 kg including 3.8 kg of cosmogenic tritium. That is in agreement with the total atmospheric input of tritium from nuclear bomb tests and the natural inventory at steady-state estimated from natural production rates in the literature (27.8–29.3 kg in the Earth). We confirm the global increase in tritium according to latitude observed in the Northern hemisphere since 1967 with a maximum in the Arctic Ocean. The minimum tritium concentrations observed in the Southern Ocean were close to steady-state with known natural tritium deposition.
We focused on the temporal evolution of surface (0 to 500 m) tritium concentrations in a selected area of the North Atlantic Ocean (30°N–60°N) where we found the 2016 concentration to be 0.60 ± 0.10 TU (1σ). Results showed that in that area, between 1988 and 2013, tritium concentrations: i) decreased faster than the sole radioactive decay, due to a mixing with lower and lateral less concentrated waters, and ii) decreased towards an apparent steady state concentration. The half-time mixing rate of surface waters and the steady state concentration were respectively calculated to be 23 ± 5 years (1σ) and 0.38 ± 0.07 TU (1σ). This apparent steady-state concentration in the North Atlantic Ocean implies a mean tritium deposition of 1870 ± 345 Bq·m−2 (1σ), five folds higher than the known inputs (natural, nuclear tests fallout and industrial releases, ~367 Bq·m−2) in this area.
Display omitted
•A database of tritium concentrations available in oceans has been gathered.•Inventory of oceanic tritium has been calculated.•North-Atlantic Ocean tritium background concentration has been estimated.•The Antarctic Ocean is at an apparent steady-state with the natural production rate.•The North-Atlantic Ocean concentration exhibits an underestimation of known sources.
► Discussion on theoretical wave-current interaction models. ► Comparison of set of equations based on the Lagrangian velocity with set of equations based on the quasi-Eulerian velocity. ► Adiabatic ...tests.
Many theoretical approaches and implementations have been proposed for the coupling of the three-dimensional ocean circulation with waves. The theoretical models are reviewed and it is shown that the formulation in terms of the quasi-Eulerian velocity circumvents the essential difficulty of alternative formulations for the Lagrangian mean velocity. Namely, models based on this Lagrangian velocity require an estimation of wave-induced motions to first order in the horizontal gradients of the wave field in order to estimate the vertical flux of wave pseudo-momentum. So far, only three-dimensional wave models have been able to provide these estimates, and all published theories based on the simpler Airy theory are not consistent at the leading order, because they ignore or incorrectly estimate the vertical momentum flux. With an adiabatic example on a sloping bottom it is shown that this inconsistency produces very large spurious velocities. These errors are independent of the slope for the inviscid case, and are still significant when a realistic vertical mixing is applied. A quick diagnostic of the potential accuracy of a theoretical model is the vertical profile of the wave-induced forcing terms: if it is not uniform over depth in adiabatic conditions then it will produce spurious artificial flow patterns in conditions with shoaling waves. Although conceptually more challenging, the quasi-Eulerian velocity theories only introduce minor modifications of the solution procedure for the standard primitive equations: a modification of the surface boundary condition for the mass conservation, the addition of the Stokes drift in the tracer advection equations, and sources of momentum and turbulent kinetic energy with associated surface and bottom fluxes. All the necessary modifications of primitive equation models are given in detail. This implementation is illustrated with the MARS3D model, which passes the test of the adiabatic shoaling waves.
In this study a three-dimensional sediment transport model was developed. The model accounts for both current and wave forcing on the sediment and was implemented over the Gulf of Lions. A two-way ...nesting technique was used to focus on the Rhone River prodelta which is considered as a sink for riverine sediment. In addition, to understand the resuspension of trapped sediment over the Rhone prodelta, an in situ experiment, called SCOPE, was conducted during the winter 2007–2008. The experiment consisted of measuring hydro-sedimentary parameters using a mooring station comprising a current profiler (ADCP) and an altimeter (acoustic transducer) located in the eastern part of the Rhone prodelta. The three-dimensional transport model was validated using these data, and used to investigate the effect of sediment dynamics at the prodelta and Gulf of Lions scale. Both modelling and data analysis highlighted the impact of the two strong storms from the south-east which characterised the experimental period. Erosion of bed material (about 2cm) and an increase in suspended material (up to about 50–100mg/l) in the water were the result of each storm as recorded at the mooring station. The erosion capacity due to waves, combined with a strong current, due to both wind and wave forcing, resulted in strong south-westward export over the whole prodelta. Each storm was responsible for an off-prodelta export estimated at around 2.1Mt. This study demonstrates that the Rhone River sediments trapped over the Rhone prodelta are subject to strong resuspension during episodic events.
•A 3D sediment transport model is developed.•The hydro-sedimentary in situ SCOPE experiment is presented.•The model is validated using the various in situ data.•The impact of winter storms over the Rhone prodelta is assessed.
Motivated by the recurrent formation of mesoscale anticyclones in the southeast of Crete, we investigated with a high resolution model the response of the ocean to orographic wind jets driven by the ...Cretean mountain range. As shown in the dynamical process study of Ioannou et al. (
2020
) which uses a simplified shallow-water model, we confirm here, using the CROCO (Coastal and Regional Ocean COmmunity) model, that the main oceanic response to the Etesian wind forcing is the formation of mesoscale anticyclones. Moreover, we found that the intensity of the wind-induced Ekman pumping acting on the eddies, once they are formed, modulates their intensity. Among the various coastal anticyclones formed during summer and fall 2015, only one of them will correspond to a long lived structure (M_IE15) which is similar to the Ierapetra Eddy detected in 2015 (O_IE15) on the AVISO/DUACS products. Thanks to the DYNED-Atlas data base, we were able to perform a quantitative comparison of the vertical structure of such long-lived anticyclone between the numerical model and the
in-situ
measurements of the various Argo profilers trapped inside the eddy core. Even without assimilation or any nudging, the numerical model was able to reproduce correctly the formation period, the seasonal evolution and the vertical structure of the O_IE15. The main discrepancy between the model and the altimetry observations is the dynamical intensity of the anticyclone. The characteristic eddy velocity derived from the AVISO/DUACS product for the O_IE15 is much lower than in the numerical model. This is probably due to the spatio temporal interpolation of the AVISO/DUACS altimetry products. More surprisingly, several coastal anticyclones were also formed in the model in the lee of Crete area during summer 2015 when the Etesian winds reach strong values. However, these coastal anticyclones respond differently to the wind forcing since they remain close to the coast, in shallow-waters, unlike the M_IE15 which propagates offshore in deep water. The impact of the bottom friction or the coastal dissipation seems to limit the wind amplification of these coastal anticyclones.
The hydrology of the Bay of Biscay was investigated using the regional ocean model MARS3D (Model for Application at Regional Scale). The simulated hydrology is compared to a set of various data ...encompassing monthly climatology, remote sensing SST, CTD casts, and coastal salinity measurements. Special focus was put on the validation over the continental shelf. This paper reports that despite some misfits, the climatological hydrology and its seasonal variability are correctly simulated. Various statistics computed over the period from 1999–2004 highlight different aspects of the hydrology. The biases and root mean square errors (RMSE) remain very weak at all depths when comparing salinity (<0.1 and <0.6
psu respectively). The predicted temperature shows a global overestimation of temperature (bias of around 0.8
°C) and the maximum errors are located near the thermocline (rmse of 1
°C at 20–40
m). The model is shown to properly reproduce the annual dynamics of sea surface temperature, as well as the dynamics of large river plumes observed by high frequency time series from coastal salinity gauges. The misfits highlighted by these various comparisons between model and observations are attributed to heat fluxes and mixing parameterisation.
Significant amounts of anthropogenic radionuclides were introduced in ocean waters following nuclear atmospheric tests and development of the nuclear industry. Dispersion of artificial dissolved ...radionuclides has been extensively measured for decades over the North-European continental shelf. In this area, the radionuclide measurement and release fluxes databases provided here between 1982 and 2016 represent an exceptional opportunity to validate dispersion hydrodynamic models. This work gives accessibility to these data in a comprehensive database. The MARS hydrodynamic model has been applied at different scales to reproduce the measured dispersion in realistic conditions. Specific methods have been developed to obtain qualitative and quantitative results and perform model/measurement comparisons. Model validation concerns short to large scales with dedicated surveys following the dispersion: it was performed within a two- and three-dimensional framework and from minutes and hours following a release up to several years. Results are presented concerning the dispersion of radionuclides in marine systems deduced from standalone measurements, or according to model comparisons. It allows characterizing dispersion over the continental shelf, pathways, transit times, budgets and source terms. This review presents the main approaches developed and types of information derived from studies of artificial radiotracers using observations, hydrodynamic models or a combination of the two, based primarily on the new featured datasets.
The seasonal variability of the carbonate system in the eastern Mediterranean Sea (EMed) was investigated based on discrete total alkalinity (A
T
), total dissolved inorganic carbon (C
T
), and pH ...measurements collected during three cruises around Crete between June 2018 and March 2019. This study presents a detailed description of this new carbonate chemistry dataset in the eastern Mediterranean Sea. We show that the North Western Levantine Basin (NWLB) is unique in terms of range of A
T
variation vs. C
T
variation in the upper water column over an annual cycle. The reasons for this singularity of the NWLB can be explained by the interplay between strong evaporation and the concomitant consumption of C
T
by autotrophic processes. The high range of A
T
variations, combined to temperature changes, has a strong impact on the variability of the seawater
p
CO
2
(
p
CO
2
S
W
). Based on Argo float data, an entire annual cycle for
p
CO
2
S
W
in the NWLB has been reconstructed in order to estimate the temporal sequence of the potential “source” and “sink” of atmospheric CO
2
. By combining this dataset with previous observations in the NWLB, this study shows a significant ocean acidification and a decrease in the oceanic surface pH
T
25
of −0.0024 ± 0.0004 pH
T
25
units.a
–1
. The changes in the carbonate system are driven by the increase of atmospheric CO
2
but also by unexplained temporal changes in the surface A
T
content. If we consider that the EMed will, in the future, encounter longer, more intense and warmer summer seasons, this study proposes some perspectives on the carbonate system functioning of the “future” EMed.