Reversible variables with hysteresis
The global mean surface temperature (GMST) increased at different rates in the forward simulations of the historical and four future scenarios. At the end of the 21st century (2081–2100), the GMST of the L, M1, M2, and H simulations was 2.9°C, 3.4°C, 3.9°C, and 4.4°C higher, respectively, than pre-industrial levels. In the corresponding time-reversed simulations, the GMST decreased in response to the time-reversed forcing, especially due to the decline of GHG concentrations (Fig.2). The GMST did not decrease immediately as the GHG forcing decreased, it continued to increase slowly for the first 10 years of the reversed simulations. The higher the CO2 concentration, the longer the GMST lag time. In addition, the GMST decreased at a slower rate than it increased in the forward simulations. Therefore, the GMST remained about 0.6°C and 1.4°C above its original value at the end points of the reversed historical and future simulations, respectively. Furthermore, the GMST returned to the same value in all of the future simulations, even though they led to different degrees of warming in 2100 because of the differences in CO2 concentrations. This indicates that the reversed GMST state depended mainly on the initial GHG concentrations. Similar surface temperature inter-annual variability was also seen in the historical and future reversed experiments.
The upper ocean heat content (OHC, 0-700 m) also showed hysteresis over different time scales in the reversed simulations. The OHC increased by about 30 × 1022 J in the forward historical simulations, and it continued to increase in the first 20 years of the corresponding reversed simulations (Fig.3a), then began to slowly decline. The OHC was only restored by 8 × 1022 J when all of the forces were returned to their 1850 states. The OHC in the future simulations showed a nearly linear growth trend. The OHC growth was 0.78 × 1022 J yr−1, 0.97 × 1022 J yr−1, 1.0 × 1022 J yr−1, and 1.28 × 1022 J yr−1 in the L, M1, M2, and H simulations, respectively. At the end of the 21st century, the OHC had increased by 70–120 × 1022 J in the four forward simulations (Fig. 3b). The OHC continued to increase in the early stages of the reversed simulations, and showed a small decline and convergence in the final stage. However, the final OHC in three of the forward simulations (those other than the H-R simulation), was still higher than that at the end of the 21st century in the forward simulations. Moreover, the OHC did not return to its initial states, and there was a disparity of about 10 × 1022 J between the two adjacent reversed simulations.
As the OHC is mainly determined by the potential temperature of the ocean, we can infer that the potential temperature is irreversible. The spatial patterns of the potential temperature (fig.4) show that the CO2 scenario affects the recovered ocean-temperature states. These patterns indicate that the higher the CO2 concentration, the more difficult it is for ocean temperature to return to its initial value, especially in high-latitude oceans where the positive sea ice feedback is relatively large. The maximum difference between the start of the H and the end of the H-R experiments exceeded 2°C, while that between the L and L-R experiments was less than 1.5°C.
Reversible variables without hysteresis
Climate variables affect various physical, biological and chemical processes, and have different roles in controlling global energy transfer and the carbon and hydrological cycles. Hence, the degree of reversibility of these variables differs. We investigated the reversibility of variables such as total runoff, soil moisture, and convective precipitation.
As CO2 concentration increased and consequently temperatures increased, there was a significant increase in globally averaged total runoff and convective precipitation and a decrease in soil moisture, ground evaporation, sea ice extent, and Atlantic Meridional Overturning Circulation (AMOC) (Fig.5). The higher the CO2 concentration, the larger the increase (or decrease) of these variables in the future simulations (Fig. 6). Generally, the total runoff was reversible, and the magnitude of the change in the reversed simulation was close to that in the forward historical and future simulations (Fig.5b and 6b). A similar phenomenon was observed for the convective precipitation in the historical simulation (Fig.5a). However, the convective precipitation did not return completely to its original states in the four future simulations, and there remained a small disparity of about 0.07 mm/day−1. The ground evaporation recovered almost completely in the historical and future simulations (Fig.5c and 6c). The soil moisture and sea ice extent in the northern hemisphere is more sensitive to increases in CO2 concentration, and both exhibited weak reversibility. The returned sea ice extent was about 1.0 million and 3.0 million km2 less than its original state in the historical and future simulations, respectively. This suggests it is difficult for the soil moisture and sea ice extent to recover from reductions caused by global warming. The AMOC showed greater reversibility in the historical and future simulations (Fig.5f and 6f). In contrast with the GMST, these variables reversed in direct response to the decline of the GHGs.
Furthermore, as with the GMST, the higher the CO2 concentration, the greater the rate of increase or decrease of the climate variables in the forward future simulations (Fig.6), which led to different climate states in 2100. However, these variables returned to the same value in all four of the reversed future simulations. This indicates that for most of the variables, the returned states mainly depended on the original states of the natural and anthropogenic forces rather than their final states.
Irreversible variables
The steric sea level increased in the forward historical and future simulations. However, unlike other variables in the climate system, it did not appear to be reversible and continued to increase after all of the forces were reversed (Fig.7). This phenomenon was consistent with other work focused on sea level reversibility under the stabilization of CO2 concentration in the atmosphere or zero CO2 emissions (Matthews, 2010; Mikolajewicz et al., 2007; Solomon et al., 2009). The sea level had increased by 0.26, 0.28, 0.30, and 0.32 m by the end of the 21st century in the four forward simulations. It demonstrated a further increase of 0.2 m during the four corresponding reversed simulations. Hence, the returned states of the sea level were about 0.5 m higher than the original states in the future simulations, and about 0.27 m higher in the historical simulations (Fig.7). However, as the sea level continued to rise in the four reversed simulations, it showed the same disparity (about 0.02 m) in the final states of the four reversed simulations. This was the most significant difference between the changes in global sea level and other variables in the climate systems.