## 1. Introduction

The hydrological cycle and dynamics in the Earth system are driven by energy received from the sun. Once averaged over a year and over the entire globe, about 71% of the solar irradiance is absorbed by Earth (e.g., Stephens et al. 2012). While dynamics in the atmosphere and ocean distributes energy absorbed by Earth, energy is converted into different forms. A part of shortwave irradiance absorbed by ocean and land is used to evaporate water vapor. Water vapor enters the atmosphere and heats the atmosphere when it condenses. Energy transport is also associated with heating and cooling by radiation, dynamics, and water vapor phase change, which in turn alter entropy of the Earth system.

Entropy is produced by heating and cooling by irreversible processes. In addition, entropy is carried by radiation. Entropy produced by a blackbody is, therefore, the sum of entropy produced by radiative cooling and entropy carried by the blackbody radiation (Planck 1913). For the Earth system, entropy is imported by radiation emitted by the sun. Heating due to absorption of solar radiation and scattering (Wu and Liu 2010) within the Earth system produce entropy. In addition to diabatic heating by irreversible processes occurring within the system, radiative cooling and heating by longwave radiation contribute to entropy production. Longwave radiation emitted to space then exports entropy. Earlier studies refer to entropy produced by irreversible processes within the Earth system as material entropy (e.g., Bannon 2015; Bannon and Lee 2017).

Global entropy budget was first estimated by Peixoto et al. (1991). Their estimate was primarily based on in situ surface and upper-atmosphere data taken between 1963 and 1973 (Peixoto and Oort 1992). Stephens and O’Brien (1993) used satellites’ broadband radiance measurements to estimate entropy production by longwave radiation emitted to space. Observational estimates of entropy production within the Earth system by absorption of shortwave radiation, however, have not been updated since the estimate by Peixoto et al. (1991) was made. Goody (2000) extended work by Peixoto et al. (1991) and estimated entropy produced by irreversible processes within the Earth system. Lucarini et al. (2011) used climate models and estimated entropy production by nonradiative irreversible processes. The analysis of thermodynamics of the climate system is also used to understand how general circulation changes under a warming climate (Kjellsson 2015; Laliberté et al. 2015).

In this study, we revise entropy budget estimated by earlier studies using satellite observations that have been taken since March 2000. We use cloud and aerosol properties derived from satellites to estimate entropy production by the absorption of shortwave irradiance and by the emission of longwave irradiance to space. Radiation also carries entropy, but we primarily focus on entropy produced by heating and cooling within the atmosphere and at Earth’s surface in this study.

Before we present entropy production estimates in section 4, we use a 1D model to illustrate the relationship between entropy produced by radiation and entropy produced by irreversible processes within the system in section 2. Section 3 explains data products used in this study and how entropy production is computed. Section 5 utilizes the results with the concept illustrated by the 1D model to discuss how entropy production changes when the energy input to Earth changes.

## 2. Conceptional 1D model

We use a 1D isothermal atmosphere model discussed by Bannon (2015) to describe the relationship between energy budget and entropy budget of the Earth climate system (Fig. 1). The 1D model also helps to understand how entropy produced by shortwave absorption and longwave emission is related to irreversible processes within the Earth system. In addition, we extend the model to help relating entropy produced by nonradiative irreversible processes and entropy produced by heating and cooling by radiation within the Earth system.

We let the incoming solar irradiance at top of the atmosphere (TOA) be *F*_{0}. The shortwave irradiance absorbed by the atmosphere is *βF*_{0}, where *β* is the absorptivity of the atmosphere. We assume that the atmosphere does not reflect shortwave irradiance, but the surface reflects shortwave irradiance *αF*_{0}, where *α* is the reflectivity of the surface. The rest of the shortwave irradiance (1 − *α* − *β*)*F*_{0} is absorbed by the surface. Excluding scattering by the atmosphere does not affect establishing the relationship between entropy produced by radiation and nonradiative processes. The surface is assumed to be black to longwave radiation and emits the longwave irradiance proportional to the fourth power of the temperature *T*_{sfc} is the surface temperature. The emissivity of the atmosphere is *ε* so that the atmosphere emits the irradiance *T*_{atm} is the temperature of the isothermal atmosphere. The atmosphere absorbs the longwave irradiance *εF*_{sfc} emitted by the surface and the rest (1 − *ε*)*F*_{sfc} is transmitted to space. Energy is also transported by turbulence from the surface to the atmosphere at the rate of *F*_{tur}, which represents a nonradiative irreversible process in the model.

At the top of the atmosphere, the system receives entropy of radiation *J*_{sun} = (4/3)(*F*_{0}/*T*_{sun}) by a blackbody radiation from the sun. A 4/3 factor appears in the expression for the entropy carried by radiation (denoted by *J*). Feistel (2011) shows that when entropy exchange between two black bodies is considered, (4/3)*σT*^{3} is the only mathematical expression for entropy carried by blackbody radiation to be consistent with the second law of thermodynamics. In addition, Yourgrau et al. (1982, sections 3–9) show, using a photon gas, that the 4/3 factor is necessary for the entropy carried by radiation. A simple example treating both processes, entropy carried by radiation and entropy produced by heating and cooling, is discussed in section 102 of Planck (1913).

*J*

_{ref}(Wu and Liu 2010). Emission of the longwave irradiance to space carries entropy,

*J*

_{atm}= (4/3)(

*F*

_{atm}/

*T*

_{atm}) and (1 −

*ε*)

*J*

_{sfc}= (1 −

*ε*)(4/3)(

*F*

_{sfc}/

*T*

_{sfc}). Therefore, entropy export to space by radiation is

*dS*/

*dt*= 0. Because of heating of ocean and land and melting ice and glacier, the global annual mean TOA net irradiance is estimated to be 0.71 ± 0.1 W m

^{−2}(Johnson et al. 2016; Loeb et al. 2018a). When entropy production due to a positive net TOA irradiance can be handled by adding entropy storage term in Eq. (14) and subsumed into entropy production within the system (Bannon and Lee 2017), then, at a steady sate, the difference of the entropy produced by longwave emission to space and shortwave absorption balances with entropy produced within the system is

*Q*

_{e}/

*T*

_{e}−

*Q*

_{a}/

*T*

_{a}equals

*Q*

_{e}/

*T*

_{e}−

*Q*

_{a}/

*T*

_{a}.

*dS*/

*dt*= 0, entropy production by nonradiative irreversible process is

## 3. Data products

We use Clouds and the Earth’s Radiant Energy System (CERES) Edition 4.1 SYN1deg-Month data product (Rutan et al. 2015a; Kato et al. 2018) in this study. The SYN1deg-Month product includes shortwave and longwave irradiance profiles computed hourly for every equal-area grid. The angle size of the equal-area grids is 1° × 1° between 45°N and 45°S and increases toward poles (Doelling et al. 2013). Cloud and aerosol properties are derived from Moderate Resolution Imaging Spectroradiometer (MODIS, collection 5 March 2000 through February 2016 and collection 6.1 from March 2016 onward for clouds and collection 6.1 throughout for aerosols) on *Terra* and *Aqua* except for the period from March 2000 through June 2002 when only MODIS on *Terra* is used. In addition, cloud properties are also derived from geostationary satellites hourly over regions between 60°N and 60°S to account for the diurnal cycle of clouds.

Temperature and humidity profiles used in producing SYN1deg-Month are taken from the Goddard Earth Observing System Data Assimilation System reanalysis product (Rienecker et al. 2008), version 5.4.1 (GEOS-5.4.1). Temperature and relative humidity profiles are 6 hourly in the equal-area grid. The temporal resolution of temperature and humidity profiles is increased to hourly below 10 m from the surface. For surface properties, we primarily use MODIS-derived spectral surface albedos except over ocean. Modeled spectral surface albedos are used for ocean (Jin et al. 2004) and other surface types when the MODIS spectral albedo is not available. Surface spectral emissivities are based on Wilber et al. (1999). Further details of inputs and algorithm to produce SYN1deg-Month is given in Rutan et al. (2015a) and Kato et al. (2018).

### Computation of entropy production

We compute shortwave and longwave upward and downward irradiances hourly at 35 levels in the atmosphere (plus cloud-top and cloud-base heights, depending on retrieved heights) for every equal-area grid with a radiative transfer model (Fu and Liou 1993; Rose et al. 2006, 2013). For shortwave irradiance computations, a gamma distribution is used to express the distribution of hourly cloud optical thickness in an equal-area grid separated by 4 different cloud types. We use two parameters to define a gamma distribution, the mean optical thickness and shape factor that is the mean optical thickness divided by the standard deviation of optical thicknesses. We infer the shape factor from linear and logarithmic means of cloud optical thickness (Wilks 1995). When the shape factor is less than 10, we use the algorithm that analytically integrates two-stream solutions of an integro-deferential equation of radiative transfer over a gamma distribution (Barker 1996; Kato et al. 2005). When the shape factor is greater than or equal to 10, we use a four-stream model (Liou et al. 1988). The four-stream model is also used for clear-sky irradiance computations. For longwave computations, we use a combination of two- and four-stream models (Fu et al. 1997).

*i*th level is defined as the downward irradiance minus upward irradiance, where the top of the atmosphere is defined as level = 0 and the level number increases downward. The top layer between level = 0 and level = 1 is layer 1 and the layer number also increases downward. The absorbed irradiance by the

*i*th layer

*F*

_{i,abs}is the net irradiance at

*i*− 1th level minus the net irradiance at

*i*th level,

*x*is either shortwave SW or longwave LW. The top-of-atmosphere shortwave + longwave net irradiance is therefore

*T*

_{i}is the temperature of the

*i*th layer,

*T*

_{sfc}is the surface skin temperature. To compute the entropy production by longwave irradiance emitted to space, Stephens and O’Brien (1993) divide the upward longwave irradiance at TOA by the effective emission temperature. In this study, we utilize the two-stream source function technique discussed in Toon et al. (1989) and the appendix. Toon et al. (1989) replace the radiance in computing the source function by the radiance computed with a two-stream approximation. With this approximation, the source function can be written explicitly. When the source function is specified, the azimuthally averaged radiance at a two-stream angle can also be written explicitly. We compute the downward irradiance starting from TOA moving downward toward the surface with this expression. We subsequently compute the upward irradiance starting at the surface moving upward toward TOA with the specified surface emissivity. We compute the entropy flux in the same way that the irradiance is computed except that the source function is divided by the temperature of the layer such that

*t*

_{i}is the transmission of the irradiance from

*i*th level to TOA,

*i*th level and

*T*

_{i}is the temperature of the

*i*th layer, and

*N*is the number of layers (~35). Once net irradiances are computed, entropy produced by nonradiative irreversible processes is computed by

## 4. Results

In this section, we discuss the result of global annual mean entropy production followed by spatial distribution of the production.

### a. Global annual mean entropy production

Table 1 summarizes computed global mean shortwave and longwave irradiances averaged from July 2005 through June 2015. Because computed TOA irradiances from SYN1deg-Month are not constrained by observed ocean heating rates averaged over a decade by the method discussed in Loeb et al. (2018a), the net TOA irradiance is 1.3 W m^{−2}, which is larger than the net TOA irradiance derived from CERES observations of 0.71 W m^{−2} (Johnson et al. 2016; Loeb et al. 2018a). Table 2 summarizes the annual global mean entropy production by shortwave and longwave irradiances computed from a 10-yr period (July 2005 through June 2015), as well as the standard deviation of annual global mean values. One year or a longer time period is needed to the system to be close to a steady state (Lucarini et al. 2011). The entropy produced by absorption of shortwave irradiance is 0.853 W m^{−2} K^{−1} and entropy produced by emission of longwave to space is −0.928 W m^{−2} K^{−1}. These give the entropy production by irreversible processes of 0.076 W m^{−2} K^{−1}. Entropy produced by longwave cooling in the atmosphere and surface is, respectively, −0.729 and −0.173 W m^{−2} K^{−1}. This gives the entropy produced by nonradiative irreversible process of 0.049 W m^{−2} K^{−1} [−(0.853–0.729–0.173)]. The difference between 0.076 and 0.049 W m^{−2} K^{−1} is 0.027 W m^{−2} K^{−1}, which is the entropy produced by longwave radiation exchange within the Earth system. Uncertainty in these estimates is difficult to quantify, but it is estimated based on uncertainties in irradiances estimated in earlier studies. Relative uncertainties in the global annual mean TOA shortwave and longwave irradiances are, respectively, 1% and 0.75% (Loeb et al. 2009). Based on uncertainties given in Kato et al. (2018), relative uncertainties in the global annual mean surface net shortwave and longwave irradiances are, respectively, 3% and 10%. If we assume that the uncertainty in global mean temperature is less than 1 K and ignore the correlation of errors in irradiances and temperatures, the uncertainty in the global annual mean entropy production estimate is about 10% (1*σ*). This uncertainty estimate also ignores the uncertainty in radiative heating rate and cloud vertical profiles discussed in Ham et al. (2017) and Kato et al. (2019).

Global annual mean and standard deviation of TOA and surface net irradiances (in parentheses) for July 2005 to June 2015 in W m^{−2}.

Global annual mean radiative entropy production in mW m^{−2} K^{−1} and standard deviation of annual means (in parentheses) computed from a 10-yr period.

For comparisons, entropy production estimated by Peixoto et al. (1991) is also included in Table 2. Our values are different from those estimated by Peixoto et al. (1991). The difference of entropy production by the longwave irradiance emitted to space is largely caused by partitioning the entropy production by the surface and atmosphere. However, entropy produced by shortwave absorption are significantly different, especially at the surface.

Nonradiative irreversible processes that increase entropy considered by Goody (2000) include turbulent enthalpy transport of 0.0024 W m^{−2} K^{−1}, frictional dissipation of 0.0113 W m^{−2} K^{−1}, and irreversible phase change, transport, and precipitation of water of 0.0188 W m^{−2} K^{−1}. Pauluis and Dias (2012) suggest that the atmospheric heating by the frictional dissipation of falling raindrops is nearly equivalent to the frictional dissipation of turbulence in the atmosphere. We therefore replace Goody’s estimate of entropy production by raindrops of 0.0013 m W m^{−2} K^{−2} by 0.010 W m^{−2} K^{−1}, which gives the frictional dissipation of 0.020 W^{−2} K^{−1}. Entropy produced by ocean estimated by Bannon and Najjar (2018) is 0.0016 W m^{−2} K^{−2}. When these productions are added, entropy production by nonradiative irreversible processes is 0.043 W m^{−2} K^{−1}. Our estimate of entropy produced by nonradiative irreversible process of 0.049 W m^{−2} K^{−1} agrees reasonably well to within 15% with the sum of the components. Also, our estimate agrees with the estimate of 0.055 W m^{−2} K^{−1} by Lucarini et al. (2011) using climate models, although their estimate is based on a simple expression [i.e. Eq. (23)] that ignores the vertical irradiance and temperature profiles in the atmosphere. The result of Lembo et al. (2019), who estimated entropy produced by nonradiative irreversible processes using CMIP5 models with the simple expression, shows a substantial spread among models and dependence on the method used for the estimate (i.e., direct and indirect methods). Nevertheless, our estimate is close to the mean value from seven CMIP5 models estimated by the direct method of 0.045 W m^{−2} K^{−2} (Lembo et al. 2019).

*σ*is the Stefan–Boltzmann constant. The annual global mean TOA upward longwave irradiance is 240 W m

^{−2}, which gives the entropy production of 0.941 W m

^{−2}K

^{−1}by Eq. (28), while the global annual mean entropy production from SYN1deg-Month is 0.928 W m

^{−2}K

^{−1}. The difference is presumably caused by the translucent atmosphere combined with decreasing temperature with height in the atmosphere. Longwave radiation emitted by the lower atmosphere transmits through the atmosphere in the wavelength regions where the atmosphere is not entirely opaque. Because of this difference, the effective emission temperature computed from dividing longwave irradiance emitted to space by longwave entropy production using SYN1deg-Month values is 259 K and larger than 255 K computed by Eq. (29) (Table 3). The result of Stephens and O’Brien (1993) also shows a smaller effective emission temperature when Eq. (29) is used. Therefore, the reason that our estimate shown in Table 2 is slightly larger than the estimate by Stephens and O’Brien (1993), who used Eq. (29), is probably due to calibration differences between ERBE and CERES instruments.

Global annual mean values that appear in Eq. (2) derived from SYN1deg-Month and the ideal energy balanced model.

^{−2}by entropy production by shortwave of 0.853 W m

^{−2}K

^{−1}is 282 K. While different definitions of the thermodynamic efficiency are proposed (e.g., Pauluis and Held 2002; Lucarini et al. 2011; Novak and Teilleux 2018), the efficiency that can be computed with these effective temperatures is the Carnot efficiency

*η*of the Earth system defined by Bannon and Lee (2017). Based on our results, the Carnot efficiency is

### b. Regional entropy production

Figure 2 shows regional annual mean entropy production by absorption of solar irradiance and by emission of longwave irradiance to space. Entropy produced by shortwave absorption is larger in the tropics and decreases toward poles. The spatial distribution of entropy produced by shortwave absorption resembles the spatial pattern of shortwave absorption, indicating that the spatial distribution is driven by shortwave absorption. The production is larger in the tropics and decreases toward poles. Absorption of shortwave irradiance produces more entropy at the surface than in the atmosphere. Similarly, entropy production by longwave emission is large over tropics and decreases toward poles. The effect of clouds is evident over the western Pacific, central Africa, and the Amazon. Similar to entropy production by shortwave, the spatial distribution is driven by the spatial distribution of emitted longwave irradiance to space. A larger contribution, however, comes from the atmosphere than from the surface. Once entropy production by shortwave absorption is subtracted from entropy production by longwave emission to space, the resulting entropy production increases toward poles (Fig. 3). This is equal to entropy produced by irreversible processes [Eq. (20)]. The regional difference of entropy produced by longwave emission to space and shortwave absorption is the net effect of entropy produced by frictional dissipation, entropy produced by diabatic heating by water phase change associated with water vapor transport (i.e., evaporation and condensation temperature differs), frictional dissipation of falling raindrops, horizontal and vertical enthalpy transport by irreversible processes, and energy transport by radiation exchange. As indicated by Eq. (13) and unlike regional energy balance in the column, vertical energy transports by irreversible processes and radiation exchange affect entropy production in the column.

Figure 4 shows the standard deviation of annual regional mean entropy produced by shortwave absorption and longwave emission to space. The variability of tropical convections causes large variability of both shortwave and longwave entropy production (Fig. 4 top). While both atmosphere and surface contribute the variability of longwave entropy production, the atmospheric contribution to the shortwave entropy production is smaller than the surface contribution. Clouds and water vapor variabilities are responsible for the variability of entropy production by longwave. The variability of clouds also causes the variability of entropy production by shortwave at the surface, but water vapor and aerosol variabilities primarily drive the variability of shortwave entropy production in the atmosphere. Once entropy production by shortwave absorption is subtracted from entropy production by longwave emission to space, the variability due to tropical convection is small (Fig. 5), which is analogous to a small net TOA irradiance over convective cloud regions due to partial cancellation of their shortwave and longwave radiative effects. As a result, large variabilities of the difference of entropy produced by shortwave absorption and longwave emission occur mostly over stratocumulus regions. Similar to Fig. 3, the variability shown in Fig. 5 is the result of variability of frictional dissipation, radiation exchange, and heating and cooling by water phase change in the column.

Figure 6 shows the mean and standard deviation of annual mean entropy produced by nonradiative irreversible process, which is equal to entropy produced by longwave cooling in the atmosphere and surface minus entropy produced by shortwave absorption in the atmosphere and surface [Eq. (23)]. The spatial pattern is similar to Figs. 3 and 5.

## 5. Discussion

At a steady state, the sum of entropy produced by absorption of solar irradiance and emission of longwave irradiance to space and entropy produced by irreversible processes in the Earth system is zero. In addition, the 1D model suggests that entropy produced by the net shortwave and longwave irradiances in the atmosphere and at the surface and entropy produced by nonradiative irreversible process is zero at a steady state. The nonradiative irreversible processes include frictional dissipation of turbulence, diabatic heating by water phase change associated with water vapor transport (i.e., evaporation and condensation temperatures differ), and frictional dissipation of falling raindrops. Although the 1D model provides insight into how entropy is produced, the 1D model can be misleading when it is used to understand the sensitivity of entropy production to energy input to Earth. Once observations are combined, however, the 1D model can provide an insight of how entropy production changes when energy input to the Earth system is changed. This point is demonstrated in the remaining of this section.

To understand the sensitivity of entropy production to energy input to Earth, we start with showing the relationship between global annual TOA emitted longwave irradiance and net shortwave irradiance anomalies derived from Edition 4.1 Energy Balanced and Filled (EBAF) product (Loeb et al. 2018a) in Fig. 7. The slope of the line derived by regressing net longwave irradiance anomalies versus shortwave irradiance anomalies is 0.39. The TOA emitted longwave irradiance, therefore, increases (decreases) with increasing (decreasing) shortwave absorption. At an annual global scale, however, the magnitude of emitted longwave irradiance anomalies is generally smaller than the magnitude of shortwave absorption anomalies. Standard deviations of global annual anomalies of the absorbed shortwave irradiance and emitted longwave irradiance are, respectively, 0.40 and 0.25 W m^{−2}. Similarly, when we plot anomalies of the annual global mean shortwave and longwave entropy production as a function of anomalies of the annual global mean TOA absorptivity *a* (TOA net shortwave irradiance divided by *F*_{0}) (Fig. 8), the slope is 0.948 W m^{−2} K^{−1} for ^{−2} K^{−1} for *T*_{a} or *T*_{e} is smaller than the relative increase of

At a global annual scale, most energy is used for heating ocean so that interannual variability of TOA net irradiance agrees with the variability of ocean heating rate (Johnson et al. 2016). When absorbed shortwave irradiance is perturbed, shortwave anomalies are damped by ocean so that the magnitude of emitted longwave anomalies is smaller than the magnitude of absorbed shortwave anomalies. These smoothing or damping processes seem to be responsible for a negative ^{−2} K^{−1}) derived by a linear regression applied to the *a* relationship (Fig. 8) is consistent with the difference of the slopes of

*T*

_{a}with shortwave absorption so that observations are needed to estimate ∂

*T*

_{a}/∂

*a*. If we set ∂

*T*

_{a}/∂

*a*= 0, the derivative

^{−2}K

^{−1}. The estimate of

^{−2}K

^{−1}. The reason for the difference between our value of −0.201 W m

^{−2}K

^{−1}and their value is that we use

*T*

_{a}= 283 K instead of 275 K. When we use the value of ∂

*T*

_{a}/∂

*a*estimated from the CERES data products, the derivative is 0.172 W m

^{−2}K

^{−1}. The difference of the sensitivity derived from the CERES data products and 1D model arises because of the difference of

We test whether Eq. (31) can also be used to estimate the critical absorptivity to make ^{−2} K^{−1} derived by directly fitting *a*, and use *T*_{a}, and ∂*T*_{a}/∂*a* derived from the CERES data products (the second column of Table 3), *a* ≤ 1. Therefore, entropy production by irreversible processes always decreases with increasing absorption of shortwave irradiance for the possible range of the planetary albedo.

## 6. Summary and conclusions

We used satellite observations of cloud and aerosol properties, and temperature and humidity profiles from reanalysis to compute entropy production within the Earth system. We focused entropy production by heating and cooling by radiation that, at a steady state, balances with entropy production by irreversible processes within the system.

Entropy produced by longwave emission to space estimated in this study, 0.929 W m^{−2} K^{−1}, is close to that of 0.925 W m^{−2} K^{−1} estimated by Peixoto et al. (1991). In addition, entropy carried by longwave irradiance emitted to space estimated in this study is in a close agreement with that estimated by Stephens and O’Brien (1993) using ERBE data. However, entropy produced by absorption of shortwave irradiance at the surface and by the atmosphere are significantly different from those estimated by Peixoto et al. (1991). With a steady-state assumption, entropy production by irreversible processes within the Earth system is estimated to be 0.076 W m^{−2} K^{−1} and by nonradiative irreversible processes to be 0.049 W m^{−2} K^{−1}.

Spatial distribution of entropy production is similar to the spatial distribution of net shortwave irradiance and longwave irradiance emitted to space. The larger standard deviation of annual mean of regional entropy productions by shortwave irradiance absorption and longwave irradiance emission to space is caused by the variability due to convections in the tropics. Once the entropy production by absorption of shortwave irradiance is subtracted from the entropy production by longwave irradiance emission to space, the variability of the difference is larger over stratocumulus regions, similar to the spatial distribution of the variability of net TOA irradiances.

Entropic absorption and emission temperatures derived from annual global mean TOA shortwave and longwave irradiance and entropy productions are, respectively, 282 and 259 K. These entropic temperature gives the Carnot efficiency of the Earth system of 8.5%.

CERES observations suggest that, at a global annual scale, entropy produced by emission of longwave irradiance to space and by absorption of shortwave irradiance increases with increasing energy input to (absorptivity of) the Earth system. Increasing the entropy produced by the shortwave absorption is larger than increasing the entropy produced by the longwave emission, which makes entropy production by irreversible processes decrease when absorption of shortwave irradiance increases.

## Acknowledgments

We thank two anonymous reviewers who provided constructive comments that improved the manuscript. This research was supported by the NASA CERES project. The Ed4.1 CERES EBAF–TOA (Loeb et al. 2018b) and SYN1deg-Month (Rutan et al. 2015b) datasets were downloaded from https://ceres.larc.nasa.gov/order_data.php.

## APPENDIX

### Computation of Material Entropy Production by Longwave Irradiance Emitted to Space

*B*

_{0n}is the Planck function evaluated at the temperature of the top of the

*n*th computation layer with the optical thickness of

*τ*

_{n}and

*B*(

*T*

_{bn}) is the Planck function evaluated at the temperature of the bottom of the

*n*th layer

*T*

_{bn}. Second, we use the radiance computed with a two-stream approximation instead of the true radiance in computing source function. With these approximations, the source function

*S*

_{νt}of longwave irradiances propagating upward and downward though the atmosphere are [Toon et al. 1989, their Eqs. (53) and (54)], respectively,

*G*,

*H*,

*J*,

*K*,

*α*, and

*σ*are defined in Table 3 of Toon et al. (1989) for the hemispheric mean two-stream approximation and superscripts + and − indicate upward (propagating from the surface to TOA) and downward (propagating from TOA to the surface). We also assume that the temperature changes linearly within the

*n*th computational layer with the optical thickness of

*τ*

_{n}so that

*T*

_{bn}and

*T*

_{0n}are the temperature at the bottom and top of the

*n*th layer. The entropy produced by the upward longwave irradiance by the

*n*th layer is

*C*

^{±}are replaced by

*ζ*,

*η*, and

*β*are defined in Fu et al. (1997). We only keep terms with up to second power of

*T*

_{0n}in the denominator because high-order terms are smaller in computing entropy production by Eqs. (A18) and (A19). These expressions compute entropy produced by longwave emission by the layer. These expressions therefore accurately compute entropy produced by longwave absorption and emission by the layer. These expressions, however, do not explicitly compute the transmission from each layer to space and hence do not explicitly keep track of the emission temperature of radiation emitted to space.

## REFERENCES

Bannon, P. R., 2015: Entropy production and climate efficiency.

,*J. Atmos. Sci.***72**, 3268–3280, https://doi.org/10.1175/JAS-D-14-0361.1.Bannon, P. R., and S. Lee, 2017: Toward quantifying the climate heat engine: Solar absorption and terrestrial emission temperatures and material entropy production.

,*J. Atmos. Sci.***74**, 1721–1734, https://doi.org/10.1175/JAS-D-16-0240.1.Bannon, P. R., and R. G. Najjar, 2018: Heat-engine and entropy-production analyses of the world ocean.

,*J. Geophys. Res. Oceans***123**, 8532–8547, https://doi.org/10.1029/2018JC014261.Barker, H. W., 1996: A parameterization for computing grid-averaged solar fluxes for inhomogeneous marine boundary layer clouds. Part I: Methodology and homogeneous biases.

,*J. Atmos. Sci.***53**, 2289–2303, https://doi.org/10.1175/1520-0469(1996)053<2289:APFCGA>2.0.CO;2.Doelling, D. R., N. G. Loeb, D. F. Keyes, M. L. Nordeen, D. Morstad, B. A. Wielicki, D. F. Young, and M. Sun, 2013: Geostationary enhanced temporal interpolation for CERES flux products.

,*J. Atmos. Oceanic Technol.***30**, 1072–1090, https://doi.org/10.1175/JTECH-D-12-00136.1.Feistel, R., 2011: Entropy flux and entropy production of stationary black-body radiation.

,*J. Non-Equilib. Thermodyn.***36**, 131–139, https://doi.org/10.1515/jnetdy.2011.009.Fu, Q., and K. N. Liou, 1993: Parameterization of the radiative properties of cirrus clouds.

,*J. Atmos. Sci.***50**, 2008–2025, https://doi.org/10.1175/1520-0469(1993)050<2008:POTRPO>2.0.CO;2.Fu, Q., K. N. Liou, M. C. Cribb, T. P. Charlock, and A. Grossman, 1997: Multiple scattering parameterization in thermal infrared radiative transfer.

,*J. Atmos. Sci.***54**, 2799–2812, https://doi.org/10.1175/1520-0469(1997)054<2799:MSPITI>2.0.CO;2.Goody, R., 2000: Sources and sinks of climate entropy.

,*Quart. J. Roy. Meteor. Soc.***126**, 1953–1970, https://doi.org/10.1002/qj.49712656619.Ham, S.-H., S. Kato, H. W. Barker, F. G. Rose, and S. Sun-Mack, 2014: Effects of 3-D clouds on atmospheric transmission of solar radiation: Cloud type dependencies inferred from A-train satellite data.

,*J. Geophys. Res. Atmos.***119**, 943–963, https://doi.org/10.1002/2013JD020683.Ham, S.-H., and et al. , 2017: Cloud occurrences and cloud radiative effects (CREs) from CERES-CALIPSO-CloudSat-MODIS (CCCM) and CloudSat radar-lidar (RL) products.

,*J. Geophys. Res. Atmos.***122**, 8852–8884, https://doi.org/10.1002/2017JD026725.Jin, Z., T. P. Charlock, W. L. Smith Jr., and K. Rutledge, 2004: A parameterization ocean surface albedo.

,*Geophys. Res. Lett.***31**, L22301, https://doi.org/10.1029/2004GL021180.Johnson, G. C., J. M. Lyman, and N. G. Loeb, 2016: Improving estimates of Earth’s energy imbalance.

,*Nat. Climate Change***6**, 639–640, https://doi.org/10.1038/NCLIMATE3043.Kato, S., F. G. Rose, and T. P. Charlock, 2005: Computation of domain averaged irradiance using satellite-derived cloud properties.

,*J. Atmos. Oceanic Technol.***22**, 146–164, https://doi.org/10.1175/JTECH-1694.1.Kato, S., and et al. , 2018: Surface irradiances of Edition 4.0 Clouds and the Earth’s Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) data product.

,*J. Climate***31**, 895–918, https://doi.org/10.1175/JCLI-D-17-0208.1.Kato, S., and et al. , 2019: Radiative heating rates computed with clouds derived from satellite-based passive and active sensors and their effects on generation of available potential energy.

,*J. Geophys. Res. Atmos.***124**, 1720–1740, https://doi.org/10.1029/2018JD028878.Kjellsson, J., 2015: Weakening of the global atmospheric circulation with global warming.

,*Climate Dyn.***45**, 975–988, https://doi.org/10.1007/s00382-014-2337-8.Laliberté, F., J. Zika, L. Mudryk, P. J. Kushner, J. Kjellsson, and K. Döös, 2015: Constrained work output of the moist atmospheric heat engine in a warming climate.

,*Science***347**, 540–543, https://doi.org/10.1126/science.1257103.Lembo, V., F. Lunkeit, and V. Lucarini, 2019: TheDiaTo (v1.0)—A new diagnostic tool for water, energy and entropy budgets in climate models.

,*Geosci. Model Dev.***12**, 3805–3834, https://doi.org/10.5194/gmd-12-3805-2019.Liou, K.-N., Q. Fu, and T. P. Ackerman, 1988: A simple formulation of the delta-four-stream approximation for radiative transfer parameterizations.

,*J. Atmos. Sci.***45**, 1940–1948, https://doi.org/10.1175/1520-0469(1988)045<1940:ASFOTD>2.0.CO;2.Loeb, N. G., B. A. Wielicki, D. R. Doelling, G. L. Smith, D. F. Keyes, S. Kato, N. Manalo-Smith, and T. Wong, 2009: Toward optimal closure of the Earth’s top-of-atmosphere radiation budget.

,*J. Climate***22**, 748–766, https://doi.org/10.1175/2008JCLI2637.1.Loeb, N. G., and et al. , 2018a: Clouds and the Earth’s Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) top-of-atmosphere (TOA) Edition-4.0 data product.

,*J. Climate***31**, 895–918, https://doi.org/10.1175/JCLI-D-17-0208.1.Loeb, N. G., and et al. , 2018b, EBAF Edition 4.1. CERES, accessed April 2019, https://ceres-tool.larc.nasa.gov/ord-tool/jsp/EBAFTOA41Selection.jsp.

Lucarini, V., K. Fraedrich, and F. Ragone, 2011: New results on the thermodynamic properties of the climate system.

,*J. Atmos. Sci.***68**, 2438–2458, https://doi.org/10.1175/2011JAS3713.1.Lucarini, V., R. Blender, C. Herbert, F. Ragone, S. Pascale, and J. Wouters, 2014: Mathematical and physical ideas for climate science.

,*Rev. Geophys.***52**, 809–859, https://doi.org/10.1002/2013RG000446.Novak, L., and R. Tailleux, 2018: On the local view of atmospheric available potential energy.

,*J. Atmos. Sci.***75**, 1891–1907, https://doi.org/10.1175/JAS-D-17-0330.1.Pauluis, O. and I. M. Held, 2002: Entropy budget of an atmosphere in radiative–convective equilibrium. Part I: Maximum work and frictional dissipation.

,*J. Atmos. Sci.***59**, 125–139, https://doi.org/10.1175/1520-0469(2002)059<0125:EBOAAI>2.0.CO;2.Pauluis, O., and J. Dias, 2012: Satellite estimates of precipitation-induced dissipation in the atmosphere.

,*Science***335**, 953–956, https://doi.org/10.1126/science.1215869.Peixoto, J. P., and A. H. Oort, 1992:

*Physics of Climate*. American Institute of Physics, 520 pp.Peixoto, J. P., A. H. Oort, M. De Almeida, and A. Tome, 1991: Entropy budget of the atmosphere.

,*J. Geophys. Res.***96**, 10 981–10 988, https://doi.org/10.1029/91JD00721.Planck, M., 1913:

*The Theory of Heat Radiation*. History of Modern Physics and Astronomy, Vol. 11, American Institute of Physics, 256 pp.Rienecker, M. M., and et al. , 2008: The GEOS-5 Data Assimilation System—Documentation of versions 5.0.1, 5.1.0, and 5.2.0. NASA Tech. Memo. NASA/TM-2008-105606, Vol. 27, 97 pp, http://gmao.gsfc.nasa.gov/pubs/docs/Rienecker369.pdf.

Rose, F., T. Charlock, Q. Fu, S. Kato, D. Rutan, and Z. Jin, 2006: CERES Proto-Edition 3 radiative transfer: Tests and radiative closure over surface validation sites.

*12th Conf. on Atmospheric Radiation*, Madison, WI, Amer. Meteor. Soc., P2.4, https://ams.confex.com/ams/pdfpapers/112358.pdf.Rose, F., D. A. Rutan, T. P. Charlock, G. L. Smith, and S. Kato, 2013: An algorithm for the constraining of radiative transfer calculations to CERES-observed broadband top-of-atmosphere irradiance.

,*J. Atmos. Oceanic Technol.***30**, 1091–1106, https://doi.org/10.1175/JTECH-D-12-00058.1.Rutan, D. A., S. Kato, D. R. Doelling, F. G. Rose, L. T. Nguyen, T. E. Caldwell, and N. G. Loeb, 2015a: CERES synoptic product: Methodology and validation of surface radiant flux.

,*J. Atmos. Oceanic Technol.***32**, 1121–1143, https://doi.org/10.1175/JTECH-D-14-00165.1.Rutan, D. A., and et al. , 2015b, SYN1-deg-Month Edition 4. CERES, accessed April 2019, https://ceres-tool.larc.nasa.gov/ord-tool/jsp/SYN1degEd41Selection.jsp.

Stephens, G. L., and D. M. O’Brien, 1993: Entropy and climate. I: ERBE observations of the entropy production of the earth.

,*Quart. J. Roy. Meteor. Soc.***119**, 121–152, https://doi.org/10.1002/qj.49711950906.Stephens, G. L., and et al. , 2012: An update on Earth’s energy balance in light of the latest global observations.

,*Nat. Geosci.***5**, 691–696, https://doi.org/10.1038/ngeo1580.Toon, O. B., C. P. Mckay, and T. P. Ackerman, 1989: Rapid calculation of radiative heating rates and photodissociation rates in inhomogeneous multiple scattering atmospheres.

,*J. Geophys. Res.***94**, 16 287–16 301, https://doi.org/10.1029/JD094iD13p16287.Wilber, A. C., D. P. Kratz, and S. K. Gupta, 1999: Surface emissivity maps for use in satellite retrievals of longwave radiation. NASA Tech. Memo. TP-1999-209362, 30 pp., https://eosweb.larc.nasa.gov/sites/default/files/project/calipso/Wilber.NASATchNote99.pdf.

Wilks, D. S., 1995:

*Statistical Methods in the Atmospheric Sciences*, Academic Press, 467 pp.Wu, W., and Y. Liu, 2010: Radiation entropy flux and entropy production of the Earth system.

,*Rev. Geophys.***48**, RG2003, https://doi.org/10.1029/2008RG000275.Yourgrau, W., A. van der Merwe, and G. Raw, 1982: Treatise on Irreversible and Statistical Thermophysics: An Introduction to Nonclassical Thermodynamics. Dover Publishing, 268 pp.