Roles of climate feedback and ocean vertical mixing in modulating global warming rate

Despite the rapid increase of greenhouse gases (GHGs) in the atmosphere during the past 50 years, observed global mean surface temperature (GMST) showed a pause in the warming trend during the first decade of the twenty-first century. This is referred to as the global warming “hiatus”. A dominant hypothesis emphasizes that the superimposition of the cold phase of the Pacific decadal variability and the global warming trend can lead to the hiatus. Using simply energy balance models, we explore two potential mechanisms that may supress the GMST warming trend: enhanced negative climate feedback and downward heat mixing. Forced by linearly increasing heating, a stronger negative climate feedback can reduce the GMST warming rate, but cannot result in a warming hiatus. Downward mixing of heat can cause a short-lived hiatus of surface warming rate due to enhanced nonlinear ocean heat uptake by the lower ocean, but the surface warming would be accelerated in the long run due to the decline of downward heat mixing rate. This study provides further evidence, both theoretically and numerically, that in the long run, the only route to contain the global warming effectively is to reduce GHG emissions.


Introduction
The increase of greenhouse gases (GHGs) in the atmosphere is commonly thought as the deterministic cause for the centennial warming trend of global mean surface temperature (GMST) since the Industrial Revolution (IPCC 2013). With the steep increase of carbon dioxide (CO 2 ) in the atmosphere since 1960, the GMST, however, showed a warming pause for about 10 years during the first decade of the twentieth-one century ( Fig. 1) (Easterling and Wehner 2009;Knight et al. 2009;Trenberth and Fasullo 2013). This phenomenon is termed as the global warming "hiatus" (IPCC 2013; Kosaka and Xie 2013;Yan et al. 2016). The observed CO 2 concentration at Mauna Loa of Hawaii has exceeded 400 ppm since 2015 ( Fig. 1), with a growth rate of annual mean CO 2 at Mauna Loa being about 2 ppm/year since 2000, well above any preceding period. This inconsistency between the warming hiatus and the intensification of anthropogenic forcing has set off fierce scientific and political debates on the authenticity of global warming.
A global warming hiatus can appear when the warming trend encounters natural climate variability. This mechanism was first proposed by Kosaka and Xie (2013), and recognized by many scientists (Clement and DiNezio 2014;England et al. 2014;Dai et al. 2015). In Kosaka and Xie (2013), the global warming hiatus was well simulated after nudging the model sea-surface temperature (SST) to the observed SST in the central to eastern tropical Pacific in a historical run. Xie et al. (2016) further indicated that the temporal variation of SST can be driven by radiative forcing and periodic natural variability. When the natural variability happens to be in a downward phase, the warming hiatus would occur. The natural variability is widely believed to be the Pacific decadal oscillation (PDO) (Clement and DiNezio 2014;Dai et al. 2015). The recent hiatus appeared because the historical SST happened to be in the cold phase of the PDO during the first decade of the twentieth-one century (Meehl et al. 2013;Trenberth and Fasullo 2013;Brown et al. 2015). This kind of hiatus also occurred during1943-1976 ( Fig. 1) (Trenberth and Fasullo 2013;England et al. 2014), when the PDO was in its cold phase. Although about 30% of heat uptake increase occurs in the eastern tropical Pacific associated with the cold phase of the PDO, studies also suggested that the Pacific was not the only contributor to the hiatus, and other ocean basins, in particular, the southern and subtropical Indian Ocean and subpolar North Atlantic, contributed substantially to the global warming hiatus via reduced heat loss to the atmosphere (Drijfhout et al. 2014).
A more intuitive explanation to the recent GMST hiatus comes from the view of Earth's energy budget: either a reduction in radiative forcing from the top of the atmosphere (TOA) or an enhanced heat uptake by the deeper ocean ) could lead to a hiatus. Weakened solar radiation during a period of low sunspot activity (Trenberth 2009), or due to reduction in stratospheric water vapor (Solomon et al. 2010), and increasing aerosols related to anthropogenic activities (Smith et al. 2016;Kaufmann et al. 2011) and 17 small volcano eruptions since 1999 (Santer et al. 2014;Smith et al. 2016) could have contributed to the recent hiatus by weakening the net TOA incoming energy. Trenberth and Fasullo (2013) showed that forcing reductions from these factors likely contributed no more than 20% of the global warming slowdown. Schmidt et al. (2014) showed that recent moderate volcano eruptions and anthropogenic pollutions accounted for half of the divergence between observation and model simulations, while the solar minimum explained one seventh. Kühn et al. (2014) found that Asian pollution from coal burning has very little effect globally or regionally, because of the cancelation between the warming effect of black carbon and the cooling effect of sulphate aerosols.
The net TOA energy gain during the recent hiatus is estimated to be 0.5 ± 0.43 W∕m 2 (Loeb et al. 2012), showing no apparent reduction. The deep ocean is thus emerged as a likely candidate to absorb the additional energy, buffering the global warming and leading to the GMST hiatus (Meehl et al. 2011;Balmaseda et al. 2013a, b;Rhein et al. 2013;Liu et al. 2016). Both observations and model simulations showed that the vertical heat redistribution, especially in the North Atlantic and Southern Ocean, is remarkably different between hiatus period and fast-warming period (Meehl et al. 2011;Chen and Tung 2014;Drijfhout et al. 2014;Liu et al. 2016). The subsurface warming in the Southern Ocean is attributed to the southward displacement and intensification of the circumpolar jet, while the North Atlantic's role is believed to be related to the Atlantic meridional overturning circulation (AMOC). Meehl et al. (2011) simulated a negative trend of the AMOC in a historical run, that is, a weakening of deep convection in the North Atlantic, which can lead to subsurface warming. In contrast, Chen and Tung (2018) showed an acceleration of the AMOC during 1993-2005 based on various independent proxies. They believed that the enhanced deep convection transports more heat downward to the deep ocean, which neutralizes the warming effect caused by human activities and slows down the surface warming.
In this work we explore theoretically whether the Earth's warming trend as a whole can be controlled under a linear increase of GHGs. Using simple energy balance models (EBMs), we examine two candidates: the overall global climate feedback and ocean vertical mixing. From the viewpoint of Earth energy balance, a stronger negative climate feedback would reduce the net energy absorbed by the Earth surface, which might have the potential to control the GMST warming. However, in reality the global climate feedback is getting even positively stronger (Armour et al. 2013;Gregory et al. 2015;Rugenstein et al. 2020), which would potentially result in a disastrous runaway warming. Enhanced ocean vertical mixing can reduce the global warming rate over a short term, by mixing heat downward to the deeper ocean; however, it would eventually accelerate the surface warming in the long run, once the deeper ocean is saturated. A stronger overturning circulation can also transport more surface warm water downward, contributing to the hiatus over a short time. However, the overturning circulation itself depends strongly on the stratification, and can be shut down by a reduced stratification. This study provides further  Santa Maria in 1902, Agung in 1963, El Chichon in 1982, and Pinatubo in 1991 Dashed purple line shows the linear trend of CO 2 concentration. The surface temperature data is from the HadCRUT4 Morice et al. 2012), downloaded at https:// cruda ta. uea. ac. uk/ cru/ data/ tempe rature/. The CO 2 data is from the station at Mauna Loa of Hawaii, downloaded at https:// gml. noaa. gov/ ccgg/ trends/ evidence that in the long run, the only route to contain the global warming effectively is to reduce GHGs, namely, to change the behaviour of human beings. This paper is organized as follows. In Sect. 2, we introduce one-box EBMs, and investigate the role of climate feedback in the hiatus of global warming rate. Two-box model is introduced in Sect. 3, and the role of vertical ocean mixing is examined. Summary and discussion are given in Sect. 4. Detailed derivation of theoretical equations is provided in an Appendix.

One-box model with climate feedback
First, we use a zero-dimensional EBM, i.e., a 1-box ocean model (Fig. 2a), to investigate the impact of climate feedback on surface warming trend. The transient response of the GMST in the 1-box model is proportional to the net energy flux at the TOA (Gregory et al. 2004), that is, where C = c p D is ocean heat capacity, = 1026 kg • m −3 is seawater density, c p = 3900J • (kg • C) −1 is seawater specific heat at constant pressure, and D is ocean depth, set to 400 m here. Symbol Δ denotes the anomaly from the Usually, F should be zero in an equilibrium climate state, and temperature stays unchanged. However, if the GHGs increase, the Earth gains extra energy that causes ΔT to increase. Here, we assume the concentration of GHGs increases linearly with time, and the anomalous radiative forcing can be simply expressed by: where is set to 0.016W • m −2 yr −1 , which is equivalent to the greenhouse effect by a gradual doubling of CO 2 concentration over 200 years.

Constant climate feedback
The situation with constant climate feedback is examined first. Under a linearly increasing radiative forcing given in Eq.
(2), the theoretical solution to Eq. (1) can be easily written as follows, For simplicity, ΔT is written as T from now on. Thus, the warming rate is, When t → ∞ , the ocean warming goes toward infinity, and the warming rate reaches a constant B , that is, Here, we can define the e-folding timescale = C B , to represent the fast transient timescale of the surface warming (Hansen et al. 1981), which is about 34 years here. When t > , the surface warming goes up roughly linearly and the warming rate approaches a constant. Figure 3a, b show the theoretical solutions of T and dT dt , respectively. T increases roughly linearly as expected (Fig. 3a). During the first several decades, the warming rate ( dT dt ) itself increases roughly linearly, and then reaches a constant value in about 100 years (Fig. 3b). The ocean temperature eventually increases linearly after 100 years (Fig. 3a), and the system will never reach an equilibrium. Of course, The ocean surface heat flux is denoted by solid blue arrows. In (b), yellow arrow represents vertical heat transport in the ocean interior stronger negative feedback will lead to a weaker warming rate, and vice versa (color curves, Fig. 3a, b). The temperature change is straightforward under the constant climate feedback. One important implication here is that to bend the warming rate curve downward, i.e., to reduce the warming rate to zero, or in other words, to realize a global warming hiatus in this simple box model, an extremely strong negative feedback is required. In fact, Eq. (5) shows that only when B → ∞, dT dt = B → 0 . This situation is practically impossible.

Varying climate feedback
Here, we explore the possibility of a global warming hiatus by assuming a temperature-dependent climate feedback, under a strong positive . This situation was studied previously (e.g., Winton et al. 2010;Armour et al. 2013;Bloch-Johnson et al. 2015;Gregory et al. 2015).
The theoretical solution to system (1), (2) and (6) is written as follows, is the Airy function (Abramowitz and Stegun 1964). There are four kinds of Airy function.
. Details on the derivation of Eq. (7) and parameters K, K 1 , K 2 , a and b can be found in Appendix.
For the Airy functions, we have, Therefore, The global warming will slow down and the warming rate will reduce substantially if the negative climate feedback increases with the global warming (i.e., > 0 ) (red curves, Fig. 3c, d). However, a global warming hiatus, i.e., the zero warming rate, appears to be unlikely for a reasonable range of the climate feedback. In contrast, the global warming will be easily out of control if the negative feedback is weakened by the global warming (i.e., < 0 ) (blue curves, Fig. 3c, d). This is also called "quadratic runaway" (Bloch-Johnson et al. 2015). To view clearly the temperature change, we examined the situations with = ±0.5, ±1.0 W • m −2 K −2 . These values are extremely strong when compared to more practical values derived from observations and coupled models. For example, based on estimates from various coupled model studies, is suggested to be in the range of ±0.06 W • m −2 K −2 (Roe and Armour 2011). Previous studies suggested that the Planck feedback would become more negative under the global warming, at a rate of about -Johnson et al. 2015). The water vapour feedback would become more positive, possibly offsetting the negative lapse rate feedback (Bony et al. 2006;Soden and Held 2006). The surface albedo feedback would be weakened with the melting of snow and ice under global warming, that is, this negative feedback will, in fact, become weaker, at a rate of about 0.1 W • m −2 K −2 (Manabe and Bryan 1985). In Fig. 3c, d, the black and red curves are plotted based on the theoretical solution (7), while the blue curves are plotted based on numerical solutions to Eq. (A1) (Appendix), since the Airy function does not work very well when < 0.
In fact, Eq. (9) suggests that, when the time goes to infinity, the warming magnitude itself goes to infinity, while the global warming rate could approach zero. To show this situation more clearly, we replot the temperature and its rate of change in the case of B = 1.5, = 1.0 in Fig. 4. We can see that in about 10 5 years, the warming rate can approach zero (it would never be equal to zero based on Eq. (A1)), the global ocean warming can be as strong as 40 °C. Of course, this situation is far beyond the reality, even for the Earth's climate over the past billion years. Therefore, we can conclude safely that it is unlikely to obtain a global warming hiatus by merely enhancing the negative feedback, in the case of linearly increasing external forcing. We can expect a substantial slowdown of the global warming under strong negative feedback, but never the zero warming rate (hiatus). Moreover, even under an extremely strong negative feedback, the mitigation effect of climate feedback works slowly at inter-decadal to centennial timescales (Fig. 3d). Figure 3 suggests that it is unlikely that the climate feedback change alone has contributed seriously to the GMST hiatus observed in the first decade of the twenty-first century century.

Two-box model with vertical heat mixing
As the 1-box EBM does not take into account ocean interior process, a 2-layer ocean EBM is introduced here to study the effect of lower ocean on the surface warming (Fig. 2b). Due to its enormous heat capacity, the interior ocean could play a role in GMST hiatus, by taking up heat from the layer above and redistributing the heat in the vertical (Gregory 2000). The 2-layer box model is widely used in the study of transient climate responses (e.g., Held et al. 2010;Geoffroy et al. 2013;Zhou and Chen 2015;Yoshimori et al. 2016), climate sensitivity (Gregory et al. 2015), and so on. The model is effective in capturing the minimal physics of transient surface temperature change (Gregory et al. 2015).
The 2-layer model consists of upper and lower oceans (Fig. 2b), and the layers have heat capacities C 1 and C 2 , respectively. Here, we set C 2 = 10C 1 , corresponding to a 400-m upper ocean and a 4000-m lower ocean, respectively. The upper ocean participates in the surface energy budget, and transports heat flux H = (T 1 − T 2 ) downward to the deep ocean through processes such as vertical advection, mixing and diffusion. We simply call them "vertical heat transport" or ocean heat uptake (OHU) by the lower ocean in this work. The model equations are, where T 1 and T 2 are the temperature anomalies of upper and lower layers, respectively, from an unperturbed climate; and T v is the change in vertical temperature gradient. represents the efficiency of vertical heat transport (units: W • m −2 K −1 ). It is also treated as the heat uptake efficiency by the deep ocean (Gregory and Mitchell 1997;Gregory et al. 2015). can vary empirically over time (Raper et al. 2002). Watanabe et al. (2013) reported a range of of about 0.5 − 1.5 W • m −2 K −1 based on outputs from 16 CMIP5 models. Gregory et al. (2015)

Constant efficiency of vertical heat transport
The analytical solutions to the 2-box system of (10)-(12) under constant are written as follows, The warming rates for the 2-layer ocean are, Here, we have defined the fast transient timescale f and the slow timescale s , respectively, as follows, Details on the derivation of Eqs. (13,14,15,18), and the parameters therein can be found in Appendix.
Based on (13-15), we have Under constant , the warming rate for the 2-box ocean will eventually approach that of the 1-box ocean (Eq. (5)).
The transient timescales f and s given in Eq. (18) are for the upper and lower oceans, respectively, which are about 25 and 1360 years, respectively, given the parameter = 0.5 , B = 1.5 and the depths of the 2-box ocean of 400 and 4000 m, respectively. Since we consider a linearly increasing external forcing in this paper, it is the warming rate, instead of warming itself, that can reach the equilibrium. Therefore, f and s can be thought as the e-folding timescales of the warming rate for the upper and deep oceans, respectively. This is qualitatively consistent with fast and slow timescales for the temperature changes in a 2-layer ocean, under sudden constant external forcing (Gregory et al. 2015). It is expected that for the lower ocean, the quasi-equilibrium timescale of the warming rate can be as long as several thousand years, under linearly enhancing external forcing. Figure 5a-c show the results of theoretical solutions (13)-(17). The OHU by the lower ocean can slow down the surface warming. This is straightforward, and can be easily deduced from Eqs. (1) and (10). Under the same external forcing and climate feedback, the surface temperature increase in the 2-box model is smaller in magnitude than that in the 1-box model. For example, at the end of the 200 th year, T 1 in the 2-box model is about 1.4 °C (Fig. 5a), 20% smaller than that in the 1-box model (~ 1.8 °C) (Fig. 3a). Here, = 0.5 and B = 1.5 . Under the linearly increasing forcing, in the long run (more than 3000 years) the 2-layer system can approach an "equilibrium" state, in terms of warming rate ( dT dt ) and vertical temperature gradient ( T v ) (Fig. 5b, c). The 2-layer ocean will eventually warm up at the same rate B as indicated in Eq. (19) and shown in Fig. 5b, which is independent of . T v will eventually reach a constant B * C 2 (~ 10.8 °C) (Eq. (19), Fig. 5c), which is inversely proportional to . It is obvious that the efficiency of vertical heat transport affects the magnitude and timescale of temperature change in each layer. It does not affect those of the whole ocean, since it only redistributes the heat in the vertical. A larger when t → ∞ will lead to a quicker and stronger warming in the deep ocean, at the cost of a slower and weaker warming in the upper ocean (Fig. 5a, b), which leads to a weaker vertical temperature gradient (dashed black, Fig. 5c). The theoretical solution (Eq. (19)) shows that the equilibrium warming rates for both the surface and deep oceans depend only on the external heating rate ( ) and the overall climate feedback ( B ) of the system, regardless of the interior parameters of the ocean. The overall warming magnitude and warming rate for the whole ocean (figure not shown) would eventually approach those in the 1-box model, given the same total ocean depth, the climate feedback in the two models and long enough time. Here, we see that under a constant (Fig. 5b), a decrease of the surface warming rate, i.e., a surface ocean hiatus, will never occur.

Varying efficiency of vertical heat transport
In reality, the surface ocean is usually warmed up faster than the deep ocean, leading to an enhanced vertical stratification at least during the transient period of global warming (Yang and Zhang 2008;Yang and Wang 2009). Although there is a lack of observational evidence supporting the proportional relationship between stratification and the efficiency of vertical heat transport, we can study a varying efficiency in a simple conceptual model, to gain insight to its role in surface warming hiatus. Here, we assume a linear relationship between the efficiency of vertical heat transport and vertical temperature gradient as follows, where 0 = 0.5 W • m −2 K −1 is the reference efficiency from an unperturbed state, and is the rate of change, related to T v . A positive denotes a more (less) efficient vertical heat transport in response to an enhanced (weakened) stratification; and a negative denotes a less (more) efficient vertical heat transport in response to an enhanced (weakened) stratification. Combining Eq. (20) with Eqs. (10, 11, 12), the 2-layer system becomes a cubic system. We can expect a third-order spline curve of the warming rate, with two turning points, when we solve this system numerically.
It is expected that an enhanced efficiency of vertical heat transport would slow down the surface warming to some extent, which is shown in Fig. 5d, e. For a strong positive (1 W • m −2 K −2 ) , the surface ocean warming rate will reduce slightly after the fast transient period (Fig. 5e). However, the warming rate will restore and increase to the control level during the slow evolution stage of the deep ocean (red curve, Fig. 5e). To understand the behaviors of surface ocean warming rate, it is useful to examine the time differential equation of Eq. (10): Equation (21) indicates that the change of surface ocean warming rate is determined by external forcing , and by changes in climate feedback −B . The contributions of these physical processes (right-hand side terms in Eq. (21)) to the warming rate change can be easily quantified.
In the beginning when the temperature change in the deep ocean is small, the stratification is enhanced with time, leading to a higher efficiency of downward heat transport, and thus a slightly reduced warming rate in the upper ocean. Figure 6 shows the time series of the terms in Eq. (21). In Fig. 5e, dT 1 dt reaches a maximum in about 50 years; hence, its change rate d 2 T 1 dt 2 is zero (black curve, Fig. 6a). This is determined by the strong negative climate feedback (solid red curve) and the heat uptake by the lower ocean (solid and dashed blue curves, Fig. 6a). In particular, the reduced warming rate during years 50 and 100, that is, the slightly negative d 2 T 1 dt 2 , is mainly due to the enhanced nonlinear OHU related to the effect (dashed blue curve, Fig. 6a), because the other two terms are weakened slightly in this period. During this period, although the vertical stratification is still increasing (red curve, Fig. 5f), its rate of change is weakened. Therefore, the rate of linear OHU becomes weaker (solid blue curve, Fig. 6a), while the rate of quadratic OHU becomes stronger (dashed blue curve, Fig. 6a), leading to a negative d 2 T 1 dt 2 (solid black curve, Fig. 6a). This state of affair cannot last very long because as dT 1 dt decreases, it begins to close the gap between dT 1 dt and dT 2 dt ; specifically, it reduces the rate of linear OHU further (solid blue curve, Fig. 6a), which soon halts the decline of dT 1 dt and makes it positive again. Note that it is the decrease of linear OHU after year 100 that makes the surface warming rate rise again, since the quadratic OHU continues to increase until about year 200. To better understand the role of the linear and nonlinear OHU in the surface warming, we also plot the terms of Eq. (21) with = 0 in Fig. 6b. It is clear that without the quadratic OHU by the lower ocean, the surface ocean warming rate would never decline.
In the late stage of the surface ocean warming, the warming rate approaches the equilibrium in about 1000 years, mainly due to the enhanced climate feedback (solid red  Fig. 6a, b). The upper ocean is still transporting heat downward to the lower ocean, but at a very weak rate (blue curve, Fig. 6a, b). There is no upward heat release from the lower ocean because both the linear and quadratic OHU are negative (positive) for the upper (lower) ocean.
Note that the contributions of processes on the right-hand side of Eq. (21) to the warming rate change depend on the parameters of the system, such as B , 0 , , C 1 , C 2 , and so on. We examined large ranges of the parameters and their effects on the change of surface ocean warming rate (figures not shown). The common features we obtained are as follows: First, it is the enhanced nonlinear OHU (if it occurs in the ocean) that leads to a short-lived decline in the surface warming rate (i.e., negative d 2 T 1 dt 2 ) during the first several decades of global warming, and second, it is the changing climate feedback that makes the warming rate approach the equilibrium in the long run.
The consequence for a negative , i.e., a weaker efficiency for a stronger stratification, is straightforward. The surface warming rate, vertical temperature gradient and surface climate will explode, that is, the runaway climate would become inevitable, which would occur much faster than expected. This situation is shown by the blue curve in Fig. 5d-f, in which is set to −0.5 W • m −2 K −2 . The quadratic OHU under this situation plays a dominant role in transporting heat upward from the lower ocean, contributing greatly to the exploding rise of surface warming rate (dashed blue curve, Fig. 6c), regardless of the strong stabilizing effects of the linear OHU and climate feedback (solid blue and red curves, Fig. 6c). We are not sure whether a negative would be possible in reality. Considering the vast ocean and the complex bio-geo-chemical and physical processes in the real ocean, one cannot exclude the possibility that a negative could occur.

Varying climate feedback and efficiency of vertical heat transport
If we consider both the enhanced negative feedback (Eq. (6)) and enhanced downward heat transport efficiency (Eq. (20)), the upper-ocean temperature Eq. (10) can be rewritten as follows, And the time differential of Eq. (22) is, Equations (22,23) can be easily solved numerically. The surface warming can be reduced substantially, as expected (blue curve, Fig. 7a), under both the enhanced negative climate feedback and downward heat transport. Still, the surface warming hiatus (i.e., zero warming rate) would be unlikely to happen under a linearly increasing external forcing. For the surface ocean warming rate (Fig. 7b), the combined effect of enhanced climate feedback and downward heat transport (blue curve) is not the simple superimposition of enhanced climate feedback and downward heat transport (red and purple curves), due to strong nonlinear effects (Eq. (23)). In fact, the enhanced climate feedback (purple curve) appears to be more effective than the enhanced downward heat transport (red curve) to reduce the surface warming rate (Fig. 7b), since the former is determined by surface temperature (whose equilibrium rate of change is a constant), while the latter is determined Purple curve is for the case with only enhanced negative feedback ( = 1 , = 0 ). Blue curve is for the case with both enhanced negative feedback and vertical heat transport ( = 1, = 1) by vertical temperature gradient (whose equilibrium rate of change is zero), as expressed in Eq. (23), under the condition of comparable parameters ∼ . Therefore, the magnitude of surface ocean warming is mainly determined by the climate feedback (purple and blue curves, Fig. 7a), which also plays a major role in determining vertical temperature gradient by controlling surface temperature (blue curve, Fig. 7c). Different from the discussion in Sect. 3.2, in which the enhanced climate feedback is not considered, here the enhanced quadratic climate feedback can lead to a relatively long-lived decline in surface warming rate (blue curve in Fig. 7b, dashed red curve in Fig. 8). It will take a much longer time for the surface warming rate to approach the equilibrium due to the quadratic climate feedback. The effects from the other terms, in the long run, tend to always increase the warming rate (Fig. 8). In the presence of the quadratic climate feedback (dashed red curve, Fig. 8), the quadratic OHU plays a minor role in reducing the surface warming rate (dashed blue curve, Fig. 8), because the former depends only on surface warming, while the latter depends on vertical stratification, whose change will eventually approach zero.
Some studies suggested that the downward heat transport has contributed a great deal to the recent global warming hiatus. Watanabe et al. (2013) found ocean heat transport efficiency changes oppositely in general circulation models (GCMs) and in nature, namely, it is weakened in GCMs and strengthened in nature. The weakening tendency of in GCMs was seen in the concomitant transient experiments, in which CO 2 is increased at 1% per year, and also in individual models. These authors believed this can be the reason that GCMs tend to overestimate the surface warming trend. This is qualitatively consistent with our results with negative (blue curve, Fig. 5e). England et al. (2014) also pointed out that anomalous winds can cause ocean energy redistribution, and thus induce a warming slowdown; however, rapid warming is expected to resume once the anomalous wind trend abates. This is in agreement with our results with positive (red curve, Fig. 5). Although the triggering mechanisms are different in different studies and people tend to agree that the downward heat redistribution can slow down the surface warming rate, the hiatus induced by enhanced subsurface heat uptake cannot last long and surface temperature will eventually climb up. The 2-box model is quite simple; nevertheless, it captures the essential mechanism, which is qualitatively consistent with that from more complex coupled models and observational studies.

Conclusions and discussion
Inspired by the heated debate on the recent warming hiatus, we investigated two factors that could potentially mitigate GMST warming: the climate feedback of Earth system and ocean vertical heat mixing. The role of the climate feedback was studied using a 1-box EBM. It is shown that under the enhanced heating of GHGs, the GMST warming rate can be reduced if the global climate feedback becomes stronger. However, the warming hiatus appears to be unlikely, unless the global overall climate feedback goes to negative infinity. We realize that even under an extremely negative feedback, its damping effect on the global warming takes a long time to come into play and the magnitude is limited. Regardless of its simplicity, the 1-box EBM tells us that fundamentally, the climate feedback alone would never lead to a short-time hiatus, as observed in the first decades of the twenty-first century.
Under current global warming, the deep ocean, as an energy reservoir, is widely expected to curb the surface warming to some extent. Using a 2-box EBM, we found that the downward heat redistribution can indeed slow down the surface warming rate in the beginning; however, the reduction induced by enhanced subsurface OHU cannot last long. In the long run, the surface warming will recover to be the same speed as that without the deep ocean heat uptake. The 2-box model is simple; nevertheless, it exhibits the essential processes of the vertical heat transport and deep ocean's roles in global warming. Results from this simple model are qualitatively consistent with those from observational studies and complex coupled modelling studies.
We need to be aware that in the real world, GMST change would never be as stable as shown in these simple models. It can be affected by volcano eruptions and other important aerosols, sunspot activities, cloud amount and distribution. Decadal temperature change can be even stronger than the long-term trend (Liebmann et al. 2010), so is the decadal climate feedback. Therefore, it is likely that the recent hiatus had some link to a more negative decadal-changed climate feedback, as suggested in Zhou et al. (2016). SST patterninduced low-cloud anomalies are substantially more negative than the long-term cloud feedback, and could have contributed to the period of reduced warming between 1998 and 2013 (Zhou et al. 2016).
In the box models, several assumptions are used, which are simple but without loss of generality. For example, in the 1-box model we assumed a linear increase of the negative feedback with the increase of GMST, to test whether a stiffening-up of damping effect of the climate feedback could possibly lead to a hiatus. We found that a zero warming rate is not plausible under the linearly increasing external heating. In the 2-box model, the vertical heat transfer efficiency was assumed to change linearly with vertical temperature gradient. Thus, we found that a slowdown of the surface ocean warming will be surely followed by an acceleration of surface warming. The box models contain several important parameters, such as the climate feedback, vertical heat transfer coefficient, the thicknesses of the upper and lower oceans, etc. We tested wide ranges of values of these parameters, and they did not affect our conclusions qualitatively.
Results from simple box models provide further evidence that we should not be blinded by the short-term warming slowdown, as shown in the observations. The long-term trend of global warming cannot be undermined by the occurrence of a one-time, 15-year hiatus. The severity of climate change is closely tied to the total emissions due to anthropogenic activities since the Industrial Revolution, because the carbon dioxide is long lived in the atmosphere. Reducing the GHG emission is the most feasible solution to prevent the Earth from becoming too warm for humans to live on. And a quicker action is the better action.

One-box model with varying climate feedback
Based on Eqs. (1), (2) and (6), the combined equation is: This is Riccati equation, a nonlinear differential equation with constant coefficients of the first order. The parameters in (A1) are defined as follows, To obtain the theoretical solution to (A1), we redefine T = 1 a F � F . Then, we have: Therefore, (A1) can be rewritten as, Eq. (A2) is a linear differential equation with variable coefficients of the second order.
If we define F(t) = Z(t)e − b 2 t , we have: We then define M = m b 2 4 + aKt , where m = (aK) − 2 3 . We have: Eq. (A4) is Airy's Equation, which has two specific solutions: where A(0, M) and A(2, M) represent two kinds of Airy function.
The general solution to Eq. (A4) is: Therefore, F(t) can be expressed as follows: Substituting F(t) into T , we have the expression of T: Since the initial condition is T(0) = 0 , we have F � (0) = 0 , that is, Therefore, The definitions of K 1 and K 2 are given as follows:

Two-box model with constant efficiency of vertical heat transport
Equations (10, 11) give the first-order nonhomogeneous linear differential equation system, which can be rewritten in the matrix form as follows: where.
The associated homogeneous system to (A8) is The eigen equation of the homogeneous system is: where is the eigenvalue of A . The equation for can be simplified as follows, . Eq. (A10) has real solutions, Therefore, the associated eigenvectors can be expressed as follows, where, The e-folding timescale is defined by, Therefore, we have the fast timescale f and slow timescale s , as follows, The fundamental matrix based on (A12) is: The general solution to this homogeneous system is where K = K 1 K 2 is a matrix with undetermined coefficients.
A particular solution is given by the expression: Substituting with Φ(t) and b(t) , we obtain Therefore, the expression of T p can be written as Here, Therefore, T p can be written as: where the fast and slow parameters of a f and a s are defined, respectively, as follows, a f = s f C 1 ( s − f ) B and a s = − f s C 1 ( s − f ) B Consider the following identities, Therefore, T p can be finally re-written as follows: The general solution to the nonhomogeneous equation consists of the general solution to the homogeneous equation and a particular solution to the nonhomogeneous equation. Therefore, the expression of the solution T(t) can be written as Consider the initial conditions T = 0 0 , att = 0 , we have.
Therefore, Substituting (A17) into (A16), we have Finally, we obtain the solutions to T 1 and T 2 , And the vertical temperature difference T v is:  Data Availability Not applicable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. (A21)