Anthropogenic climate change is projected to lead to ocean warming, acidification, deoxygenation, reductions in near-surface nutrients, and changes to primary production, all of which are expected to ...affect marine ecosystems. Here we assess projections of these drivers of environmental change over the twenty-first century from Earth system models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) that were forced under the CMIP6 Shared Socioeconomic Pathways (SSPs). Projections are compared to those from the previous generation (CMIP5) forced under the Representative Concentration Pathways (RCPs). A total of 10 CMIP5 and 13 CMIP6 models are used in the two multi-model ensembles. Under the high-emission scenario SSP5-8.5, the multi-model global mean change (2080–2099 mean values relative to 1870–1899) ± the inter-model SD in sea surface temperature, surface pH, subsurface (100–600 m) oxygen concentration, euphotic (0–100 m) nitrate concentration, and depth-integrated primary production is +3.47±0.78 ∘C, -0.44±0.005, -13.27±5.28, -1.06±0.45 mmol m−3 and -2.99±9.11 %, respectively. Under the low-emission, high-mitigation scenario SSP1-2.6, the corresponding global changes are +1.42±0.32 ∘C, -0.16±0.002, -6.36±2.92, -0.52±0.23 mmol m−3, and -0.56±4.12 %. Projected exposure of the marine ecosystem to these drivers of ocean change depends largely on the extent of future emissions, consistent with previous studies. The ESMs in CMIP6 generally project greater warming, acidification, deoxygenation, and nitrate reductions but lesser primary production declines than those from CMIP5 under comparable radiative forcing. The increased projected ocean warming results from a general increase in the climate sensitivity of CMIP6 models relative to those of CMIP5. This enhanced warming increases upper-ocean stratification in CMIP6 projections, which contributes to greater reductions in upper-ocean nitrate and subsurface oxygen ventilation. The greater surface acidification in CMIP6 is primarily a consequence of the SSPs having higher associated atmospheric CO2 concentrations than their RCP analogues for the same radiative forcing. We find no consistent reduction in inter-model uncertainties, and even an increase in net primary production inter-model uncertainties in CMIP6, as compared to CMIP5.
The Australian Earth System Model: ACCESS-ESM1.5 Ziehn, Tilo; Chamberlain, Matthew A.; Law, Rachel M. ...
Journal of Southern Hemisphere earth systems science,
01/2020, Volume:
70, Issue:
1
Journal Article
Peer reviewed
Open access
The Australian Community Climate and Earth System Simulator (ACCESS) has been extended to include land and ocean carbon cycle components to form an Earth System Model (ESM). The current version, ...ACCESS-ESM1.5, has been mainly developed to enable Australia to participate in the Coupled Model Intercomparison Project Phase 6 (CMIP6) with an ESM version. Here we describe the model components and changes to the previous version, ACCESS-ESM1. We use the 500-year pre-industrial control run to highlight the stability of the physical climate and the carbon cycle. The long spin-up, negligible drift in temperature and small pre-industrial net carbon fluxes (0.02 and 0.08 PgC year−1 for land and ocean respectively) highlight the suitability of ACCESS-ESM1.5 to explore modes of variability in the climate system and coupling to the carbon cycle. The physical climate and carbon cycle for the present day have been evaluated using the CMIP6 historical simulation by comparing against observations and ACCESS-ESM1. Although there is generally little change in the climate simulation from the earlier model, many aspects of the carbon simulation are improved. An assessment of the climate response to CO2 forcing indicates that ACCESS-ESM1.5 has an equilibrium climate sensitivity of 3.87°C.
Decadal trends in the ocean carbon sink DeVries, Tim; LeQuéré, Corinne; Andrews, Oliver ...
Proceedings of the National Academy of Sciences - PNAS,
06/2019, Volume:
116, Issue:
24
Journal Article
Peer reviewed
Open access
Measurements show large decadal variability in the rate of CO₂ accumulation in the atmosphere that is not driven by CO₂ emissions. The decade of the 1990s experienced enhanced carbon accumulation in ...the atmosphere relative to emissions, while in the 2000s, the atmospheric growth rate slowed, even though emissions grew rapidly. These variations are driven by natural sources and sinks of CO₂ due to the ocean and the terrestrial biosphere. In this study, we compare three independent methods for estimating oceanic CO₂ uptake and find that the ocean carbon sink could be responsible for up to 40% of the observed decadal variability in atmospheric CO₂ accumulation. Data-based estimates of the ocean carbon sink from pCO₂ mapping methods and decadal ocean inverse models generally agree on the magnitude and sign of decadal variability in the ocean CO₂ sink at both global and regional scales. Simulations with ocean biogeochemical models confirm that climate variability drove the observed decadal trends in ocean CO₂ uptake, but also demonstrate that the sensitivity of ocean CO₂ uptake to climate variability may be too weak in models. Furthermore, all estimates point toward coherent decadal variability in the oceanic and terrestrial CO₂ sinks, and this variability is not well-matched by current global vegetation models. Reconciling these differences will help to constrain the sensitivity of oceanic and terrestrial CO₂ uptake to climate variability and lead to improved climate projections and decadal climate predictions.
Grazing dynamics are one of the most poorly constrained components of the marine carbon cycle. We use inverse modeling to infer the distribution of community‐integrated zooplankton grazing dynamics ...based on the ability of different grazing formulations to recreate the satellite‐observed seasonal cycle in phytoplankton biomass after controlling for physical and bottom‐up controls. We find large spatial variability in the optimal community‐integrated half saturation concentration for grazing (K1/2), with lower (higher) values required in more oligotrophic (eutrophic) biomes. This leads to a strong sigmoidal relationship between observed mean‐annual phytoplankton biomass and the optimally inferred grazing parameterization. This relationship can be used to help constrain, validate and/or parameterize next‐generation biogeochemical models.
Plain Language Summary
To improve predictions of the ocean's ability to feed a growing human population and buffer a changing climate, we need to improve our understanding of what happens to carbon once it is absorbed into the surface ocean. One of the largest knowledge gaps in marine carbon cycling is the role of zooplankton grazing. The rate at which zooplankton graze phytoplankton modifies the size and seasonal evolution of phytoplankton populations and in turn, the associated rates of net primary production at the base of the food‐web, secondary production of grazers (an indicator of fisheries potential) and export production (the biological sequestration of carbon). However, regional differences in grazing, which are difficult to measure outside of the laboratory, remain poorly constrained by observations and thus difficult to model. Here, we run a suite of model simulations, which each simulate grazing differently, then compare the results to infer which grazing dynamics best match observations. We find that there is dramatic spatial variability in how zooplankton, as a community, appear to be grazing and that this variability maps well onto observed phytoplankton abundance, suggesting that the type of zooplankton present may be determined by the amount of prey available.
Key Points
Oligotrophic (eutrophic) biomes exhibit more (less) efficient community‐integrated grazing, characteristic of micro‐ (meso‐) zooplankton
We find a strong link between observed mean‐annual phytoplankton biomass and the grazing dynamics required to recreate its seasonal cycle
A type III functional response typically does a better job recreating observed phytoplankton seasonal cycles than a type II response
Abstract
The Great Barrier Reef (GBR) is a globally significant coral reef system supporting productive and diverse ecosystems. The GBR is under increasing threat from climate change and local ...anthropogenic stressors, with its general condition degrading over recent decades. In response to this, a number of techniques have been proposed to offset or ameliorate environmental changes. In this study, we use a coupled hydrodynamic-biogeochemical model of the GBR and surrounding ocean to simulate artificial ocean alkalinisation (AOA) as a means to reverse the impact of global ocean acidification on GBR reefs. Our results demonstrate that a continuous release of 90 000 t of alkalinity every 3 d over one year along the entire length of the GBR, following the Gladstone-Weipa bulk carrier route, increases the mean aragonite saturation state (
Ω
a
r
) across the GBR’s 3860 reefs by 0.05. This change offsets just over 4 years (∼4.2) of ocean acidification under the present rate of anthropogenic carbon emissions. The injection raises
Ω
a
r
in the 250 reefs closest to the route by
⩾
0.15
, reversing further projected Ocean Acidification. Following cessation of alkalinity injection
Ω
a
r
returns to the value of the waters in the absence of AOA over a 6 month period, primarily due to transport of additional alkalinity into the Coral Sea. Significantly, our study provides for the first time a model of AOA applied along existing shipping infrastructure that has been used to investigate shelf scale impacts. Thus, amelioration of decades of OA on the GBR is feasible using existing infrastructure, but is likely to be extremely expensive, include as yet unquantified risks, and would need to be undertaken continuously until such time, probably centuries in the future, when atmospheric CO
2
concentrations have returned to today’s values.
Based on the 2019 assessment of the Global Carbon Project, the ocean took up on average, 2.5+/-0.6PgCyr-1 or 23+/-5% of the total anthropogenic CO2 emissions over the decade 2009-2018. This sink ...estimate is based on global ocean biogeochemical models (GOBMs) and is compared to data-products based on surface ocean pCO2 (partial pressure of CO2) observations accounting for the outgassing of river-derived CO2. Here we evaluate the GOBM simulations by comparing the simulated pCO2 to observations. The simulations are well suited for quantifying the global ocean carbon sink on the time-scale of the annual mean and its multi-decadal trend (RMSE <20 μatm), as well as on the time-scale of multi-year variability (RMSE <10 μatm), despite the large model-data mismatch on the seasonal time-scale (RMSE of 20-80 μatm). Biases in GOBMs have a small effect on the global mean ocean sink (0.05 PgC yr−1), but need to be addressed to improve the regional budgets and model-data comparison. Accounting for non-mapped areas in the data-products reduces their spread as measured by the standard deviation by a third. There is growing evidence and consistency among methods with regard to the patterns of the multi-year variability of the ocean carbon sink, with a global stagnation in the 1990s and an extra-tropical strengthening in the 2000s. GOBMs and data-products point consistently to a shift from a tropical CO2 source to a CO2 sink in recent years. On average, the GOBMs reveal less variations in the sink than the data-based products. Despite the reasonable simulation of surface ocean pCO2 by the GOBMs, there are discrepancies between the resulting sink estimate from GOBMs and data-products. These discrepancies are within the uncertainty of the river flux adjustment, increase over time, and largely stem from the Southern Ocean. Progress in our understanding of the global ocean carbon sink necessitates significant advancement in modelling and observing the Southern Ocean including (i) a game-changing increase in high-quality pCO2 observations, and (ii) a critical re-evaluation of the regional river flux adjustment.
Earth system models (ESMs) that incorporate carbon–climate feedbacks represent the present state of the art in climate modelling. Here, we describe the Australian Community Climate and Earth System ...Simulator (ACCESS)-ESM1, which comprises atmosphere (UM7.3), land (CABLE), ocean (MOM4p1), and sea-ice (CICE4.1) components with OASIS-MCT coupling, to which ocean and land carbon modules have been added. The land carbon model (as part of CABLE) can optionally include both nitrogen and phosphorous limitation on the land carbon uptake. The ocean carbon model (WOMBAT, added to MOM) simulates the evolution of phosphate, oxygen, dissolved inorganic carbon, alkalinity and iron with one class of phytoplankton and zooplankton. We perform multi-centennial pre-industrial simulations with a fixed atmospheric CO2 concentration and different land carbon model configurations (prescribed or prognostic leaf area index). We evaluate the equilibration of the carbon cycle and present the spatial and temporal variability in key carbon exchanges. Simulating leaf area index results in a slight warming of the atmosphere relative to the prescribed leaf area index case. Seasonal and interannual variations in land carbon exchange are sensitive to whether leaf area index is simulated, with interannual variations driven by variability in precipitation and temperature. We find that the response of the ocean carbon cycle shows reasonable agreement with observations. While our model overestimates surface phosphate values, the global primary productivity agrees well with observations. Our analysis highlights some deficiencies inherent in the carbon models and where the carbon simulation is negatively impacted by known biases in the underlying physical model and consequent limits on the applicability of this model version. We conclude the study with a brief discussion of key developments required to further improve the realism of our model simulation.
A biogeochemical ocean general circulation model, driven with NCEP‐R1 and observed atmospheric CO2 history, is used to investigate and quantify the role that the Southern Annular Mode (SAM), ...identified as the leading mode of climate variability, has in driving interannual variability in Southern Ocean air‐sea CO2 fluxes between 1980 and 2000. Our simulations show the Southern Ocean to be a region of decreased CO2 uptake during the positive SAM phase. The SAM induces changes in Southern Ocean CO2 uptake with a 2‐month time lag explaining 42% of the variance in the total interannual variability in air‐sea CO2 fluxes. Our analysis shows that the response of the Southern Ocean to the SAM is primarily governed by changes in ΔpCO2 (67%), and that this response is driven by changes in ocean physics that control the supply of nutrients to the upper ocean, primarily Dissolved Inorganic Carbon (DIC). The SAM is predicted to become stronger and more positive in response to climate change and our results suggest this will decrease the Southern Ocean CO2 uptake by 0.1PgC/yr per unit change in the SAM.
The Great Barrier Reef (GBR) is founded on reef-building corals. Corals build their exoskeleton with aragonite, but ocean acidification is lowering the aragonite saturation state of seawater (Ωa). ...The downscaling of ocean acidification projections from global to GBR scales requires the set of regional drivers controlling Ωa to be resolved. Here we use a regional coupled circulation-biogeochemical model and observations to estimate the Ωa experienced by the 3,581 reefs of the GBR, and to apportion the contributions of the hydrological cycle, regional hydrodynamics and metabolism on Ωa variability. We find more detail, and a greater range (1.43), than previously compiled coarse maps of Ωa of the region (0.4), or in observations (1.0). Most of the variability in Ωa is due to processes upstream of the reef in question. As a result, future decline in Ωa is likely to be steeper on the GBR than currently projected by the IPCC assessment report.