ALBERT

Description: 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.

Description: Results from the fully and biogeochemically coupled simulations in which CO2 increases at a rate of 1 % yr−1 (1pctCO2) from its preindustrial value are analyzed to quantify the magnitude of carbon–concentration and carbon–climate feedback parameters which measure the response of ocean and terrestrial carbon pools to changes in atmospheric CO2 concentration and the resulting change in global climate, respectively. The results are based on 11 comprehensive Earth system models from the most recent (sixth) Coupled Model Intercomparison Project (CMIP6) and compared with eight models from the fifth CMIP (CMIP5). The strength of the carbon–concentration feedback is of comparable magnitudes over land (mean ± standard deviation = 0.97 ± 0.40 PgC ppm−1) and ocean (0.79 ± 0.07 PgC ppm−1), while the carbon–climate feedback over land (−45.1 ± 50.6 PgC ∘C−1) is about 3 times larger than over ocean (−17.2 ± 5.0 PgC ∘C−1). The strength of both feedbacks is an order of magnitude more uncertain over land than over ocean as has been seen in existing studies. These values and their spread from 11 CMIP6 models have not changed significantly compared to CMIP5 models. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The transient climate response to cumulative emissions (TCRE) from the 11 CMIP6 models considered here is 1.77 ± 0.37 ∘C EgC−1 and is similar to that found in CMIP5 models (1.63 ± 0.48 ∘C EgC−1) but with somewhat reduced model spread. The expressions for feedback parameters based on the fully and biogeochemically coupled configurations of the 1pctCO2 simulation are simplified when the small temperature change in the biogeochemically coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters is used to gain insight into the reasons for differing responses among ocean and land carbon cycle models.

Description: Here we show results from Earth system model simulations from the marine cloud brightening experiment G4cdnc of the Geoengineering Model Intercomparison Project (GeoMIP). The nine contributing models prescribe a 50 % increase in the cloud droplet number concentration (CDNC) of low clouds over the global oceans in an experiment dubbed G4cdnc, with the purpose of counteracting the radiative forcing due to anthropogenic greenhouse gases under the RCP4.5 scenario. The model ensemble median effective radiative forcing (ERF) amounts to −1.9 W m−2, with a substantial inter-model spread of −0.6 to −2.5 W m−2. The large spread is partly related to the considerable differences in clouds and their representation between the models, with an underestimation of low clouds in several of the models. All models predict a statistically significant temperature decrease with a median of (for years 2020–2069) −0.96 [−0.17 to −1.21] K relative to the RCP4.5 scenario, with particularly strong cooling over low-latitude continents. Globally averaged there is a weak but significant precipitation decrease of −2.35 [−0.57 to −2.96] % due to a colder climate, but at low latitudes there is a 1.19 % increase over land. This increase is part of a circulation change where a strong negative top-of-atmosphere (TOA) shortwave forcing over subtropical oceans, caused by increased albedo associated with the increasing CDNC, is compensated for by rising motion and positive TOA longwave signals over adjacent land regions.

Description: Geoengineering, or climate intervention, describes methods of deliberately altering the climate system to offset anthropogenic climate change. As an idealized representation of near-surface solar geoengineering over the ocean, such as marine cloud brightening, this paper discusses experiment G1ocean-albedo of the Geoengineering Model Intercomparison Project (GeoMIP), involving an abrupt quadrupling of the CO2 concentration and an instantaneous increase in ocean albedo to maintain approximate net top-of-atmosphere radiative flux balance. A total of 11 Earth system models are relatively consistent in their temperature, radiative flux, and hydrological cycle responses to this experiment. Due to the imposed forcing, air over the land surface warms by a model average of 1.14 K, while air over most of the ocean cools. Some parts of the near-surface air temperature over ocean warm due to heat transport from land to ocean. These changes generally resolve within a few years, indicating that changes in ocean heat content play at most a small role in the warming over the oceans. The hydrological cycle response is a general slowing down, with high heterogeneity in the response, particularly in the tropics. While idealized, these results have important implications for marine cloud brightening, or other methods of geoengineering involving spatially heterogeneous forcing, or other general forcings with a strong land–ocean contrast. It also reinforces previous findings that keeping top-of-atmosphere net radiative flux constant is not sufficient for preventing changes in global mean temperature.

Description: We examine extreme temperature and precipitation under two potential geoengineering methods forming part of the Geoengineering Model Intercomparison Project (GeoMIP). The solar dimming experiment G1 is designed to completely offset the global mean radiative forcing due to a CO2-quadrupling experiment (abrupt4 × CO2), while in GeoMIP experiment G4, the radiative forcing due to the representative concentration pathway 4.5 (RCP4.5) scenario is partly offset by a simulated layer of aerosols in the stratosphere. Both G1 and G4 geoengineering simulations lead to lower minimum temperatures (TNn) at higher latitudes and on land, primarily through feedback effects involving high-latitude processes such as snow cover, sea ice and soil moisture. There is larger cooling of TNn and maximum temperatures (TXx) over land compared with oceans, and the land–sea cooling contrast is larger for TXx than TNn. Maximum 5-day precipitation (Rx5day) increases over subtropical oceans, whereas warm spells (WSDI) decrease markedly in the tropics, and the number of consecutive dry days (CDDs) decreases in most deserts. The precipitation during the tropical cyclone (hurricane) seasons becomes less intense, whilst the remainder of the year becomes wetter. Stratospheric aerosol injection is more effective than solar dimming in moderating extreme precipitation (and flooding). Despite the magnitude of the radiative forcing applied in G1 being ∼ 7.7 times larger than in G4 and despite differences in the aerosol chemistry and transport schemes amongst the models, the two types of geoengineering show similar spatial patterns in normalized differences in extreme temperatures changes. Large differences mainly occur at northern high latitudes, where stratospheric aerosol injection more effectively reduces TNn and TXx. While the pattern of normalized differences in extreme precipitation is more complex than that of extreme temperatures, generally stratospheric aerosol injection is more effective in reducing tropical Rx5day, while solar dimming is more effective over extra-tropical regions.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the "global carbon budget" – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of our imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2007–2016), EFF was 9.4 ± 0.5 GtC yr−1, ELUC 1.3 ± 0.7 GtC yr−1, GATM 4.7 ± 0.1 GtC yr−1, SOCEAN 2.4 ± 0.5 GtC yr−1, and SLAND 3.0 ± 0.8 GtC yr−1, with a budget imbalance BIM of 0.6 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1. Also for 2016, ELUC was 1.3 ± 0.7 GtC yr−1, GATM was 6.1 ± 0.2 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1 and SLAND was 2.7 ± 1.0 GtC yr−1, with a small BIM of −0.3 GtC. GATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the higher fossil emissions and smaller SLAND for that year consistent with El Niño conditions. The global atmospheric CO2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data indicate a renewed growth in EFF of +2.0 % (range of 0.8 % to 3.0 %) based on national emissions projections for China, USA, and India, and projections of Gross Domestic Product corrected for recent changes in the carbon intensity of the economy for the rest of the world. For 2017, initial data indicate an increase in atmospheric CO2 concentration of around 5.3 GtC (2.5 ppm), attributed to a combination of increasing emissions and receding El Niño conditions. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set (Le Quéré et al., 2016; 2015b; 2015a; 2014; 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2017.

Description: 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.

Description: Over the last decade many climate models have evolved into Earth system models (ESMs), which are able to simulate both physical and biogeochemical processes through the inclusion of additional components such as the carbon cycle. The Australian Community Climate and Earth System Simulator (ACCESS) has been recently extended to include land and ocean carbon cycle components in its ACCESS-ESM1 version. A detailed description of ACCESS-ESM1 components including results from pre-industrial simulations is provided in Part 1. Here, we focus on the evaluation of ACCESS-ESM1 over the historical period (1850–2005) in terms of its capability to reproduce climate and carbon-related variables. Comparisons are performed with observations, if available, but also with other ESMs to highlight common weaknesses. We find that climate variables controlling the exchange of carbon are well reproduced. However, the aerosol forcing in ACCESS-ESM1 is somewhat larger than in other models, which leads to an overly strong cooling response in the land from about 1960 onwards. The land carbon cycle is evaluated for two scenarios: running with a prescribed leaf area index (LAI) and running with a prognostic LAI. We overestimate the seasonal mean (1.7 vs. 1.4) and peak amplitude (2.0 vs. 1.8) of the prognostic LAI at the global scale, which is common amongst CMIP5 ESMs. However, the prognostic LAI is our preferred choice, because it allows for the vegetation feedback through the coupling between LAI and the leaf carbon pool. Our globally integrated land–atmosphere flux over the historical period is 98 PgC for prescribed LAI and 137 PgC for prognostic LAI, which is in line with estimates of land use emissions (ACCESS-ESM1 does not include land use change). The integrated ocean–atmosphere flux is 83 PgC, which is in agreement with a recent estimate of 82 PgC from the Global Carbon Project for the period 1959–2005. The seasonal cycle of simulated atmospheric CO2 is close to the observed seasonal cycle (up to 1 ppm difference for the station at Mace Head and up to 2 ppm for the station at Mauna Loa), but shows a larger amplitude (up to 6 ppm) in the high northern latitudes. Overall, ACCESS-ESM1 performs well over the historical period, making it a useful tool to explore the change in land and oceanic carbon uptake in the future.

Description: Atmospheric CO2 levels continue to rise, increasing the risk of severe impacts on the Earth system, and on the ecosystem services that it provides. Artificial Ocean Alkalization (AOA) is capable of reducing atmospheric CO2 concentrations, surface warming and addressing ocean acidification. Here we simulate global and regional responses to alkalinity addition (0.25 PmolAlk/year) using the CSIRO-Mk3L-COAL Earth System Model in the period 2020–2100, under high (RCP8.5) and low (RCP2.6) emissions. While regionally there are large changes associated with locations of AOA, globally we see only a very weak dependence on where and when AOA is applied. We see that under RCP2.6, while the carbon uptake associated with AOA is only ~ 60 % of the total under RCP8.5, the relative changes in temperature are larger, as are the changes in pH (1.4×) and aragonite saturation (1.7×). The results of this modelling study are significant as they demonstrate that AOA is more effective under lower emissions, and the higher the emissions the more AOA required to achieve the same reduction in global warming and ocean acidification. Finally, our simulations show AOA in the period 2020–2100 is capable of offsetting global warming and ameliorating ocean acidification increases due to low emissions, but regionally the response is more variable.

Description: Carbon-climate feedbacks have the potential to significantly impact the future climate by altering atmospheric CO2 concentrations (Zaehle et al., 2010). By modifying the future atmospheric CO2 concentrations, the carbon-climate feedbacks will also influence the future trajectory for ocean acidification. Here, we use the CO2 emissions scenarios from 4 Representative Concentration Pathways (RCPs) with an Earth System Model to project the future trajectories of ocean acidification with the inclusion of carbon-climate feedbacks. We show that simulated carbon-climate feedbacks can significantly impact the onset of under-saturated aragonite conditions in the Southern and Arctic Oceans, the suitable habitat for tropical coral and the deepwater saturation states. Under higher emission scenarios (RCP8.5 and RCP6.0), the carbon-climate feedbacks advance the onset of under-saturation conditions and the reduction in suitable coral reef habitat by a decade or more. The impact of the carbon-climate feedback is most significant for the medium (RCP4.5) and low emission (RCP2.6) scenarios. For RCP4.5 scenario by 2100, the carbon-climate feedbacks nearly double the area of surface water under-saturated respect to aragonite and reduce by 50 % the surface water suitable for coral reefs. For RCP2.6 scenario by 2100, the carbon-climate feedbacks reduce the area suitable for coral reefs by 40 % and increase the area of under-saturated surface water by 20 %. The high sensitivity of the impact of ocean acidification to the carbon-climate feedbacks in the low to medium emissions scenarios is important because our recent commitments to reduce CO2 emissions are trying to move us on to such an emissions scenario. The study highlights the need to better characterise the carbon-climate feedbacks to ensure we do not excessively stress the oceans by under-estimating the future impact of ocean acidification.