Abstract

Projected rates of emissions reductions are unlikely to keep global temperatures from crossing the Paris Agreement temperature targets. Large-scale carbon dioxide removal (CDR) could help recover a target temperature after it has been exceeded, producing an overshoot scenario. Solar radiation modification (SRM) is the idea to cool the planet by increasing the reflection of incoming solar radiation. SRM could be used in an overshoot scenario for ‘peak shaving’, where SRM is deployed to maintain a temperature target during the overshoot. Here, we quantify the effect of SRM peak shaving on the duration of the overshoot using an adapted extension of the SSP2-4.5 scenario and an ensemble of variants of the FaIR simple climate model. We find a substantial reduction in overshoot duration, which ranges from approximately 5% for multi-decade overshoots up to approximately 20% for multi-century overshoots. The shortening is predominantly driven by the ocean response to peak shaving. Peak shaving results in lower ocean temperatures relative to the overshoot scenario, inducing a stronger surface temperature response to decreasing and negative emissions, driving overshoot shortening. Thus, SRM, when deployed as a complement to emissions reductions and CDR, could end overshoot decades earlier than otherwise.

Lay Summary

Recently, scientists have been researching new ways to help tackle climate change. One idea is ‘solar radiation modification’ (SRM), which involves cooling the planet by increasing the reflection of incoming sunlight. ‘Peak shaving’ is a proposed scenario involving SRM, where SRM is used to limit global mean temperature to a target level until it is no longer needed, due to large-scale emissions removals gradually reducing the greenhouse gas concentration in the atmosphere during a period of net negative emissions. This study investigates the effect of peak shaving on the time taken for temperatures to drop below the target temperature without the need for continued SRM. We find that peak shaving reduces this duration by up to approximately 20%, which we term ‘overshoot shortening’. We find that this is driven by avoided increases in ocean temperatures, which impact surface temperatures.

Introduction

Through the Paris Agreement in 2015, almost every nation agreed to limit global average temperature increases to well below 2°C above pre-industrial levels, while pursuing efforts to not exceed 1.5°C [1]. By the end of 2024, global surface temperatures have risen approximately 1.31°C above the 1850 to 1900 baseline, according to the real-time global warming index from Haustein et al. [2]. The planet is projected to surpass the 1.5°C target in less than a decade [3], and temperatures are expected to warm further, reaching an estimated 2.7°C by the end of this century under a moderate emissions (SSP2-4.5) scenario [4]. This is significantly above the upper limit of the Paris Agreement, motivating research into novel climate policy options, including carbon dioxide removal and solar radiation modification [5].

Carbon dioxide removal (CDR) is the anthropogenic removal of CO2 from the atmosphere to long-term storage [6]. CDR in some form is required to drive the reduction in temperature under an overshoot scenario. Several approaches are being explored. Direct air carbon capture and storage (DACCS), which involves extracting CO2 from ambient air using chemical sorbents [7], is arguably the leading technology for large-scale CDR [8]. DACCS is primarily limited by the costs associated with its high energy demand [9], with cost estimates ranging from $100 to $1,000 per ton of carbon removed [10]. CDR technologies such as DACCS currently remove 1.35 MtCO2 per year [11], less than 0.01% of current annual global CO2 emissions [12]. However, large-scale CDR is expected to play a more significant role in future climate strategies [13]. Two of the five illustrative scenarios from the latest IPCC Working Group 1 Assessment Report [4] feature net negative emissions (requiring CDR): SSP1-1.9 and SSP1-2.6, with roughly –14 and –9 GtCO2 per year by 2100, respectively. The maximum potential of global CDR has been estimated to be about 35 GtCO2 per year [14]. The global capacity for underground CO2 storage for CDR, while deeply uncertain [15], has been estimated at between 5,000 and 25,000 GtCO2, which is greater than the expected long-term global demand [16].

Solar radiation modification (SRM), also known as solar geoengineering, is the idea to cool the planet by increasing the reflection of incoming solar radiation back out to space [17]. Stratospheric aerosol injection (SAI) is a leading approach for SRM deployment [18]. It would aim to mimic the cooling effect of large volcanic eruptions by creating a layer of reflective aerosols in the stratosphere [19]. The leading engineering approach for SAI deployment suggests using a fleet of specially designed aircraft to loft aerosols (or their precursor gases) into the stratosphere [20]. The cost of this has been estimated at about $18 billion per year for each 1°C of cooling, which is inexpensive relative to mitigation, adaptation, and CDR technologies [21].

SAI has the potential to offset temperature changes caused by global warming [22], but it would not affect temperature in isolation. It would cause a reduction in global mean precipitation [23], potentially causing worsened droughts in some areas already particularly vulnerable to climate hazards [24]. There would also be several other side effects, including worsened acid rain in some regions [25] and delayed ozone hole recovery [26]. There are several other issues related to the advancement and potential deployment of SRM in general, including concerns over its governability [27], mitigation deterrence or ‘moral hazard’ [28], and termination shock—the harms associated with its sudden cessation [29]. These issues pose significant challenges relating to the feasibility and risks of SRM.

Net negative emissions through CDR would enable temperature overshoot scenarios, in which a temperature target is exceeded but later recovered [30]. Temporary overshoot came to relevance given the apparent inevitability of exceeding the 1.5°C temperature target of the Paris Agreement; most scenarios assessed by the IPCC which achieve 1.5°C by the end of the century initially overshoot that threshold [31]. Even if temperatures return gradually to a given target level, the higher temperatures during the overshoot period would have lasting effects and irreversible consequences, such as extinctions [31]. ‘Peak shaving’ refers to the deployment of SRM to cap temperatures at a certain threshold while greenhouse gas concentrations are reduced [32]. This is commonly framed as a method of ‘buying time’ for mitigation and adaptation [33]. It has been argued to be the only acceptable form of SRM deployment [34], since indefinite deployment would, among other issues, substantially increase the risk of termination shock [35]. However, the plausibility of peak shaving is not guaranteed [36]—it requires the implementation of large-scale and long-term CDR and SRM in tandem, a prospect some consider a high-risk speculation [37]. Acknowledging this, it remains important to explore scenarios involving these emerging technologies to allow for a more holistic analysis of potential pathways to mitigate climate risks.

Previous studies [38, 39] modelling a peak-shaving scenario have shown, but not discussed, that the SRM deployment can be ended before the threshold temperature would have been recovered without peak shaving. We term this phenomenon ‘overshoot shortening’, and this study aims to both quantify this effect and analyse the mechanisms behind it. We consider two potential drivers of overshoot shortening, both related to the broader climatic effects of avoided increases in temperature under peak shaving.

Firstly, the carbon cycle response to peak shaving is proposed as a potential driver. Most models agree that there is a positive carbon–climate feedback due to climate change reducing the net efficiency of carbon sinks [40]. The land and ocean have each removed about a quarter of anthropogenic CO2 emissions from the atmosphere [41], but the efficacy of these sinks is changing as the climate changes. The response of the land carbon sink to climate change is obscured by uncertain compensatory effects [42], while global warming causes a clear reduction in the efficiency of the ocean carbon sink due to strengthened ocean stratification and reduced thermohaline circulation [43]. The avoided increases in temperature following SRM deployment would strengthen carbon sinks through several mechanisms [44, 45]. For example, the relatively lower temperatures would reduce the rate of heterotrophic respiration (i.e. decomposition), increasing soil carbon retention [46]; suppress permafrost thaw, limiting the associated carbon releases [47]; and cool the ocean surface, increasing ocean CO2 solubility [48]. As a result, atmospheric carbon dioxide would likely be reduced under a peak-shaving scenario compared the overshoot scenario, which would drive overshoot shortening.

Secondly, a relative reduction in ocean heat content as a result of peak shaving is another potential driver of overshoot shortening. The ocean plays a key role in regulating climate changes, dampening surface warming by absorbing and storing vast amounts of both heat and CO2 [49]. Ocean warming has accounted for 91% of global heat content increases since 1971 [4]. Ocean heat uptake is controlled by the interaction of several processes, including ocean circulation (currents driven by wind and buoyancy), stratification (layering due to density differences), and ventilation (downwelling from near-surface to deep ocean) [50]. Although literature on the effect of peak shaving on ocean heat content is lacking, large volcanic eruptions have been found to have long-term effects on the climate through their impact on ocean heat uptake [51]. We propose that increases in ocean temperatures would be partially avoided under a peak-shaving scenario, which could drive overshoot shortening.

Methods

The model

This study uses version 2.1.4 of the Finite-amplitude Impulse Response (FaIR) model [52]. It is a reduced complexity climate model, comprising a set of six equations, capable of capturing the global mean temperature response to greenhouse gas and aerosol emissions. FaIR is a zero-dimensional model of globally averaged variables, designed to emulate the behaviour of more complex Earth system models, while avoiding their large computational expense. The model features a carbon cycle but does not represent permafrost thaw. The energy balance model in FaIR comprises three layers, corresponding to a surface, upper ocean, and deep ocean.

This study uses a calibrated, constrained ensemble of model parameters [53] that reproduces both historic observations and assessed climate metrics consistent with the IPCC WG1 Sixth Assessment Report [4]. The ensemble comprises an 841-member posterior created from a 1.6 million-member prior ensemble. 199 members were filtered out because their projections either did not exceed the temperature threshold for overshoot by at least 0.1°C or did not return below it by 2500 (within the time frame of this study). A further 39 members were filtered out due to irregular behaviour of the feedback controller (discussed in Supplementary Text S1), leaving the final ensemble for this study with 603 members. For the entire ensemble, the median of the maximum temperature anomaly without SRM under our adjusted SSP2-4.5 scenario was 2.60°C with a 5th to 95th percentile range of 1.87°C to 3.81°C. The corresponding values for the final ensemble are 2.69°C and 2.16°C to 3.51°C.

A key model parameter for this study is the carbon cycle sensitivity to temperature change, which modulates the strength of the carbon sink based on changes in temperature. The value of this parameter varies across the ensemble, with a median of 2.1 K1, and 5th to 95th percentile range of –1.1 to 5.9 K1, where positive values mean a decrease in strength of carbon sink with increasing temperature. See Leach et al. [52] equations 1–4 for the definition of this parameter, labelled rT. The upper ocean heat transfer coefficient is a second key parameter in this study. It defines the efficacy of heat transfer at the ocean–atmosphere interface, which is between the first and second layers of our three-layer energy balance model—see equations 1–4 of Cummins et al. [54]. This value also varies across the ensemble, with a median of 1.2 Wm2 K1 and a 5th to 95th percentile range of 0.9 to 1.7 Wm2 K1.

The scenarios

In order to analyse the effect of peak shaving, a scenario with a pronounced overshoot is required. We select the SSP2-4.5 scenario as a moderate scenario, consistent with current climate policies [55]. In order to provide an overshoot, we adapted an extension of this scenario from Meinshausen et al. [56], as shown in Supplementary Fig. S1. We adjust their idealised extension so fossil and industry CO2 emissions maintain the trajectory of the standard SSP2-4.5 scenario, reaching net zero in 2130, rather than 2250. This divergence is attributed in part to the implementation of more drastic CDR. Beyond 2130, emissions continue to decline at the same rate, reaching –10 GtCO2 per year in 2155, remaining constant thereafter. This is consistent with the projections for CDR and storage potential detailed in the introduction. These adjustments create an overshoot scenario, where temperatures exceed and recover the upper temperature target of the Paris Agreement [1] within the scenario time frame (up to 2500). Although substantial uncertainties arise when modelling over such an extended time frame [57], this approach allows for a quantitative analysis of peak shaving under one potentially plausible scenario.

Peak shaving requires a threshold temperature at which global temperature anomalies are capped through the use of SRM. For this study, it was set at 2°C of warming relative to the 1850–1900 baseline, corresponding to the upper limit of the Paris Agreement temperature targets [1]. The year-by-year amount of SRM forcing required to keep to the threshold temperature was computed using a feedback controller—the proportional-integral-derivative (PID) algorithm [58]—which is discussed in Supplementary Text S1. SRM was implemented in the model through the volcanic forcing input, which serves as an analogue for SAI forcing. This study compares two scenarios: an overshoot scenario without the deployment of any SRM (‘OS’ from here on), and a peak-shaving scenario where SRM is deployed to maintain the threshold temperature of 2°C during the overshoot (‘PS’ from here on).

Results

Overshoot shortening

Figure 1 shows the temperature responses for a single model configuration to illustrate the concept of this study. The OS temperature anomaly can be seen to exceed the threshold temperature (at Line 1) and later recover it (at Line 3). The PS temperatures are capped until SRM is no longer required in order to keep below the threshold temperature (at Line 2). This occurs before the threshold temperature is recovered under the OS scenario, demonstrating overshoot shortening, which is quantified by the difference between the years indicated by Lines 2 and 3. For this illustrative case, the overshoot of 340 years is shortened by 60 years—an 18% reduction. The SRM forcing, which is implemented through volcanic forcing, deviates from zero in 2042 when the threshold temperature is reached and PS temperatures separate from OS. After reaching a maximum of approximately –5.4 Wm2, it returns to zero when SRM deployment is ended in 2322, aligning with when PS temperatures decline from the threshold temperature.

Surface temperature anomalies (K), relative to the 1850–1900 baseline, for the OS (overshoot without SRM) and PS (peak shaving with SRM) scenarios under one example model configuration. The SRM forcing is shown on the secondary y-axis. The threshold temperature anomaly of 2 K is marked with a dashed line. Vertical dotted lines mark the beginning and ends of the overshoots.
Figure 1.

Surface temperature anomalies (K), relative to the 1850–1900 baseline, for the OS (overshoot without SRM) and PS (peak shaving with SRM) scenarios under one example model configuration. The SRM forcing is shown on the secondary y-axis. The threshold temperature anomaly of 2 K is marked with a dashed line. Vertical dotted lines mark the beginning and ends of the overshoots.

Alternative illustrative temperature profiles are presented in Supplementary Fig. S2 for scenarios where net negative emissions do not continue beyond the end of the overshoot. For the case of net zero emissions from the point of the threshold temperature being recovered under OS, temperatures rebound by approximately 0.05°C for the alternative OS scenario, and approximately 0.12°C for the alternative PS scenario.

Figure 2 shows the temperature profiles for the ensemble. The median year for exceeding the 2°C threshold is 2059, which is consistent with the results for the SSP2-4.5 scenario from the Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations [4]. The OS temperatures exhibit a broad range of overshoots, and the PS temperatures display effective peak shaving. However, the threshold temperature is exceeded by up to 0.2 K for some ensemble members at the beginning of SRM deployment. This is due to a slight deficiency in the feedback controller, as discussed in Supplementary Text S1. For the ensemble (median values provided with the 5th to 95th percentile range), the overshoot of 231 years (91–417 years) is shortened by 40 years (7–91 years)—a 17% (8%–22%) reduction.

Surface temperature anomalies (K) for the 603-member ensemble, relative to the 1850–1900 baseline, for the OS and PS scenarios. The plotted lines are the median temperature anomalies, and the shaded regions cover the 5th to 95th percentile ranges across the ensemble. The threshold temperature is marked with a dashed line.
Figure 2.

Surface temperature anomalies (K) for the 603-member ensemble, relative to the 1850–1900 baseline, for the OS and PS scenarios. The plotted lines are the median temperature anomalies, and the shaded regions cover the 5th to 95th percentile ranges across the ensemble. The threshold temperature is marked with a dashed line.

Figure 3 presents the percentage overshoot shortening under the PS scenario plotted against the OS overshoot duration. Every member of the ensemble has its overshoot shortened under PS compared to OS. Although there is considerable spread, an increasing trend in percentage shortening with overshoot duration can be seen from roughly 5% to 20%. The increasing trend spans from OS overshoot durations of about 60 to 300 years. Judging from this Figure, a 100-year overshoot would be shortened by approximately 8% (2–15%, 5th to 95th percentile range), equivalent to 8 years (2–15 years), and a 300-year overshoot would be shortened by approximately 20% (13–26%), equivalent to 60 years (39–78 years). The colour gradient across the plots displays a horizontal trend, indicating a positive correlation between OS temperature overshoot and overshoot duration. The OS temperature overshoot accounts for approximately 3% of the variance in the vertical spread of data points about the trend line. This suggests the relationship between OS temperature overshoot and PS overshoot shortening is largely indirect, mediated through OS overshoot duration. The PS and OS overshoot durations are compared more directly in Supplementary Fig. S3. OS overshoot durations span a range of 91–417 years (5th to 95th percentile) with median 231 years, which shorten to 84–325 years with median 192 years under the PS scenario.

Percentage reduction in overshoot duration from the OS to the PS scenario versus the duration of the OS overshoot (years). A quadratic trend line has been fitted to the data, and the shaded region covers the 5th to 95th percentile range. Each data point represents an individual member of the ensemble. The colour scale shows the maximum difference between the OS temperature anomaly and the threshold temperature.
Figure 3.

Percentage reduction in overshoot duration from the OS to the PS scenario versus the duration of the OS overshoot (years). A quadratic trend line has been fitted to the data, and the shaded region covers the 5th to 95th percentile range. Each data point represents an individual member of the ensemble. The colour scale shows the maximum difference between the OS temperature anomaly and the threshold temperature.

Carbon cycle response

The atmospheric CO2 concentration profiles under each scenario are shown in Fig. 4. PS and OS follow closely aligned trends, with PS remaining slightly below OS from the onset of SRM deployment. The average difference between the two scenarios increases up to a maximum of 2.5 ppm in 2143 (shortly after peak OS temperatures), and steadily decreases over time thereafter. By 2500, the average difference between the scenarios is minimal, at approximately 0.1 ppm. This difference between the scenarios is due to the lower PS temperatures causing a relative increase in the efficiency of carbon sinks. This response is controlled by the carbon cycle sensitivity to temperature change parameter. This parameter varies across the ensemble and its influence on CO2 concentrations is shown in Supplementary Fig. S4. We find no evidence for a significant correlation between the percentage overshoot shortening and the carbon cycle sensitivity to temperature change parameter (P-value: 0.45). A scatter plot is shown in Supplementary Fig. S5. The variation in overshoot shortening across the ensemble in our study therefore does not depend on the carbon cycle response to peak shaving.

Atmospheric CO2 concentrations (ppm) under the PS and OS scenarios. The plotted lines are the median concentrations, and the shaded regions cover the 5th to 95th percentile ranges. The secondary y-axis shows the mean difference in CO2 concentrations between the scenarios.
Figure 4.

Atmospheric CO2 concentrations (ppm) under the PS and OS scenarios. The plotted lines are the median concentrations, and the shaded regions cover the 5th to 95th percentile ranges. The secondary y-axis shows the mean difference in CO2 concentrations between the scenarios.

Ocean response

Figure 5a shows a negative correlation between the OS overshoot duration and the upper ocean heat transfer coefficient. Figure 5b shows a negative trend between the overshoot shortening and this coefficient, although the spread is fairly broad. The 5th percentile coefficient value of 0.9 Wm2 K1 corresponds approximately to an overshoot duration of 373 years and a shortening of 22%, while the 95th percentile coefficient value of 1.7 Wm2 K1 corresponds approximately to a duration of 115 years and a 10% shortening.

Dependence of OS overshoot duration (a) and percentage reduction in overshoot duration (b) on the upper ocean heat transfer coefficient (Wm−2K−1). The colour scale shows the maximum difference between the OS temperature anomaly and the threshold temperature.
Figure 5.

Dependence of OS overshoot duration (a) and percentage reduction in overshoot duration (b) on the upper ocean heat transfer coefficient (Wm2K1). The colour scale shows the maximum difference between the OS temperature anomaly and the threshold temperature.

The median temperature responses of the three model layers (surface, upper ocean, and deep ocean) are compared in Fig. 6. In a similar format to Fig. 2, the ensemble temperature profiles for each model layer are shown in Supplementary Fig. S6. For OS temperatures, the ocean layers can be seen to follow similar trends as the surface, but with delayed and diminished temperature responses, particularly for the deep ocean. The delay in reaching peak OS temperatures, measured from when surface OS temperatures peak, is 43 years for the upper ocean and 163 years for the deep ocean. PS temperatures remain lower than OS temperatures from the onset of SRM deployment for all three model layers. After 2300, this disparity is more distinct for the ocean layers, particularly the deep ocean, compared to the surface. The ocean layers are also slower to cool than the surface for both OS and PS scenarios. In 2500, the deep ocean layer has greater OS and PS temperature anomalies than the upper ocean layer, which in turn has greater anomalies than the surface layer. Once OS and PS ocean temperatures start to decline, the rate of cooling is greater under OS compared to PS, indicating that the ocean is losing heat less quickly under PS and thus the heat flux to the surface is smaller.

Median temperature profiles (K) relative to average temperatures from 1850–1900 for the OS and PS scenarios and each model layer.
Figure 6.

Median temperature profiles (K) relative to average temperatures from 1850–1900 for the OS and PS scenarios and each model layer.

Figure 7a shows the median cooling under PS compared to the OS scenario for each model layer. The magnitude and timing of peak cooling varies between the layers. For the surface, upper and deep ocean layers respectively, the peak cooling is roughly 0.7 K in 2130, 0.6 K in 2140, and 0.4 K in 2210. At the median year of ending SRM deployment (2250), PS cooling is roughly 0.2 K, 0.3 K and 0.35 K. After this point, greater cooling continues to be seen in the ocean layers, particularly the deep ocean, compared to the surface. This is indicated in panel b by the negative cooling difference after the ‘SRM end’ point. The initially positive values here indicate the ocean layers being cooled less than the surface, while the negative values indicate the opposite.

Panel (a) shows the median temperature difference (‘cooling’) between the OS and the PS scenarios over time for each model layer. The negative values indicate cooler temperatures under the PS scenario. Panel (b) shows the differences in these temperature changes for the upper and deep ocean layers compared to the surface. Negative values correspond to a greater cooling in the layer under the PS scenario compared to the surface. The vertical dotted lines correspond to the median year that SRM deployment is ended.
Figure 7.

Panel (a) shows the median temperature difference (‘cooling’) between the OS and the PS scenarios over time for each model layer. The negative values indicate cooler temperatures under the PS scenario. Panel (b) shows the differences in these temperature changes for the upper and deep ocean layers compared to the surface. Negative values correspond to a greater cooling in the layer under the PS scenario compared to the surface. The vertical dotted lines correspond to the median year that SRM deployment is ended.

Discussion

This study analysed the effect of peak shaving with SRM on the duration of the overshoot with an ensemble of variants of the FaIR model, using an adapted extension of the SSP2-4.5 scenario. Every member of the ensemble had its overshoot shortened by peak shaving, with longer overshoots having a greater proportional shortening up to about 20% (Fig. 3). Overshoot shortening as a result of peak shaving is exhibited in other studies [38, 39]. While these studies did not quantify the shortening, their findings are qualitatively similar to the results of this study (Fig. 3), despite using different scenarios and climate models.

Overshoot shortening implies a shorter SRM deployment period than would be assumed based on the background emission scenario’s projected overshoot duration. A shorter deployment would, compared to a longer one, reduce the exposure to side effects of SRM, such as the impacts on precipitation, acid rain, and ozone hole recovery, as discussed in the introduction. It would also slightly lessen the risk of termination shock, since there would be less time for disruptions to arise [29]. However, in line with the findings of Baur et al. [59], we stress that peak shaving at 2°C under our extension of SSP2-4.5 could still require a multi-century commitment to SRM, despite overshoot shortening.

With peak shaving ending (in the median case) decades before the overshoot, the cumulative net negative emissions by the end of peak shaving would be hundreds of GtCO2 lower than by the end of the overshoot. However, as shown by Supplementary Fig. S2, net negative emissions would need to continue beyond the end of SRM deployment in order to stay below the threshold temperature. Peak shaving would bring temperatures to within the threshold in advance of some of the net negative emissions effort and the associated costs.

Carbon cycle response

The carbon cycle response to peak shaving is explored as a potential driver of overshoot shortening. There is a marginal difference between the atmospheric CO2 concentration profiles of the OS and PS scenarios (Fig. 4). The peak average difference of 2.5 ppm corresponds to emissions of roughly 19.5 GtCO2 [60], slightly over half current annual CO2 emissions [12]. According to the IPCC 6th Assessment Report [4], each 1,000 GtCO2 of cumulative CO2 emissions induces an estimated warming of 0.45°C, meaning this CO2 difference corresponds to less than 0.01°C of temperature difference. The cooling effect of the carbon cycle response to peak shaving is therefore being outweighed by another factor in driving the overshoot shortening. This is supported by the lack of correlation between the carbon cycle sensitivity to temperature change parameter and the overshoot shortening (Supplementary Fig. S5), indicating the carbon cycle response has minimal influence on the overshoot shortening in our results.

However, the FaIR model does not capture key processes relevant to the carbon cycle response to peak shaving. For example, permafrost thaw is not accounted for since the majority of the models it was calibrated on do not include it. The latest IPCC report estimates that each 1°C of global warming would release approximately 11–150 GtCO2 from permafrost thaw [61]. The carbon cycle response to peak shaving would therefore be underestimated, as reasonably significant avoided CO2 emissions are not accounted for. Judging from the peak PS cooling for the median overshoot and the peak mean difference in CO2 concentrations, accounting for permafrost thaw would increase the OS–PS atmospheric CO2 difference by between approximately 40% and 520%. Zhao et al. [45] investigated the carbon cycle response to SAI under a similar scenario (stabilising temperatures under SSP2-4.5 at 1.5°C until 2069) using an earth system model with a permafrost component. They found that an SAI cooling of approximately 0.8°C (0.1°C more than the median of this study) increases cumulative carbon uptake by 92 GtCO2. This value is significantly larger than the corresponding 19.5 GtCO2 figure from this study but falls within the above-stated range when accounting for permafrost thaw.

Ocean response

The ocean response to peak shaving appears to be the primary factor causing the overshoot shortening. There is a clear correlation between the overshoot shortening and the upper ocean heat transfer coefficient (Fig. 5b). However, there also is a clear correlation between the coefficient and the OS overshoot duration (Fig. 5a), which is in turn correlated with the overshoot shortening (Fig. 3). The coefficient may therefore only have an impact on overshoot shortening through its relationship with the OS overshoot duration.

The ocean response can also be investigated by inspecting the temperature profiles of the ocean layers (Fig. 6). The temperature effects of peak shaving extend beyond the surface to the ocean, which could be driving overshoot shortening through the ocean’s relationship with surface temperatures. Under the OS scenario, ocean temperatures increase and then decrease, implying the ocean changes from a heat sink (energy transfers from the atmosphere to the ocean) to a heat source (energy transfers from the ocean to the atmosphere). While acting as a heat source, the ocean inhibits the decrease in surface temperatures in response to net negative emissions. Under the PS scenario, ocean temperatures remain lower, peak later, and, once they decrease, do so at a slower rate compared to OS. The ocean under PS therefore acts as a heat sink for longer and subsequently as a weaker heat source. This induces a stronger surface temperature response to net negative emissions, as ocean heat uptake is prolonged and a weaker ocean heat source results in a weaker inhibition of decreasing global temperatures. This causes faster cooling, causing overshoot shortening. After the end of SRM deployment, PS surface temperature remains markedly lower than OS for centuries due to the reduced ocean temperatures (Fig. 7).

This can also be understood through the Held et al. [62] framing of fast and slow components of global warming. In the case of SRM deployment, there is a fast component of cooling as the surface quickly responds to the negative forcing. The slow component relates to the ocean heat content, whereby the negative forcing causes an anomaly in ocean heat storage, which affects surface temperatures over a longer time frame. A similar effect is explored by Stenchikov et al. [51], who note that changes in ocean heat uptake caused by large volcanic eruptions affect the climate over decades, separately to the short-term surface temperature response to volcanic forcing. This presents a parallel to SRM deployment, and our results indicate a slow cooling component in addition to the fast surface response.

Our ocean response results are limited by the simplicity of the FaIR model. The energy balance model in FaIR comprises just three layers, representing a surface, upper ocean, and deep ocean. This broad simplification fails to capture many important complexities related to the ocean’s response in this study. For example, changes in ocean circulations, such as the Atlantic Meridional Overturning Circulation (AMOC), are not adequately represented. The AMOC is sensitive to temperature changes and, if disturbed, would impact ocean heat uptake [63] and ocean carbon uptake [64]. The AMOC plays a significant role in the climate’s nonlinear response to an overshoot scenario [65], which would also be relevant to this study.

Limitations and further research

The complexity of the FaIR model is low, even compared to other simple climate models [52]. While it is able to emulate the behaviour of more complex models, there are drawbacks to its simplicity. It is highly parameterised and simulates globally averaged systems, meaning many processes and regional differences are not captured in the model. Also, the scenario, as shown in Supplementary Fig. S1, follows a single emissions pathway. We expect our findings to remain qualitatively similar across different overshoot scenarios, given the physical mechanisms underlying the response. However, different pathways would lead to overshoots with different OS temperature profiles, likely impacting the quantitative response to peak shaving.

The modelled scenarios depend on large-scale CDR and SRM, neither of which are certain prospects, and both remain distant from large-scale deployment. Furthermore, our scenario spans a multi-century time frame, in contrast to most climate projections, which finish at 2100 [66]. While this is necessary for an analysis of multi-century overshoots, the scenario uncertainty increases substantially over such an extended time frame [57].

Further research could conduct the analysis of this study using more and more complex climate models. This could be carried out in future phases of the Reduced Complexity Model Intercomparison Project (RCMIP) [67] or the Geoengineering Model Intercomparison Project (GeoMIP) [68], which could be achieved by adapting and extending the newly proposed G6-1.5K-SAI experiment [69]. Future studies using earth system models [70] could address many of the limitations of the FaIR model. Future work could also model a broader range of emissions scenarios and threshold temperatures to further explore overshoot shortening. Lastly, further insights regarding the ocean response to peak shaving could be gained from the analysis of energy fluxes into and out of the ocean.

Conclusions

This study has analysed the shortening of global temperature overshoot under a peak-shaving scenario with SRM deployment. We used an ensemble of variants of the FaIR simple climate model and an adapted extension of the SSP2-4.5 scenario. In the absence of SRM, temperature overshoots 2°C in this scenario for 91–417 years (5th to 95th percentile) with median 231 years. With SRM deployed to limit warming to 2°C, this is shortened to 84–325 years with median 192 years. We find that the overshoot shortening ranges from approximately 5% for multi-decade overshoots up to approximately 20% for multi-century overshoots. Our results thus suggest that deploying SRM alongside emissions cuts and CDR would substantially accelerate the end of temperature overshoot.

This study has also explored the mechanisms behind overshoot shortening through peak shaving. The carbon cycle response to peak shaving is relatively minor in this study, and it was determined that its role in driving the overshoot shortening is minimal. Instead, we attribute the overshoot shortening effect primarily to the ocean temperature response. This is due to peak shaving resulting in comparatively lower ocean temperatures, which induces a stronger surface temperature response to decreasing and negative emissions. This causes expedited cooling due to peak shaving, resulting in overshoot shortening. This study provides an initial analysis of overshoot shortening through peak shaving with SRM and serves as a basis for future studies to explore the subject further.

Acknowledgments

The authors would like to thank Chris Smith (Vrije Universiteit Brussel) for his timely and helpful support with FaIR model functionalities and troubleshooting.

Author contributions

Linus Boselius (Data curation [lead], Formal analysis [lead], Methodology [equal], Software [lead], Visualization [lead], Writing—original draft [lead], Writing—review & editing [equal]), Alistair Duffey (Methodology [equal], Project Administration [supporting], Software [supporting], Supervision [supporting], Writing—review & editing [equal]), and Peter J. Irvine (Conceptualization [lead], Methodology [equal], Project administration [lead], Supervision [lead], Writing—review & editing [equal])

Supplementary data

Supplementary data is available at Oxford Open Climate Change online.

Conflict of interest: None declared.

Funding

LB’s contribution was in part funded by the University College London MAPS Summer Research Internship scheme. AD’s contribution was funded by the London Natural Environment Research Council (NERC) Doctoral Training Partnership (DTP) Grant NE/S007229/1.

Data availability

The FaIR model and instructions for its use are publicly available at https://docs.fairmodel.net. SSP scenario data are publicly available at https://tntcat.iiasa.ac.at/SspDb. The Meinshausen et al. (2020) scenario extensions are publicly available at https://greenhousegases.science.unimelb.edu.au

References

1

UNFCCC
. Adoption of the Paris agreement. In 21st Conference of the Parties. Paris: United Nations Framework Convention on Climate Change,
2015
.

2

Haustein
K
,
Allen
MR
,
Forster
PM
 et al.  
A real-time global warming index
.
Sci Rep
 
2017
;
7
:
15417
.

3

Matthews
HD
,
Wynes
S.
 
Current global efforts are insufficient to limit warming to 1.5 C
.
Science
 
2022
;
376
:
1404
9
.

4

IPCC
. Summary for policymakers. In: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK: Cambridge University Press,
2021
.

5

Fawzy
S
,
Osman
AI
,
Doran
J
,
Rooney
DW.
 
Strategies for mitigation of climate change: a review
.
Environ Chem Lett
 
2020
;
18
:
2069
94
.

6

IPCC
. Annex VII: glossary. In: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK: Cambridge University Press,
2021
, pp.
2215
56
.

7

Gambhir
A
,
Tavoni
M.
 
Direct air carbon capture and sequestration: how it works and how it could contribute to climate-change mitigation
.
One Earth
 
2019
;
1
:
405
9
.

8

Cobo
S
,
Negri
V
,
Valente
A
 et al.  
Sustainable scale-up of negative emissions technologies and practices: where to focus
.
Environ Res Lett
 
2023
;
18
:
023001
.

9

Fuss
S
,
Lamb
WF
,
Callaghan
MW
 et al.  
Negative emissions—Part 2: costs, potentials and side effects
.
Environ Res Lett
 
2018
;
13
:
063002
.

10

Meckling
J
,
Biber
E.
 
A policy roadmap for negative emissions using direct air capture
.
Nat Commun
 
2021
;
12
:
2051
.

11

Smith
S
,
Geden
O
,
Gidden
M
 et al. The State of Carbon Dioxide Removal. Smith School of Enterprise and the Environment, Univ. Oxford,
2024
.

12

IEA
. CO2 Emissions in 2023. Paris: International Energy Agency,
2024
. https://www.iea.org/reports/co2-emissions-in-2023

13

Nemet
GF
,
Callaghan
MW
,
Creutzig
F
 et al.  
Negative emissions—part 3: innovation and upscaling
.
Environ Res Lett
 
2018
;
13
:
063003
.

14

Psarras
P
,
Krutka
H
,
Fajardy
M
 et al.  
Slicing the pie: how big could carbon dioxide removal be?
 
WIREs Energy Environment
 
2017
;
6
:
e253
.

15

Lane
J
,
Greig
C
,
Garnett
A.
 
Uncertain storage prospects create a conundrum for carbon capture and storage ambitions
.
Nat Clim Chang
 
2021
;
11
:
925
36
.

16

Kelemen
P
,
Benson
SM
,
Pilorgé
H
 et al.  
An overview of the status and challenges of CO2 storage in minerals and geological formations
.
Front Clim
 
2019
;
1
:
9
.

17

Irvine
PJ
,
Kravitz
B
,
Lawrence
MG
,
Muri
H.
 
An overview of the Earth system science of solar geoengineering
.
WIREs Clim Change
 
2016
;
7
:
815
33
.

18

Kravitz
B
,
MacMartin
DG.
 
Uncertainty and the basis for confidence in solar geoengineering research
.
Nat Rev Earth Environ
 
2020
;
1
:
64
75
.

19

Hulme
M.
 
Climate change: climate engineering through stratospheric aerosol injection
.
Progr Phys Geograph
 
2012
;
36
:
694
705
.

20

Lockley
A
,
MacMartin
D
,
Hunt
H.
 
An update on engineering issues concerning stratospheric aerosol injection for geoengineering
.
Environ Res Commun
 
2020
;
2
:
082001
.

21

Smith
W.
 
The cost of stratospheric aerosol injection through 2100
.
Environ Res Lett
 
2020
;
15
:
114004
.

22

Kravitz
B
,
MacMartin
DG
,
Tilmes
S
 et al.  
Comparing surface and stratospheric impacts of geoengineering with different SO2 injection strategies
.
Jgr Atmospheres
 
2019
;
124
:
7900
18
.

23

Ferraro
AJ
,
Griffiths
HG.
 
Quantifying the temperature-independent effect of stratospheric aerosol geoengineering on global-mean precipitation in a multi-model ensemble
.
Environ Res Lett
 
2016
;
11
:
034012
.

24

Abiodun
BJ
,
Odoulami
RC
,
Sawadogo
W
 et al.  
Potential impacts of stratospheric aerosol injection on drought risk managements over major river basins in Africa
.
Climatic Change
 
2021
;
169
:
1
19
.

25

Visioni
D
,
Slessarev
E
,
MacMartin
DG
 et al.  
What goes up must come down: impacts of deposition in a sulfate geoengineering scenario
.
Environ Res Lett
 
2020
;
15
:
094063
.

26

Tilmes
S
,
Richter
JH
,
Kravitz
B
 et al.  
Sensitivity of total column ozone to stratospheric sulfur injection strategies
.
Geophys Res Lett
 
2021
;
48
:e2021GL094058.

27

Reynolds
JL.
 
Is solar geoengineering ungovernable? A critical assessment of governance challenges identified by the Intergovernmental Panel on Climate Change
.
WIREs Clim Change
 
2021
;
12
:
e690
.

28

McLaren
D.
 
Mitigation deterrence and the “moral hazard” of solar radiation management
.
Earth Future
 
2016
;
4
:
596
602
.

29

Parker
A
,
Irvine
PJ.
 
The risk of termination shock from solar geoengineering
.
Earth Future
 
2018
;
6
:
456
67
.

30

Huntingford
C
,
Lowe
J.
 
“Overshoot” scenarios and climate change
.
Science
 
2007
;
316
:
829
829
.

31

Reisinger
A
,
Geden
O.
 
Temporary overshoot: origins, prospects, and a long path ahead
.
One Earth
 
2023
;
6
:
1631
7
.

32

Horton
JB.
 
The emergency framing of solar geoengineering: time for a different approach
.
Anthropocene Rev
 
2015
;
2
:
147
51
.

33

Neuber
F
,
Ott
K.
 
The buying time argument within the solar radiation management discourse
.
Appl Sci
 
2020
;
10
:
4637
.

34

Ott
KK.
 
On the political economy of solar radiation management
.
Front Environ Sci
 
2018
;
6
:
43
.

35

Low
S
,
Honegger
M.
 
A precautionary assessment of systemic projections and promises from sunlight reflection and carbon removal modeling
.
Risk Anal
 
2022
;
42
:
1965
79
.

36

Schleussner
C-F
,
Ganti
G
,
Lejeune
Q
 et al.  
Overconfidence in climate overshoot
.
Nature
 
2024
;
634
:
366
73
.

37

Asayama
S
,
Hulme
M.
 
Engineering climate debt: temperature overshoot and peak-shaving as risky subprime mortgage lending
.
Clim Policy
 
2019
;
19
:
937
46
.

38

Cao
L
,
Jiang
J.
 
Simulated effect of carbon cycle feedback on climate response to solar geoengineering
.
Geophys Res Lett
 
2017
;
44
:
12
484
.

39

MacMartin
DG
,
Ricke
KL
,
Keith
DW.
 
Solar geoengineering as part of an overall strategy for meeting the 1.5 C Paris target
.
Phil Trans R Soc A
 
2018
;
376
:
20160454
.

40

Friedlingstein
P
,
Prentice
IC.
 
Carbon–climate feedbacks: a review of model and observation based estimates
.
Curr Opin Environ Sustain
 
2010
;
2
:
251
7
.

41

Friedlingstein
P
,
O'Sullivan
M
,
Jones
MW
 et al.  
Global carbon budget 2023
.
Earth Syst Sci Data
 
2023
;
15
:
5301
69
.

42

Huntzinger
DN
,
Michalak
AM
,
Schwalm
C
 et al.  
Uncertainty in the response of terrestrial carbon sink to environmental drivers undermines carbon–climate feedback predictions
.
Sci Rep
 
2017
;
7
:
4765
.

43

Friedlingstein
P.
 
Carbon cycle feedbacks and future climate change
.
Phil Trans R Soc A
 
2015
;
373
:
20140421
.

44

Plazzotta
M
,
Seferian
R
,
Douville
H.
 
Impact of solar radiation modification on allowable CO2 emissions: What can we learn from multimodel simulations?
 
Earth Future
 
2019
;
7
:
664
76
.

45

Zhao
M
,
Cao
L
,
Visioni
D
,
MacMartin
DG.
 
Carbon cycle response to stratospheric aerosol injection with multiple temperature stabilization targets and strategies
.
Earth Future
 
2024
;
12
:e2024EF004474.

46

Tjiputra
JF
,
Grini
A
,
Lee
H.
 
Impact of idealized future stratospheric aerosol injection on the large-scale ocean and land carbon cycles
.
Jgr Biogeosciences
 
2016
;
121
:
2
27
.

47

Chen
Y
,
Moore
JC
,
Ji
D.
 
Simulated responses and feedbacks of permafrost carbon under future emissions pathways and idealized solar geoengineering scenarios
.
Environ Res Lett
 
2024
;
19
:
024050
.

48

Cao
L.
 
The effects of solar radiation management on the carbon cycle
.
Curr Clim Change Rep
 
2018
;
4
:
41
50
.

49

Reid
PC
,
Fischer
AC
,
Lewis-Brown
E
 et al.  
Impacts of the oceans on climate change
.
Adv Mar Biol
 
2009
;
56
:
1
150
.

50

Marshall
DP
,
Zanna
L.
 
A conceptual model of ocean heat uptake under climate change
.
Journal of Climate
 
2014
;
27
:
8444
65
.

51

Stenchikov
G
,
Delworth
TL
,
Ramaswamy
V
 et al.  
Volcanic signals in oceans
.
J Geophys Res
 
2009
;
114
:D16104.

52

Nicholas
J
,
Leach
S
,
Jenkins
Z
,
Nicholls
 et al.  
A generalized impulse response model for climate uncertainty and future scenario exploration
.
Geosci Model Dev
 
2021
;
14
:
3007
36
.

53

Smith
C.
 Fair calibration data version 1.4.1.
Zenodo
.
2024
. .

54

Cummins
DP
,
Stephenson
DB
,
Stott
PA.
 
Optimal estimation of stochastic energy balance model parameters
.
Journal of Climate
 
2020
;
33
:
7909
26
.

55

Burgess
MG
,
Becker
SL
,
Fredston
A
,
Brooks
CM.
Perspectives on most plausible climate futures, and recommendations for using scenarios in fisheries and aquatic conservation research. Technical report, Center for Open Science,
2022
.

56

Meinshausen
M
,
Nicholls
ZRJ
,
Lewis
J
 et al.  
The shared socio-economic pathway (SSP) greenhouse gas concentrations and their extensions to 2500
.
Geosci Model Dev
 
2020
;
13
:
3571
605
.

57

Latif
M.
 Uncertainty in climate change projections. In:
Reframing the Problem of Climate Change
.
Routledge
,
2013
,
17
34
.

58

Mj
W.
 
Proportional-integral-derivative control
.
Dept. Chem Process Eng Univ Newcastle
 
1999
;
6
.

59

Baur
S
,
Nauels
A
,
Nicholls
Z
 et al.  
The deployment length of solar radiation modification: an interplay of mitigation, net-negative emissions and climate uncertainty
.
Earth Syst Dynam
 
2023
;
14
:
367
81
.

60

O’Hara
FM.
 
Glossary: Carbon Dioxide and Climate
. Oak Ridge, Tennessee:
Oak Ridge National Laboratory
,
1990
.

61

IPCC
. Global carbon and other biogeochemical cycles and feedbacks. In: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom: Cambridge University Press,
2021
, pp.
673
816
.

62

Held
IM
,
Winton
M
,
Takahashi
K
 et al.  
Probing the fast and slow components of global warming by returning abruptly to preindustrial forcing
.
J Clim
 
2010
;
23
:
2418
27
.

63

Kostov
Y
,
Armour
KC
,
Marshall
J.
 
Impact of the Atlantic meridional overturning circulation on ocean heat storage and transient climate change
.
Geophys Res Lett
 
2014
;
41
:
2108
16
.

64

Liu
Y
,
Moore
JK
,
Primeau
F
,
Wang
WL.
 
Reduced CO2 uptake and growing nutrient sequestration from slowing overturning circulation
.
Nat Clim Chang
 
2023
;
13
:
83
90
.

65

Palter
JB
,
Frölicher
TL
,
Paynter
D
,
John
JG.
 
Climate, ocean circulation, and sea level changes under stabilization and overshoot pathways to 1.5 K warming
.
Earth Syst Dynam
 
2018
;
9
:
817
28
.

66

Boulton
CA.
 
Looking to the (far) future of climate projection
.
Glob Chang Biol
 
2022
;
28
:
346
8
.

67

Nicholls
Z
,
Meinshausen
M
,
Lewis
J
 et al.  
Reduced complexity model intercomparison project phase 2: synthesizing earth system knowledge for probabilistic climate projections
.
Earth Future
 
2021
;
9
:
e2020
001900
.

68

Kravitz
B
,
Robock
A
,
Tilmes
S
 et al.  
The geoengineering model intercomparison project phase 6 (GeoMIP6): Simulation design and preliminary results
.
Geosci Model Dev
 
2015
;
8
:
3379
92
.

69

Visioni
D
,
Robock
A
,
Haywood
J
 et al.  
G6-1.5 K-SAI: a new Geoengineering Model Intercomparison Project (GeoMIP) experiment integrating recent advances in solar radiation modification studies
.
Geosci Model Dev
 
2024
;
17
:
2583
96
.

70

Flato
GM.
 
Earth system models: an overview
.
WIREs Clim Change
 
2011
;
2
:
783
800
.

71

Nishikawa
Y
,
Sannomiya
N
,
Ohta
T
,
Tanaka
H.
 
A method for auto-tuning of PID control parameters
.
Automatica
 
1984
;
20
:
321
32
.

72

Somefun
OA
,
Akingbade
K
,
Dahunsi
F.
 
The dilemma of PID tuning
.
Ann Rev Control
 
2021
;
52
:
65
74
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Supplementary data