Main

Delivering on the many national and corporate net-zero emission pledges will probably require the gross removal of atmospheric carbon dioxide (CO2) on top of conventional emission reductions6,7. To achieve the Paris Agreement, global gross CO2 removals will need to exceed gross residual emissions4,8 after the middle of the century1,9. The resultant net-negative emissions compensate for the carbon debt3 accrued by CO2 emissions that overshoot the remaining carbon budget10,11. Carbon debt is projected to amount to roughly the equivalent of 9 years of global emissions before the COVID-19 pandemic according to the 1.5 °C ‘middle of the road’ scenario P3/S2 of the Intergovernmental Panel on Climate Change (IPCC)1 (Extended Data Table 1). Such large-scale deployment of CDR is controversial mainly for the implied economic and technological risks12,13,14,15,16 and environmental effects17,18; and because reliance on CDR in mitigation scenarios often goes hand-in-hand with a substantial shift of the mitigation burden to future generations19.

Here we would like to highlight a fundamental economic problem associated with the existing assessments of climate mitigation scenarios, aiming to inform international climate negotiations. Existing economic policy instruments for emission control are inadequate to incentivize a global transformation towards a net-negative carbon economy without imposing excessive fiscal burden from 2050 onwards. Currently envisaged carbon tax schemes would turn into public subsidies under net-negative emissions with potentially prohibitive fiscal implications5. Emission trading schemes (ETS), on the other hand, are presently designed to handle only positive emission caps. Negative emissions are merely treated as offsets, suggesting that CO2 emissions from one point in time cannot be compensated by an equivalent quantity of negative emissions at another point in time, as required by most mitigation scenarios. Notably, we observe that pricing the depletion of the remaining carbon budget is fundamentally different to pricing overshot emissions after the depletion of the budget, which has profound implications for the consistent earmarking of accrued revenues from a price on CO2.

We argue that establishing the responsibility of emitters for carbon debt is a prerequisite to ensuring viable net-negative carbon futures. Carbon debt could therefore be treated similar to financial debt, including interest payments on physical liabilities (that is, as a CRO) to internalize the inherent risks. On the basis of this idealized global carbon policy proposal motivated by the IPCC’s mitigation scenarios, our numerical results address the shortcomings of the existing climate mitigation literature20. Despite the conceptual character of this study, we establish profound implications for national carbon policies, which are strongly influenced by the IPCC’s global mitigation pathways in many high-emission countries21.

Carbon pricing for net-negative emissions

Integrated assessment models (IAMs) provide global carbon price paths that serve as a proxy for a wider range of cost-effective climate policy options to achieve specified greenhouse gas mitigation goals9. Such carbon prices typically increase exponentially with the interest rate as a consequence of the Hotelling rule, which defines the intertemporally optimal extraction schedule and price of an non-renewable resource22,23, such as the carbon budget. If understood as a global common, revenues generated from pricing its depletion should consistently add to public budgets, for instance to compensate for the associated welfare effects, which may be unfairly distributed across society. However, in scenarios in which the carbon budget is overshot and subsequently replenished, the budget can no longer be regarded as a non-renewable resource. In this case, the Hotelling rule lends itself to an ‘intertemporal interpretation’ for carbon policy: revenues from carbon pricing after the depletion of the budget can be invested at the market interest rate to finance net carbon removal later in the century. Because marginal abatement costs increase at the market interest rate, this calculation is exact under perfect foresight conditions—as assumed in most IAMs—if the retained funds purchase net-negative emissions at marginal costs later on. Because emitters pay for future net CDR through the carbon price, this intertemporal interpretation is compatible with the ‘polluter pays principle’. The resultant intertemporal financial transfer thereby addresses concerns of intergenerational equity because public budgets in the near-term no longer spuriously benefit from pricing an already depleted resource, while future generations thereafter are forced to replenish the carbon budget through other sources, such as income, sales or payroll taxes. According to the ‘conventional interpretation’ of the Hotelling rule, revenues from carbon pricing are merely treated as contemporaneous additions to public budgets, with no clear earmarking of accrued funds. Notably, as both approaches are simply interpretations of the same underlying carbon price paths, emitters also pay the discounted future costs of net emission removal in case of the conventional interpretation. However, in the absence of consistent earmarking, the financial viability of net CDR in the second half of the twenty-first century is highly doubtful5, and intergenerational equity remains unaccounted for.

To operationalize a future net-negative carbon economy, carbon tax revenues could be partially retained and transferred over generations to finance net CDR in the style of a nuclear decommissioning trust fund or a sovereign wealth fund. The value of such a global net carbon removal fund is potentially enormous, yet in the range of comparable funds, peaking at roughly 100% of global gross domestic product (GDP) in the median of the Shared Socioeconomic Pathway (SSP) scenarios that are compatible with Representative Concentration Pathway 1.9 (RCP 1.9)24 (Fig. 1). For comparison, Norway’s large sovereign wealth fund has passed 250% of national GDP25. Given this order of magnitude, intermediate investment portfolios could be a game changer to lift CDR out of the pilot phase even before pay-out of the fund. However, protecting financial resources from diversion for other purposes as political environments change, or as public finances become stressed, will surely be extremely challenging. For instance, sovereign borrowing to cushion the effects of the COVID-19 pandemic meant that by the end of 2020 the debt-to-GDP ratio of governments according to the Organisation for Economic Co-operation and Development had increased by about 13.4 percentage points26. Severe crises in the future could induce considerable pressure for governments to appropriate savings originally reserved for net CDR.

Fig. 1: Idealized global tax scheme with net carbon removal fund.
figure 1

a, Bottom, public income and expenditure from a tax on net emissions expressed as a percentage of GDP. Hotelling-compatible (exponential) carbon prices from SSP–RCP 1.9 scenarios are multiplied by net emissions and divided by GDP (grey dashed lines). An idealized income/expenditure curve (black solid line) was derived from these scenarios using a strictly exponential median carbon price, median net emissions and GDP. Instead of reserving 100% of tax revenues after depletion of the carbon budget, we assert that a fraction ϕ = 0.76 of revenues is earmarked for net carbon removal, from 2020 onwards. This share of income (green area) would need to be accrued into a net carbon removal fund invested at the market rate of interest to account for later expenditure when net emissions turn negative. See Methods for a definition of ϕ. b, Top, cumulative payments into the net carbon removal fund (green) and interest (orange) in theory pay exactly for cumulative tax expenditure (blue), such that the net value of the fund (brown solid line) gets exhausted as the warming target is achieved in 2100.

The success of a net CDR fund also depends on the appropriate choice of several inherently uncertain parameters, including future abatement costs. If costs and other socioeconomic parameters are not estimated in line with the precautionary principle, or if regulators are reluctant to adequately reflect future carbon removal in near-term price instruments, insufficient financial resources would be collected as observed for nuclear decommissioning27. Because the carbon debt and associated risks would be mutualized by a net CDR fund, missing financial resources would need to be replenished by public budgets.

Dynamic emission trading

Emission trading with fully liberalized banking and borrowing of allowances can be regarded as a response to these concerns. Decentralized decision-making and price determination in a competitive market is believed to improve efficiency by leveraging the ability of carbon markets to determine cost-effective time paths of mitigation28. In an idealized global scheme, the remaining carbon budget would be distributed over time resulting in positive emission caps for consecutive auctioning periods. Emitters would decide in each period what fraction of their CO2 emissions to compensate for by allowances and how much carbon debt to generate for compensation by future allowances—or future CDR in the absence of a positive emission cap. Effectively, emitters generating carbon debt would remain liable for the timing and delivery of net-negative emissions (Fig. 2) and can therefore balance present against future abatement based on individual expectations, such as those concerning technological breakthroughs. Stranded assets can be avoided by harmonizing abatement investments with natural renewal cycles of capital; and fluctuations in the business cycle can be addressed. Fixed price schedules under a carbon tax suggest lower costs for hedging risks related to the long-run costs of negative emissions and low-carbon investments. However, increased intertemporal flexibility in emission trading stabilizes the price—which reflects discounted future marginal abatement costs—compared with currently implemented ETS with no intertemporal trade of allowances29,30. At least in principle, this ETS arrangement enables emitters to develop optimal investments over longer time horizons, increasing the dynamic efficiency of emission trading. Although emission caps can be overshot, the quantity of cumulative emissions remains exactly controlled under an ETS with intertemporal trade of carbon debt, which is, more generally, the main advantage of cap-and-trade schemes compared with carbon taxes. If caps no longer directly control emission reductions, they can be set to equitably distribute ETS revenues over time. However, as the carbon budget diminishes rapidly—the 1.5 °C compatible budget is projected to become depleted roughly within the next 10 years1—the importance of carbon debt management increasingly outweighs the requirement of an adequate temporal distribution of the remaining carbon budget.

Fig. 2: Idealized ETS with intertemporal trade of carbon debt.
figure 2

Illustrative 2 °C pathway with gross carbon emissions from FFI, LUC and non-specified sources of CDR including the schematic architecture of idealized global intertemporal emission trading. ETS emission caps bt are obtained by distributing the carbon budget in tranches over consecutive periods. The amount by which emission caps bt are exceeded by net emissions is conceptualized as ‘carbon debt’ (dt). In this idealized illustration, dt is compensated later by corresponding net-negative emissions (eNN,t) such that dt = −eNN,t. In a conventional ETS, emission caps would be set to bt + dt, and eNN,t would have to be incentivized by public subsidies. dt, bt and eNN,t (in which t indicates 1, 2 and so on) are simplified discrete analogues of the continuous variables d(t), b(t) and eNN(t), respectively, which are described in the Methods. Historical emissions are from a previously published study39.

Privately managed carbon debt within an ETS also has considerable drawbacks: the enforcement of carbon debt, assessment of creditworthiness of emitters, the potential for speculation on future softening of emission targets and subsequent deferral of mitigation (time inconsistency)—which is stronger the lower the solvency of emitters (adverse selection)—and the resultant incentive to lobby for cancellation of carbon debt (moral hazard) are crucial obstacles that explain why such intertemporal mechanisms are severely restricted in currently implemented ETS28. Moreover, intertemporal trade of carbon debt by means of forward and future markets trading negative emissions over potentially long periods at a fixed price is perceived as infeasible, given the deep uncertainty in the parameters guiding a large-scale CDR rollout31.

Carbon removal obligations

Intertemporal emission trading would necessarily come at the cost of considerable regulation to address these drawbacks. We argue, however, that practices from the financial industry and monetary policy could be leveraged to reduce risks and adaptively balance potentially competing interests of economic development and climate mitigation by treating carbon debt in a similar manner to a financial debt obligation, and thereby invoking an interest on carbon debt. Economic growth, aggregate demand for carbon debt and individual financial ratings of debtors would define a general base rate, individual mark-ups, term structures and debt maturities. To assure its physical conservation and exert control over its aggregate level, carbon debt would initially be issued at the base rate by managing authorities—for example, Central Banks—to which commercial banks would be held liable in case of insolvent debtors. Commercial banks, or their equivalents, would issue debt to emitters and, assisted by rating agencies, assess and hedge their insolvency risk by determining individual mark-ups on the base rate. Carbon debt would enter the balance sheets of firms as a physical liability in tonnes (t) CO2—a carbon removal obligation, for which interest payments would be due (Extended Data Fig. 7a). This chain of legal liabilities across layers of public and private actors reduces the moral hazard that governments would ultimately pick up the bill for net emission removal, and limit the issuance of CROs to debtors who are reluctant to fulfil their (interest) obligations. Individual interest mark-ups would also balance the push of the market for adverse selection and incentivize a debt transfer from agents losing ground under stringent climate policy to low-risk agents; or lead to more near-term abatement (see below) if risks are deemed non-insurable. The rate controls the volatility of the carbon price (Extended Data Fig. 1) and therefore directly affects the price–risk costs of scheduled abatement investments. More generally, interest and debt maturities would need to reflect the speculative nature of CDR, leading to short—but potentially renewable—repayment terms and elevated rates in the near-term. A concrete phase-in scenario of CROs in the ETS of the European Union and beyond is described in Box 1.

For intertemporal emission trading to work efficiently—for instance to reduce issues of time inconsistency and price volatility—emission caps would need to be credibly announced as early as possible. As a consequence, regulators would lose the flexibility of adapting caps as new knowledge concerning the Earth system becomes available. In an idealized global scheme, emission caps need to exactly reflect the remaining carbon budget. Budget uncertainties related to the issuing of carbon debt, similar to those of permafrost thaw after a temperature overshoot2, could be hedged by collecting risk funds through base rate payments and by incentivizing more-ambitious emission reductions to minimize the risk of climate feedback effects (see below). Such uncertainties should remain manageable by risk reserves, allowing for the budget to be replenished by drawing on risk funds rather than requiring a downwards correction of scheduled emissions caps. In the best case, uncertainties and base rates would decrease over time as updated estimates of the carbon budget converge to a value within the expected range of the previously announced budget. However, new findings might realistically also lead to exceeding of the abilities of risk management, requiring a combined effort of future generations to counter potentially abrupt climate change. Management of physical risks therefore remains limited to what is presently perceivable and realistically quantifiable.

Climate mitigation under carbon debt

In IAMs, abatement costs are discounted at the market interest rate, implying a cost advantage for abatement in the distant future compared with near-term decarbonization in terms of net present value. The interest rate is therefore a key driver of carbon debt accrual in IAMs32,33. This ‘discounting effect’ is balanced by imposing interest on carbon debt. Longer CRO maturities indicate lower net present costs for CDR. Simultaneously, carbon debt interest is paid over a longer period, compensating for these gains. When the market rate of interest and the carbon debt interest rate (rd) coincide, the gains from discounting are balanced exactly, as we analytically show in the Methods. In Fig. 3, we illustrate the sensitivity of 2 °C-compatible global mitigation pathways to interest on carbon debt, with rates constant over the 2020–2100 period ranging from rd = 0 to rd = 0.08. For each rate, 13 scenarios are computed based on different SSPs and IAMs that are used to calibrate the marginal abatement cost curves of our model (Extended Data Table 2).

Fig. 3: The 2 °C (RCP 2.6) mitigation scenarios for a range of interest rates on carbon debt.
figure 3

a, Net CO2 emissions of all scenarios with rd = 0 (turquoise) and rd = 0.08 (yellow), including geometric median paths (bold solid lines) and minimum to maximum ranges (shaded areas). b, Marginal abatement costs of scenarios with rd = 0 (turquoise) and rd = 0.08 (yellow). Bold solid lines indicate geometric medians, shaded areas indicate 25–75% interquartile ranges. c, d, Geometric median net emissions as in a, including gross emissions from FFI, BECCS and LUC. c, rd = 0. d, rd = 0.08. e, Total discounted abatement costs (net present value, including interest costs) expressed as a percentage cost increase compared with the baseline (rd = 0) are shown as function of total carbon debt D. The boxes indicate the 25–75% interquartile ranges around the median values of the costs and D. Symbols linked by grey solid lines indicate the medians grouped by SSP. The entire dataset is shown in the top right corner, in which each scenario is reflected by a symbol, grouped by SSP (symbol type) and rd (colour).

For comparison, only the two extreme cases—rd = 0 and rd = 0.08—are illustrated in Fig. 3a–d. Notably, when rd = 0.08, the cumulative emission target is achieved without the accrual of carbon debt in the median path (Fig. 3d), suggesting that emissions remain at the net-zero level once achieved. This is accomplished by the contemporaneous compensation of residual CO2 from fossil fuels and industry (FFI) with negative emissions from bioenergy with carbon capture and storage (BECCS) and land-use change (LUC). Complete decarbonization of FFI emissions is, however, not cost-effective owing to the high marginal costs of emission reductions from hard-to-abate sectors. Notably, net-negative emissions of individual scenarios in Fig. 3a turn back to zero before 2100, thereby minimizing the ‘problem of phasedown’34. With reduced reliance on net-negative emissions, marginal costs are higher in the near-term due to the more-rapid reduction in FFI emissions and increase in BECCS, but considerably lower in 2100 (Fig. 3b). Figure 3e shows a reduction in the total carbon debt D as rd is gradually increased. Carbon debt risks are therefore greatly reduced at a moderate cost increase of below 12.5% in more than 75% of scenarios in which rd > 0.02.

A similar analysis was performed for the 1.5 °C global warming target; however, direct air capture and storage (DACS) is added to the mitigation technology mix, represented by six different DACS-specific marginal abatement cost curves with low, medium and high costs as well as low- and high-capacity limits. This results in a set of 78 scenarios for each rate rd. Not surprisingly, the higher the potential for DACS to be deployed, the larger the level of D when rd = 0. By contrast, when interest is invoked, this discounting effect is reversed and scenarios with large-capacity low-cost DACS simultaneously exhibit the lowest levels of D (Extended Data Fig. 2). The pathways in Fig. 4 show baseline (Fig. 4a, b) and reduced D (Fig. 4c, d) scenarios for those scenarios that achieve a reduction in D of at least 30% compared with their associated baselines (see Extended Data Figs. 35 for reductions of 5%, 15% and 45%, respectively). For illustration, we interpret the CRO-ETS baseline scenarios, in which rd = 0, as conventional ETS scenarios because both schemes are theoretically equivalent in terms of the resultant emission profiles while they imply a qualitatively different timing of financial flows.

Fig. 4: The 1.5 °C (RCP 1.9) pathways under a conventional ETS or a CRO-ETS.
figure 4

a, b, A conventional ETS is used. c, d, A CRO-ETS is used. a, c, Geometric median net emissions (solid line) and gross emissions from FFI, BECCS, LUC and DACS. Net emissions from a are also displayed in c (dashed line) and vice versa. The total carbon debt D is shown as a box-and-whiskers plot. Boxes indicate the 25–75% interquartile range around the median values (bold line), whiskers indicate minimum to maximum ranges, points mark the outliers. b, d, Annual mitigation costs as a percentage of GDP, including the share of average abatement costs attributed to emission reductions (ABM), to the compensation of residual emissions by CDR (RES) and to net-negative emissions (NNE), as well as expenditures for allowances (ETS) and interest costs (INT). Total mitigation costs (that is, ABM + RES + NNE + ETS + INT) from d are also displayed in b (dashed line) and vice versa. Box-and-whiskers plots show the total discounted abatement costs (that is, ABM + RES + NNE) as a percentage of GDP, the number above the chart indicates out-of-range outliers. Pie charts in d summarize the properties of the underlying set of scenarios (see Methods). The distribution of rd in CRO-ETS scenarios is depicted in c.

Despite the earlier increase in DACS in Fig. 4c, causing emissions to turn net zero around 2050, an emission overshoot appears to be inevitable if warming is to be limited to 1.5 °C. Remaining net-negative emissions might cause problems of phasedown in 2100, unless CRO maturities are further extended to enable a smooth transition to net-zero emissions; or more net-negative emissions are needed to stabilize the climate in the twenty-second century35. Therefore, the median D is equal to roughly 7 years of global net emissions in 2019 in Fig. 4c. Yet, the role of CDR changes considerably: without considering risks, CDR seems to justify late-century compensation of carbon debt. In this case, the median D is equivalent to about 11 years of 2019 global net emissions, with compensation starting roughly 10 years later in Fig. 4a. However, when risks are accounted for by imposing interest, CDR supports a rapid decrease in net emissions by balancing the residual emissions. Controversially, the availability of cheap and large-scale CDR options, such as DACS, is key in 1.5 °C scenarios with reduced reliance on the accrual of carbon debt. As illustrated by the pie charts in Fig. 4d, the share of high-capacity DACS scenarios among feasible scenarios with respect to the 30% reduction requirement grows to 81% (50% in the underlying set) and the share of low-cost DACS scenarios to 54% (33% in the underlying set). Should CDR not become readily available as asserted in IAMs36,37,38, this would be reflected in an elevated carbon debt interest rate, incentivizing emission reductions provided by other sources, such as the replacement of fossil fuel with renewable energy sources in hard-to-abate sectors—even if this leads to much higher costs.

Effect on financial flows over time

Figure 4, moreover, illustrates the distribution of annual mitigation cost shares, including investments in emission reductions and negative emissions and the financial flows associated with ETS allowances and interest for CROs. The share of abatement costs for emission reductions (ABM), negative emissions compensating for residual emissions (RES) and net-negative emissions (NNE) incurred in the near- versus the long-term increases with larger levels of rd (compare Fig. 4d to Extended Data Figs. 3d–5d). Here, CROs with interest induce a more equitable temporal distribution of these cost items, in sum peaking at 2.4% in Fig. 4d compared to 4.5% of GDP in Fig. 4b. This is partly because the CRO-ETS requires carbon debtors to reserve financial resources early in the century, and such funds earn interest until they are spent for net-negative emissions. By contrast, net-negative emissions expenditures in Fig. 4b are incurred at the time of net carbon removal and would need to be funded by public sources in the absence of intertemporal financial transfers. Note that here we show average abatement costs. If marginal costs are paid by incentivizing net CDR on a market, public expenditures are much higher (for comparison, see Fig. 1). Pricing overshot emissions under a conventional ETS, moreover, implies much larger revenues (‘ETS’) than under the CRO-ETS, where emission caps reflect exactly the remaining carbon budget. Median total discounted abatement costs, excluding ETS costs and interest costs (‘INT’), increase from 1.6% to 2.0% of GDP when interest is invoked within the CRO-ETS. Median interest costs in these scenarios are substantial, peaking at above 1.3% of GDP, and 0.4% to 1.5% in Extended Data Figs. 35. These numbers are, however, highly uncertain and will need to be determined considering the viability and scalability of near-term CDR options and other emission reduction technologies.

Enlarging IAM CDR portfolios to reduce technological risks and environmental effects would probably lead to further burdening of future generations in scenarios if CDR remains primarily a motivation for reducing net present costs by accrual of carbon debt. This is especially problematic if such results trickle down through the IPCC and international climate negotiations into national target setting because no viable mechanisms for the repayment of carbon debt have entered the policy debate at the moment. Simultaneously, mitigation pathways with reduced carbon debt heavily rely on CDR, requiring that risks be appropriately managed. Similar pathways result from lowering the market interest rate in IAMs32 or from adequately setting intermediate climate targets or constraints on net emissions20. However, such measures would individually not resolve the more profound issue of finance of net-negative emissions discussed here.

Conclusion

In view of the rapid depletion of the global carbon budget, CROs seem to be indispensable for any robust climate mitigation framework. CROs imply a paradigm shift from pricing the permanent to pricing the temporary storage of CO2 in the atmosphere, with carbon debtors being responsible for delivering net CDR. The implied flexibility for emitters also bears the largest drawback of intertemporal emission trading, if public bailout of carbon debtors becomes necessary. To minimize such risks, the ‘conservation of carbon debt’ needs to take top priority by controlling the total amount of carbon debt and by establishing liability across several layers of actors. Risk management under a CRO-ETS relies on imposing interest on carbon debt. For higher and risk-adjusted carbon debt interest rates, net-negative emission investments no longer benefit from net present cost gains when mitigation is deferred to the distant future. By implication, CDR under a CRO-ETS will need to prove its viability compared with conventional options for the reduction of emissions already in the near-term. This will promote bottom-up CDR market development with the accompanying benefits of price discovery, earlier technological learning, testing of scalability and identification of socio-environmental co-benefits and hazards, and ultimately, eliminating the uncertainties surrounding CDR.

Methods

Basic analytical setup

Emission reductions induced by a CRO-ETS are quantified using a Hotelling-type optimization problem (see ref. 32 for an analytical solution of the model). A global social planner is tasked to implement emission reductions at minimum costs to meet a cumulative emission target—that is, the remaining carbon budget B—by T = 2100 (t0 = 2020). Exogenously given baseline emissions Ebase—that is, future emission paths based on ‘business-as-usual’ climate policy assumptions, are reduced by a fraction a to obtain net emissions e:

$$e(t)={E}_{{\rm{base}}}(t)(1-a(t)).$$
(1)

Total abatement costs ctot are discounted at the market interest rate r to obtain the net present value of total abatement costs, which is minimized:

$$\mathop{{\rm{\min }}}\limits_{a(t),\,{c}_{{\rm{tot}}}(t,a(t))\,}{\int }_{{t}_{0}}^{T}{c}_{{\rm{tot}}}(t,a(t))\exp (-r(t-{t}_{0})){\rm{d}}t,$$
(2)

subject to:

$${\int }_{{t}_{0}}^{T}e(t){\rm{d}}t=B.$$
(3)

Integrating over marginal abatement costs MAC(a) gives the cost per tonne CO2 for an instantaneous emission reduction of a compared to the baseline. Consequently, total abatement costs ctot are defined as:

$${c}_{{\rm{tot}}}(t,a(t))={E}_{{\rm{base}}}(t)\,{\int }_{0}^{a(t)}{\rm{MAC}}(t,\tilde{a}){\rm{d}}\tilde{a}.$$
(4)

Assume that under an idealized CRO-ETS, a constant fraction 1 − ϕ of net-positive emissions eNP is equivalent to a (continuous) emission cap b (that is, the amount of conventional emission allowances issued over time) and ϕ < 1 of eNP equals carbon debt d (that is, the quantity of CROs issued). Then ϕ is defined as the ratio of cumulative net-negative emissions eNN to cumulative net-positive emissions eNP:

$$\varphi =\frac{{\int }_{{t}_{0}}^{T}{e}_{{\rm{NN}}}(t){\rm{d}}t}{{\int }_{{t}_{0}}^{T}{e}_{{\rm{NP}}}(t){\rm{d}}t},$$
(5)

and net-negative emissions and net-positive emissions equal the negative and positive parts of net emissions (eNPeNN > 0)

$${e}_{{\rm{NN}}}(t)=\{\begin{array}{cc}-e(t) & {\rm{for}}\,e(t) < 0\\ 0 & {\rm{for}}\,e(t)\ge 0\end{array},$$
(6)
$${e}_{{\rm{NP}}}(t)=\{\begin{array}{cc}e(t) & {\rm{for}}\,e(t) > 0\\ 0 & {\rm{for}}\,e(t)\le 0\end{array}.$$
(7)

Carbon debt d and the continuous emission cap b are defined as:

$$d(t)=\varphi {e}_{{\rm{NP}}}(t),\,b(t)=(1-\varphi ){e}_{{\rm{NP}}}(t),$$
(8)

and total carbon debt D is obtained by integration over the planning horizon T (combining equations (5) and (8)):

$$D={\int }_{{t}_{0}}^{T}d(t){\rm{d}}t={\int }_{{t}_{0}}^{T}{e}_{{\rm{NN}}}(t){\rm{d}}t.$$
(9)

Consequently, we can write ϕ as:

$$\varphi =\frac{D}{B+D}.$$
(10)

By implication, a fraction ϕ of cumulative net-positive emissions overshoots B and thereby generates D, and a fraction 1 − ϕ depletes the budget B.

Instead of exogenously imposing ETS emission caps, ϕ allows us to endogenously compute caps b and carbon debt d to conceptualize the intertemporal allocation of carbon debt such that debt is solely compensated by net-negative emissions eNN. On this basis, we can compute a ‘physical repayment term’ TR linking the timing of net-positive to net-negative emissions (equation (12)). CROs in this idealized ETS therefore represent a long-term intertemporal net transaction for financing net-negative emissions. This aggregate can be regarded as a proxy for a multitude of smaller carbon debt transfers over shorter timeframes that are possible in real ETS implementations in which CROs can be compensated by (gross) carbon removal, issuance of new CROs or allowances at a later point in time.

Average abatement costs are obtained from total abatement costs by dividing by the abated quantity of CO2:

$${c}_{{\rm{avg}}}(t,a(t))=\frac{{c}_{{\rm{tot}}}(t,a(t))}{a(t){E}_{{\rm{base}}}(t)}.\,$$
(11)

Next, we introduce interest payments that are due for carbon debt d over the repayment term t → t + TR(t), that is, from issuance of the CRO until its retirement (Extended Data Fig. 7b–e). TR is implicitly defined as:

$${\int }_{{t}_{0}}^{t}d(\tau ){\rm{d}}\tau ={\int }_{{t}_{0}}^{t+{T}_{{\rm{R}}}(t)}{e}_{{\rm{NN}}}(\tau ){\rm{d}}\tau ,$$
(12)

and instantaneous interest payments are obtained by multiplication of the quantity of CO2 for which CROs have been issued (d = ϕeNP) and the average abatement costs, cavg, the moment of retirement of the CRO, t + TR(t), with the interest rate on carbon debt rd:

$${i}_{{\rm{inst}}}(t)=d(t){c}_{{\rm{avg}}}(t+{T}_{{\rm{R}}}(t),\,a(t+{T}_{{\rm{R}}}(t))){r}_{{\rm{d}}}.$$
(13)

Integrating and discounting instantaneous interest payments over the repayment term TR gives the total net present interest costs at t for carbon debt d(t):

$${i}_{{\rm{tot}}}(t)={i}_{{\rm{inst}}}(t){\int }_{t}^{t+{T}_{{\rm{R}}}(t)}\exp (-r(\tau -{t}_{0})){\rm{d}}\tau .$$
(14)

Now we add interest costs to the standard objective function (equation (2)) to obtain the optimization problem for a CRO-ETS:

$$\mathop{{\rm{\min }}}\limits_{a(t),\,{c}_{{\rm{tot}}}(t,a(t)),{i}_{{\rm{tot}}}(t)\,}{\int }_{{t}_{0}}^{T}({c}_{{\rm{tot}}}(t,a(t))\exp (-r(t-{t}_{0}))+{i}_{{\rm{tot}}}(t)){\rm{d}}t.$$
(15)

Mitigation cost discounting

If we set r = rd, the objective function can be written as (Supplementary Information section 3):

$$\begin{array}{c}\mathop{{\rm{\min }}}\limits_{a(t),\,{e}_{{\rm{NP}}}(t),d(t),{c}_{{\rm{avg}}}(t,a(t))}{\int }_{{t}_{0}}^{T}(d(t){c}_{{\rm{avg}}}(t+{T}_{{\rm{R}}}(t),\,a(t+{T}_{{\rm{R}}}(t)))\\ +({E}_{{\rm{base}}}(t)-{e}_{{\rm{NP}}}(t)){c}_{{\rm{avg}}}(t,a(t)))\exp (-r(t-{t}_{0})){\rm{d}}t.\end{array}$$
(16)

Notably, instead of pricing eNN, carbon debt d in this new formulation is paid for the moment it is created; however, it is paid for at the average (undiscounted) future costs during removal at t + TR, which is due to our definition of interest costs in equation (13). This is because when we set r = rd, interest payments exactly compensate for the cost reduction in the net present value terms from discounting. The second term, Ebase(t) − eNP(t) equals emission reductions in the net-positive/net-zero domain. These reductions can be achieved by a mix of CDR, low-carbon and zero-carbon technologies; however, CDR is deployed only to offset contemporaneous emissions and not to recapture previously released CO2.

Intuitively, rd therefore controls to what extent cost discounting becomes a driver for accruing carbon debt. If rd equals the market interest rate, future costs cavg at t + TR—which depend on technological learning and the aggregate demand for abatement a in t + TR—determine whether the carbon debt route (d) proves competitive compared to instantaneous emission reductions (Ebase − eNP). However, if d is reduced, eNP needs to be reduced simultaneously to meet the emission target (less carbon debt leads to higher demand for near-term emission reductions and, therefore, an increase in near-term marginal costs). Because near-term emission reductions potentially include CDR, technological and socio-environmental learning associated with CDR is induced earlier, leading to a reduction in the uncertainty, which is key for operating in the net-negative domain later in the century. In this Article, we provide some intuition about the dynamic effects of invoking an interest on carbon debt, but do not determine optimal risk-reducing rates, which could—but do not necessarily need to—coincide with market interest rates. However, we expect, under circumstances in which physical and financial risks associated with carbon debt are managed by appropriately setting an interest rate on carbon debt, that rd is driven by the market interest rate. In our model, an increase in the market interest rate induces deferral of mitigation due to discounting, leading to higher quantities of D, thereby also to an increase in risks and, finally, the necessity to correct rd upwards to account for the increased risks.

Assessing the value of carbon debt

Net present cost gains from discounting are only cancelled exactly if the market interest rate is invoked on abatement costs at t + TR. Costs are known in our model, but are potentially impossible to determine in the context of real emission control policy. Therefore, given their liability for issued debt, managing authorities and financial institutions need to estimate the financial value of carbon debt as a basis for interest payments and CRO maturities. The incentive to correctly value debt has a societal benefit of gradually reducing uncertainty with respect to CDR and other technologies relied on at large scales in mitigation scenarios. Notably, by prudently valuing debt, issuing bodies assure the quality of price signals on carbon markets, instead of relying on the carbon price to value debt. In fact, carbon prices on their own are insufficient benchmarks for valuing debt. For instance, a large demand for carbon debt would lead to a lower near-term carbon price if this is not balanced by an increase in rd, which in turn would lead to an undervaluation of risks.

Supply and demand of CROs

The supply of allowances under a pure ETS is completely inelastic, whereas the supply in a tax system is infinitely elastic. By contrast, the supply of CROs (adding to the supply of allowances) is finitely elastic. Generally, the supply curve is increasing because the larger the demand for CROs, the more abatement is required in the future, making future abatement and thus CROs more expensive. Because the total discounted interest costs (itot) are reflected in the supply curve, by valuing carbon debt and setting the rate rd accordingly, debt-issuing bodies can partly control its slope. The slope, however, determines the level of price volatility, for example, resulting from a demand shock, as depicted in Extended Data Fig. 1. By implication, price volatility is the largest in a pure ETS, and zero in a tax system. By increasing the interest costs in a CRO-ETS, the potential for volatile prices increases, and vice versa. On the other hand, net emissions are fixed in an idealized ETS and subject to demand fluctuations in a tax system. In a CRO-ETS, the cumulative quantity of net emissions is fixed, however, only if default risks are adequately managed.

Numerical solution of the model

The model used to solve the CRO optimization problem (equation (15)) is based on marginal abatement cost curves (MACCs; equation (4)) that are derived from scenarios reported in the SSP scenario database24,54. MACCs are derived for each IAM and SSP by combining and fitting a curve to carbon prices from different RCPs in each time step. For instance, a MACC in 2040 for a specific IAM–SSP configuration is composed of the carbon prices reported for RCP 1.9–RCP 6.0. The set of parameters of each IAM–SSP configuration of our model is therefore composed of net emissions from the baseline scenario (Ebase), the interest rate r (derived from the slope of log-transformed carbon prices), the carbon budget B derived from the sum of net emissions compatible with specific climate targets and eight MACCs for the period from 2030 to 2100, that is, one per decade. We fix abatement rates a during optimization in which no MACCs could be derived because the reported prices pi are (close to) zero over the whole range of ai—that is, in 2020 for all configurations; for IMAGE–SSP 2 in 2030; for IMAGE–SSP 3 in 2030 and 2040; for IMAGE–SSP 5 in 2030. In decades in which abatement is fixed, costs are set to zero. Moreover, for the 1.5 °C and 2 °C case studies, carbon budgets were corrected using historical emission data39 (scenarios reported in the SSP database start in 2005 or 2010 and were exceeded by estimated net emissions in the past decades). Baseline emissions in 2020 were replaced by the projection for 2019 in ref. 39.

We fit the inverse of the generalized logistic function55 to log-transformed prices pi as reported in the SSP database. Abatement rates ai are computed by subtracting net emissions in a scenario with a climate target from the net emissions in the baseline scenario and dividing by the baseline. The index i denotes the different RCPs within the same IAM–SSP configuration and the same year (see Supplementary Information section 2.1 for the cost curves of all IAM–SSP configurations):

$$\mathrm{ln}({p}_{i})=P+\frac{1}{k}\,\mathrm{ln}\left(\frac{1}{\nu }\left({\left(\frac{L-A}{{a}_{i}-A}\right)}^{\nu }-1\right)\right).$$
(17)

MACs are therefore a power law defined for the interval \(\{a\,\in {\mathbb{R}}|L < a < A\}\):

$${\rm{MAC}}({a}_{i})=b{\left(\frac{1}{\nu }\left({\left(\frac{L-A}{{a}_{i}-A}\right)}^{\nu }-1\right)\right)}^{c},$$
(18)

where b = exp(P) and \(c=\frac{1}{k}\). An interpretation of the parameters is provided in Extended Data Fig. 8a. L ≈ 0 (subject to model fitting) and A = max(ai) + ε (such that a can become max(ai) without MAC(a) becoming ∞)—that is, A is set to the maximum abatement (plus ε = 0.01) observed in each decade for each IAM–SSP configuration because this level cannot be exceeded. For numerical reasons, however, a is also constrained by A during optimization such that a < A.

In most IAMs, the carbon prices are either imposed exogenously as driver of mitigation (for example, in the recursive dynamic models AIM–CGE or GCAM4) or prices are derived after optimization from Lagrange multipliers of emission caps (for example, in the intertemporal optimization models MESSAGE–GLOBIOM, WITCH–GLOBIOM or REMIND–MAgPIE). In these cases, the carbon prices typically increase exponentially with the interest rate, as explained by the Hotelling rule22,23. In heavily constrained, detailed process-based IAMs, intertemporally optimal carbon prices are a good proxy of marginal costs; however, they do not necessarily reflect MACs exactly in each point of time, as a consequence of growth constraints or caps on total deployment levels of specific mitigation technologies. This is also the case here because we limit a < A with an additional constraint and we fix a in cases in which no MACCs could be derived, suggesting that MAC and carbon prices (derived from the Lagrange multiplier of the budget constraint (equation (3))) do not necessarily coincide.

We compare MAC and carbon prices from our model with carbon prices as reported for the individual scenarios in the SSP database for all models, SSPs and RCPs (Supplementary Information section 2.2). Reported carbon prices in the database for AIM–CGE do not follow an exact exponential curve because these prices reflect marginal costs from the SSP 1 scenario, which was initially constrained by emission caps to obtain the climate target. Then, prices were manually scaled and imposed on other SSP scenarios to achieve the respective climate targets56. Therefore, AIM–CGE prices are better replicated by the MAC of our model than by carbon prices. The same is true for the IMAGE framework, which contains simulation as well as optimization components and does not report Hotelling-type carbon prices.

Furthermore, we show abatement costs for all SSPs, IAMs and RCPs computed with our model (Supplementary Information section 2.3). Because abatement costs are not explicitly reported in the database, we compare costs from our model with close proxies—that is, GDP loss and consumption loss in SSP scenarios. For GCAM4 and IMAGE, GDP loss and consumption loss are either not reported or losses are close to zero. For GCAM4, we therefore added abatement costs for some scenarios as reported in the supplementary information of ref. 57, which are well replicated by our model. No comparable data could be retrieved for the IMAGE model. For the other IAMs, abatement costs of our model mainly coincide with consumption loss. Net CO2 emissions are also compared for all SSPs, IAMs and RCPs (Supplementary Information section 2.4).

Abatement rates a cannot exceed A, hence only IAM–SSP configurations with RCP 1.9 data are used for our 2 °C case studies, because more-ambitious mitigation under a CRO-ETS requires our model to partly operate in the 1.5 °C abatement domain to achieve 2 °C. Therefore, thirteen IAM–SSP parameter sets of our model are used for the case studies: AIM–CGE (SSP 1 and SSP 2), GCAM4 (SSP 1, SSP 2 and SSP 5), IMAGE (SSP 1), MESSAGE–GLOBIOM (SSP 1 and SSP 2), REMIND–MAgPIE (SSP 1, SSP 2 and SSP 5), WITCH–GLOBIOM (SSP 1 and SSP 4); that is, six parameter sets for SSP 1, four for SSP 2, one for SSP 4 and two for SSP 5 (Extended Data Table 2). All 2 °C scenarios are shown graphically in Supplementary Information section 1.1 and numerically in Supplementary Information section 1.2.

Scenarios for DACS

For our 1.5 °C (RCP 1.9) case study, additional sources of abatement are required to assess compatible pathways of more-ambitious mitigation than suggested by RCP 1.9 scenarios. We therefore add DACS to the mitigation portfolio; however, we treat this technology in a stylized manner as completely stand-alone and independent of other abatement technologies (for example, the energy needs for DACS are assumed to be met by additional local renewable sources that do not interfere with the ramp-up of renewable energy as part of conventional abatement). DACS is less controversial than BECCS with respect to land use and has potentially limited environmental effects compared to other large-scale CDR options16, making it more independently scalable. However, capital and energy requirements are uncertain and potentially enormous. Costs range between US$20 and US$1,000 per t CO2 (refs. 16,18,58,59,60) and potentials for CDR range from 0.5–5 Gt CO2 yr−1 in 2050 to 15–40 Gt CO2 yr−1 in 210018; however, these potentials are mainly constrained by cost considerations rather than biophysical limits61. Here we derive six idealized MACCs for DACS covering three cost ranges and two maximum abatement rates (Extended Data Fig. 8b). Instead of modifying the MACCs derived from SSP scenarios to account for DACS, we add aDACS to equation (1):

$$e(t)={E}_{{\rm{base}}}(t)(1-a(t)-{a}_{{\rm{DACS}}}(t)).$$
(19)

and change the total costs in equation (4) to:

$${c}_{{\rm{tot}}}(t,\,a(t))={E}_{{\rm{base}}}(t)({\int }_{0}^{a(t)}{\rm{MAC}}(t,\,\tilde{a}){\rm{d}}\tilde{a}+{\int }_{0}^{{a}_{{\rm{DACS}}}(t)}{{\rm{MAC}}}_{{\rm{DACS}}}(\bar{a}){\rm{d}}\bar{a}).$$
(20)

Moreover, MACs are always required to be equal:

$${\rm{MAC}}(t,a(t))={{\rm{MAC}}}_{{\rm{DACS}}}({a}_{{\rm{DACS}}}(t))$$
(21)

To obtain a detailed technology downscaling of sources and sinks of CO2 (fossil fuels and industry, including residual emissions from carbon capture and storage; BECCS and land-use emissions) we interpolate linearly between the closest abatement levels reported in the SSP database—that is, ai < a < ai+1 (again, i denotes different RCPs within the same IAM–SSP configuration and the same year)—and add DACS after the interpolation.

The 1.5 °C scenarios used

For Fig. 4, the set of all 468 scenarios (13 IAM–SSP parameter sets, 6 rates rd and 6 DACS parameters sets) is filtered for scenarios that achieved at least a 30% reduction in D compared with their baselines (that is, where rd = 0). For scenarios depicted in Extended Data Figs. 35 this reduction in D needs to be at least 5%, 15% and 45%, respectively. From scenarios with different rates rd but otherwise identical parameters, only the lowest rate is kept, resulting in a potential set of 78 scenarios, of which 26 are feasible regarding the 30% carbon debt reduction requirement in Fig. 4. Hence, in Fig. 4c, d the 26 scenarios for which rd > 0 are compared to the associated 26 baselines in Fig. 4a, b in which rd = 0. Baselines are interpreted as ‘conventional ETS’ scenarios, which are—in terms of emission paths—equivalent to CRO-ETS scenarios with rd = 0. All underlying scenarios are shown graphically in Supplementary Information section 1.3 and numerically in Supplementary Information section 1.4, abatement and interest costs are illustrated in Extended Data Fig. 6 for all scenarios.

A note on technological learning

Technological learning in most IAMs is either exogenous—that is, purely time dependent—or induced by learning-by-doing, which is strongly backed by empirical evidence. However, learning is best perceived as a complex interplay between research and development, learning-by-doing and different types of spillovers62, which only few models attempt to fully address. The MACCs derived here from SSP scenario results reflect learning rates in the IAMs that are used to generate these scenarios, resulting in typically decreasing marginal costs over time for similar abatement rates. Therefore, learning in our model is exogenous (purely time dependent), which is one of the main caveats of this model, because fixed learning rates over time imply an incentive to wait until abatement becomes cheaper. More-ambitious near-term mitigation under a CRO-ETS, however, would probably lead to earlier cost reductions than reflected in the model. Owing to the complexity of learning and the simplicity of our model, we disregard DACS-related technological change.

Software and solver

The model is solved using the CONOPT solver in GAMS v.26.1. CONOPT is based on the generalized reduced gradient algorithm, one of the most robust and commonly applied methods for solving models with highly nonlinear objective functions or constraints63.