Radiative Energy Flux Variation from 2001–2020

: Radiative energy flux data, downloaded from CERES, are evaluated with respect to their variations from 2001 to 2020. We found the declining outgoing shortwave radiation to be the most important contributor for a positive TOA (top of the atmosphere) net flux of 0.8 W/m 2 in this time frame. We compare clear sky with cloudy areas and find that changes in the cloud structure should be the root cause for the shortwave trend. The radiative flux data are compared with ocean heat content data and analyzed in the context of a longer ‐ term climate system enthalpy estimation going back to the year 1750. We also report differences in the trends for the Northern and Southern hemi ‐ sphere. The radiative data indicate more variability in the North and higher stability in the South. The drop of cloudiness around the millennium by about 1.5% has certainly fostered the positive net radiative flux. The declining TOA SW (out) is the major heating cause (+1.42 W/m 2 from 2001 to 2020). It is almost compensated by the growing chilling TOA LW (out) ( − 1.1 W/m 2 ). This leads to ‐ gether with a reduced incoming solar of − 0.17 W/m 2 to a small growth of imbalance of 0.15 W/m 2 . We further present surface flux data which support the strong influence of the cloud cover on the radiative budget.


Introduction
In the big picture, climate variations originate from variations of the radiative balance at the top of atmosphere (TOA). Surpluses of the EEI (Earth energy imbalance) or net radiative energy fluxes, as measured by satellite mounted radiometers, lead to an increase of the climate system enthalpy and vice versa [1][2][3][4]. For about two decades, the CERES Energy Balanced and Filled (EBAF) Ed4.1 [5,6] offers datasets for a variety of radiative fluxes, and, thus, provides a basis to scrutinize the radiative climate driving forces and shine light on the cause-and-effect relation between radiation and temperature change.
As an independent but less direct source of information of the climate system enthalpy change, there are several studies and reconstructions [7][8][9][10][11][12] of the ocean heat content (OHC) which represents the bulk of the climate system enthalpy, estimated to be about 90%. Assuming this fraction of 90% were a longer-term constant, one can trace back the time-development of the climate system enthalpy. Von Schuckmann et al. [13] have combined radiative, ocean heat and other data to reconstruct an enthalpy curve back to 1960 and have found an accelerated heating since 2010. Recently, Loeb et al. [14] found a good agreement between radiative (CERES) and OHC data for the period mid-2005 to mid-2019. These authors have further studied the influencing factors for the shortwave (SW) and longwave (LW) radiative fluxes and concluded that cloud changes have fostered the downwelling shortwave radiation.
Dewitte et al. [15] have analyzed CERES datasets for the period from 2000 to 2018 and found an EEI value of about 0.9 W/m 2 but with a declining trend going in line with a declining time-derivative of the latest OHC data obtained from Cheng et al. [9]. Based upon recent CERES data, Loeb et al. [16], Wong et al. [17] and Ollila [18] reported an increasing downwelling shortwave (SW) radiation. Loeb et al. reported a decreasing TOA SW trend, mainly caused by a reduction in low cloud cover, and Ollila concluded that this increasing downwelling SW, which is particularly strong since 2014, may be responsible for a new wave of heating after the hiatus. This finding is in conflict with the assumption that further global warming originates mainly from the LW radiation capture caused by greenhouse gases, i.e., a decline of outgoing LW.
The obvious and substantial, if not overwhelming role of clouds for the radiation budget and climate system enthalpy and, hence, for the question about the root cause of the further development of global warming, is nowadays still a vaguely known factor. The cloud-albedo feedback is deemed to be essential for climate modeling [19] but is still poorly understood. According to a cloudiness dataset from EUMETSAT/CM SAF [20] there was a significant drop in global cloudiness around the year 2000, which has not yet fully recovered, and which certainly has affected the radiative net flux in the time-period considered here. In this paper, we report radiative flux data and trends in cloudy and cloud-free regions, obtained from CERES and other sources and relate them to the TOA and surface radiative budgets and climate system enthalpy. We further attribute the differences between the Northern and the Southern hemisphere.
Finally, we discuss these results in a longer-term context and suggest a possible correlation of cloud cover shifts such as the one around the millennium with the AMO (Atlantic Multidecadal Oscillation).

Datasets and Evaluation
We have downloaded the datasets presented and analyzed here as monthly averages from Jan-2001 to Dec-2020, obtained from the NASA Langley Research Centre CERES [5,6] ordering tool at (http://ceres.larc.nasa.gov/), (accessed on 29 July 2021), spatial resolution 1° X 1°, unless marked otherwise. The changes from Edition 4.0 [5,6] to Edition 4.1 are summarized in: https://ceres.larc.nasa.gov/documents/DQ_summaries/CERES_EBAF_Ed4.1_DQS.pdf (accessed on 29 July 2021). We transformed them into annual averages and used these annual averages without further changes. Only the datasets marked as "Cloudy Areas" were calculated by the authors from the "All Sky" and the "Clear Sky (for cloud-free areas of region)" data according to Equation (1), year by year, making use of the annual averages of CERES "cloud area fraction" (abbreviated as c) and with obvious notation for the fluxes (F) involved.
To interpret the results, it is of some importance to realize the wavelength channels in CERES, which are:  Hence, the infrared active molecular vibrations and rotations of the greenhouse gases fall into different categories: fundamental bending and deformation vibrations (H2O, CO2, CH4, etc.) as well as rotations (H2O) are in the LW-range and the OH (oxygen-hydrogen bond) and CH (carbon-hydrogen bond) stretching vibrations of H2O and CH4 including their overtone and combination band spectra fall into SW. Even the asymmetric stretching vibration of carbon dioxide falls into the SW. That means that not only the LW but also the SW fluxes are influenced by infrared absorption.
The effective TOA ("Top of Atmosphere") focal height of CERES is ca. 20 km. Surface fluxes as well as "Cloud Effective Temperature" and "Cloud Area Fraction" are provided by CERES, transformed into annual averages without further changes with the exception of the "Cloudy Areas" as explained above.
Enthalpy data were calculated from the TOA radiative fluxes by means of Equation (2). Statistical error analysis is made by means of linear regression and summarized in the Appendix A. It includes the slope (trend) and errors, p-value, R 2 and confidence intervals. In the main text, we are just referring to a statistical uncertainty in the form x ± ∆x. Details can be seen in the Appendix A.
In the text, we use the absolute values of the fluxes (in W/m 2 ), respectively, their actual changes from 2001 to 2020 as well as the statistical trend data (in W/m 2 per decade). The further are used to discuss the actual changes, whereas the latter are linear estimations and specify the average tendencies.

Results and Discussion: Top of Atmosphere (TOA) Radiative Balance
Figures 1-4 show the incoming solar, the outgoing shortwave (SW), the outgoing longwave (LW) and the net radiative fluxes for "All Sky" conditions, which includes "Cloudy Areas" and "Clear Sky" areas. The figures include the global, Northern (NH) and Southern (SH) hemispheric fluxes. The incoming solar irradiation is slightly decreasing in the total window of 0.3 to 200 μm, whereas as the outgoing LW flux became stronger during the last two decades. These both effects alone would have diminished the climate system's enthalpy, if not the reduction of the outgoing shortwave flux would have overcompensated these two effects. This SW effect caused an increasing net flux during these two decades. Loeb et al. [16] and Ollila [18] had already noticed the declining SW (out) which corresponds to an increasing downwelling SW trend, because the incoming solar flux was almost constant during the period of time considered.   Dewitte et al. [15] have analyzed CERES data until 2017, and found a positive TOA net flux but with a declining trend. In the data presented here, the TOA net flux has a similar average positive value but is slightly rising with time. The TOA net flux is a tiny difference of large numbers and is therefore prone to a rather large relative uncertainty. Therefore, the authors are unable to decide whether the trend is, actually, going up or down. Figure 4. TOA net flux for global, NH, and SH. The Southern hemisphere obtains a higher flux (+ 1.7 W/m 2 ) than the North which is due to the aphelion/perihelion of the Earth's rotation around the sun, resulting in a stronger radiative flux during the SH summer, not compensated by the weaker SH winter. Nevertheless, the NH has warmed up more strongly during the given period of time.
The TOA net flux was +0.75 W/m 2 in 2020. The data shown in Figures 1-3 suggest that the root cause for the positive TOA net flux and, hence, for a further accumulation of energy during the last two decades was a declining outgoing shortwave flux and not a retained LW flux. The dominating influence of a declining TOA outgoing SW is not in line with the assumption that the global warming in the last 20 years is primarily caused by holding back the LW due to the greenhouse effect. As discussed below, the greenhouse effect plays a role, however, for the period 2001-2020 at TOA it was certainly not dominating.
Some residual uncertainty remains regarding the sign of the TOA net flux. By simply adding the uncertainties of the gross in-and outgoing fluxes, the uncertainty of the difference could easily be in the order of 1 W/m 2 or more, which is larger than the 0.75 W/m 2 reported for the year 2020. However, as shown below, the independently observed OHC data comply well with a positive radiative net flux, and it can be rather safely assumed that there was indeed a positive net flux during the last two decades.
The declining outgoing SW was mainly caused by less reflected solar light which points towards the role of the clouds. Only a minor SW portion may be absorbed by the increasing amount of water vapor by means of their OH stretching vibrations (ca. 3 μm) and overtone, respectively combination band absorption (ca. 1.9 μm, 1.5 μm, 1.0 μm, etc.). Figure 5 compares outgoing SW under "Clear Sky" conditions (left side) with the "Cloudy Areas" (right side) using identical scale ranges. The decline is pronounced in the NH and the absolute effect is greater over the "Cloudy Areas". The "Clear Sky" SH trend is almost flat. Over the clouds, we find a declining trend in NH and SH, still stronger in the North but with a visible decline in the South also. It seems as if the trend accelerated around 2010-2014. In fact, for the global "Cloudy Areas" data the slope is −0.043 W/m 2 a for 2001-2010 and −0.16 W/m 2 a for 2011-2020, i.e., four times larger. Shortwave outgoing radiative flux (TOA) over "Clear Sky" (a) and "Cloudy Areas" (b) for global average, NH, and SH. Please note, there is an apparent inconsistency when comparing these plots with Figure 3, in which there is no significant difference between NH and SH. The reason is a different average cloud area fraction of about 70% in the South and 64% in the North. Therefore, the cloud-area-weighted "All Sky" data are very similar.
The changes in the cloud albedo are overwhelming the changes of the land/ocean albedo. If we define surface albedo as the ratio of the TOA outgoing SW and the incoming solar fluxes, we find a reduction from 15.80% to 15.56% (−0.23% absolute) in the "Clear Sky" from 2001 to 2020 stemming from changes in the land and ocean surface. In the "Cloudy Areas" we find a much stronger reduction from 35.85% to 35.28% (−0.58% absolute). The "All Sky" values changed from 29.33% to 28.80% (−0.53% absolute). The area-weighted effect of the albedo change is −0.08% (absolute) for the land/ocean part and −0.39% (absolute) for the clouds. The cloud effect on the albedo change is 5.2 times larger than the land/ocean effect. In other terms, about 15% of the observed albedo change is due to the diminishing land/ocean albedo, and 85% is due to the diminishing cloud albedo. Hartmann and Ceppi [21] as well as Dewitte et al. [15] have studied the influence of sea ice extend (SIE) in the Northern polar region on the reflected SW and found a good correlation of the "Clear Sky" reflected solar light with SIE for that region. We also find such a pronounced regional effect, however, it is limited to a too small area to explain the globally declining "All Sky" outgoing SW presented here.
The outgoing TOA LW flux shows a different trend. Figure 6 reveals the increasing LW flux which is much steeper over the clouds (0.35 ± 0.13 W/m 2 per decade) than over the clear sky (0.04 ± 0.1 W/m 2 per decade) which is essentially constant. Furthermore, there is a striking difference between North and South. The difference between NH and SH is smaller for the "Clear Sky" areas, and the North emits a higher LW flux than the South. For the "Cloudy Areas" the difference between NH and SH is twice as large but the South emits a higher LW flux than the North. The latter goes in line with a higher cloud effective temperature ( Figure 7) in the South and the former corresponds to a higher surface temperature on the Northern hemisphere. The smaller absolute value of the flux over the "Cloudy Areas" is due to a lower cloud temperature compared with the surface. Only the much steeper trend for the "Cloudy Areas" is puzzling and points towards a change in the cloud structure.  The TOA net fluxes over the "Clear Sky" and "Cloudy Areas" exhibit-as expecteda striking difference ( Figure 8). The former has a strongly positive ("heating") net flux whereas the latter has a large negative ("cooling") value. About 2/3rds of the sky is covered with clouds of all types, heights, thicknesses and other properties, and the mix of ca. The difference between NH and SH is much more pronounced over the "Clear Sky", but the rising trend is somewhat steeper over the "Cloudy Areas". The enormous influence of the clouds for the energy balance and, hence, for the climate system enthalpy is obvious. This is also visible from Figure 9 in which the cloudiness and the outgoing LW (here: OLR) are presented for a longer time span. The (inverse) correlation coefficient between OLR and cloudiness is 74%. The cloud window that opened around the year 2000 has certainly fostered the heating. The CERES data agree quite well with the HIRS OLR data, but the comparison also shows that some methodical uncertainties in the order of 1 W/m 2 must be taken into account.  The terms "cloudiness" and the CERES term "cloud area fraction" have different definition and methodology but both are associated with the cloud cover. The cloudiness is currently ca. 61.3% whereas the cloud area fraction from CERES is ca. 67.5%. Nevertheless, from Figure 8 it follows, that a change of cloud cover by 1% would cause a change of the TOA global net flux of 0.26 W/m 2 so that the drop of 1.86% of cloudiness in Figure 9 could have caused an effect of ca. 0.5 W/m 2 , a large portion of the observed average value for 2001-2020 (0.8 W/m 2 ). Actually, adding the 0.5 W/m 2 to the 20th century average (ca. 0.3 W/m 2 ) gives the current value.

Results and Discussion: Surface Fluxes
CERES data provide an opportunity to differentiate between "Clear Sky" and "Cloudy Areas" not only at TOA, but also for the surface fluxes and this enables us to analyze the influence of clouds on the surface radiation budget as well. At first, under idealized balanced conditions, the TOA radiative budget is exactly zero. At the surface, however, the radiation budget under idealized balanced conditions is far away from zero.
It is ca. 110 W/m 2 and this radiative imbalance is compensated by enthalpy changes, mainly by the evaporation of water so that the sum of the overall balance of radiative flux and the time derivative of the enthalpy at the surface strives towards zero.
The surface fluxes show striking differences between "Clear Sky" and "Cloudy Areas" as well. Figure 10 shows SW up-and downwelling surface fluxes for different cloud coverages. While in the "Clear Sky" regions the downwelling SW (Figure 10a (Figure 10a-c), with decreasing trend. The surface SW net flux is rising more pronounced in the "Cloudy Areas" by +0.93 ± 0.18 W/m 2 per decade compared with the "Clear Sky" regions with only +0.24 ± 0.10 W/m 2 per decade. We conclude that the clouds have increased their transmittance for the SW radiation. This is possible by means of changes of the cloud type, height, particle size, phase mix, geographical location, interaction, etc. The parameter "Cloud Optical Density", which is available from CERES, however, could neither support nor deny this interpretation because the rather high variance prevents to identify a clear trend. The surface downwelling LW flux ( Figure 11) has a different trend than the SW. The "Clear Sky" areas ( Figure 11a) show a strongly rising trend, compared with a flat or slightly declining trend in the "Cloudy Areas". In the "All Sky" average, the downwelling LW has a rising trend half as strong as the downwelling SW and with high variance, namely +0.27 ± 0.34 W/m 2 per decade. The LW upwelling fluxes (Figure 11d-f) show strong trends around 1 W/m 2 per decade and are almost independent of the cloud coverage in size and trend. Here, the surface temperature is the main driver, and the radiative flux increases as a function of the surface temperature. The LW net surface fluxes show a strongly declining trend that is accelerating after 2016 with the exception of the "Clear Sky" regions. As a consequence of the rising SW and the falling LW net fluxes, the "All Sky" surface net flux has no clear trend but a significant fluctuation ( Figure 12).  At this point, we describe the direct observation of the greenhouse effect, namely by comparing the "Clear Sky" LW upwelling radiation (Figure 11d) and the outgoing TOA LW flux (Figure 6a, left side). The first is rising with +1.22 ± 0.22 W/m 2 per decade, the latter trend is literally flat with +0.04 ± 0.10 W/m 2 per decade. Hence, in the absence of clouds, a large portion of this additional upwelling LW radiation is absorbed by the increasing greenhouse gas concentration (CO2 has increased from 371 to 414 ppm obtained from https://gml.noaa.gov/webdata/ccgg/trends/co2/co2_annmean_mlo.txt (accessed on 29 July 2021), water from 10.44 to 10.71 g/kg at 1000 mbar from 2001 to 2020 obtained from https://psl.noaa.gov/cgi-bin/data/timeseries/timeseries1.pl, (accessed on 29 July 2021). We have applied the Lambert-Beer's law and correlated the logarithmic ratio of the outgoing TOA LW flux and the "Clear Sky" LW upwelling flux with the rising greenhouse gas concentrations and found a very good linear correlations for CO2 (R 2 = 0.92) and water vapor (R 2 = 0.72).
Another observation is that these correlations vanish for the "Cloudy Areas". The clouds absorb most of the LW flux and are emitting a weaker thermal flux in all directions. The part designated to reach the TOA point has little chance to be absorbed by greenhouse gases because the pathway is shorter, the pressure is lower, the water vapor is sparse, and the initial flux is weaker. All this works towards much lower infrared absorption above the clouds and the TOA outgoing LW rises with 0.35 ± 0.13 W/m 2 per decade.
In absolute terms, the LW surface upwelling flux is about 398 W/m 2 for "clear sky" (Figure 11d) and only ca. 268 W/m 2 at TOA (Figure 6a) which is an attenuation of 130 W/m 2 or 33%, caused by greenhouse gases, by scattering, and by small clouds below the measurement threshold. In the "Cloudy Areas" we start with about 399 W/m 2 LW surface emission (Figure 11f), which is absorbed by the clouds and re-emitted at a lower temperature of about 258 K (Figure 7). This cloud average temperature corresponds to about 260 W/m 2 emission (estimated from Stefan-Boltzmann's law), and ca. 227 W/m 2 of it reaches the TOA (Figure 6b). The attenuation here is only 33 W/m 2 or 12.6%. There again, the corresponding downwelling back-radiation is absorbed by the clouds and contributes to elevating their temperature. Only a minor share of this smaller part reaches the Earth's surface.
Hence, the rise of the greenhouse gas concentration from 2001 to 2020 had a measurable effect on the LW flux in the "Clear Sky", covering about 1/3rd of the Earth surface. In the cloudy part, about 2/3rd, this effect was much smaller, if significant at all.

Effect on the Climate System's Enthalpy
With the one and only exception of radiation, the Earth is an almost perfect adiabatic thermodynamical system. Material exchange, convection and heat diffusion is prevented by gravity and all other sources or mechanism of energy exchange with the environment, the space, such as starlight and cosmic particle flux are negligible in terms of enthalpy. The cosmic particle flux could, however, have a significant indirect effect by its influence on the cloud formation which again provides a very strong leverage for the radiation budget. This potentially very important indirect effect was described by Svensmark [22] and by Shaviv [23].
Radiative datasets are available for only a few decades, however, ocean heat content (OHC) data were reconstructed for much longer and, assuming that constantly ca. 90% of the climate system's enthalpy ended up as ocean heat, provide an opportunity to reconstruct the enthalpy development for a longer period of time. OHC reconstructions were reported by Levitus et al. [7], Roemmich et al. [24], Cheng et al. [10] and Gebbie and Huybers [11]. Von Schuckmann et al. [13] have developed an enthalpy summary for 1971-2018. They also reported the distribution over the various degrees of freedom such as the heat of oceans, air and land, evaporation of water and melting of ice and volume work. decade over that same time-period. The data analyzed in this paper resulted in a net flux of 0.8 W/m 2 for the period 2001-2020 and this corresponds to an enthalpy increase of about 240 ZJ during this period. However, the root cause of this enthalpy gain is in dispute. This becomes clear from the strong effects of the shortwave radiation changes, the cloud cover changes and, also, from a longer-term development of the enthalpy. Figure 13 presents a longer-term enthalpy curve, constructed from various sources, radiative and oceanic, the latter under the assumption of a constant 90% OHC share. Figure 14 shows a shorter-term view of the same data. In both figures, the distinct zero-level of each of the enthalpy data was shifted to fit. The current enthalpy is still ca. 600 ZJ below the medieval maximum 1000 years ago grounded on the 2000 years OHC reconstruction reported by Gebbie and Huybers [11]. This would be compensated in only 50 years with the presently observed +0.8 W/m 2 net flux and in about 100-200 years for the average rate as during the 20th century. With a possible interception by another phase of a negative radiative net flow, it would take centuries until we reach the medieval OHC maximum again.
Hence, the decisive question is; whether the currently observed strong heating phase C is temporary like the phases A and B in Figure 13 or the beginning of a steadily rising warming period. A third possibility is that phase C is a combination of the natural induced A/B type and steadily warming. The enthalpy-enhancing effect of the greenhouse gases (including water vapor) is clearly visible in our analysis, however, we found, in line with other authors [14,15,18], that this effect is only one out of several other important factors. In the period 2001-2020, the greenhouse gas effect is prominent in the "Clear Sky" areas but overcompensated by the SW fluxes and cloud effects. This is also evident from the bridge chart in Figure 15 below. As an alternative, the phases A, B and C since 2001 could be interpreted as a kind of "global brightening" phases followed by "global dimming" (for a review of global brightening/dimming see [27]).
The rather short period of only 20 years prevents us from answering this key question at this time. Furthermore, there is no validated physical mechanism that led to the declining TOA shortwave emission and to the cloudiness change around the year 2000. Figure 13. Climate system enthalpy since 1750, reconstructed from radiative and ocean heat data taken from the publications of Schuckmann et al. [13], Ollila [18], Dewitte et al. [15], Levitus et al. [7] and NODC [7], Cheng et al. [10], Gebbie and Huybers [11] and this work. The OHC-Data were divided by 0.9. The two heating impulses A and B had a span of about 25 years, and similarly, high net fluxes as the presently observed phase C. The zero levels of the enthalpy datasets was set to match. The AMO values were downloaded from http://www.psl.noaa.gov/data/timeseries/AMO/ (accessed on 29 July 2021) and transformed into annual averages. A sine-fit is added to the AMOdata (right scale). The arrows at ca. 1860, 1925 and 2000 indicate the sign change of the AMO-Index from negative to positive. Figure 14. Expanded enthalpy plot. The cloudiness [20] is included (right scale). The rise of the enthalpy after the millennium coincides with the drop of cloudiness. The latest OHC data from NODC exhibit such a flatter tendency for about one year. Figure 15. Bridge-chart for the actual changes at TOA from 2001 to 2020 ("start-to-end change"), dominated by the change in the outgoing shortwave radiation over "Cloudy Areas". The declining SW (out) is the major heating cause. It is almost compensated by the chilling LW (out). The incremental effects are weighted by area (see the example in the caption of Figure 16).

Conclusions
Radiative energy flux data from CERES were analyzed and showed in accordance with OHC data a further increasing climate system enthalpy during the period 2001-2020. The total enthalpy rise amounted to about 240 ZJ in these two decades. As Figure 15 shows, the major driving effect was the declining shortwave TOA emission. The TOA outgoing longwave emission has increased and therefore reduced the TOA net flux.
Generating the CERES data is a demanding task and requires sophisticated technology and models which are vulnerable and prone to uncertainties. Liu et al. [28] have pointed out significant uncertainties of satellite datasets and discussed the role of the lateral energy flow. Su et al. [29] have recently pointed out the importance of maintaining the consistency among the components of the measuring system. On the other hand, Loeb et al. [14], Johnson et al. [30] and, also, Dewitte et al. [15] have shown that the CERES net flux agrees well with the independently observed OHC data. This good agreement, also confirmed by our analysis, justifies some confidence in the CERES datasets used.
We could identify the effect of the anthropogenic greenhouse gas emissions from 2001 to 2020 in the "Clear Sky" LW part but not in the "Cloudy Areas" and not in the SW. At the same time, we find, in accordance with the analysis of Loeb et al. [14] and Ollila [18], that the major changes for the TOA energy budget during this period of time stemmed from the clouds for SW and LW, as well as the ground temperature in the LW.
Loeb et al. [14] pointed out, that the direct aerosol effect is rather small, but the indirect effect via the cloud formation may be larger. The shift from a negative to a positive PDO (Pacific Decadal Oscillation) index as an additional factor for the net TOA flux is mentioned in [14].
Our analysis, which differentiates between clear sky and cloudy areas, support that in view of the historical heating steps shown in Figure 13, that the currently observed high radiative net flux has a large intrinsic component. As shown in Figure 13, the heating phases coincide with the AMO change from negative to positive. It has been shown by several authors that AMO may be an important intrinsic climate factor [31][32][33]. A similar discussion was held about 20 years ago in papers of Chen et al. [34] and Wielicki et al. [35], both of them emphasized the underestimated decadal natural variabilities in the tropical regions.
The start-to-end bridge charts in Figures 15 (TOA) and 16 (surface) are highlighting the effect of the "Cloudy Areas" and the shortwave radiation in a slightly different view. In these figures, we look at the actual start and end data, their actual differences from 2001 to 2020, and the actual incremental contribution of each of the radiative categories. Please note, that the increments are weighted by area. In Figure 15 (TOA), the biggest changes originate from the "Cloudy Areas". In Figure 16 (surface) the largest heating impulse stemmed from the increasing downwelling and decreasing upwelling shortwave fluxes in the "Cloudy Areas" (together +1.23 W/m 2 ), the strongest cooling effect was the decreased longwave downwelling flux in the "Cloudy Areas" (−1.48 W/m 2 ), followed by increasing LW upwelling fluxes in the "Cloudy Areas" (−1.34 W/m 2 ) and the "Clear Sky" (−0.93 W/m 2 ). The change of the longwave downwelling radiation can be interpreted in part as the additional effect of the increased greenhouse gas concentration. For the "Clear Sky" it is +1.20 W/m 2 . In the "Cloudy Areas", this effect is negative (−1.48 W/m 2 ) so that the sum of these values is −0.14 W/m 2 . The −0.93 W/m 2 of the "Clear Sky" upwelling longwave should be caused by the increased thermal emission due to the higher surface temperature.
There are distinct differences between the Northern and the Southern hemisphere. Generally speaking, the South was more stable than the North in trends and variances of almost all radiative quantities. This could be due to the larger ocean share of the surface in the South.
Finally, the key issue, i.e., whether the current heating phase is a temporary phase or a permanent phenomenon, can be judged only on the basis of a longer observation time.
In the latter case, the physical mechanism behind the "shortwave heating" [18] or a possible "cloud thinning", as discussed by several other authors [36][37][38] should be understood, because it could accelerate the warming trend. In the former case, the strong net flux of +0.8 W/m 2 should decrease naturally.

Acknowledgments:
We are grateful to Nic Lewis and Frank Bosse, whose reviews helped to improve the manuscript.

Conflicts of Interest:
H-R D is a physical chemist with many years of research experience in the chemical industry. His scientific training enables him to work effectively with interdisciplinary datasets and interpret their respective implications. The project did not receive any funding, nor was it commissioned or supported by industry. FV has been for many years Professor at the University of Hamburg, Department of Chemistry and was employed in the Renewables Energy sector. He has many years of academic experience and a solid publication track record in chemistry and climate science. Please note that the first year 2001 was set to zero so that "Intercept" is equal to the linear estimate for the year 2001. "Error Intercept" denotes the statistical standard error. "Slope" is the result of a linear regression. "Error of Slope" denotes the statistical standard error. "Trend" is the slope multiplied by 10 and the "Rel. Error" is the Ratio of "Error of Slope" and "Slope", multiplied by 100. The "p-value" for the slope is the result of a linear regression. A p-value below 0.05 rejects the Null hypothesis (i.e., slope = 0) with 95% confidence. p-values below 1E-4 are cut off after four digits. "CI low" is the lower boundary and "CI high" is the upper boundary of the 95% confidence interval for the slope. "∆ CI" is the width of the confidence interval. "R 2 " is coefficient of determination and indicates the strength of the correlation. An R 2 near zero means that the respective flux is approximately a constant for the time-period considered. For convenience, " Fig." gives the reference to the figures in this paper. Please note that the first year 2001 was set to zero so that "Intercept" is equal to the linear estimate for the year 2001. "Error Intercept" denotes the statistical standard error. "Slope" is the result of a linear regression. "Error of Slope" denotes the statistical standard error. "Trend" is the slope multiplied by 10 and the "Rel. Error" is the Ratio of "Error of Slope" and "Slope", multiplied by 100. The "p-value" for the slope is the result of a linear regression. A p-value below 0.05 rejects the Null hypothesis (i.e., slope = 0) with 95% confidence. p-values below 1E-4 are cut off after four digits. "CI low" is the lower boundary and "CI high" is the upper boundary of the 95% confidence interval for the slope. "∆ CI" is the width of the confidence interval. "R 2 " is coefficient of determination and indicates the strength of the correlation. An R 2 near zero means that the respective flux is approximately a constant for the time period considered. For convenience, " Fig." gives the reference to the figures in this paper.