Search

Driven by climate change, marine biodiversity is undergoing a phase of rapid change that has proven to be even faster than changes observed in terrestrial ecosystems. Understanding how these changes in species composition will affect future marine life is crucial for conservation management, especially due to increasing demands for marine natural resources. Here, we analyse predictions of a multiparameter habitat suitability model covering the global projected ranges of >33,500 marine species from climate model projections under three CO2 emission scenarios (RCP2.6, RCP4.5, RCP8.5) up to the year 2100. Our results show that the core habitat area will decline for many species, resulting in a net loss of 50% of the core habitat area for almost half of all marine species in 2100 under the high-emission scenario RCP8.5. As an additional consequence of the continuing distributional reorganization of marine life, gaps around the equator will appear for 8% (RCP2.6), 24% (RCP4.5), and 88% (RCP8.5) of marine species with cross-equatorial ranges. For many more species, continuous distributional ranges will be disrupted, thus reducing effective population size. In addition, high invasion rates in higher latitudes and polar regions will lead to substantial changes in the ecosystem and food web structure, particularly regarding the introduction of new predators. Overall, our study highlights that the degree of spatial and structural reorganization of marine life with ensued consequences for ecosystem functionality and conservation efforts will critically depend on the realized greenhouse gas emission pathway.

Decreasing concentrations of dissolved oxygen in the ocean are considered one of the main threats to marine ecosystems as they jeopardize the growth of higher organisms. They also alter the marine nitrogen cycle, which is strongly bound to the carbon cycle and climate. While higher organisms in general start to suffer from oxygen concentrations < ∼ 63 µM (hypoxia), the marine nitrogen cycle responds to oxygen concentration below a threshold of about 20 µM (microbial hypoxia), whereas anoxic processes dominate the nitrogen cycle at oxygen concentrations of < ∼ 0.05 µM (functional anoxia). The Arabian Sea and the Bay of Bengal are home to approximately 21 % of the total volume of ocean waters revealing microbial hypoxia. While in the Arabian Sea this oxygen minimum zone (OMZ) is also functionally anoxic, the Bay of Bengal OMZ seems to be on the verge of becoming so. Even though there are a few isolated reports on the occurrence of anoxia prior to 1960, anoxic events have so far not been reported from the open northern Indian Ocean (i.e., other than on shelves) during the last 60 years. Maintenance of functional anoxia in the Arabian Sea OMZ with oxygen concentrations ranging between > 0 and ∼ 0.05 µM is highly extraordinary considering that the monsoon reverses the surface ocean circulation twice a year and turns vast areas of the Arabian Sea from an oligotrophic oceanic desert into one of the most productive regions of the oceans within a few weeks. Thus, the comparably low variability of oxygen concentration in the OMZ implies stable balances between the physical oxygen supply and the biological oxygen consumption, which includes negative feedback mechanisms such as reducing oxygen consumption at decreasing oxygen concentrations (e.g., reduced respiration). Lower biological oxygen consumption is also assumed to be responsible for a less intense OMZ in the Bay of Bengal. According to numerical model results, a decreasing physical oxygen supply via the inflow of water

Climate and marine biogeochemistry changes over the Holocene are investigated based on transient global climate and biogeochemistry model simulations over the last 9500 years. The simulations are forced by accelerated and non-accelerated orbital parameters, respectively, and atmospheric PCO2, CH4, and N2O. The analysis focusses on key climatic parameters of relevance to the marine biogeochemistry, and on the physical and biogeochemical processes that drive atmosphere-ocean carbon fluxes and changes in the oxygen minimum zones (OMZs). The simulated global mean ocean temperature is characterized by a mid-Holocene cooling and a late Holocene warming, a common feature among Holocene climate simulations which, however, contradicts a proxy-derived mid-Holocene climate optimum. As the most significant result, and only in the non-accelerated simulation, we find a substantial increase in volume of the OMZ in the eastern equatorial Pacific (EEP) continuing into the late Holocene. The concurrent increase in apparent oxygen utilization (AOU) and age of the water mass within the EEP OMZ can be attributed to a weakening of the deep northward inflow into the Pacific. This results in a large-scale mid-to-late Holocene increase in AOU in most of the Pacific and hence the source regions of the EEP OMZ waters. The simulated expansion of the EEP OMZ raises the question of whether the deoxygenation that has been observed over the last 5 decades could be a - perhaps accelerated - continuation of an orbitally driven decline in oxygen. Changes in global mean biological production and export of detritus remain of the order of 10 %, with generally lower values in the mid-Holocene. The simulated atmosphere-ocean CO2 flux would result in atmospheric pCO2 changes of similar magnitudes to those observed for the Holocene, but with different timing. More technically, as the increase in EEP OMZ volume can only be simulated with the non-accelerated model simulation, non-accelerated model simulations

© The Author(s), 2016. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Geoscientific Model Development 9 (2016): 1827-1851, doi:10.5194/gmd-9-1827-2016., During the fifth phase of the Coupled Model Intercomparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hindcasts were compared with available modern biogeochemical observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to-model uncertainty. This requires more attention in future model intercomparison exercises in order to provide quantitatively more correct ESM results on marine biogeochemistry and carbon cycle feedbacks., This work was supported by H2020 project CRESCENDO “Coordinated Research in Earth Systems and Climate: Experiments, kNowledge, Dissemination and Outreach”, which received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 641816 and by the EU FP7 project CARBOCHANGE “Changes in carbon uptake and emissions by oceans in a changing climate” which received funding from the European community’s Seventh Framework Programme under grant agreement no. 264879. Jerry Tjiputra acknowledges ORGANIC project (239965/F20) funded by the Research Council of Norway. Christoph Heinze and Jerry Tjiputra are grateful for support through project EVA – Earth system modelling of climate variations in the Anthropocene (229771/E10) funded by the Research Council of Norway, as well as CPU-time and mass storage provided through NOTUR project NN2345K as well as NorStore project NS2345K. Keith Lindsay and Scott C. Doney acknowledge support from the National Science Foundation.

We performed simulations with a global model of ocean biogeochemistry forced with orbitally driven anomalies of oceanic conditions for the mid-Holocene, known as Holocene climate optimum, to investigate natural variability in the eastern equatorial Pacific oxygen minimum zone (EEP OMZ). While the global mean temperature during the mid-Holocene was likely slightly higher than the 1961–1990 mean, the sea surface temperature in the EEP was slightly lower. Mid-Holocene oxygen concentrations in the EEP OMZ are generally increased, locally by up to 50%, and the EEP OMZ volume was, depending on definition of the OMZ threshold, at least 6% lower. These higher oxygen levels are the combined result of competing physical and biogeochemical processes. Our results imply that mechanisms for past changes in the EEP OMZ intensity and extension can differ from the global warming driven decline in oxygen levels observed for the recent decades and predicted for the future.

We estimate regional long-term surface ocean pCO(2) growth rates using all available underway and bottled biogeochemistry data collected over the past four decades. These observed regional trends are compared with those simulated by five state-of-the-art Earth system models over the historical period. Oceanic pCO(2) growth rates faster than the atmospheric growth rates indicate decreasing atmospheric CO2 uptake, while ocean pCO(2) growth rates slower than the atmospheric growth rates indicate increasing atmospheric CO2 uptake. Aside from the western subpolar North Pacific and the subtropical North Atlantic, our analysis indicates that the current observation-based basin-scale trends may be underestimated, indicating that more observations are needed to determine the trends in these regions. Encouragingly, good agreement between the simulated and observed pCO(2) trends is found when the simulated fields are subsampled with the observational coverage. In agreement with observations, we see that the simulated pCO(2) trends are primarily associated with the increase in surface dissolved inorganic carbon (DIC) associated with atmospheric carbon uptake, and in part by warming of the sea surface. Under the RCP8.5 future scenario, DIC continues to be the dominant driver of pCO(2) trends, with little change in the relative contribution of SST. However, the changes in the hydrological cycle play an increasingly important role. For the contemporary (1970-2011) period, the simulated regional pCO(2) trends are lower than the atmospheric growth rate over 90% of the ocean. However, by year 2100 more than 40% of the surface ocean area has a higher oceanic pCO(2) trend than the atmosphere, implying a reduction in the atmospheric CO2 uptake rate. The fastest pCO(2) growth rates are projected for the subpolar North Atlantic, while the high-latitude Southern Ocean and eastern equatorial Pacific have the weakest growth rates, remaining below the atmospheric pCO(2) growth rate. Our work

The marine CaCO3 cycle is an important component of the oceanic carbon system and directly affects the cycling of natural and the uptake of anthropogenic carbon. In numerical models of the marine carbon cycle, the CaCO3 cycle component is often evaluated against the observed distribution of alkalinity. Alkalinity varies in response to the formation and remineralization of CaCO3 and organic matter. However, it also has a large conservative component, which may strongly be affected by a deficient representation of ocean physics (circulation, evaporation, and precipitation) in models. Here we apply a global ocean biogeochemical model run into preindustrial steady state featuring a number of idealized tracers, explicitly capturing the model's CaCO3 dissolution, organic matter remineralization, and various preformed properties (alkalinity, oxygen, phosphate). We compare the suitability of a variety of measures related to the CaCO3 cycle, including alkalinity (TA), potential alkalinity and TA*, the latter being a measure of the time-integrated imprint of CaCO3 dissolution in the ocean. TA* can be diagnosed from any data set of TA, temperature, salinity, oxygen and phosphate. We demonstrate the sensitivity of total and potential alkalinity to the differences in model and ocean physics, which disqualifies them as accurate measures of biogeochemical processes. We show that an explicit treatment of preformed alkalinity (TA0) is necessary and possible. In our model simulations we implement explicit model tracers of TA0 and TA*. We find that the difference between modelled true TA* and diagnosed TA* was below 10% (25%) in 73% (81%) of the ocean's volume. In the Pacific (and Indian) Oceans the RMSE of A* is below 3 (4) mmol TA m−3, even when using a global rather than regional algorithms to estimate preformed alkalinity. Errors in the Atlantic Ocean are significantly larger and potential improvements of TA0 estimation are discussed. Applying the TA* approach to the output of thr

Decadal-to-century scale trends for a range of marine environmental variables in the upper mesopelagic layer (UML, 100–600 m) are investigated using results from seven Earth System Models forced by a high greenhouse gas emission scenario. The models as a class represent the observation-based distribution of oxygen (O2) and carbon dioxide (CO2), albeit major mismatches between observation-based and simulated values remain for individual models. By year 2100 all models project an increase in SST between 2 °C and 3 °C, and a decrease in the pH and in the saturation state of water with respect to calcium carbonate minerals in the UML. A decrease in the total ocean inventory of dissolved oxygen by 2% to 4% is projected by the range of models. Projected O2 changes in the UML show a complex pattern with both increasing and decreasing trends reflecting the subtle balance of different competing factors such as circulation, production, remineralization, and temperature changes. Projected changes in the total volume of hypoxic and suboxic waters remain relatively small in all models. A widespread increase of CO2 in the UML is projected. The median of the CO2 distribution between 100 and 600m shifts from 0.1–0.2 mol m−3 in year 1990 to 0.2–0.4 mol m−3 in year 2100, primarily as a result of the invasion of anthropogenic carbon from the atmosphere. The co-occurrence of changes in a range of environmental variables indicates the need to further investigate their synergistic impacts on marine ecosystems and Earth System feedbacks.

Climate modes such as the North Atlantic Oscillation (NAO), representing internal variability of the climate system, influence the ocean carbon cycle and may mask trends in the sink of anthropogenic carbon. Here, utilising control runs of six fully coupled Earth System Models, the response of the ocean carbon cycle to the NAO is quantified. The dominating response, a seesaw pattern between the subtropical gyre and the subpolar Northern Atlantic, is instantaneous (<3 months) and dynamically consistent over all models and with observations for a range of physical and biogeochemical variables. All models show asymmetric responses to NAO+ and NAO− forcing, implying non-linearity in the connection between NAO and the ocean carbon cycle. However, model differences in regional expression and magnitude and conflicting results with regard to air–sea flux and CO2 partial pressure remain. Typical NAO-driven variations are ±10 mmol/m3 in the surface concentration of dissolved inorganic carbon and alkalinity and ±8 ppm in the air–sea partial pressure difference. The effect on the basin-wide air–sea CO2 flux is small due to compensating fluxes on the sub-basin scale. Two models show a reduced carbon sink in the north-eastern North Atlantic during negative NAO phases, qualitatively in accordance with the observed decline during a phase of predominantly negative NAO. The results indicate that wind-driven dynamics are the main driver of the response to the NAO, which – via vertical mixing, upwelling and the associated entrainment of dissolved inorganic carbon and nutrients – leave an imprint on surface pCO2 and the air–sea CO2 flux as well as on biological export production, pH and the calcium carbonate saturation state. The biogeochemical response to the NAO is predominantly governed by vertical exchange between the surface and the thermocline; large-scale horizontal transport mechanisms are of minor importance.

Phosphate distributions simulated by seven state-of-the-art biogeochemical ocean circulation models are evaluated against observations of global ocean nutrient distributions. The biogeochemical models exhibit different structural complexities, ranging from simple nutrient-restoring to multi-nutrient NPZD type models. We evaluate the simulations using the observed volume distribution of phosphate. The errors in these simulated volume class distributions are significantly larger when preformed phosphate (or regenerated phosphate) rather than total phosphate is considered. Our analysis reveals that models can achieve similarly good fits to observed total phosphate distributions for a very different partitioning into preformed and regenerated nutrient components. This has implications for the strength and potential climate sensitivity of the simulated biological carbon pump. We suggest complementing the use of total nutrient distributions for assessing model skill by an evaluation of the respective preformed and regenerated nutrient components.