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.
•PO4 availability in the sub-surface layer of the EMB significantly lowered in the last decades.•The mean DOC concentrations in the surface layer has increased by 20% on average over the EMB.•A new ...potential explanation is provided for the shift between the top of N and P nutriclines in the EMB .•A decline in abundance and the reduction of size of copepods has been modelled over the last decades.
The Mediterranean Sea (MS) is a semi-enclosed sea characterized by a zonal west-east gradient of oligotrophy, where microbial growth is controlled by phosphate availability in most situations. External inputs of nutrients including Gibraltar inputs, river inputs and atmospheric deposition are therefore of major importance for the biogeochemistry of the MS. The latter has long been considered to be driven mainly by nutrient exchanges at Gibraltar. However, recent studies indicate that river inputs significantly affect nutrients concentrations in the Mediterranean Sea, although their resulting impact on its biogeochemistry remains poorly understood.
In this study, our aim was to help fill this knowledge gap by addressing the large-scale and long-term impact of variations in river inputs on the biogeochemistry of the Mediterranean Sea over the last decades, using a coupled physical-biogeochemical 3D model (NEMO-MED12/Eco3M-Med).
As a first result, it has been shown by the model that the strong diminution (60%) of phosphate (PO4) in river inputs into the Mediterranean Sea since the end of the 1980s induced a significant lowering of PO4 availability in the sub-surface layer of the Eastern Mediterranean Basin (EMB). One of the main consequences of PO4 diminution is the rise, never previously documented, of dissolved organic carbon (DOC) concentrations in the surface layer (by 20% on average over the EMB). Another main result concerns the gradual deepening of the top of the phosphacline during the period studied, thus generating a shift between the top of the nitracline and the top of the phosphacline in the EMB. This shift has already been observed in situ and documented in literature, but we propose here a new explanation for its occurrence in the EMB. The last main result is the evidence of the decline in abundance and the reduction of size of copepods calculated by the model over the years 1985–2010, that could partially explain the reduction in size of anchovy and sardine recently recorded in the MS.
In this study, it is shown for the first time that the variations in river inputs that occurred in the last decades may have significantly altered the biogeochemical cycles of two key elements (P and C), in particular in the EMB. To conclude, the magnitude of the biogeochemical changes induced by river inputs and runoff alone over the last thirty years clearly calls for the use of realistic scenarios of river inputs along with climate scenarios for coupled physical-biogeochemical forecasts in the MS.
For the first time, modeling of the neodymium (Nd) isotopic composition (expressed as
ε
Nd) is performed in an Ocean Global Circulation Model (OPA-ORCA2.0 model). In the interest of simplicity, Nd ...concentration is not explicitly included and
ε
Nd is modeled as a passive tracer. Recent studies suggest that dissolved/particulate exchanges between continental margin sediments and seawater (termed boundary exchange, BE), could be the dominant source-sink terms that determine the distribution of neodymium isotopes in the global ocean. To test this hypothesis, we only included this BE source-sink term in the model. BE was parameterized as a relaxing term that forces
ε
Nd of seawater to the isotopic composition of the continental margin over a characteristic exchange time. Sensitivity tests were made using different exchange times and parameterizations. The model reproduces most of the general trends of
ε
Nd composition in seawater at global scale (inter-basin gradient, composition of some of the main water masses). These results suggest that BE process may account for the major sources and sinks of oceanic Nd, while indicating that the characteristic time of interaction between seawater and margins varies from about half a year in surface waters to several years at depth. The disparities between the model and observations might reflect the typical shortcomings of utilizing a coarse resolution ocean model circulation fields or neglected sources of neodymium such as dissolution of aeolian inputs or dissolved river loads.
The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model ...Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948-2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF subscript 6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.
A model of aluminium has been developed and implemented in an Ocean General Circulation Model (NEMO-PISCES). In the model, aluminium enters the ocean by means of dust deposition. The internal oceanic ...processes are described by advection, mixing and reversible scavenging. The model has been evaluated against a number of selected high-quality datasets covering much of the world ocean, especially those from the West Atlantic Geotraces cruises of 2010 and 2011. Generally, the model results are in fair agreement with the observations. However, the model does not describe well the vertical distribution of dissolved Al in the North Atlantic Ocean. The model may require changes in the physical forcing and the vertical dependence of the sinking velocity of biogenic silica to account for other discrepancies. To explore the model behaviour, sensitivity experiments have been performed, in which we changed the key parameters of the scavenging process as well as the input of aluminium into the ocean. This resulted in a better understanding of aluminium in the ocean, and it is now clear which parameter has what effect on the dissolved aluminium distribution and which processes might be missing in the model, among which boundary scavenging and biological incorporation of aluminium into diatoms.
A suite of standard ocean hydrographic and circulation metrics are applied to the equilibrium physical solutions from 13 global carbon models participating in phase 2 of the Ocean Carbon‐cycle Model ...Intercomparison Project (OCMIP‐2). Model‐data comparisons are presented for sea surface temperature and salinity, seasonal mixed layer depth, meridional heat and freshwater transport, 3‐D hydrographic fields, and meridional overturning. Considerable variation exists among the OCMIP‐2 simulations, with some of the solutions falling noticeably outside available observational constraints. For some cases, model‐model and model‐data differences can be related to variations in surface forcing, subgrid‐scale parameterizations, and model architecture. These errors in the physical metrics point to significant problems in the underlying model representations of ocean transport and dynamics, problems that directly affect the OCMIP predicted ocean tracer and carbon cycle variables (e.g., air‐sea CO2 flux, chlorofluorocarbon and anthropogenic CO2 uptake, and export production). A substantial fraction of the large model‐model ranges in OCMIP‐2 biogeochemical fields (±25–40%) represents the propagation of known errors in model physics. Therefore the model‐model spread likely overstates the uncertainty in our current understanding of the ocean carbon system, particularly for transport‐dominated fields such as the historical uptake of anthropogenic CO2. A full error assessment, however, would need to account for additional sources of uncertainty such as more complex biological‐chemical‐physical interactions, biases arising from poorly resolved or neglected physical processes, and climate change.
Daily modeled fields of phosphate deposition to the Mediterranean from natural dust, anthropogenic combustion and wildfires were used to assess the effect of this external nutrient on marine ...biogeochemistry. The ocean model used is a high-resolution (1∕12°) regional coupled dynamical–biogeochemical model of the Mediterranean Sea (NEMO-MED12/PISCES). The input fields of phosphorus are for 2005, which are the only available daily resolved deposition fields from the global atmospheric chemical transport model LMDz-INCA. Traditionally, dust has been suggested to be the main atmospheric source of phosphorus, but the LMDz-INCA model suggests that combustion is dominant over natural dust as an atmospheric source of phosphate (PO4, the bioavailable form of phosphorus in seawater) for the Mediterranean Sea. According to the atmospheric transport model, phosphate deposition from combustion (Pcomb) brings on average 40.5×10−6 mol PO4 m−2 yr−1 over the entire Mediterranean Sea for the year 2005 and is the primary source over the northern part (e.g., 101×10−6 mol PO4 m−2 yr−1 from combustion deposited in 2005 over the north Adriatic against 12.4×10−6 from dust). Lithogenic dust brings 17.2×10−6 mol PO4 m−2 yr−1 on average over the Mediterranean Sea in 2005 and is the primary source of atmospheric phosphate to the southern Mediterranean Basin in our simulations (e.g., 31.8×10−6 mol PO4 m−2 yr−1 from dust deposited in 2005 on average over the south Ionian basin against 12.4×10−6 from combustion). The evaluation of monthly averaged deposition flux variability of Pdust and Pcomb for the 1997–2012 period indicates that these conclusions may hold true for different years. We examine separately the two atmospheric phosphate sources and their respective flux variability and evaluate their impacts on marine surface biogeochemistry (phosphate concentration, chlorophyll a, primary production). The impacts of the different phosphate deposition sources on the biogeochemistry of the Mediterranean are found localized, seasonally varying and small, but yet statistically significant. Differences in the geographical deposition patterns between phosphate from dust and from combustion will cause contrasted and significant changes in the biogeochemistry of the basin. We contrast the effects of combustion in the northern basin (Pcomb deposition effects are found to be 10 times more important in the northern Adriatic, close to the main source region) to the effects of dust in the southern basin. These different phosphorus sources should therefore be accounted for in modeling studies.
The distribution of dissolved aluminium in the West Atlantic Ocean shows a mirror image with that of dissolved silicic acid, hinting at intricate interactions between the ocean cycling of Al and Si. ...The marine biogeochemistry of Al is of interest because of its potential impact on diatom opal remineralisation, hence Si availability. Furthermore, the dissolved Al concentration at the surface ocean has been used as a tracer for dust input, dust being the most important source of the bio-essential trace element iron to the ocean. Previously, the dissolved concentration of Al was simulated reasonably well with only a dust source, and scavenging by adsorption on settling biogenic debris as the only removal process. Here we explore the impacts of (i) a sediment source of Al in the Northern Hemisphere (especially north of ~ 40° N), (ii) the imposed velocity field, and (iii) biological incorporation of Al on the modelled Al distribution in the ocean. The sediment source clearly improves the model results, and using a different velocity field shows the importance of advection on the simulated Al distribution. Biological incorporation appears to be a potentially important removal process. However, conclusive independent data to constrain the Al / Si incorporation ratio by growing diatoms are missing. Therefore, this study does not provide a definitive answer to the question of the relative importance of Al removal by incorporation compared to removal by adsorptive scavenging.
Dissolved manganese (Mn) is a biologically essential element. Moreover, its oxidised form is involved in removing itself and several other trace elements from ocean waters. Here we report the longest ...thus far (17 500 km length) full-depth ocean section of dissolved Mn in the west Atlantic Ocean, comprising 1320 data values of high accuracy. This is the GA02 transect that is part of the GEOTRACES programme, which aims to understand trace element distributions. The goal of this study is to combine these new observations with new, state-of-the-art, modelling to give a first assessment of the main sources and redistribution of Mn throughout the ocean. To this end, we simulate the distribution of dissolved Mn using a global-scale circulation model. This first model includes simple parameterisations to account for the sources, processes and sinks of Mn in the ocean. Oxidation and (photo)reduction, aggregation and settling, as well as biological uptake and remineralisation by plankton are included in the model. Our model provides, together with the observations, the following insights: – The high surface concentrations of manganese are caused by the combination of photoreduction and sources contributing to the upper ocean. The most important sources are sediments, dust, and, more locally, rivers. – Observations and model simulations suggest that surface Mn in the Atlantic Ocean moves downwards into the southward-flowing North Atlantic Deep Water (NADW), but because of strong removal rates there is no elevated concentration of Mn visible any more in the NADW south of 40° N. – The model predicts lower dissolved Mn in surface waters of the Pacific Ocean than the observed concentrations. The intense oxygen minimum zone (OMZ) in subsurface waters is deemed to be a major source of dissolved Mn also mixing upwards into surface waters, but the OMZ is not well represented by the model. Improved high-resolution simulation of the OMZ may solve this problem. – There is a mainly homogeneous background concentration of dissolved Mn of about 0.10–0.15 nM throughout most of the deep ocean. The model reproduces this by means of a threshold on particulate manganese oxides of 25 pM, suggesting that a minimal concentration of particulate Mn is needed before aggregation and removal become efficient. – The observed distinct hydrothermal signals are produced by assuming both a strong source and a strong removal of Mn near hydrothermal vents.
We use a state-of-the-art ocean general circulation and biogeochemistry model to examine the impact of changes in ocean circulation and biogeochemistry in governing the change in ocean carbon-13 and ...atmospheric CO2 at the last glacial maximum (LGM). We examine 5 different realisations of the ocean's overturning circulation produced by a fully coupled atmosphere-ocean model under LGM forcing and suggested changes in the atmospheric deposition of iron and phytoplankton physiology at the LGM. Measured changes in carbon-13 and carbon-14, as well as a qualitative reconstruction of the change in ocean carbon export are used to evaluate the results. Overall, we find that while a reduction in ocean ventilation at the LGM is necessary to reproduce carbon-13 and carbon-14 observations, this circulation results in a low net sink for atmospheric CO2. In contrast, while biogeochemical processes contribute little to carbon isotopes, we propose that most of the change in atmospheric CO2 was due to such factors. However, the lesser role for circulation means that when all plausible factors are accounted for, most of the necessary CO2 change remains to be explained. This presents a serious challenge to our understanding of the mechanisms behind changes in the global carbon cycle during the geologic past.