Rapid emergence of climate change in environmental drivers of marine ecosystems

Climate change is expected to modify ecological responses in the ocean, with the potential for important effects on the ecosystem services provided to humankind. Here we address the question of how rapidly multiple drivers of marine ecosystem change develop in the future ocean. By analysing an ensemble of models we find that, within the next 15 years, the climate change-driven trends in multiple ecosystem drivers emerge from the background of natural variability in 55% of the ocean and propagate rapidly to encompass 86% of the ocean by 2050 under a ‘business-as-usual’ scenario. However, we also demonstrate that the exposure of marine ecosystems to climate change-induced stress can be drastically reduced via climate mitigation measures; with mitigation, the proportion of ocean susceptible to multiple drivers within the next 15 years is reduced to 34%. Mitigation slows the pace at which multiple drivers emerge, allowing an additional 20 years for adaptation in marine ecological and socio-economic systems alike

Decadal trends in the ocean carbon sink

Measurements show large decadal variability in the rate of [Formula: see text] accumulation in the atmosphere that is not driven by [Formula: see text] 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 [Formula: see text] due to the ocean and the terrestrial biosphere. In this study, we compare three independent methods for estimating oceanic [Formula: see text] uptake and find that the ocean carbon sink could be responsible for up to 40% of the observed decadal variability in atmospheric [Formula: see text] accumulation. Data-based estimates of the ocean carbon sink from [Formula: see text] mapping methods and decadal ocean inverse models generally agree on the magnitude and sign of decadal variability in the ocean [Formula: see text] sink at both global and regional scales. Simulations with ocean biogeochemical models confirm that climate variability drove the observed decadal trends in ocean [Formula: see text] uptake, but also demonstrate that the sensitivity of ocean [Formula: see text] uptake to climate variability may be too weak in models. Furthermore, all estimates point toward coherent decadal variability in the oceanic and terrestrial [Formula: see text] 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 [Formula: see text] uptake to climate variability and lead to improved climate projections and decadal climate predictions

Irreversible loss in marine ecosystem habitability after a temperature overshoot

Anthropogenic warming of the oceans and associated deoxygenation are altering marine ecosystems. Current knowledge suggests these changes may be reversible on a centennial timescale at the ocean surface but irreversible at deeper depths even if global warming were to ameliorate. In contrast, the marine ecosystem’s response to these persistent changes remains poorly elucidated. Here we explore to what extent global warming may drive alterations in marine habitats by exploring the evolution of a metabolic index that captures marine organisms’ ecophysiological response to both temperature and oxygen changes, throughout an idealised ramp-up/ramp-down atmospheric carbon dioxide concentration and an overshoot scenarios. Using a multi-model approach; we find that changes in ocean temperature and oxygen drive a centuries-long irreversible loss in the habitable volume of the upper 1000 m of the world ocean. These results suggest that the combined effect of warming and deoxygenation will have profound and long-lasting impacts on the viability of marine ecosystems, well after global temperatures have peaked.publishedVersio

Evaluation of global fire simulations in CMIP6 Earth system models

Fire is the primary form of terrestrial ecosystem disturbance on a global scale and an important Earth system process. Most Earth system models (ESMs) have incorporated fire modeling, with 19 out of them submitting model outputs of fire-related variables to the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study provides the first comprehensive evaluation of CMIP6 historical fire simulations by comparing them with multiple satellite-based products and charcoal-based historical reconstructions. Our results show that most CMIP6 models simulate the present-day global burned area and fire carbon emissions within the range of satellite-based products. They also capture the major features of observed spatial patterns and seasonal cycles, the relationship of fires with precipitation and population density, and the influence of El Niño-Southern Oscillation (ENSO) on the interannual variability of tropical fires. Regional fire carbon emissions simulated by the CMIP6 models from 1850 to 2010 generally align with the charcoal-based reconstructions, although there are regional mismatches, such as in southern South America and eastern temperate North America prior to the 1910s and in temperate North America, eastern boreal North America, Europe, and boreal Asia since the 1980s. The CMIP6 simulations have addressed three critical issues identified in the CMIP5: (1) the simulated global burned area less than half of the observations, (2) the failure to reproduce the high burned area fraction observed in Africa, and (3) the weak fire seasonal variability. Furthermore, the CMIP6 models exhibit improved accuracy in capturing the observed relationship between fires and both climatic and socioeconomic drivers, and better align with the historical long-term trends indicated by charcoal-based reconstructions in most regions worldwide. However, the CMIP6 models still fail to reproduce the decline in global burned area and fire carbon emissions observed over the past two decades, mainly attributed to an underestimation of anthropogenic fire suppression, and the spring peak in fires in the Northern Hemisphere mid-latitudes, mainly due to an underestimation of crop fires. In addition, the model underestimates the fire sensitivity to wet-dry conditions, indicating the need to improve fuel wetness estimation. Based on these findings, we present specific guidance for fire scheme development and suggest the post-processing methodology for using CMIP6 multi-model outputs to generate reliable fire projection products

Consistency and Challenges in the Ocean Carbon Sink Estimate for the Global Carbon Budget

Based on the 2019 assessment of the Global Carbon Project, the ocean took up on average, 2.5 ± 0.6 PgC yr−1 or 23 ± 5% of the total anthropogenic CO2 emissions over the decade 2009–2018. This sink estimate is based on simulation results from global ocean biogeochemical models (GOBMs) and is compared to data-products based on observations of surface ocean pCO2 (partial pressure of CO2) accounting for the outgassing of river-derived CO2. Here we evaluate the GOBM simulations by comparing the simulated surface ocean pCO2 to observations. Based on this comparison, 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 modeling and observing the Southern Ocean carbon sink including (i) a game-changing increase in high-quality pCO2 observations, and (ii) a critical re-evaluation of the regional river flux adjustment

The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) contribution to C4MIP: Quantifying committed climate changes following zero carbon emissions

The amount of additional future temperature change following a complete cessation of CO2 emissions is a measure of the unrealized warming to which we are committed due to CO2 already emitted to the atmosphere. This “zero emissions commitment” (ZEC) is also an important quantity when estimating the remaining carbon budget – a limit on the total amount of CO2 emissions consistent with limiting global mean temperature at a particular level. In the recent IPCC Special Report on Global Warming of 1.5 ∘C, the carbon budget framework used to calculate the remaining carbon budget for 1.5 ∘C included the assumption that the ZEC due to CO2 emissions is negligible and close to zero. Previous research has shown significant uncertainty even in the sign of the ZEC. To close this knowledge gap, we propose the Zero Emissions Commitment Model Intercomparison Project (ZECMIP), which will quantify the amount of unrealized temperature change that occurs after CO2 emissions cease and investigate the geophysical drivers behind this climate response. Quantitative information on ZEC is a key gap in our knowledge, and one that will not be addressed by currently planned CMIP6 simulations, yet it is crucial for verifying whether carbon budgets need to be adjusted to account for any unrealized temperature change resulting from past CO2 emissions. We request only one top-priority simulation from comprehensive general circulation Earth system models (ESMs) and Earth system models of intermediate complexity (EMICs) – a branch from the 1 % CO2 run with CO2 emissions set to zero at the point of 1000 PgC of total CO2 emissions in the simulation – with the possibility for additional simulations, if resources allow. ZECMIP is part of CMIP6, under joint sponsorship by C4MIP and CDRMIP, with associated experiment names to enable data submissions to the Earth System Grid Federation. All data will be published and made freely available

Inconsistent strategies to spin up models in CMIP5: implications for ocean biogeochemical model performance assessment

International audienceDuring the fifth phase of the Coupled Model Inter-comparison 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 hind-casts were compared with available modern biogeochemi-cal 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

The Zero Emissions Commitment and climate stabilization

How do we halt global warming? Reaching net zero carbon dioxide (CO2) emissions is understood to be a key milestone on the path to a safer planet. But how confident are we that when we stop carbon emissions, we also stop global warming? The Zero Emissions Commitment (ZEC) quantifies how much warming or cooling we can expect following a complete cessation of anthropogenic CO2 emissions. To date, the best estimate by the Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report is zero change, though with substantial uncertainty. In this article, we present an overview of the changes expected in major Earth system processes after net zero and their potential impact on global surface temperature, providing an outlook toward building a more confident assessment of ZEC in the decades to come. We propose a structure to guide research into ZEC and associated changes in the climate, separating the impacts expected over decades, centuries, and millennia. As we look ahead at the century billed to mark the end of net anthropogenic CO2 emissions, we ask: what is the prospect of a stable climate in a post-net zero world

Evaluation of ocean dimethylsulfide concentration and emission in CMIP6 models

Characteristics and trends of surface ocean dimethylsulfide (DMS) concentrations and fluxes into the atmosphere of four Earth system models (ESMs: CNRM-ESM2-1, MIROC-ES2L, NorESM2-LM, and UKESM1-0-LL) are analysed over the recent past (1980–2009) and into the future, using Coupled Model Intercomparison Project 6 (CMIP6) simulations. The DMS concentrations in historical simulations systematically underestimate the most widely used observed climatology but compare more favourably against two recent observation-based datasets. The models better reproduce observations in mid to high latitudes, as well as in polar and westerlies marine biomes. The resulting multi-model estimate of contemporary global ocean DMS emissions is 16–24 Tg S yr−1, which is narrower than the observational-derived range of 16 to 28 Tg S yr−1. The four models disagree on the sign of the trend of the global DMS flux from 1980 onwards, with two models showing an increase and two models a decrease. At the global scale, these trends are dominated by changes in surface DMS concentrations in all models, irrespective of the air–sea flux parameterisation used. In turn, three models consistently show that changes in DMS concentrations are correlated with changes in marine productivity; however, marine productivity is poorly constrained in the current generation of ESMs, thus limiting the predictive ability of this relationship. In contrast, a consensus is found among all models over polar latitudes where an increasing trend is predominantly driven by the retreating sea-ice extent. However, the magnitude of this trend between models differs by a factor of 3, from 2.9 to 9.2 Gg S decade−1 over the period 1980–2014, which is at the low end of a recent satellite-derived analysis. Similar increasing trends are found in climate projections over the 21st century

Assessing Model Predictions of Carbon Dynamics in Global Drylands

Drylands cover ca. 40% of the land surface and are hypothesised to play a major role in the global carbon cycle, controlling both long-term trends and interannual variation. These insights originate from land surface models (LSMs) that have not been extensively calibrated and evaluated for water-limited ecosystems. We need to learn more about dryland carbon dynamics, particularly as the transitory response and rapid turnover rates of semi-arid systems may limit their function as a carbon sink over multi-decadal scales. We quantified aboveground biomass carbon (AGC; inferred from SMOS L-band vegetation optical depth) and gross primary productivity (GPP; from PML-v2 inferred from MODIS observations) and tested their spatial and temporal correspondence with estimates from the TRENDY ensemble of LSMs. We found strong correspondence in GPP between LSMs and PML-v2 both in spatial patterns (Pearson’s r = 0.9 for TRENDY-mean) and in inter-annual variability, but not in trends. Conversely, for AGC we found lesser correspondence in space (Pearson’s r = 0.75 for TRENDY-mean, strong biases for individual models) and in the magnitude of inter-annual variability compared to satellite retrievals. These disagreements likely arise from limited representation of ecosystem responses to plant water availability, fire, and photodegradation that drive dryland carbon dynamics. We assessed inter-model agreement and drivers of long-term change in carbon stocks over centennial timescales. This analysis suggested that the simulated trend of increasing carbon stocks in drylands is in soils and primarily driven by increased productivity due to CO$_2$ enrichment. However, there is limited empirical evidence of this 50-year sink in dryland soils. Our findings highlight important uncertainties in simulations of dryland ecosystems by current LSMs, suggesting a need for continued model refinements and for greater caution when interpreting LSM estimates with regards to current and future carbon dynamics in drylands and by extension the global carbon cycle