Declining Temporal Effectiveness of Carbon Sequestration: Implications for Compliance with the United National Framework Convention on Climate Change

Carbon sequestration is increasingly being promoted as a potential response to the risks of unrestrained emissions of CO₂, either in place of or as a complement to reductions in the use of fossil fuels. However, the potential role of carbon sequestration as an (at-least partial) substitute for reductions in fossil fuel use can be properly evaluated only in the context of a long-term acceptable limit (or range of limits) to the increase in atmospheric CO₂ concentration, taking into account the response of the entire carbon cycle to artificial sequestration. Under highly stringent emission-reduction scenarios for non-CO₂ greenhouse gases, 450 ppmv CO₂ is the equivalent, in terms of radiative forcing of climate,to a doubling of the pre-industrial concentration of CO₂. It is argued in this paper that compliance with the United Nations Framework Convention on Climate Change (henceforth, the UNFCCC) implies that atmospheric CO₂ concentration should be limited, or quickly returned to, a concentration somewhere below 450 ppmv. A quasi-one-dimensional coupled climate-carbon cycle model is used to assess the response of the carbon cycle to idealized carbon sequestration scenarios. The impact on atmospheric CO₂ concentration of sequestering a given amount of CO₂ that would otherwise be emitted to the atmosphere, either in deep geological formations or in the deep ocean, rapidly decreases over time. This occurs as a result of a reduction in the rate of absorption of atmospheric CO₂ by the natural carbon sinks (the terrestrial biosphere and oceans) in response to the slower buildup of atmospheric CO₂ resulting from carbon sequestration. For 100 years of continuous carbon sequestration, the sequestration fraction (defined as the reduction in atmospheric CO₂ divided by the cumulative sequestration) decreases to 14% 1000 years after the beginning of sequestration in geological formations with no leakage, and to 6% 1000 years after the beginning of sequestration in the deep oceans. The difference (8% of cumulative sequestration) is due to an eflux from the ocean to the atmosphere of some of the carbon injected into the deep ocean.The coupled climate-carbon cycle model is also used to assess the amount of sequestration needed to limit or return the atmospheric CO₂ concentration to 350–400 ppmv after phasing out all use of fossil fuels by no later than 2100. Under such circumstances, sequestration of 1–2 Gt C/yr by the latter part of this century could limit the peak CO₂ concentration to 420–460 ppmv, depending on how rapidly use of fossilfuels is terminated and the strength of positive climate-carbon cycle feedbacks. To draw down the atmospheric CO₂ concentration requires creating negative emissions through sequestration of CO₂ released as a byproduct of the production of gaseous fuels from biomass primary energy. Even if fossil fuel emissions fall to zero by 2100, it will be difficult to create a large enough negative emission using biomass energy to return atmospheric CO₂ to 350 ppmv within 100 years of its peak. However, building up soil carbon could help in returning CO₂ to 350 ppmv within 100 years of its peak. In any case, a 100-year period of climate corresponding to the equivalent of a doubled-CO₂ concentration would occur before temperatures decreased. Nevertheless, returning the atmospheric CO₂concentration to 350 ppmv would reduce longterm sea level rise due to thermal expansion and might be sufficient to prevent the irreversible total melting of the Greenland ice sheet, collapse of the West Antarctic ice sheet, and abrupt changes in ocean circulation that might otherwise occur given a prolonged doubled-CO₂ climate. Recovery of coral reef ecosystems, if not already driven to extinction, could begin.     

A simulation study for the global carbon cycle,including man's impact on the biosphere

The simulation model accounts for four major compartments in the global carbon cycle: atmosphere, ocean, terrestrial biosphere and fossil carbon reservoir. The ocean is further compartmentalized into a high and a low latitude surface layer, and into 10 deep sea strata. The oceanic carbon fluxes are caused by massflow of descending and upwelling water, by precipitation of organic material and by diffusion exchange.The biosphere is horizontally subdivided into six ecosystems and vertically into leaves, branches, stemwood, roots, litter, young humus and stable soil carbon. Deforestation, slash and burn agriculture, rangeland burning and shifts in land use have been included. The atmosphere is treated as one well mixed reservoir. Fossil fuel consumption is simulated with historic data, and with IIASA scenario's for the future. Using the low IIASA scenario an atmospheric CO₂ concentration of 431 ppmv is simulated for 2030 AD. A sensitivity analysis shows the importance of different parameters and of human behaviour. Notwithstanding the large size of the biosphere fluxes, the atmospheric CO₂ concentration in the next century will be predominantly determined by the growth rate of fossil fuel consumption.     

Millennial timescale carbon cycle and climate change in an efficient Earth system model

A new Earth system model, GENIE-1, is presented which comprises a 3-D frictional geostrophic ocean, phosphate-restoring marine biogeochemistry, dynamic and thermodynamic sea-ice, land surface physics and carbon cycling, and a seasonal 2-D energy-moisture balance atmosphere. Three sets of model climate parameters are used to explore the robustness of the results and for traceability to earlier work. The model versions have climate sensitivity of 2.8–3.3°C and predict atmospheric CO₂ close to present observations. Six idealized total fossil fuel CO₂ emissions scenarios are used to explore a range of 1,100–15,000 GtC total emissions and the effect of rate of emissions. Atmospheric CO₂ approaches equilibrium in year 3000 at 420–5,660 ppmv, giving 1.5–12.5°C global warming. The ocean is a robust carbon sink of up to 6.5 GtC year⁻¹. Under ‘business as usual’, the land becomes a carbon source around year 2100 which peaks at up to 2.5 GtC year⁻¹. Soil carbon is lost globally, boreal vegetation generally increases, whilst under extreme forcing, dieback of some tropical and sub-tropical vegetation occurs. Average ocean surface pH drops by up to 1.15 units. A Greenland ice sheet melt threshold of 2.6°C local warming is only briefly exceeded if total emissions are limited to 1,100 GtC, whilst 15,000 GtC emissions cause complete Greenland melt by year 3000, contributing 7 m to sea level rise. Total sea-level rise, including thermal expansion, is 0.4–10 m in year 3000 and ongoing. The Atlantic meridional overturning circulation shuts down in two out of three model versions, but only under extreme emissions including exotic fossil fuel resources.     

Managing atmospheric CO2

A coupled carbon cycle-climate model is used to compute global atmospheric CO₂ and temperature variation that would result from several future CO₂ emission scenarios. The model includes temperature and CO₂ feedbacks on the terrestrial biosphere, and temperature feedback on the oceanic uptake of CO₂. The scenarios used include cases in which fossil fuel CO₂ emissions are held constant at the 1986 value or increase by 1% yr^–1 until either 2000 or 2020, followed by a gradual transition to a rate of decrease of 1 or 2% yr^–1. The climatic effect of increases in non-CO₂ trace gases is included, and scenarios are considered in which these gases increase until 2075 or are stabilized once CO₂ emission reductions begin. Low and high deforestation scenarios are also considered. In all cases, results are computed for equilibrium climatic sensitivities to CO₂ doubling of 2.0 and 4.0 °C.Peak atmospheric CO₂ concentrations of 400–500 ppmv and global mean warming after 1980 of 0.6–3.2 °C occur, with maximum rates of global mean warming of 0.2–0.3 °C decade^–1. The peak CO₂ concentrations in these scenarios are significantly below that commonly regarded as unavoidable; further sensitivity analyses suggest that limiting atmospheric CO₂ to as little as 400 ppmv is a credible option.Two factors in the model are important in limiting atmospheric CO₂: (1) the airborne fraction falls rapidly once emissions begin to decrease, so that total emissions (fossil fuel + land use-induced) need initially fall to only about half their present value in order to stabilize atmospheric CO₂, and (2) changes in rates of deforestation have an immediate and proportional effect on gross emissions from the biosphere, whereas the CO₂ sink due to regrowth of forests responds more slowly, so that decreases in the rate of deforestation have a disproportionately large effect on net emission.If fossil fuel emissions were to decrease at 1–2% yr^–1 beginning early in the next century, emissions could decrease to the rate of CO₂ uptake by the predominantly oceanic sink within 50–100 yrs. Simulation results suggest that if subsequent emission reductions were tied to the rate of CO₂ uptake by natural CO₂ sinks, these reductions could proceed more slowly than initially while preventing further CO₂ increases, since the natural CO₂ sink strength decreases on time scales of one to several centuries. The model used here does not account for the possible effect on atmospheric CO₂ concentration of possible changes in oceanic circulation. Based on past rates of atmospheric CO₂ variation determined from polar ice cores, it appears that the largest plausible perturbation in ocean-air CO₂ flux due to changes of oceanic circulation is substantially smaller than the permitted fossil fuel CO₂ emissions under the above strategy, so tieing fossil fuel emissions to the total sink strength could provide adequate flexibility for responding to unexpected changes in oceanic CO₂ uptake caused by climatic warming-induced changes of oceanic circulation.     

Effects of physical changes in the ocean on the atmospheric pCO2: glacial-interglacial cycles

Based on LGM experiments with an atmosphere–ocean general circulation model, we systematically investigated the effects of physical changes in the ocean and induced biological effects as well on the low atmospheric CO₂ concentration (pCO₂) at the last glacial maximum (LGM). Numerical experiments with an oceanic carbon-cycle model showed that pCO₂ was lowered by ~30 ppm in the LGM ocean. Most of the pCO₂ reduction was explained by the change in CO₂ solubility in the ocean due to lower sea surface temperature (SST) during the LGM. Moreover, we found that SST changes in the high-latitude Northern Atlantic could explain more than one-third of the overall change in pCO₂ induced by global SST change, suggesting an important feedback between the Laurentide ice sheet and pCO₂.     

Ocean dynamics determine the response of oceanic CO2 uptake to climate change

The increase of atmospheric CO₂ concentrations due to anthropogenic activities is substantially damped by the ocean, whose CO₂ uptake is determined by the state of the ocean, which in turn is influenced by climate change. We investigate the mechanisms of the ocean’s carbon uptake within the feedback loop of atmospheric CO₂ concentration, climate change and atmosphere/ocean CO₂ flux. We evaluate two transient simulations from 1860 until 2100, performed with a version of the Max Planck Institute Earth System Model (MPI-ESM) with the carbon cycle included. In both experiments observed anthropogenic CO₂ emissions were prescribed until 2000, followed by the emissions according to the IPCC Scenario A2. In one simulation the radiative forcing of changing atmospheric CO₂ is taken into account (coupled), in the other it is suppressed (uncoupled). In both simulations, the oceanic carbon uptake increases from 1 GT C/year in 1960 to 4.5 GT C/year in 2070. Afterwards, this trend weakens in the coupled simulation, leading to a reduced uptake rate of 10% in 2100 compared to the uncoupled simulation. This includes a partial offset due to higher atmospheric CO₂ concentrations in the coupled simulation owing to reduced carbon uptake by the terrestrial biosphere. The difference of the oceanic carbon uptake between both simulations is primarily due to partial pressure difference and secondary to solubility changes. These contributions are widely offset by changes of gas transfer velocity due to sea ice melting and wind changes. The major differences appear in the Southern Ocean (?45%) and in the North Atlantic (?30%), related to reduced vertical mixing and North Atlantic meridional overturning circulation, respectively. In the polar areas, sea ice melting induces additional CO₂ uptake (+20%).     

Adjustment of the natural ocean carbon cycle to negative emission rates

Carbon dioxide removal (CDR) is the only geoengineering technique that allows negative emissions and the reduction of anthropogenic carbon in the atmosphere. Since the time scales of the global carbon cycle are largely driven by the exchanges with the natural oceanic stocks, the implementation of CDR actions is anticipated to create outgassing from the ocean that may reduce their efficiency. The adjustment of the natural carbon cycle to CDR was studied with a numerical Earth System Model, focusing on the oceanic component and considering two idealized families of CDR policies, one based on a target atmospheric concentration and one based on planned negative emissions. Results show that both actions are anticipated to release the anthropogenic carbon stored in the surface ocean, effectively increasing the required removal effort. The additional negative emissions are expected to be lower when the CDR policy is driven by planned removal rates without prescribing a target atmospheric CO₂ concentration.     

Wind-driven changes in Southern Ocean residual circulation, ocean carbon reservoirs and atmospheric CO2

The effect of idealized wind-driven circulation changes in the Southern Ocean on atmospheric CO₂ and the ocean carbon inventory is investigated using a suite of coarse-resolution, global coupled ocean circulation and biogeochemistry experiments with parameterized eddy activity and only modest changes in surface buoyancy forcing, each experiment integrated for 5,000 years. A positive correlation is obtained between the meridional overturning or residual circulation in the Southern Ocean and atmospheric CO₂: stronger or northward-shifted westerly winds in the Southern Hemisphere result in increased residual circulation, greater upwelling of carbon-rich deep waters and oceanic outgassing, which increases atmospheric pCO₂ by ～20 μatm; weaker or southward-shifted winds lead to the opposing result. The ocean carbon inventory in our model varies through contrasting changes in the saturated, disequilibrium and biogenic (soft-tissue and carbonate) reservoirs, each varying by O(10–100) PgC, all of which contribute to the net anomaly in atmospheric CO₂. Increased residual overturning deepens the global pycnocline, warming the upper ocean and decreasing the saturated carbon reservoir. Increased upwelling of carbon- and nutrient-rich deep waters and inefficient biological activity results in subduction of unutilized nutrients into the ocean interior, decreasing the biogenic carbon reservoir of intermediate and mode waters ventilating the Northern Hemisphere, and making the disequilibrium carbon reservoir more positive in the mode waters due to the reduced residence time at the surface. Wind-induced changes in the model carbon inventory are dominated by the response of the global pycnocline, although there is an additional abyssal response when the peak westerly winds change their latitude, altering their proximity to Drake Passage and changing the depth extent of the southward return flow of the overturning: a northward shift of the westerly winds isolates dense isopycnals, allowing biogenic carbon to accumulate in the deep ocean of the Southern Hemisphere, while a southward shift shoals dense isopycnals that outcrop in the Southern Ocean and reduces the biogenic carbon store in the deep ocean.     

Global transport and inter-reservoir exchange of carbon dioxide with particular reference to stable isotopic distributions

A two-dimensional model of global atmospheric transport is used to relate estimated air-to-surface exchanges of carbon dioxide (CO₂) to spatial and temporal variations of atmospheric CO₂ concentrations and isotopic composition. The atmospheric model coupled with models of the biosphere and mixed layer of the ocean describes the gross features of the global carbon cycle. In particular this paper considers the change in isotopic composition due to interreservoir exchanges and thus the potential application and measurement requirements of new isotopic observational programs.A comparison is made between the model-generated CO₂ concentration variation and those observed on secular, interannual and seasonal time scales and spatially through the depth of the troposphere and meridionally from pole-to-pole.The relationship between isotopic and concentration variation on a seasonal time-scale is discussed and it is shown how this can be used to quantitatively estimate relative contributions of biospheric and oceanic CO₂ exchange. Further, it is shown that the interhemispheric gradient of concentration and isotopic ratio results primarily from the redistribution of fossil fuel CO₂. Both isotopic and concentration data indicate that tropical deforestation contributes less than 2 Gt yr^-1 of carbon to the atmosphere.The study suggests that changes in the rate of change of the ratio of ¹³C to ¹²C in the atmosphere of less than 0.03 yr^-1 might be expected if net exchanges with the biosphere are the cause of interannual variations of CO₂ concentrations.     

Glacial-Interglacial Atmospheric CO2 Change—The Glacial Burial Hypothesis

Organic carbon buried under the great ice sheets of the Northern Hemisphere is suggested to be the missing link in the atmospheric CO2 change over the glacial-interglacial cycles. At glaciation, the advancement of continental ice sheets buries vegetation and soil carbon accumulated during warmer pe-riods. At deglaciation, this burial carbon is released back into the atmosphere. In a simulation over two glacial-interglacial cycles using a synchronously coupled atmosphere-land-ocean carbon model forced by reconstructed climate change, it is found that there is a 547-Gt terrestrial carbon release from glacial maximum to interglacial, resulting in a 60-Gt (about 30-ppmv) increase in the atmospheric CO2, with the remainder absorbed by the ocean in a scenario in which ocean acts as a passive buffer. This is in contrast to previous estimates of a land uptake at deglaciation. This carbon source originates from glacial burial,continental shelf, and other land areas in response to changes in ice cover, sea level, and climate. The input of light isotope enriched terrestrial carbon causes atmospheric δ^13C to drop by about 0.3‰ at deglaciation,followed by a rapid rise towards a high interglacial value in response to oceanic warming and regrowth on land. Together with other ocean based mechanisms such as change in ocean temperature, the glacial burial hypothesis may offer a full explanation of the observed 80-100-ppmv atmospheric CO2 change.     

Historical and future anthropogenic emission pathways derived from coupled climate�Ccarbon cycle simulations

Using a coupled climate?Ccarbon cycle model, fossil fuel carbon dioxide (CO₂) emissions are derived through a reverse approach of prescribing atmospheric CO₂ concentrations according to observations and future projections, respectively. In the second half of the twentieth century, the implied fossil fuel emissions, and also the carbon uptake by land and ocean, are within the range of observational estimates. Larger discrepancies exist in the earlier period (1860?C1960), with small fossil fuel emissions and uncertain emissions from anthropogenic land cover change. In the IPCC SRES A1B scenario, the simulated fossil fuel emissions more than double until 2050 (17 GtC/year) and then decrease to 12 GtC/year by 2100. In addition to A1B, an aggressive mitigation scenario was employed, developed within the European ENSEMBLES project, that peaks at 530 ppm CO₂(equiv) around 2050 and then decreases to approach 450 ppm during the twenty-second century. Consistent with the prescribed pathway of atmospheric CO₂ in E1, the implied fossil fuel emissions increase from currently 8 GtC/year to about 10 by 2015 and decrease thereafter. In the 2050s (2090s) the emissions decrease to 3.4 (0.5) GtC/year, respectively. As in previous studies, our model simulates a positive climate?Ccarbon cycle feedback which tends to reduce the implied emissions by roughly 1 GtC/year per degree global warming. Further, our results suggest that the 450 ppm stabilization scenario may not be sufficient to fulfill the European Union climate policy goal of limiting the global temperature increase to a maximum of 2°C compared to pre-industrial levels.     

Transport and storage of CO2 in the ocean ——an inorganic ocean-circulation carbon cycle model

Inorganic carbon in the ocean is modelled as a passive tracer advected by a three-dimensional current field computed from a dynamical global ocean circulation model. The carbon exchange between the ocean and atmosphere is determined directly from the (temperature-dependent) chemical interaction rates in the mixed layer, using a standard CO₂ flux relation at the air-sea interface. The carbon cycle is closed by coupling the ocean to a one-layer, horizontally diffusive atmosphere. Biological sources and sinks are not included. In this form the ocean carbon model contains essentially no free tuning parameters. The model may be regarded as a reference for interpreting numerical experiments with extended versions of the model including biological processes in the ocean (Bacastow R and Maier-Reimer E in prep.) and on land (Esser G et al in prep.). Qualitatively, the model reproduces the principal features of the observed CO₂ distribution bution in the surface ocean. However, the amplitudes of surface pCO₂ are underestimated in upwelling regions by a factor of the order of 1.5 due to the missing biological pump. The model without biota may, nevertheless, be applied to compute the storage capacity of the ocean to first order for anthropogenic CO₂ emissions. In the linear regime, the response of the model may be represented by an impulse response function which can be approximated by a superposition of exponentials with different amplitudes and time constants. This provides a simple reference for comparison with box models. The largest-amplitude (0.35) exponential has a time constant of 300 years. The effective storage capacity of the oceans is strongly dependent on the time history of the anthropogenic input, as found also in earlier box model studies.     

Global and regional ocean carbon uptake and climate change: sensitivity to a substantial mitigation scenario

Under future scenarios of business-as-usual emissions, the ocean storage of anthropogenic carbon is anticipated to decrease because of ocean chemistry constraints and positive feedbacks in the carbon-climate dynamics, whereas it is still unknown how the oceanic carbon cycle will respond to more substantial mitigation scenarios. To evaluate the natural system response to prescribed atmospheric ??target?? concentrations and assess the response of the ocean carbon pool to these values, 2 centennial projection simulations have been performed with an Earth System Model that includes a fully coupled carbon cycle, forced in one case with a mitigation scenario and the other with the SRES A1B scenario. End of century ocean uptake with the mitigation scenario is projected to return to the same magnitude of carbon fluxes as simulated in 1960 in the Pacific Ocean and to lower values in the Atlantic. With A1B, the major ocean basins are instead projected to decrease the capacity for carbon uptake globally as found with simpler carbon cycle models, while at the regional level the response is contrasting. The model indicates that the equatorial Pacific may increase the carbon uptake rates in both scenarios, owing to enhancement of the biological carbon pump evidenced by an increase in Net Community Production (NCP) following changes in the subsurface equatorial circulation and enhanced iron availability from extratropical regions. NCP is a proxy of the bulk organic carbon made available to the higher trophic levels and potentially exportable from the surface layers. The model results indicate that, besides the localized increase in the equatorial Pacific, the NCP of lower trophic levels in the northern Pacific and Atlantic oceans is projected to be halved with respect to the current climate under a substantial mitigation scenario at the end of the twenty-first century. It is thus suggested that changes due to cumulative carbon emissions up to present and the projected concentration pathways of aerosol in the next decades control the evolution of surface ocean biogeochemistry in the second half of this century more than the specific pathways of atmospheric CO₂ concentrations.     

Future ocean uptake of CO2: interaction between ocean circulation and biology

We discuss the potential variations of the biological pump that can be expected from a change in the oceanic circulation in the ongoing global warming. The biogeochemical model is based on the assumption of a perfect stoichiometric composition (Redfield ratios) of organic material. Upwelling nutrients are transformed into organic particles that sink to the deep ocean according to observed profiles. The physical circulation model is driven by the warming pattern as derived from scenario computations of a fully coupled ocean-atmosphere model. The amplitude of the warming is determined from the varying concentration of atmospheric CO₂. The model predicts a pronounced weakening of the thermohaline overturning. This is connected with a reduction of the transient uptake capacity of the ocean. It yields also a more effective removal of organic material from the surface which partly compensates the physical effects of solubility. Both effects are rather marginal for the evolution of atmospheric pCO₂. Running climate models and carbon cycle models separately seems to be justified. Received: 9 August 1995 / Accepted: 22 April 1996     

Role of the ocean in controlling atmospheric CO2 concentration in the course of global glaciations

Responses of ocean circulation and ocean carbon cycle in the course of a global glaciation from the present Earth conditions are investigated by using a coupled climate-biogeochemical model. We investigate steady states of the climate system under colder conditions induced by a reduction of solar constant from the present condition. A globally ice-covered solution is obtained under the solar constant of 92.2% of the present value. We found that because almost all of sea water reaches the frozen point, the ocean stratification is maintained not by temperature but by salinity just before the global glaciation (at the solar constant of 92.3%). It is demonstrated that the ocean circulation is driven not by the surface cooling but by the surface freshwater forcing associated with formation and melting of sea ice. As a result, the deep ocean is ventilated exclusively by deep water formation in southern high latitudes where sea ice production takes place much more massively than northern high latitudes. We also found that atmospheric CO₂ concentration decreases through the ocean carbon cycle. This reduction is explained primarily by an increase of solubility of CO₂ due to a decrease of sea surface temperature, whereas the export production weakens by 30% just before the global glaciation. In order to investigate the conditions for the atmospheric CO₂ reduction to cause global glaciations, we also conduct a series of simulations in which the total amount of carbon in the atmosphere?Cocean system is reduced from the present condition. Under the present solar constant, the results show that the global glaciation takes place when the total carbon decreases to be 70% of the present-day value. Just before the glaciation, weathering rate becomes very small (almost 10% of the present value) and the organic carbon burial declines due to weakened biological productivity. Therefore, outgoing carbon flux from the atmosphere?Cocean system significantly decreases. This suggests the atmosphere?Cocean system has strong negative feedback loops against decline of the total carbon content. The results obtained here imply that some processes outside the atmosphere?Cocean feedback loops may be required to cause global glaciations.     

A numerical study of the influence of an air temperature-inversion layer and a seawater density-jump layer on the structure of interacting boundary layers

The effects of an air-temperature inversion in the atmosphere and a seawater density jump in the ocean on the structure of the atmospheric and oceanic boundary layers are studied by use of a coupled model. The numerical model consists of a closed system of equations for velocities, turbulent kinetic energy, turbulent exchange coefficient, local turbulent length scale, and stratification expressions for both air and sea boundary layers. The effects of the temperature inversion and the density jump are incorporated into the equations of turbulent kinetic energy of the atmosphere and ocean by a parameterization. A series of numerical experiments was conducted to determine the effects of various strengths of the inversion layer and surface heat fluxes in the atmosphere and of the density-jump layer in the ocean on the structure of the interacting boundary layers.The numerical results show that the temperature inversion in the atmosphere and density jump in the ocean have strong influences on turbulent structure [especially on the turbulent exchange coefficient (TEC) and turbulent kinetic energy (TKE)] and on air-sea interaction characteristics. Maxima of TKE and TEC strongly decrease with increasing strength of the inversion layer, and they disappear for strong inversions in the atmosphere. Certain strengths (density differences between the upper and the lower layers) of the density-jump layer in the ocean (₂ 0.1 g/cm³) produce double maxima in TEC-profiles and TKE-profiles in the ocean. The magnitudes of air-sea interaction characteristics such as geostrophic drag coefficient, and surface drift current increase with increasing strength of the density-jump layer in the ocean. The density-jump layer plays the role of a barrier that limits vertical mixing in the ocean. The numerical results agree well with available observed data and accepted quantitive understanding of the influences of a temperature inversion layer and a density-jump layer on the interacting atmospheric and oceanic boundary layers.     

The Variability of Air-sea O2 Flux in CMIP6: Implications for Estimating Terrestrial and Oceanic Carbon Sinks

The measurement of atmospheric O₂ concentrations and related oxygen budget have been used to estimate terrestrial and oceanic carbon uptake. However, a discrepancy remains in assessments of O₂ exchange between ocean and atmosphere (i.e. air-sea O₂ flux), which is one of the major contributors to uncertainties in the O₂-based estimations of the carbon uptake. Here, we explore the variability of air-sea O₂ flux with the use of outputs from Coupled Model Intercomparison Project phase 6 (CMIP6). The simulated air-sea O₂ flux exhibits an obvious warming-induced upward trend (~1.49 Tmol yr^?2) since the mid-1980s, accompanied by a strong decadal variability dominated by oceanic climate modes. We subsequently revise the O₂-based carbon uptakes in response to this changing air-sea O₂ flux. Our results show that, for the 1990?2000 period, the averaged net ocean and land sinks are 2.10±0.43 and 1.14±0.52 GtC yr^?1 respectively, overall consistent with estimates derived by the Global Carbon Project (GCP). An enhanced carbon uptake is found in both land and ocean after year 2000, reflecting the modification of carbon cycle under human activities. Results derived from CMIP5 simulations also investigated in the study allow for comparisons from which we can see the vital importance of oxygen dataset on carbon uptake estimations.     

The millennial atmospheric lifetime of anthropogenic CO<Subscript>2</Subscript>

The notion is pervasive in the climate science community and in the public at large that the climate impacts of fossil fuel CO₂ release will only persist for a few centuries. This conclusion has no basis in theory or models of the atmosphere/ocean carbon cycle, which we review here. The largest fraction of the CO₂ recovery will take place on time scales of centuries, as CO₂ invades the ocean, but a significant fraction of the fossil fuel CO₂, ranging in published models in the literature from 20–60%, remains airborne for a thousand years or longer. Ultimate recovery takes place on time scales of hundreds of thousands of years, a geologic longevity typically associated in public perceptions with nuclear waste. The glacial/interglacial climate cycles demonstrate that ice sheets and sea level respond dramatically to millennial-timescale changes in climate forcing. There are also potential positive feedbacks in the carbon cycle, including methane hydrates in the ocean, and peat frozen in permafrost, that are most sensitive to the long tail of the fossil fuel CO₂ in the atmosphere.     

Dangerous anthropogenic interference, dangerous climatic change, and harmful climatic change: non-trivial distinctions with significant policy implications

Article 2 of the United Nations Framework Convention on Climate Change (UNFCCC) calls for stabilization of greenhouse gas (GHG) concentrations at levels that prevent dangerous anthropogenic interference (DAI) in the climate system. However, some of the recent policy literature has focused on dangerous climatic change (DCC) rather than on DAI. DAI is a set of increases in GHGs concentrations that has a non-negligible possibility of provoking changes in climate that in turn have a non-negligible possibility of causing unacceptable harm, including harm to one or more of ecosystems, food production systems, and sustainable socio-economic systems, whereas DCC is a change of climate that has actually occurred or is assumed to occur and that has a non-negligible possibility of causing unacceptable harm. If the goal of climate policy is to prevent DAI, then the determination of allowable GHG concentrations requires three inputs: the probability distribution function (pdf) for climate sensitivity, the pdf for the temperature change at which significant harm occurs, and the allowed probability (“risk”) of incurring harm previously deemed to be unacceptable. If the goal of climate policy is to prevent DCC, then one must know what the correct climate sensitivity is (along with the harm pdf and risk tolerance) in order to determine allowable GHG concentrations. DAI from elevated atmospheric CO₂ also arises through its impact on ocean chemistry as the ocean absorbs CO₂. The primary chemical impact is a reduction in the degree of supersaturation of ocean water with respect to calcium carbonate, the structural building material for coral and for calcareous phytoplankton at the base of the marine food chain. Here, the probability of significant harm (in particular, impacts violating the subsidiary conditions in Article 2 of the UNFCCC) is computed as a function of the ratio of total GHG radiative forcing to the radiative forcing for a CO₂ doubling, using two alternative pdfs for climate sensitivity and three alternative pdfs for the harm temperature threshold. The allowable radiative forcing ratio depends on the probability of significant harm that is tolerated, and can be translated into allowable CO₂ concentrations given some assumption concerning the future change in total non-CO₂ GHG radiative forcing. If future non-CO₂ GHG forcing is reduced to half of the present non-CO₂ GHG forcing, then the allowable CO₂ concentration is 290–430 ppmv for a 10% risk tolerance (depending on the chosen pdfs) and 300–500 ppmv for a 25% risk tolerance (assuming a pre-industrial CO₂ concentration of 280 ppmv). For future non-CO₂ GHG forcing frozen at the present value, and for a 10% risk threshold, the allowable CO₂ concentration is 257–384 ppmv. The implications of these results are that (1) emissions of GHGs need to be reduced as quickly as possible, not in order to comply with the UNFCCC, but in order to minimize the extent and duration of non-compliance; (2) we do not have the luxury of trading off reductions in emissions of non-CO₂ GHGs against smaller reductions in CO₂ emissions, and (3) preparations should begin soon for the creation of negative CO₂ emissions through the sequestration of biomass carbon.