Victoria Earth System Climate Model Research Articles

Impact of plastic pollution on atmospheric carbon dioxide

Since the beginning of its large-scale production in the early 20th century, plastics have remained an important material in widespread use throughout modern society. Nevertheless, despite possessing many benefits, plastics are resistant to degradation and instead accumulate in the ocean and terrestrial sediments, thereby potentially affecting marine and terrestrial ecosystems. Plastics release CO2 throughout their entire lifecycle; during the extraction of materials used in their production, through plastic–carbon leaching in the marine and terrestrial environment, and during their different end-of-life scenarios, which include recycling, landfill, and incineration. Here, we use the University of Victoria earth system climate model to quantity the effects on atmospheric CO2 and the ocean carbon cycle by using upper-bound estimates of carbon emissions from marine plastic–carbon leaching or land-based incineration. Despite the suggestions of some, our results indicate that it has only a very minor influence and an insignificant effect on the earth's global climate system. This holds even if plastic contamination increases well beyond current levels. On the other hand, carbon emissions associated with plastic production and incineration have a greater impact on climate while still dwarfed by emissions associated with the combustion of fossil fuels (coal, oil, and natural gas) and other anthropogenic sources. Our results have important policy implications for ongoing United Nations Environment Programme Intergovernmental Negotiating Committee on Plastic Pollution negotiations.

Explicit silicate cycling in the Kiel Marine Biogeochemistry Model version 3 (KMBM3) embedded in the UVic ESCM version 2.9

Abstract. We describe and test a new model of biological marine silicate cycling, implemented in the Kiel Marine Biogeochemical Model version 3 (KMBM3), embedded in the University of Victoria Earth System Climate Model (UVic ESCM) version 2.9. This new model adds diatoms, which are a key component of the biological carbon pump, to an existing ecosystem model. This new model combines previously published parameterizations of a diatom functional type, opal production and export with a novel, temperature-dependent dissolution scheme. Modelled steady-state biogeochemical rates, carbon and nutrient distributions are similar to those found in previous model versions. The new model performs well against independent ocean biogeochemical indicators and captures the large-scale features of the marine silica cycle to a degree comparable to similar Earth system models. Furthermore, it is computationally efficient, allowing both fully coupled, long-timescale transient simulations and “offline” transport matrix spinups. We assess the fully coupled model against modern ocean observations, the historical record starting from 1960 and a business-as-usual atmospheric CO2 forcing to the year 2300. The model simulates a global decline in net primary production (NPP) of 1.4 % having occurred since the 1960s, with the strongest declines in the tropics, northern midlatitudes and Southern Ocean. The simulated global decline in NPP reverses after the year 2100 (forced by the extended RCP8.5 CO2 concentration scenario), and NPP returns to 98 % of the pre-industrial rate by 2300. This recovery is dominated by increasing primary production in the Southern Ocean, mostly by calcifying phytoplankton. Large increases in calcifying phytoplankton in the Southern Ocean offset a decline in the low latitudes, producing a global net calcite export in 2300 that varies only slightly from pre-industrial rates. Diatom distribution moves southward in our simulations, following the receding Antarctic ice front, but diatoms are outcompeted by calcifiers across most of their pre-industrial Southern Ocean habitat. Global opal export production thus drops to 75 % of its pre-industrial value by 2300. Model nutrients such as phosphate, silicate and nitrate build up along the Southern Ocean particle export pathway, but dissolved iron (for which ocean sources are held constant) increases in the upper ocean. This different behaviour of iron is attributed to a reduction of low-latitude NPP (and consequently, a reduction in both uptake and export and particle, including calcite scavenging), an increase in seawater temperatures (raising the solubility of particulate iron) and stratification that “traps” the iron near the surface. These results are meant to serve as a baseline for sensitivity assessments to be undertaken with this model in the future.

Optimality-based non-Redfield plankton–ecosystem model (OPEM v1.1) in UVic-ESCM 2.9 – Part 2: Sensitivity analysis and model calibration

Abstract. We analyse 400 perturbed-parameter simulations for two configurations of an optimality-based plankton–ecosystem model (OPEM), implemented in the University of Victoria Earth System Climate Model (UVic-ESCM), using a Latin hypercube sampling method for setting up the parameter ensemble. A likelihood-based metric is introduced for model assessment and selection of the model solutions closest to observed distributions of NO3-, PO43-, O2, and surface chlorophyll a concentrations. The simulations closest to the data with respect to our metric exhibit very low rates of global N2 fixation and denitrification, indicating that in order to achieve rates consistent with independent estimates, additional constraints have to be applied in the calibration process. For identifying the reference parameter sets, we therefore also consider the model's ability to represent current estimates of water-column denitrification. We employ our ensemble of model solutions in a sensitivity analysis to gain insights into the importance and role of individual model parameters as well as correlations between various biogeochemical processes and tracers, such as POC export and the NO3- inventory. Global O2 varies by a factor of 2 and NO3- by more than a factor of 6 among all simulations. Remineralisation rate is the most important parameter for O2, which is also affected by the subsistence N quota of ordinary phytoplankton (Q0,phyN) and zooplankton maximum specific ingestion rate. Q0,phyN is revealed as a major determinant of the oceanic NO3- pool. This indicates that unravelling the driving forces of variations in phytoplankton physiology and elemental stoichiometry, which are tightly linked via Q0,phyN, is a prerequisite for understanding the marine nitrogen inventory.

Evaluation of the University of Victoria Earth System Climate Model version 2.10 (UVic ESCM 2.10)

Abstract. The University of Victoria Earth System Climate Model (UVic ESCM) of intermediate complexity has been a useful tool in recent assessments of long-term climate changes, including both paleo-climate modelling and uncertainty assessments of future warming. Since the last official release of the UVic ESCM 2.9 and the two official updates during the last decade, considerable model development has taken place among multiple research groups. The new version 2.10 of the University of Victoria Earth System Climate Model presented here will be part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). More precisely it will be used in the intercomparison of Earth system models of intermediate complexity (EMIC), such as the C4MIP, the Carbon Dioxide Removal and Zero Emissions Commitment model intercomparison projects (CDR-MIP and ZECMIP, respectively). It now brings together and combines multiple model developments and new components that have come about since the last official release of the model. The main additions to the base model are (i) an improved biogeochemistry module for the ocean, (ii) a vertically resolved soil model including dynamic hydrology and soil carbon processes, and (iii) a representation of permafrost carbon. To set the foundation of its use, we here describe the UVic ESCM 2.10 and evaluate results from transient historical simulations against observational data. We find that the UVic ESCM 2.10 is capable of reproducing changes in historical temperature and carbon fluxes well. The spatial distribution of many ocean tracers, including temperature, salinity, phosphate and nitrate, also agree well with observed tracer profiles. The good performance in the ocean tracers is connected to an improved representation of ocean physical properties. For the moment, the main biases that remain are a vegetation carbon density that is too high in the tropics, a higher than observed change in the ocean heat content (OHC) and an oxygen utilization in the Southern Ocean that is too low. All of these biases will be addressed in the next updates to the model.

Irreversible ocean thermal expansion under carbon dioxide removal

Abstract. In the Paris Agreement in 2015 countries agreed on holding global mean surface air warming to well below 2 °C above pre-industrial levels, but the emission reduction pledges under that agreement are not ambitious enough to meet this target. Therefore, the question arises of whether restoring global warming to this target after exceeding it by artificially removing CO2 from the atmosphere is possible. One important aspect is the reversibility of ocean heat uptake and associated sea level rise, which have very long (centennial to millennial) response timescales. In this study the response of sea level rise due to thermal expansion to a 1 % yearly increase of atmospheric CO2 up to a quadrupling of the pre-industrial concentration followed by a 1 % yearly decline back to the pre-industrial CO2 concentration is examined using the University of Victoria Earth System Climate Model (UVic ESCM). We find that global mean thermosteric sea level (GMTSL) continues to rise for several decades after atmospheric CO2 starts to decline and does not return to pre-industrial levels for over 1000 years after atmospheric CO2 is restored to the pre-industrial concentration. This finding is independent of the strength of vertical sub-grid-scale ocean mixing implemented in the model. Furthermore, GMTSL rises faster than it declines in response to a symmetric rise and decline in atmospheric CO2 concentration partly because the deep ocean continues to warm for centuries after atmospheric CO2 returns to the pre-industrial concentration. Both GMTSL rise and decline rates increase with increasing vertical ocean mixing. Exceptions from this behaviour arise if the overturning circulations in the North Atlantic and Southern Ocean intensify beyond pre-industrial levels in model versions with lower vertical mixing, which leads to rapid cooling of the deep ocean.

Climate and health implications of future aerosol emission scenarios

Anthropogenic aerosols have a net cooling effect on climate and also cause adverse health effects by degrading air quality. In this global-scale sensitivity study, we used a combination of the aerosol-climate model ECHAM-HAMMOZ and the University of Victoria Earth System Climate Model to assess the climate and health effects of aerosols emissions from three Representative Concentration Pathways (RCP2.6, RCP4.5, and RCP8.5) and two new (LOW and HIGH) aerosol emission scenarios derived from RCP4.5, but that span a wider spectrum of possible future aerosol emissions. All simulations had CO2 emissions and greenhouse gas forcings from RCP4.5. Aerosol forcing declined similarly in the standard RCP aerosol emission scenarios: the aerosol effective radiative forcing (ERF) decreased from −1.3 W m−2 in 2005 to between −0.1 W m−2 and −0.4 W m−2 in 2100. The differences in ERF were substantially larger between LOW (−0.02 W m−2 in 2100) and HIGH (−0.8 W m−2) scenarios. The global mean temperature difference between the simulations with standard RCP aerosol emissions was less than 0.18 °C, whereas the difference between LOW and HIGH reached 0.86 °C in 2061. In LOW, the rate of warming peaked at 0.48 °C per decade in the 2030s, whereas in HIGH it was the lowest of all simulations and never exceeded 0.23 °C per decade. Using present-day population density and baseline mortality rates for all scenarios, PM2.5-induced premature mortality was 2 371 800 deaths per year in 2010 and 525 700 in 2100 with RCP4.5 aerosol emissions; in HIGH, the premature mortality reached its maximum value of 2 780 800 deaths per year in 2030, whereas in LOW the premature mortality at 2030 was below 299 900 deaths per year. Our results show potential trade-offs in aerosol mitigation with respect to climate change and public health as ambitious reduction of aerosol emissions considerably increased warming while decreasing mortality.

The importance of terrestrial weathering changes in multimillennial recovery of the global carbon cycle: a two-dimensional perspective

Abstract. In this paper, we describe the development and application of a new spatially explicit weathering scheme within the University of Victoria Earth System Climate Model (UVic ESCM). We integrated a dataset of modern-day lithology with a number of previously devised parameterizations for weathering dependency on temperature, primary productivity, and runoff. We tested the model with simulations of future carbon cycle perturbations, comparing a number of emission scenarios and model versions with each other and with zero-dimensional equivalents of each experiment. Overall, we found that our two-dimensional weathering model versions were more efficient in restoring the carbon cycle to its pre-industrial state following the pulse emissions than their zero-dimensional counterparts; however, in either case the effect of this weathering negative feedback on the global carbon cycle was small on timescales of less than 1000 years. According to model results, the largest contribution to future changes in weathering rates came from the expansion of tropical and mid-latitude vegetation in grid cells dominated by weathering-vulnerable rock types, whereas changes in temperature and river runoff had a more modest direct effect. Our results also confirmed that silicate weathering is the only mechanism that can lead to a full recovery of the carbon cycle to pre-industrial levels on multimillennial timescales.

Evaluation of the transport matrix method for simulation of ocean biogeochemical tracers

Abstract. Conventional integration of Earth system and ocean models can accrue considerable computational expenses, particularly for marine biogeochemical applications. Offline numerical schemes in which only the biogeochemical tracers are time stepped and transported using a pre-computed circulation field can substantially reduce the burden and are thus an attractive alternative. One such scheme is the transport matrix method (TMM), which represents tracer transport as a sequence of sparse matrix–vector products that can be performed efficiently on distributed-memory computers. While the TMM has been used for a variety of geochemical and biogeochemical studies, to date the resulting solutions have not been comprehensively assessed against their online counterparts. Here, we present a detailed comparison of the two. It is based on simulations of the state-of-the-art biogeochemical sub-model embedded within the widely used coarse-resolution University of Victoria Earth System Climate Model (UVic ESCM). The default, non-linear advection scheme was first replaced with a linear, third-order upwind-biased advection scheme to satisfy the linearity requirement of the TMM. Transport matrices were extracted from an equilibrium run of the physical model and subsequently used to integrate the biogeochemical model offline to equilibrium. The identical biogeochemical model was also run online. Our simulations show that offline integration introduces some bias to biogeochemical quantities through the omission of the polar filtering used in UVic ESCM and in the offline application of time-dependent forcing fields, with high latitudes showing the largest differences with respect to the online model. Differences in other regions and in the seasonality of nutrients and phytoplankton distributions are found to be relatively minor, giving confidence that the TMM is a reliable tool for offline integration of complex biogeochemical models. Moreover, while UVic ESCM is a serial code, the TMM can be run on a parallel machine with no change to the underlying biogeochemical code, thus providing orders of magnitude speed-up over the online model.

A model study of warming-induced phosphorus–oxygen feedbacks in open-ocean oxygen minimum zones on millennial timescales

Abstract. Observations indicate an expansion of oxygen minimum zones (OMZs) over the past 50 years, likely related to ongoing deoxygenation caused by reduced oxygen solubility, changes in stratification and circulation, and a potential acceleration of organic matter turnover in a warming climate. The overall area of ocean sediments that are in direct contact with low-oxygen bottom waters also increases with expanding OMZs. This leads to a release of phosphorus from ocean sediments. If anthropogenic carbon dioxide emissions continue unabated, higher temperatures will cause enhanced weathering on land, which, in turn, will increase the phosphorus and alkalinity fluxes into the ocean and therefore raise the ocean's phosphorus inventory even further. A higher availability of phosphorus enhances biological production, remineralisation and oxygen consumption, and might therefore lead to further expansions of OMZs, representing a positive feedback. A negative feedback arises from the enhanced productivity-induced drawdown of carbon and also increased uptake of CO2 due to weathering-induced alkalinity input. This feedback leads to a decrease in atmospheric CO2 and weathering rates. Here, we quantify these two competing feedbacks on millennial timescales for a high CO2 emission scenario. Using the University of Victoria (UVic) Earth System Climate Model of intermediate complexity, our model results suggest that the positive benthic phosphorus release feedback has only a minor impact on the size of OMZs in the next 1000 years. The increase in the marine phosphorus inventory under assumed business-as-usual global warming conditions originates, on millennial timescales, almost exclusively (> 80 %) from the input via terrestrial weathering and causes a 4- to 5-fold expansion of the suboxic water volume in the model.

Carbon cycle implications of terrestrial weathering changes since the last glacial maximum

We examine the importance of the rock weathering feedback mechanism during the last deglacial period (∼16 000–4000 BCE) using an Earth system model of intermediate complexity (the University of Victoria Earth System Climate Model (UVic ESCM)) with four box-model parameterizations of terrestrial weathering. The deglacial climate change is driven by changes in orbital parameters, ice core reconstructions of atmospheric CO2 variability, and prescribed removal of continental ice sheets. Over the course of the 12 000 year simulation period, increases in weathering provide a mechanism that slowly removes CO2 from the atmosphere, in opposition to the observed atmospheric CO2 increase during this period. These processes transfer both carbon and alkalinity to the ocean, the combination of which results in as much as a 1000 Pg C increase in total ocean carbon, relative to a control simulation with constant weathering. However, the rapid expansion of northern hemisphere vegetation introduces a significant uncertainty among the weathering parameterizations. Further experiments to test the individual impacts of weathering dissolved inorganic carbon and alkalinity fluxes on ocean biogeochemistry suggest that the worldwide distribution of rock types and the ratio of carbonate to silicate weathering may be crucially important in obtaining an accurate estimate of changes in global weathering rates.

The Sensitivity of the Proportionality between Temperature Change and Cumulative CO2 Emissions to Ocean Mixing

The ratio of global mean surface air temperature change to cumulative CO2 emissions, referred to as transient climate response to cumulative CO2 emissions (TCRE), has been shown to be approximately constant on centennial time scales. The mechanisms behind this constancy are not well understood, but previous studies suggest that compensating effects of ocean heat and carbon fluxes, which are governed by the same ocean mixing processes, could be one cause for this approximate constancy. This hypothesis is investigated by forcing different versions of the University of Victoria Earth System Climate Model, which differ in the ocean mixing parameterization, with an idealized scenario of 1% annually increasing atmospheric CO2 until quadrupling of the preindustrial CO2 concentration and constant concentration thereafter. The relationship between surface air warming and cumulative emissions remains close to linear, but the TCRE varies between model versions, spanning the range of 1.2°–2.1°C EgC−1 at the time of CO2 doubling. For all model versions, the TCRE is not constant over time while atmospheric CO2 concentrations increase. It is constant after atmospheric CO2 stabilizes at 1120 ppm, because of compensating changes in temperature sensitivity (temperature change per unit radiative forcing) and cumulative airborne fraction. The TCRE remains approximately constant over time even if temperature sensitivity, determined by ocean heat flux, and cumulative airborne fraction, determined by ocean carbon flux, are taken from different model versions with different ocean mixing settings. This can partially be explained with temperature sensitivity and cumulative airborne fraction following similar trajectories, which suggests ocean heat and carbon fluxes scale approximately linearly with changes in vertical mixing.

A cloud feedback emulator (CFE, version 1.0) for an intermediate complexity model

Abstract. The dominant source of inter-model differences in comprehensive global climate models (GCMs) are cloud radiative effects on Earth's energy budget. Intermediate complexity models, while able to run more efficiently, often lack cloud feedbacks. Here, we describe and evaluate a method for applying GCM-derived shortwave and longwave cloud feedbacks from 4 × CO2 and Last Glacial Maximum experiments to the University of Victoria Earth System Climate Model. The method generally captures the spread in top-of-the-atmosphere radiative feedbacks between the original GCMs, which impacts the magnitude and spatial distribution of surface temperature changes and climate sensitivity. These results suggest that the method is suitable to incorporate multi-model cloud feedback uncertainties in ensemble simulations with a single intermediate complexity model.

Quantifying the Limits of a Linear Temperature Response to Cumulative CO2 Emissions

Abstract Recent studies have shown that the transient climate response to cumulative carbon emissions (TCRE) of the global temperature can be well approximated by a constant value for cumulative emissions up to about 2 TtC. However, there has been little attention given in the literature to how the TCRE varies across the range of emissions rates represented by the current RCP emissions scenarios. The authors use an ensemble of simulations generated using the University of Victoria Earth System Climate Model to quantify how the temperature response to cumulative emissions varies as a function of both the total magnitude and the rate of CO2 emissions. This study shows that the 500-yr response to a pulse CO2 emission (1.81°C TtC−1) does not depend on the magnitude of cumulative emissions up to 3 TtC. The TCRE (1.66°C TtC−1), which relates to the short-term response, is relatively insensitive to constant-rate emissions up to 30 GtC yr−1. This experiment shows that the formal way of estimating the TCRE—that is, at the point of CO2 doubling in an idealized scenario with a 1% yr−1 increase of the atmospheric concentration—is a highly robust measure. The authors conclude that the TCRE provides a good estimate of the temperature response to CO2 emissions in RCP scenarios 2.6, 4.5, and 6, whereas a constant TCRE value significantly overestimates the temperature response to CO2 emissions in RCP8.5.

Competition between ocean carbon pumps in simulations with varying Southern Hemisphere westerly wind forcing

We analyse the impact of migration and strength of Southern Hemisphere westerly winds on the ocean carbon cycle in a systematic sensitivity study with the University of Victoria Earth System Climate Model. We find that changes in the biological pump are mainly driven by changes in ocean residence times while changes in export production are negligible. Changes in the biological and physical pumps are always of opposite sign; with the physical pump being dominant for southward shifts and the biological pump being dominant for northward shifts. Furthermore, changes in the Pacific Ocean carbon budget dictate the overall changes in global marine and atmospheric carbon. Overall, atmospheric \(\hbox {CO}_2\) increases (and \(\Delta ^{14}\hbox {C}\) decreases) for northward shifts or a strengthening in wind forcing. The opposite is true for a southward shift or a weakening in wind forcing. Combining forcings (shift and intensity change) results in a combination of their impacts with the direction of the shift being the first order forcing. The terrestrial carbon reservoir absorbs (releases) 50–70 % of the net oceanic carbon loss (increase), counterbalancing the effect on atmospheric \(\hbox {CO}_2\).

Deglacial climate, carbon cycle and ocean chemistry changes in response to a terrestrial carbon release

Researchers have proposed that a significant portion of the post-glacial rise in atmospheric CO2 could be due to the respiration of permafrost carbon stocks that formed over the course of glaciation. In this paper, we used the University of Victoria Earth System Climate Model v. 2.9 to simulate the deglacial and interglacial carbon cycle from the last glacial maximum to the present. The model’s sensitivity to mid and high latitude terrestrial carbon storage is evaluated by including a 600 Pg C carbon pool parameterized to respire in concert with decreases in ice sheet surface area. The respiration of this stored carbon during the early stages of deglaciation had a large effect on the carbon cycle in these simulations, allowing atmospheric CO2 to increase by 40 ppmv in the model, with an additional 20 ppmv increase occurring in the case of a more realistic, prescribed CO2 radiative warming. These increases occurred prior to large-scale carbon uptake due to the reestablishment of boreal forests and peatlands in the proxy record (beginning in the early Holocene). Surprisingly, the large external carbon input to the atmosphere and oceans did not increase sediment dissolution and mean ocean alkalinity relative to a control simulation without the high latitude carbon reservoir. In addition, our simulations suggest that an early deglacial terrestrial carbon release may come closer to explaining some observed deglacial changes in deep-ocean carbonate concentrations than simulations without such a release. We conclude that the respiration of glacial soil carbon stores may have been an important contributor to the deglacial CO2 rise, particularly in the early stages of deglaciation.

Explicit Planktic Calcifiers in the University of Victoria Earth System Climate Model, Version 2.9

ABSTRACTMarine calcifiers as a plankton functional type (PFT) are a crucial part of the global carbon cycle, being responsible for much of the carbon export to the deep ocean entering via biological pathways. Deep ocean carbon export through calcifiers is controlled by physiological, ecological, and biogeochemical factors. This paper describes the implementation of a calcifying phytoplankton PFT in the University of Victoria Earth System Climate Model, version 2.9 (UVic ESCM), and mechanistic improvements to the representation of model carbon export (a full calcite tracer, carbonate chemistry dependent calcite dissolution rates, and a ballasting scheme). An iterative method for stabilizing and tuning the biogeochemistry is furthermore described. The UVic ESCM now fills a niche in Earth system modelling that was previously unoccupied in that it is relatively inexpensive to run, yet resolves the complete Earth system carbon cycle including prognostic calcium carbonate and a separate phytoplankton calcifier PFT. The model is now well suited to testing feedbacks between the carbonate and carbon cycles and the climate system as transient simulations. The modifications described improve the UVic ESCM's mechanistic realism without compromising performance with respect to observed carbon and nutrient fluxes. Primary production, export production, particulate organic carbon, and calcite fluxes all fall within independently observed estimates.

A dynamic marine iron cycle module coupled to the University of Victoria Earth System Model: the Kiel Marine Biogeochemical Model 2 for UVic 2.9

Abstract. Marine biological production as well as the associated biotic uptake of carbon in many ocean regions depends on the availability of nutrients in the euphotic zone. While large areas are limited by nitrogen and/or phosphorus, the micronutrient iron is considered the main limiting nutrient in the North Pacific, equatorial Pacific and Southern Ocean. Changes in iron availability via changes in atmospheric dust input are discussed to play an important role in glacial–interglacial cycles via climate feedbacks caused by changes in biological ocean carbon sequestration. Although many aspects of the iron cycle remain unknown, its incorporation into marine biogeochemical models is needed to test our current understanding and better constrain its role in the Earth system. In the University of Victoria Earth System Climate Model (UVic) iron limitation in the ocean was, until now, simulated pragmatically with an iron concentration masking scheme that did not allow a consistent interactive response to perturbations of ocean biogeochemistry or iron cycling sensitivity studies. Here, we replace the iron masking scheme with a dynamic iron cycle and compare the results to available observations and the previous marine biogeochemical model. Sensitivity studies are also conducted with the new model to test the sensitivity of the model to parameterized iron ligand concentrations, the importance of considering the variable solubility of iron in dust deposition, the importance of considering high-resolution bathymetry for the sediment release of iron, the effect of scaling the sedimentary iron release with temperature and the sensitivity of the iron cycle to a climate change scenario.

A revised global estimate of dissolved iron fluxes from marine sediments

AbstractLiterature data on benthic dissolved iron (DFe) fluxes (µmol m−2 d−1), bottom water oxygen concentrations (O2BW, μM), and sedimentary carbon oxidation rates (COX, mmol m−2 d−1) from water depths ranging from 80 to 3700 m were assembled. The data were analyzed with a diagenetic iron model to derive an empirical function for predicting benthic DFe fluxes: where γ (= 170 µmol m−2 d−1) is the maximum flux for sediments at steady state located away from river mouths. This simple function unifies previous observations that COX and O2BW are important controls on DFe fluxes. Upscaling predicts a global DFe flux from continental margin sediments of 109 ± 55 Gmol yr−1, of which 72 Gmol yr−1 is contributed by the shelf (<200 m) and 37 Gmol yr−1 by slope sediments (200–2000 m). The predicted deep‐sea flux (>2000 m) of 41 ± 21 Gmol yr−1 is unsupported by empirical data. Previous estimates of benthic DFe fluxes derived using global iron models are far lower (approximately 10–30 Gmol yr−1). This can be attributed to (i) inadequate treatment of the role of oxygen on benthic DFe fluxes and (ii) improper consideration of continental shelf processes due to coarse spatial resolution. Globally averaged DFe concentrations in surface waters simulated with the intermediate‐complexity University of Victoria Earth System Climate Model were a factor of 2 higher with the new function. We conclude that (i) the DFe flux from marginal sediments has been underestimated in the marine iron cycle and (ii) iron scavenging in the water column is more intense than currently presumed.

If Anthropogenic CO2 Emissions Cease, Will Atmospheric CO2 Concentration Continue to Increase?

If anthropogenic CO2 emissions were to suddenly cease, the evolution of the atmospheric CO2 concentration would depend on the magnitude and sign of natural carbon sources and sinks. Experiments using Earth system models indicate that the overall carbon sinks dominate, such that upon the cessation of anthropogenic emissions, atmospheric CO2 levels decrease over time. However, these models have typically neglected the permafrost carbon pool, which has the potential to introduce an additional terrestrial source of carbon to the atmosphere. Here, the authors use the University of Victoria Earth System Climate Model (UVic ESCM), which has recently been expanded to include permafrost carbon stocks and exchanges with the atmosphere. In a scenario of zeroed CO2 and sulfate aerosol emissions, whether the warming induced by specified constant concentrations of non-CO2 greenhouse gases could slow the CO2 decline following zero emissions or even reverse this trend and cause CO2 to increase over time is assessed. It is found that a radiative forcing from non-CO2 gases of approximately 0.6 W m−2 results in a near balance of CO2 emissions from the terrestrial biosphere and uptake of CO2 by the oceans, resulting in near-constant atmospheric CO2 concentrations for at least a century after emissions are eliminated. At higher values of non-CO2 radiative forcing, CO2 concentrations increase over time, regardless of when emissions cease during the twenty-first century. Given that the present-day radiative forcing from non-CO2 greenhouse gases is about 0.95 W m−2, the results suggest that if all CO2 and aerosols emissions were eliminated without also decreasing non-CO2 greenhouse gas emissions CO2 levels would increase over time, resulting in a small increase in climate warming associated with this positive permafrost–carbon feedback.

Assessing the impact of late Pleistocene megafaunal extinctions on global vegetation and climate

Abstract. The end of the Pleistocene was a turning point for the Earth system as climate gradually emerged from millennia of severe glaciation in the Northern Hemisphere. The deglacial climate change coincided with an unprecedented decline in many species of Pleistocene megafauna, including the near-total eradication of the woolly mammoth. Due to an herbivorous diet that presumably involved large-scale tree grazing, the mammoth extinction has been associated with the rapid expansion of dwarf deciduous trees in Siberia and Beringia, thus potentially contributing to the changing climate of the period. In this study, we use the University of Victoria Earth System Climate Model (UVic ESCM) to simulate the possible effects of these extinctions on climate during the latest deglacial period. We have explored various hypothetical scenarios of forest expansion in the northern high latitudes, quantifying the biogeophysical effects in terms of changes in surface albedo and air temperature. These scenarios include a Maximum Impact Scenario (MIS) which simulates the greatest possible post-extinction reforestation in the model, and sensitivity tests which investigate the timing of extinction, the fraction of trees grazed by mammoths, and the southern extent of mammoth habitats. We also show the results of a simulation with free atmospheric CO2-carbon cycle interactions. For the MIS, we obtained a surface albedo increase and global warming of 0.006 and 0.175 °C, respectively. Less extreme scenarios produced smaller global mean temperature changes, though local warming in some locations exceeded 0.3 °C even in the more realistic extinction scenarios. In the free CO2 simulation, the biogeophysical-induced warming was amplified by a biogeochemical effect, whereby the replacement of high-latitude tundra with shrub forest led to a release of soil carbon to the atmosphere and a small atmospheric CO2 increase. Overall, our results suggest the potential for a small, though non-trivial, effect of megafaunal extinctions on Pleistocene climate.