Skip to main content

Explaining the eventual transient saturation of climate-carbon cycle feedback

Abstract

Background

Coupled climate-carbon cycle simulations generally show that climate feedbacks amplify the buildup of CO2 under respective anthropogenic emission. The effect of climate-carbon cycle feedback is characterised by the feedback gain: the relative increase in CO2 increment as compared to uncoupled simulations. According to the results of the recent Coupled Climate-Carbon Cycle Model Intercomparison Project (C4MIP), the gain is expected to increase during the 21st century. This conclusion is not supported by the climate model developed at the A.M. Obukhov Institute of Atmospheric Physics at the Russian Academy of Sciences (IAP RAS CM). The latter model shows an eventual transient saturation of the feedback gain. This saturation is manifested in a change of climate-carbon cycle feedback gain which grows initially, attains a maximum, and then decreases, eventually tending to unity.

Results

Numerical experiments with the IAP RAS CM as well as an analysis of the conceptual framework demonstrate that this eventual transient saturation results from the fact that transient climate sensitivity decreases with time.

Conclusion

One may conclude that the eventual transient saturation of the climate-carbon cycle feedback is a fundamental property of the coupled climate-carbon system that manifests itself on a relevant time scale.

Background

Starting from the works [1, 2], climate-carbon cycle interactions in the global numerical models attain a lot of scientific attention. These and subsequent works [3–12] have found that an interactive coupling between climate and carbon cycle enhances the build up of the carbon dioxide in the atmosphere, C, in comparison to the hypothetical case when the carbon cycle does not respond to the climate changes. If the changes of the atmospheric carbon dioxide in these two cases are Δ q CO 2 c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaem4yamMaemiCaaNaemiBaWgaaaaa@3648@ and Δ q CO 2 u c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaaaa@37BB@ , respectively, then the parameter of the climate-carbon cycle feedback is

f = Δ q CO 2 c p l / Δ q CO 2 u c p l . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOzayMaeyypa0JaeuiLdqKaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaem4yamMaemiCaaNaemiBaWgaaOGaei4la8IaeuiLdqKaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaOGaeiOla4caaa@46B5@
(1)

In [5], the following expression for climate-carbon cycle feedback gain g = (f - 1)/f is derived

g = -α (γ l + γ oc )/(c0 + β l + β oc ) = -αG (2)

with transient temperature sensitivity α = d T g c p l / d q CO 2 c p l , T g c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeqySdeMaeyypa0JaemizaqMaemivaq1aa0baaSqaaiabdEgaNbqaaiabdogaJjabdchaWjabdYgaSbaakiabc+caViabdsgaKjabdghaXnaaBaaaleaacqqGdbWqcqqGpbWtdaWgaaadbaGaeGOmaidabeaaliabdogaJjabdchaWjabdYgaSbqabaGccqGGSaalcqWGubavdaqhaaWcbaGaem4zaCgabaGaem4yamMaemiCaaNaemiBaWgaaaaa@499E@ is globally and annually averaged surface air temperature in the coupled simulation, and c0 = 2.1 PgC/ppmv. Coefficients β X and γ X come from linear decomposition of respective differential (i.e., per year) terrestrial and oceanic differential carbon uptakes (X = l or X = oc respectively):

F X = β X d q CO 2 d t + γ X d T g d t , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aaSbaaSqaaiabdIfaybqabaGccqGH9aqpcqaHYoGydaWgaaWcbaGaemiwaGfabeaajuaGdaWcaaqaaiabdsgaKjabdghaXnaaBaaabaGaee4qamKaee4ta80aaSbaaeaacqaIYaGmaeqaaaqabaaabaGaemizaqMaemiDaqhaaOGaey4kaSIaeq4SdC2aaSbaaSqaaiabdIfaybqabaqcfa4aaSaaaeaacqWGKbazcqWGubavdaWgaaqaaiabdEgaNbqabaaabaGaemizaqMaemiDaqhaaOGaeiilaWcaaa@4871@

where t stands for time.

For temporal variations of the parameter f in the course of a given emission scenario, [13] demonstrated a monotonic increase of f during the long-term integration of the Lawrence Livermore National Laboratory coupled model. In the framework of the Coupled Climate-Carbon Cycle Model Intercomparison Project (C4MIP) [14], time variations of the respective gains g have been studied for the 21st sentury (their Fig. 2b). In the latter case, a tendency for increase of this gain during the course of integration forced by the emission scenario SRES A2 has been exhibited as well. However, according to Fig. 2b in [14], one model exhibits decrease of g in the late 21st century and a few others show a slower decrease in the second part of this century in comparison to the first part.

Figure 1
figure 1

Parameter of climate-carbon cycle interaction as a function of t / t p , obtained in simulations with the IAP RAS CM versions FocL and FocNL (gray and black curves correspondingly) for t p = 50 yr (thin solid lines), t p = 100 yr (thin dashed lines), t p = 150 yr (thick solid lines), and t p = 250 yr (thick dashed lines). In this and subsequent Figures, the data beyond t = 7 × t p are not shown because they either out of the simulation length or start to show signs of numerical instability. In addition, data with Δ q CO 2 u c p l , c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWMaeiilaWIaem4yamMaemiCaaNaemiBaWgaaaaa@3CB4@ < 5 ppmv are eliminated from plots.

Figure 2
figure 2

a) Coefficients β l and β oc (solid and dashed curves respectively) computed for 100-year running segments from the simulations with the IAP RAS CM versions FocL and FocNL (gray and black lines correspondingly) as a function of t / t p . b) Similar to a) but for coefficients γ l and γ oc .

In turn, a non-monotonic change of f in the simulation with the IAP RAS climate model of intermediate complexity (IAP RAS CM) under the same scenario SRES A2 was obtained [11, 15]. In this simulation, parameter f grows during the most part of the run but starts to decrease late in the 21st century. Futher, based on a conceptual climate-carbon cycle model, it has been demonstarted that plausible physical reason for this eventual transient saturation is a weak, logarithmic dependence of the CO2 radiative forcing on the atmospheric concentration of carbon dioxide [16]. This leads to small influence of an additional (due to climate-carbon cycle interactions) build up of the carbon dioxide on climate state when C is large. Given a length of integration, the latter condition is fulfilled only if emissions are aggressive enough.

However, the build up of the carbon dioxide in the atmosphere in the uncoupled simulation was prescribed in [16] and it is unclear how to relate directly time scales associated with this build up with time scales associated with changes in CO2 emissions. In the present paper, ensemble simulations with an Earth system model IAP RAS CM [10–12, 15, 17] forced by idealised emission scenarios are performed. In this, it is demonstrated that an eventual transient climate-carbon cycle saturation may be exhibited also under moderate emission scenarios if respective integration proceeds for a sufficient time. The results obtained in [16] are supported and it is shown that the above-mentioned hypothesis is proved and even weak but continuing emissions lead to eventual saturation of the climate-carbon cycle feedback. In addition, the IAP RAS CM simulations are supported by integrations with a conceptual coupled model. The latter model is similar, but not identical to that used in [16]. In particular, in the model used here, emissions of CO2 are used instead of prescribing the atmospheric carbon dioxide build up in the uncoupled simulation. It is also demonstrated that eventual transient saturation of climate-carbon cycle feedback in the IAP RAS CM is consistent with that derived from the earlier version of conceptual model.

Results and Discussion

Analysis of a numerical coupled model

A series of numerical experiments with two IAP RAS CM versions were performed. In one version, oceanic uptake of anthropogenic carbon is formulated as a bilinear function of time derivatives of atmospheric concentration of CO2 and globally averaged annual mean sea surface temperature. Hereafter, this version is denoted as FocL. In other version, oceanic uptake of carbon F oc is determined employing a nonlinear model [18] but with chemical constants computed as functions of temperature in accordance to [19]. In particular, this version takes into account carbonate dissolution in the ocean. Hereafter, this version is denoted as FocNL. Terrestrial uptakes for both model versions are determined based on a zero-dimensional module taking into account direct plants fertilisation by CO2 and influence of climate (expressed via temperature anomaly from the model's preindustrial state employing Q10-relationships) on biospheric productivity and autotrophic and heterotrophic respiration. Half-saturation point for the direct plant fertilisation was set to 150 ppmv in the version FocL and to 460 ppmv in the FocNL. Performance and extensive comparison of these two model versions between each other and with observational estimates are described in [10–12, 15, 17]. There, it is shown that both versions behave realistically in the 20th century and their projections for the 21st century are basically in line with those obtained in the C4MIP project. The simulations performed here were forced by idealised emission scenarios. In these scenarios, fossil fuel emissions increase exponentially in time with a time scale t p from 25 yr to 250 yr starting from the small initial value E0 = 0.1 PgC/yr. These emissions are not meant to represent any historical or future projected emissions. However, one notes that the combined fossil fuel emissions taken as historical for the 19–20th centiries [20] and adopted from the Special Report on Emission Scenarios [21] for the 21st century (except the scenario B1 where fossil fuel emissions decline in the late 21st century) may be approximated by an exponential curve with t p changing from 50 yr to 200 yr depending on scenario. One notes that an emission intensity depends both on E0 and t p . Thus, strictly speaking, quantitative results presented below are only valid for this particular choice of E0. However, qualitative conclusions are unlikely to be changed if other sufficiently small value of E0 is selected. In some sense, larger t p may be partly compensated by larger E0. In the performed simulations, land use emissions are neglected. Every simulation starts from the model's preindustrial equilibrated state and proceeds for 1, 500 yr.

As shown in Fig. 1, climate-carbon cycle feedback parameter attains maximum and than eventualy falls down tending to unity. This behaviour is in general agreement with those obtained in the earlier IAP RAS CM simulations [11, 15] and employing a conceptual climate-carbon cycle model [16]. There is a difference in timing t m when f is at maximum between the two employed here IAP RAS CM versions. For the studied here range of t p , t m ranges from 3 × t p to 5 × t p in simulations with the FocL and from 5 × t p to 8 × t p in simulations with the FocNL. For both model versions, the larger t p the earlier (in units of t p ) this maximum of f occurs. Difference in responses of the globally and annually averaged surface air temperatures between coupled and uncoupled simulation also increases early in the simulations, then attains a maximum (which occurs later than t m ) and diminishes afterwards tending to zero (not shown). To diagnose the behavior of the model, an approach by [5] has been adopted. In this, coefficients γ X and β X (X = l, oc) entering (2) were computed for the running 100-yr segments. For both model versions, β l increases till the plants fertilisation half-saturation point is reached and then decreases to zero (see Fig. 2a as an example). This behaviour reflects the Michaelis-Menten-type dependence of terrestrial productivity on q CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aaSbaaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@30C8@ as it is implemented in the IAP RAS CM. Coefficient β oc behaves similarly in the version FocNL while changes in time only slightly in the version FocL. This is expected from the formulations of F oc in these model versions. In contrast, the magnitudes of γ l and γ oc increse during the course of integration (Fig. 2b).

This behaviour of β s and γ s obviously leads to the magnitude of G (see (2)) increasing with time if the latter is large enough (perticularly, after t m , not shown). Consequently, general decrease of gain g and feedback parameter f is due to drastic decrease of transient temperature sensitivity α in the course of integration. Being computed for 100-yr running segments, α shows almost monotonic decrease during the course of integration by an order of magnitude (Fig. 3). This decrease overcompensates an increase of G with time. Such a decrease of sensitivity is again consistent with a commonly accepted weak, logarithmic dependence of the carbon dioxide instantaneous radiative forcing on the CO2 atmospheric content. To get an impression, how large q CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aaSbaaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@30C8@ should be for eventual transient saturation, the time when f eventually becomes smaller than 1.05 is computed and the corresponding value q CO 2 , s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aaSbaaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaSGaeiilaWIaem4Camhabeaaaaa@3317@ is stored. This value ranges between 3226 ppmv and 6152 ppmv (3679 ppmv and 12799 ppmv) for the version FocL (FocNL). The smaller values in these ranges correspond to the larger values of t p . While the upper part of these ranges seems unrealistic (especially that obtained for the version FocNL), the lower part can be achieved if future CO2 emission come along the worst (e.g., SRES A2) path. It is important with this respect that this lower part of range is obtained for t p = 150 – 250 yr which is applicable for the SRES scenarios.

Figure 3
figure 3

Transient temperature sensitivity a computed for 100-year running segments from the simulations with the IAP RAS CM versions FocL and FocNL with different t p . Notations are the same as in Fig. 1.

Analysis of a conceptual model

In this section, a conceptual coupled model similar, but not identical to [16] is employed. An important distinction between the present paper and [16] is as follows. In [16], a build up of carbon dioxide in the atmosphere in the uncoupled simulation is prescribed. In contrast, in the present paper, the model is forced by idealised emission scenarios. This results in different time scales used in these two papers. Here, t p means time scale of emission growth while in the earlier paper the time scale for q CO 2 u c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaaaa@3655@ growth was used. The latter approach implicitly includes all details of emissions and direct fertilisation via prescribed curve q CO 2 u c p l ( t ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaOGaeiikaGIaemiDaqNaeiykaKcaaa@3982@ . However, the time scale associated with q CO 2 u c p l ( t ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaOGaeiikaGIaemiDaqNaeiykaKcaaa@3982@ is larger than t p by an order of magnitude. As a result, to be consistent with the above analysis of the IAP RAS CM, an emission time scale is used here. The equations for the airborne carbon storage for the coupled and uncoupled cases are, respectively,

c 0 d q CO 2 c p l / d t = E − F o c c p l − F l c p l , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4yam2aaSbaaSqaaiabicdaWaqabaGccqWGKbazcqWGXbqCdaqhaaWcbaGaee4qamKaee4ta80aaSbaaWqaaiabikdaYaqabaaaleaacqWGJbWycqWGWbaCcqWGSbaBaaGccqGGVaWlcqWGKbazcqWG0baDcqGH9aqpcqWGfbqrcqGHsislcqWGgbGrdaqhaaWcbaGaem4Ba8Maem4yamgabaGaem4yamMaemiCaaNaemiBaWgaaOGaeyOeI0IaemOray0aa0baaSqaaiabdYgaSbqaaiabdogaJjabdchaWjabdYgaSbaakiabcYcaSaaa@505A@
(3)
c 0 d q CO 2 u c p l / d t = E − F o c u c p l − F l u c p l , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4yam2aaSbaaSqaaiabicdaWaqabaGccqWGKbazcqWGXbqCdaqhaaWcbaGaee4qamKaee4ta80aaSbaaWqaaiabikdaYaqabaaaleaacqWG1bqDcqWGJbWycqWGWbaCcqWGSbaBaaGccqGGVaWlcqWGKbazcqWG0baDcqGH9aqpcqWGfbqrcqGHsislcqWGgbGrdaqhaaWcbaGaem4Ba8Maem4yamgabaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaOGaeyOeI0IaemOray0aa0baaSqaaiabdYgaSbqaaiabdwha1jabdogaJjabdchaWjabdYgaSbaakiabcYcaSaaa@54B3@
(4)

where c0 = 2.1 PgC/ppmv, E stands for emissions, and F l and F oc are terrestrial and oceanic carbon uptakes, respectively. Linearising carbon fluxes [5] and neglecting difference between the changes of globally averaged annual mean surface air and sea surface temperatures one obtains

F X c p l = β X d q CO 2 c p l / d t + γ X d T g c p l / d t , F X u c p l = β X d q CO 2 u c p l / d t , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiqaaaqaaiabdAeagnaaDaaaleaacqWGybawaeaacqWGJbWycqWGWbaCcqWGSbaBaaGccqGH9aqpcqaHYoGydaWgaaWcbaGaemiwaGfabeaakiabdsgaKjabdghaXnaaDaaaleaacqqGdbWqcqqGpbWtdaWgaaadbaGaeGOmaidabeaaaSqaaiabdogaJjabdchaWjabdYgaSbaakiabc+caViabdsgaKjabdsha0jabgUcaRiabeo7aNnaaBaaaleaacqWGybawaeqaaOGaemizaqMaemivaq1aa0baaSqaaiabdEgaNbqaaiabdogaJjabdchaWjabdYgaSbaakiabc+caViabdsgaKjabdsha0jabcYcaSaqaaiabdAeagnaaDaaaleaacqWGybawaeaacqWG1bqDcqWGJbWycqWGWbaCcqWGSbaBaaGccqGH9aqpcqaHYoGydaWgaaWcbaGaemiwaGfabeaakiabdsgaKjabdghaXnaaDaaaleaacqqGdbWqcqqGpbWtdaWgaaadbaGaeGOmaidabeaaaSqaaiabdwha1jabdogaJjabdchaWjabdYgaSbaakiabc+caViabdsgaKjabdsha0jabcYcaSaaaaaa@7226@
(5)

where X = l, oc, T g stands for globally averaged annual mean temperature, and γ X and β X were defined in Introduction. From (3)–(5) one derives

Δ q CO 2 c p l = Δ q CO 2 u c p l + κ Δ T g c p l , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaem4yamMaemiCaaNaemiBaWgaaOGaeyypa0JaeuiLdqKaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaOGaey4kaSIaeqOUdSMaeuiLdqKaemivaq1aa0baaSqaaiabdEgaNbqaaiabdogaJjabdchaWjabdYgaSbaakiabcYcaSaaa@4F48@
(6)

with κ = -G = (γ l + γ oc )/(c0 + β l + β oc ), see Eq. (2). For the C4MIP models under the emission scenario SRES A2, from [14] one derives κ = (8 – 48) ppmv/K for the 21st century. Coefficient β = β l + β oc varies from 1.7 PgC/ppmv to 3.7 PgC/ppmv with a mean value 2.5 PgC/ppmv [14]. For the the IAP RAS CM ensemble simulation, κ depends on t p . Here, at a time short before t m , κ varies from very small values for t p = 25 yr up to 25 ppmv/K (35 ppmv/K) for the version FocNL (FocL) if t p = 250 yr. Coefficient β depends on t p as well. At time of a few t p , it again attains very small values for t p = 25 yr and amounts 1.6 PgC/ppmv (2.5 PgC/ppmv) for the version FocL (FocNL) for t p = 250 yr.

For the same exponentially growing emissions as those used to force the IAP RAS CM in the previos Subsection, and for the uncoupled case, Eq. 4 may be integrated analytically to give

q CO 2 u c p l = q CO 2 , 0 + E 0 t p c 0 + β ( exp ( t / t p ) − 1 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaOGaeyypa0JaemyCae3aaSbaaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaSGaeiilaWIaeGimaadabeaakiabgUcaRKqbaoaalaaabaGaemyrau0aaSbaaeaacqaIWaamaeqaaiabdsha0naaBaaabaGaemiCaahabeaaaeaacqWGJbWydaWgaaqaaiabicdaWaqabaGaey4kaSIaeqOSdigaaOGaeiikaGIagiyzauMaeiiEaGNaeiiCaaNaeiikaGIaemiDaqNaei4la8IaemiDaq3aaSbaaSqaaiabdchaWbqabaGccqGGPaqkcqGHsislcqaIXaqmcqGGPaqkaaa@58E2@
(7)

with the initial atmospheric carbon dioxide concentration C0.

For coupled simulation, one may write [22]

F H c p l = R c p l − λ Δ T g c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aa0baaSqaaiabdIeaibqaaiabdogaJjabdchaWjabdYgaSbaakiabg2da9iabdkfasnaaCaaaleqabaGaem4yamMaemiCaaNaemiBaWgaaOGaeyOeI0Iaeq4UdWMaeuiLdqKaemivaq1aa0baaSqaaiabdEgaNbqaaiabdogaJjabdchaWjabdYgaSbaaaaa@43F9@
(8)

where F H c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aa0baaSqaaiabdIeaibqaaiabdogaJjabdchaWjabdYgaSbaaaaa@3249@ is oceanic heat uptake, Rcplis radiative forcing. The latter is substituted as R c p l = R 0 ln ( q CO 2 c p l / q CO 2 , 0 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOuai1aaWbaaSqabeaacqWGJbWycqWGWbaCcqWGSbaBaaGccqGH9aqpcqWGsbGudaWgaaWcbaGaeGimaadabeaakiGbcYgaSjabc6gaUnaabmaabaGaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaem4yamMaemiCaaNaemiBaWgaaOGaei4la8IaemyCae3aaSbaaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaSGaeiilaWIaeGimaadabeaaaOGaayjkaiaawMcaaaaa@49C0@ with R0 = 5.4 Wm-2K-1 [21, 23]. Climate feedback parameter is λ = R 0 ln 2 / Δ T 2 × CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeq4UdWMaeyypa0JaemOuai1aaSbaaSqaaiabicdaWaqabaGccyGGSbaBcqGGUbGBcqaIYaGmcqGGVaWlcqqHuoarcqWGubavdaWgaaWcbaGaeGOmaiJaey41aqRaee4qamKaee4ta80aaSbaaWqaaiabikdaYaqabaaaleqaaaaa@3EA6@ with Δ T 2 × CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemivaq1aaSbaaSqaaiabikdaYiabgEna0kabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@34FD@ standing for the equilibrium model sensitivity to doubling of the atmospheric carbon dioxide concentartion. Upon substituting F H c p l = C d Δ T g c p l d t MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aa0baaSqaaiabdIeaibqaaiabdogaJjabdchaWjabdYgaSbaakiabg2da9iabdoeadjabdsgaKjabfs5aejabdsfaunaaDaaaleaacqWGNbWzaeaacqWGJbWycqWGWbaCcqWGSbaBaaGccqWGKbazcqWG0baDaaa@40B9@ (C is heat capacity per unit area) one gets

C = d Δ T g c p l d t = R 0 ln ( 1 + Δ q C O 2 c p l q C O 2 , 0 ) − λ Δ T g c p l , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4qamKaeyypa0tcfa4aaSaaaeaacqWGKbazcqqHuoarcqWGubavdaqhaaqaaiabdEgaNbqaaiabdogaJjabdchaWjabdYgaSbaaaeaacqWGKbazcqWG0baDaaGccqGH9aqpcqWGsbGudaWgaaWcbaGaeGimaadabeaakiGbcYgaSjabc6gaUnaabmaabaGaeGymaeJaey4kaSscfa4aaSaaaeaacqqHuoarcqWGXbqCdaqhaaqaaiabdoeadjabd+eapnaaBaaabaGaeGOmaidabeaaaeaacqWGJbWycqWGWbaCcqWGSbaBaaaabaGaemyCae3aaSbaaeaacqWGdbWqcqWGpbWtdaWgaaqaaiabikdaYaqabaGaeiilaWIaeGimaadabeaaaaaakiaawIcacaGLPaaacqGHsislcqaH7oaBcqqHuoarcqWGubavdaqhaaWcbaGaem4zaCgabaGaem4yamMaemiCaaNaemiBaWgaaOGaeiilaWcaaa@61F3@
(9)

Alternatively, Eq. (9) may be treated as a simple zero-dimensional climate model [24].

Eqs. (6), (7) and (9) have been numerically integrated subject to the above-mentioned exponential emissions (expressed via (7)) and to initial condition Δ T g c p l | t = 0 = 0 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaaqGaaeaacqqHuoarcqWGubavdaqhaaWcbaGaem4zaCgabaGaem4yamMaemiCaaNaemiBaWgaaaGccaGLiWoadaWgaaWcbaGaemiDaqNaeyypa0JaeGimaadabeaakiabg2da9iabicdaWaaa@3B38@ . An ensemble of conceptual model integrations is performed varying t p in the same range as it was for IAP RAS CM and varying β in the range from 0 to 5.0 PgC/ppmv (this range is wider than the corresponding C4MIP range, see above). Equilibrium model sensitivity to doubling of the carbon dioxide in the atmosphere was varied between 1 K and 9 K roughly corresponding to combined range from [25–28] (in [29], a narrower range from 2.0 K to 4.5 K is figured). Parameter κ was varied from 8 ppmv/K to 48 ppmv/K (see above).

In all these integrations, f attains maximum and then falls down to unity (Fig. 4a). The timing of this maximum is t m = (2 – 8) × t p depending on governing parameters. The latter range is similar to that obtained in the simulations with the IAP RAS CM despite of the fixed value of Δ T 2 × CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemivaq1aaSbaaSqaaiabikdaYiabgEna0kabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@34FD@ of smaller range of variations in β in the IAP RAS CM. In particular, the latter parameter was varied only between two values, one corresponding to the version FocL and the other corresponding to the version FocNL. Moreover, if a narrower range of Δ T 2 × CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemivaq1aaSbaaSqaaiabikdaYiabgEna0kabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@34FD@ figured in [29] is considered and/or β is constrained to be in a C4MIP range, an interval for t m /t p shrinks only slightly in the integrations with a conceptual model. The obtained t m /t p is insensitive to variations in κ This is consistent with the results obtained in [16]. Typically, t m /t p increases if any of Δ T 2 × CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemivaq1aaSbaaSqaaiabikdaYiabgEna0kabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@34FD@ and β increases or t p decreases (see Fig. 4a). For variations in t p , this behaviour is consistent with that obtained employing the IAP RAS CM. For variations in β, it is consistent as well since combined β for the version FocL is larger than for the version FocNL (basically, due to larger β oc for the former model version, see [15]).

Figure 4
figure 4

Parameter of climate-carbon cycle interaction in the conceptual model (a) and transient temperature sensitivity α computed for 100-year running segments (b) as functions of t / t p . Shown are curves for t p = 50 and 150 yr (thin and thick lines respectively) with β = 1.5 and 2.5 PgC/ppmv (dashed and solid lines correspondingly) and Δ T 2 × CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemivaq1aaSbaaSqaaiabikdaYiabgEna0kabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@34FD@ = 2 and 4 K (black and gray lines respectively) for κ = 28 ppmv/K.

In a conceptual model, by construction, G = -κ does not change during the course of integration. As a result, as it was for the IAP RAS CM, eventual tending of f to unity is due to decrease in transient temperature sensitivity a during the course of integration (Fig. 4b). In turn, this is due to weak, logarithmic dependence of the CO2 radiative forcing on the atmospheric concentration of carbon dioxide. One may estimate q CO 2 , s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aaSbaaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaSGaeiilaWIaem4Camhabeaaaaa@3317@ from this conceptual model as well. For simplicity, one may assume that temperature response to carbon dioxide forcing is stationary. This is achieved by dropping F H c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aa0baaSqaaiabdIeaibqaaiabdogaJjabdchaWjabdYgaSbaaaaa@3249@ in Eq. (8). Taking into accout Eq. (6), one obtains

κ Δ T 2 × CO 2 q CO 2 , s ln 2 ln q CO 2 , s q CO 2 , 0 = 0.05 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqaH6oWAcqqHuoarcqWGubavdaWgaaqaaiabikdaYiabgEna0kabboeadjabb+eapnaaBaaabaGaeGOmaidabeaaaeqaaaqaaiabdghaXnaaBaaabaGaee4qamKaee4ta80aaSbaaeaacqaIYaGmaeqaaiabcYcaSiabdohaZbqabaGagiiBaWMaeiOBa4MaeGOmaidaaOGagiiBaWMaeiOBa4wcfa4aaSaaaeaacqWGXbqCdaWgaaqaaiabboeadjabb+eapnaaBaaabaGaeGOmaidabeaacqGGSaalcqWGZbWCaeqaaaqaaiabdghaXnaaBaaabaGaee4qamKaee4ta80aaSbaaeaacqaIYaGmaeqaaiabcYcaSiabicdaWaqabaaaaOGaeyypa0JaeGimaaJaeiOla4IaeGimaaJaeGynaudaaa@5846@
(10)

where again the threshold value for f is set to 1.05. This equation has been solved numerically appling the dichotomy method. In general, the larger κ and/or Δ T 2 × CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemivaq1aaSbaaSqaaiabikdaYiabgEna0kabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@34FD@ the larger q CO 2 , s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aaSbaaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaSGaeiilaWIaem4Camhabeaaaaa@3317@ . For Δ T 2 × CO 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemivaq1aaSbaaSqaaiabikdaYiabgEna0kabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbeaaaaa@34FD@ = 2.2 K, the values obtained for the conceptual model are similar, but slightly larger, to those estimated from the IAP RAS CM runs.

As an aside issue, one may quantify how "small" E0 should be in order not to make the approach followed here meaningless. Physically, E0 has not to result in sudden kink in q CO 2 u c p l MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyCae3aa0baaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaaWcbaGaemyDauNaem4yamMaemiCaaNaemiBaWgaaaaa@3655@ at t = 0. From Eq. (7), it needs

E 0 ≪ ( c 0 + β ) q CO 2 , 0 / t p ≤ c 0 q CO 2 , 0 / t p ≤ 2.4 P g C / y r . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyrau0aaSbaaSqaaiabicdaWaqabaGccqWIQjspcqGGOaakcqWGJbWydaWgaaWcbaGaeGimaadabeaakiabgUcaRiabek7aIjabcMcaPiabdghaXnaaBaaaleaacqqGdbWqcqqGpbWtdaWgaaadbaGaeGOmaidabeaaliabcYcaSiabicdaWaqabaGccqGGVaWlcqWG0baDdaWgaaWcbaGaemiCaahabeaakiabgsMiJkabdogaJnaaBaaaleaacqaIWaamaeqaaOGaemyCae3aaSbaaSqaaiabboeadjabb+eapnaaBaaameaacqaIYaGmaeqaaSGaeiilaWIaeGimaadabeaakiabc+caViabdsha0naaBaaaleaacqWGWbaCaeqaaOGaeyizImQaeGOmaiJaeiOla4IaeGinaqJaeeiiaaIaemiuaaLaem4zaCMaem4qamKaei4la8IaemyEaKNaemOCaiNaeiOla4caaa@5D8F@
(11)

The value of E0 selected in this paper, 0.1 PgC/yr, fulfills this condition. If (11) is satisfied then qualitative conclusion of the present paper and the respective conclusions related to the order of magnitudes of basic characteristics are expected to remain unchanged.

Discussion

An eventual transient saturation, as studied in the present framework, implies that initially (at a time ≪ t p ) and eventually (at a time ≫ t p ), the state of the coupled system is close to the state of the uncoupled one. However, during an intermediate period when t' ≃ t p , f deviates substantinally from unity and climate-carbon cycle interaction does matter. Taking into account typical time scales of CO2 emission growth for the SRES emission, t p = 50 – 200 yr, this intermediate period extends for the next several centuries.

While it is beyond the scope of the present paper to concern CO2 emissions changing non-exponentially in time, it is possible to make a note with this respect. In particular, if the future emissions would change in time slower than the exponential ones, it would be expressed in apparent increase of t p during the course of simulation. In turn, the timing t m of maximum of f and an eventual transient saturation would be delayed in comparison to the case with constant t p . On the other hand, for faster-than-exponential emissions this eventual transient saturation would come closer to the present day.

The long-term fate of anthropogenic carbon is basically governed by calcium carbonate dissolution in the ocean [30, 31]. This effect is taken into account only in the version FocNL and neglected in the version FocL and in the conceptual model. However, carbonate dissolution is important only for time scales larger than ≈ 5, 000 yr [30] which is well above the length of every individual simulation performed in the present paper (1, 500 yr). As a result, the version FocL and the conceptual climate-carbon cycle model are applicable for the problem considered here as well.

Returning to C4MIP simulations, one notes that it is possible to find some hints for this eventual transient saturation in these simulations [14] as well (see Background). However, in the C4MIP integrations, an eventual climate-carbon cycle carbon saturation is less marked. The most probable reason for this is the length of the C4MIP simulations which is too small to make this saturation visible given an emission intensity in these simulations. One expects that if these simulations would extended to the future, eventually, parameter f would converge to unity in the C4MIP ensemble as well. Some complications for this convergence could arise from a specific behaviour of a for particular models, e.g., due to behaviour of the oceanic heat uptake in the course of integration and its interrelation with equilibrium temperature sensitivity [32–35].

As a final note, one may distinguish manifestation of a climate-carbon cycle feedback eventual transient saturation in terms of feedback parameter, on one hand, and in terms of feedback gain, on the other one. In particular, as dg/dt = f-2df/dt and, generally, f > 1 in the course of simulations, this saturation is somewhat masked for g in comparison to f. This also adds to the masking of the eventual transient saturation of the climate-carbon cycle feedback in the C4MIP simulations.

Conclusion

In this paper, simulations with the IAP RAS climate model of intermediate complexity have been performed to study temporal variations of the climate-carbon cycle feedback parameter. Two model versions were considered differing between each other by the formulation of the oceanic uptake of carbon dioxide and by the governing parameters of the terrestrial carbon uptake module. Both versions were forced by idealised scenarios of fossil fuel emissions. In these scenarios, emissions grow exponentially in time (with a characteristic timescale t p = 25 – 250 yr) starting from a small initial emission value. Land use emissions were set to zero. In all simulations with both model versions, climate-carbon cycle feedback parameter f grows initially, attains maximum, and then decreases eventually tending to unity. The timing of this maximum is of the order of a few t p for the selected small value of initial emission. This general behaviour is consistent to that obtained in a earlier simulations with the same model.

In particular, coefficients of climate-carbon cycle interaction diagnosed for 100-yr running segments lead to increase of climate-carbon cycle gain g during the course of integration. This increase, however, is overcompensated by decrease of transient climate sensitivity α. The latter leads to the above-mentioned eventual transient saturation of climate-carbon cycle feedback.

The IAP RAS CM simulations are supported by an analysis of a conceptual model with linear dependence of carbon sinks from the atmosphere on changes in atmospheric carbon dioxide concentration and global temperature rise. This simple model exhibits an eventual transient saturation of the climate-carbon cycle feedback which is similar to the model of intermediate complexity.

References

  1. Cox P, Betts R, Jones C, Spall S, Totterdell I: Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model. Nature 2000, 408: 184–187. 10.1038/35041539

    Article  CAS  Google Scholar 

  2. Friedlingstein P, Bopp L, Ciais P, Dufresne JL, Fairhead L, Le Treut H, Monfray P, Orr J: Positive feedback between future climate change and the carbon cycle. Geophys Res Lett 2001,28(8):1543–1546. 10.1029/2000GL012015

    Article  CAS  Google Scholar 

  3. Brovkin V, Bendtsen J, Claussen M, Ganopolski A, Kubatzki C, Petoukhov V, Andreev A: Carbon cycle, vegetation, and climate dynamics in the Holocene: Experiments with the CLIMBER-2 model. Glob Biogeochem Cycles 2002,16(4):1139. 10.1029/2001GB001662

    Article  Google Scholar 

  4. Brovkin V, Sitch S, von Bloh W, Claussen M, Bauer E, Cramer W: Role of land cover changes for atmospheric CO 2 increase and climate change during the last 150 years. Glob Change Biol 2004, 10: 1253–1266. 10.1111/j.1365-2486.2004.00812.x

    Article  Google Scholar 

  5. Friedlingstein P, Dufresne JL, Cox P, Rayner P: How positive is the feedback between climate change and the carbon cycle? Tellus 2003,55B(2):692–700.

    Article  CAS  Google Scholar 

  6. Jones C, Cox P, Essery R, Roberts D, Woodage M: Strong carbon cycle feedbacks in a climate model with interactive CO 2 and sulphate aerosols. Geophys Res Lett 2003,30(9):1479. 10.1029/2003GL016867

    Article  Google Scholar 

  7. Matthews H, Weaver A, Meissner K, Gillett N, Eby M: Natural and anthropogenic climate change: incorporating historical land cover change, vegetation dynamics and the global carbon cycle. Clim Dyn 2004, 22: 461–479. 10.1007/s00382-004-0392-2

    Article  Google Scholar 

  8. Matthews H, Weaver A, Meissner K: Terrestrial carbon cycle dynamics under recent and future climate change. J Climate 2005,18(10):1609–1628. 10.1175/JCLI3359.1

    Article  Google Scholar 

  9. Govindasamy B, Thompson S, Mirin A, Wickett M, Caldeira K, Delire C: Increase of carbon cycle feedback with climate sensitivity: results from a coupled climate and carbon cycle model. Tellus 2005,57B(2):153–163.

    Article  CAS  Google Scholar 

  10. Mokhov I, Eliseev A, Karpenko A: Sensitivity of the IAP RAS Global Climatic Model with an interactive carbon cycle to anthropogenic influence. Doklady Earth Sci 2006,407(3):424–428. 10.1134/S1028334X06030172

    Article  Google Scholar 

  11. Eliseev A, Mokhov I, Karpenko A: Climate and carbon cycle variations in the 20th and 21st centuries in a model ofintermediate complexity. Izvestiya, Atmos Ocean Phys 2007, 43: 1–14. 10.1134/S000143380701001X

    Article  Google Scholar 

  12. Eliseev A, Mokhov I: Carbon cycle-climate feedback sensitivity to parameter changes of a zero-dimensional terrestrial carbon cycle scheme in a climate model of intermediate complexity. Theor Appl Climatol 2007,89(1–2):9–24. 10.1007/s00704-006-0260-6

    Article  Google Scholar 

  13. Bala G, Caldeira K, Mirin A, Wickett M, Delire C: Multicentury changes to the global climate and carbon cycle: Results from a coupled climate and carbon cycle model. J Climate 2005,18(21):4531–4544. 10.1175/JCLI3542.1

    Article  Google Scholar 

  14. Friedlingstein P, Cox P, Betts R, Bopp L, von Bloh W, Brovkin V, Doney S, Eby M, Fung I, Govindasamy B, John J, Jones C, Joos F, Kato T, Kawamiya M, Knorr W, Lindsay K, Matthews H, Raddatz T, Rayner P, Reick C, Roeckner E, Schnitzler KG, Schnur R, Strassmann K, Weaver A, Yoshikawa C, Zeng N: Climate-carbon cycle feedback analysis: Results from the C4MIP model intercomparison. J Climate 2006,19(22):3337–3353. 10.1175/JCLI3800.1

    Article  Google Scholar 

  15. Mokhov I, Eliseev A, Karpenko A: Decadal-to-centennial scale climate-carbon cycle interactions from global climate models simulations forced by anthropogenic emissions. In Climate Change Research Progress. Edited by: Peretz L. Hauppauge, NY: Nova Publishers; 2008:in press.

    Google Scholar 

  16. Eliseev A, Mokhov I: Eventual saturation of the climate-carbon cycle feedback studied with a conceptual model. Ecol Mod 2008, 213: 127–132. 10.1016/j.ecolmodel.2007.11.015

    Article  Google Scholar 

  17. Mokhov I, Eliseev A, Demchenko P, Khon V, Akperov M, Arzhanov M, Karpenko A, Tikhonov V, Chernokulsky A, Sigaeva E: Climate changes and their assessment based on the IAP RAS global model simulations. Doklady Earth Sci 2005,402(4):591–595.

    CAS  Google Scholar 

  18. Bacastow R: Numerical evaluation of the evasion factor. In Carbon Cycle Modelling, SCOPE-16. Edited by: Bolin B. N.Y.: J. Wiley and Sons; 1981:95–101.

    Google Scholar 

  19. Millero F: Thermodynamics of carbon dioxide system in the ocean. Geophys Cosmophys Acta 1995,59(4):661–677.

    CAS  Google Scholar 

  20. Marland G, Boden T, Andres R: Global, regional, and national CO2 emissions. In Trends: A Compendium of Data on Global Change. Oak Ridge, Tenn.: Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, U.S. Department of Energy; 2005.

    Google Scholar 

  21. Houghton J, Ding Y, Griggs D, Noguer M, Linden P, Dai X, Maskell K, Johnson C, (Eds): Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge/New York: Cambridge University Press; 2001.

    Google Scholar 

  22. Gregory J, Stouffer R, Raper S, Stott P, Rayner N: An observationally based estimate of the climate sensitivity. J Climate 2002,15(22):3117–3121. Publisher Full Text 10.1175/1520-0442(2002)015<3117:AOBEOT>2.0.CO;2

    Article  Google Scholar 

  23. Myhre G, Highwood E, Shine K, Stordal F: New estimates of radiative forcing due to well mixed greenhouse gases. Geophys Res Lett 1998,25(14):2715–2718. 10.1029/98GL01908

    Article  CAS  Google Scholar 

  24. Andreae M, Jones C, Cox P: Strong present-day aerosol cooling implies a hot future. Nature 2005, 435: 1187–1190. 10.1038/nature03671

    Article  CAS  Google Scholar 

  25. Harvey L, Kaufmann R: Simultaneously constraining climate sensitivity and aerosol radiative forcing. J Climate 2002,15(20):2837–2861. Publisher Full Text 10.1175/1520-0442(2002)015<2837:SCCSAA>2.0.CO;2

    Article  Google Scholar 

  26. Stainforth D, Aina T, Christensen C, Collins M, Faull N, Frame D, Kettleborough J, Knight S, Martin A, Murphy J, Piani C, Sexton D, Smith L, Spicer R, Thorpe A, Allen M: Uncertainty in predictions of the climate response to rising levels of greenhouse gases. Nature 2005, 433: 403–406. 10.1038/nature03301

    Article  CAS  Google Scholar 

  27. Forest C, Stone P, Sokolov A: Estimated PDFs of climate system properties including natural and anthropogenic forcings. Geophys Res Lett 2006, 33: L01705. 10.1029/2005GL023977

    Article  Google Scholar 

  28. Hegerl G, Crowley T, Hyde W, DJ F: Climate sensitivity constrained by temperature reconstructions over the past seven centuries. Nature 2006, 440: 1029–1032. 10.1038/nature04679

    Article  CAS  Google Scholar 

  29. Solomon S, Qin D, Manning M, Marquis M, Averyt K, Tignor M, LeRoy Miller H, Chen Z, (Eds): Climate Change 2007: The Physical Science Basis. Cambridge/New York: Cambridge University Press; 2007.

    Google Scholar 

  30. Archer D, Martin P, Buffett B, Brovkin V, Rahmstorf S, Ganopolski A: The importance of ocean temperature to global biogeochemistry. Earth Planet Sci Lett 2004,222(2):333–348. 10.1016/j.epsl.2004.03.011

    Article  CAS  Google Scholar 

  31. Scheffer M, Brovkin V, Cox P: Positive feedback between global warming and atmospheric CO 2 concentration inferred from past climate change. Geophys Res Lett 2006,33(10):L10702. 10.1029/2005GL025044

    Article  Google Scholar 

  32. Hansen J, Russell G, Lacis A, Fung I, Rind D, Stone P: Climate response times: dependence on climate sensitivity and ocean mixing. Science 1985,229(4716):857–859. 10.1126/science.229.4716.857

    Article  CAS  Google Scholar 

  33. Harvey L: Effect of ocean mixing on the transient climate response to a CO 2 increase: analysis of recent model results. J Geophys Res 1986, 90: 2709–2718. 10.1029/JD091iD02p02709

    Article  Google Scholar 

  34. Schlesinger M, Jiang X: Simple model representation of atmosphere-ocean GCMs and estimation of the time scale of CO 2 -induced climate change. J Climate 1990,3(12):1297–1315. Publisher Full Text 10.1175/1520-0442(1990)003<1297:SMROAO>2.0.CO;2

    Article  Google Scholar 

  35. Raper S, Gregory J, Stouffer R: The role of climate sensitivity and ocean heat uptake on AOGCM transient temperature and thermal expansion response. J Climate 2002, 15: 124–130. Publisher Full Text 10.1175/1520-0442(2002)015<0124:TROCSA>2.0.CO;2

    Article  Google Scholar 

Download references

Acknowledgements

The authors are indebted to L.L. Golubyatnikov and A.M. Tarko for useful discussions on the topic of the manuscript. Comments by anonymous reviewers has led to improved presentation of the results. This work has been supported by the Russian Foundation for Basic Research (grants 05-05-64907, 05-05-65167, and 07-05-00273), the President of Russia grant 4166.2005.5, and by the Programs of the Russian Ministry for Science and Education, Russian Federal Agency for Science and Innovations and the Russian Academy of Sciences (contracts 02.515.11.5031 and 02.515.11.5046).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alexey V Eliseev.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

Both authors equally contributed to the paper.

Authors’ original submitted files for images

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Mokhov, I.I., Eliseev, A.V. Explaining the eventual transient saturation of climate-carbon cycle feedback. Carbon Balance Manage 3, 4 (2008). https://doi.org/10.1186/1750-0680-3-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1750-0680-3-4

Keywords