The Iris Effect: A Review

This study reviews the research of the past 20-years on the role of anvil cirrus in the Earth’s climate – research initiated by Lindzen et al. (Bull. Am. Meteor. Soc. 82:417-432, 2001). The original study suggested that the anvil cirrus would shrink with warming, which was estimated to induce longwave cooling for the Earth. This is referred to as the iris effect since the areal change hypothetically resembles the light control by the human eye’s iris. If the effect is strong enough, it exerts a significant negative climate feedback which stabilizes tropical temperatures and limits climate sensitivity. Initial responses to Lindzen et al. (Bull. Am. Meteor. Soc. 82:417-432, 2001) denied the existence and effectiveness of the iris effect. Assessment of the debatable issues in these responses will be presented later in this review paper. At this point, the strong areal reduction of cirrus with warming appears very clearly in both climate models and satellite observations. Current studies found that the iris effect may not only come from the decreased cirrus outflow due to increased precipitation efficiency, but also from concentration of cumulus cores over warmer areas (the so-called aggregation effect). Yet, different opinions remain as to the radiative effect of cirrus clouds participating in the iris effect. For the iris effect to be most important, it must involve cirrus clouds that are not as opaque for visible radiation as they are for infrared radiation. However, current climate models often simulate cirrus clouds that are opaque in both visible and infrared radiation. This issue requires thorough examination as it seems to be opposed to conventional wisdom based on explicit observations. This paper was written in the hope of stimulating more effort to carefully evaluate these important issues.


Introduction
Cirrus refers to the upper-level thin clouds that consist almost entirely of ice particles. As saturation vapor pressure of ice is low, particles of cirrus clouds evaporate slowly and are readily advected over long distances by strong upper-level winds (Houze Jr. 2014). For this reason, cirrus cover large areas of the Earth with daily variability on the order of millions of square meters. Very high-altitude cirrus whose base is above 10 km are concentrated within the western tropical Pacific (Mace et al. 2009). Of those cirrus, about half appear to originate from the upper parts of cumulonimbus (anvil cirrus) and the other half appear to form in situ due to the large-scale dynamics (Luo and Rossow 2004;Reverdy et al. 2012). Such anvil cirrus clouds play a key role in moistening the upperlevel air that otherwise would have been extremely dry (the water vapor above 5 km is only 10% of the total column water amount on the average). Consequently, the upper tropospheric anvil cirrus and moist area in the vicinity of clouds may reduce outgoing longwave (or infrared) radiation escaping from the Earth. The anvil cirrus, especially the optically thick part near the cumulonimbus, also blocks incoming shortwave (or visible) radiation (Kiehl 1994). Such anvil cirrus play an important role in the Earth's radiative budget.
As tropical deep convective activities are highly affected by surface temperature, anvil cirrus from the deep convective system would also be responsive to surface temperature (a cloud feedback). This is being found quite clearly with active satellite observations such as CloudSat (Igel et al. 2014), but it came as a surprise when it was first suggested by a paper titled "Does the Earth have an adaptive infrared iris?" by Lindzen et al. (2001) (hereafter LCH01). LCH01 hypothesized that the Responsible Editor: Hyemi Kim. region of anvil cirrus could control the tropical outgoing energy by reducing their areal coverage with warming. By the hypothesis, the counteractive effect on the incoming energy was considered as well, but it turned out to be secondary. It is this behavior that was referred to as the iris effect. LCH01 observationally tested the iris effect with geostationary infrared imagery over the tropical western Pacific. Although several assumptions are required for the estimation of the impact on climate sensitivity, LCH01 argued that the iris effect would be of fundamental importance in climate projectiona consequence of the combination of several radiative feedbacks associated with cloud, water vapor, etc. This paper reviews the research in the 20-years since LCH01.
LCH01 showed that it was no longer possible to disentangle the radiative impact of changes in water vapor from that due to changes in upper-level cloud cover, especially since the upper-level cirrus, when present, rather than the water vapor, determine the infrared emission level. Instead, it is necessary when dealing with feedbacks, to consider a longwave feedback that incorporates the combination of cloud, water vapor and areal responses to temperature. This has led subsequent studies to examine the change in the top-of-the-atmosphere longwave radiation instead of cloud and water vapor individually. As we will show later, many studies found negative longwave feedback consistent with the iris effect Choi 2009, 2011;Choi et al. 2014;Mauritsen and Stevens 2015). Nevertheless, it is a subject of controversy to what extent the shortwave feedback compensates for the longwave feedback; some studies also assert that the compensation would give nearly zero net radiative feedback in the iris effect (Sherwood et al. 2020).
Recently, papers have investigated sources other than what was considered by LCH01 for the origin of the iris effect. The original idea of LCH01 was that the iris effect is associated with decreased cirrus outflow due to the increased precipitation efficiency with warming. However, subsequent studies showed the iris effect may also be the manifestation of additional atmospheric processes. Wing and Emanuel (2014) and Muller and Bony (2015) suggested that something like the iris effect (i.e., the shrinkage or tardy response of anvil cirrus area with warming) clearly appears in the increased convective aggregation with the increased temperature. Similarly, the iris effect was observed as the manifestation of the concentration of cumulonimbus over warmer regions Saint-Lu et al. 2020), and the strengthening of Hadley circulation (Su et al. 2017). Bony et al. (2016) argued that the convective anvil shrinkage is rooted in basic energetics and thermodynamic properties of the atmosphere. They said the clouds find themselves in a more stable atmosphere which reduces the convective outflow in the upper troposphere. These mechanisms, while different from LCH01 in terms of the origin, reconfirm the existence of the iris effect. The present review will consider all these studies that seemingly imply the iris effect.
In order to understand the iris effect, it may be useful to review how the identification of the iris effect came to be. To this end, we begin with the historical background in Section 2. Section 3 will deal with the comments claiming to have refuted the results of LCH01. Section 4 will discuss where things appear to stand at the moment. Section 5 points out that the current climate feedback study is confronted with signal-tonoise problems. Section 6 will briefly describe some areas in particular need of more work. Some of these things have been dealt with earlier (Lindzen 2012). As to the original idea of the iris effect, the matter is treated in greater detail in the original paper (LCH01).

Historical Background of the Iris Effect
Most people working on climate initially, assumed that relative humidity remained fixed as climate changed, and that, therefore, following Manabe and Wetherald (1975), there was a water vapor feedback (Lindzen et al. 1982). However, Lindzen (1990), in response to a review, had suggested the possibility of a mechanism for reducing the water vapor feedback. Subsequent calculations showed the proposed mechanism would be ineffective. A careful study of the water vapor budget was begun leading to the following paper (Sun and Lindzen 1993). Microwave observations of water vapor by Spencer and Braswell (1997) made it clear that there was an additional factor that was going to be crucial: namely the fact that above the boundary layer, there was a clear delineation between moist and dry regionsa result, possibly, of an instability in radiative-convective equilibrium (Nilsson and Emanuel 1999), and that the behavior of any water vapor feedback would have to take account of the relative areas of these regions. These findings motivated LCH01.
The original emphasis in LCH01 was, indeed, on the behavior of the water vapor, but an important practical consideration forced LCH01 to focus on cloud cover. Namely, in order to observe the detailed response of cloudy area to surface temperature, LCH01 had to use geostationary data; the available satellite (Japanese GMS5) measured brightness temperature which could be used to estimate cloud cover, but could not measure water vapor. It was also unable to detect thin cirrus. However, LCH01 used measurements of thicker cirrus as surrogates for total cirrus. Their first instinct was to use the relation found by Udelhofen and Hartmann (1995) to assume a nearly linear positive relation between high cloud and humidity. However, in the course of the calculations, LCH01 realized that the clouds themselves were a major factor in the longwave budget, and allowed for water vapor and clouds to be unrelated. This, of course, minimized negative feedbacks (i.e., feedbacks that diminished rather than amplified the response to radiative forcing; conversely, positive feedbacks amplify the response.), which might otherwise be overestimated by a coherent relation of humidity and clouds. Nevertheless, the feedbacks were negative because it was found that high thin cirrus cloud cover, when normalized by a measure of cumulonimbus mass flux, diminished with increases in temperature allowing enhanced cooling. Tropical cirrus originate from convective cores. The normalization was crucial to obtain cirrus cloudiness per a unit convective core; without normalization, the results would simply reflect the concentration of cumulus activity over warm regions .
Much of the rest of LCH01 was devoted to the description and implementation of simple calculations to infer the impact of the iris effect on the climate sensitivity. The limitations of LCH01 were carefully described. Although the iris effect was naturally restricted to the tropics where cumulonimbus convection is a dominant process, the radiative impact of the iris effect was averaged over the globe in order to estimate the impact on climate sensitivity. Again, the details are given in LCH01. The results suggested a very strong negative feedback that would cancel the positive feedbacks in current climate models. LCH01 also showed that the few climate models, then available, differed sufficiently from the observations to overwhelm the uncertainties in the observations. The final suggestion of LCH01 was that modelers should incorporate the iris effect into climate models to see how it would impact results. It took about 14 years for Mauritsen and Stevens (2015) to take up this modest suggestion. In the interim, LCH01 was subjected to harsh criticismsespecially between 2002 and 2005. These will be described in the next section.

Subsequent Criticisms and Rebuttals
The first response to LCH01 came as a separate article by Hartmann and Michelsen (2002a). The title of the article was "No Evidence for Iris" and the new editor of the Bulletin of the American Meteorological Society appended a subtitle "Careful analysis of data reveals no shrinkage of tropical cloud anvil area with increasing sea surface temperature (SST).". The basic assertion of Hartmann and Michelsen (2002a) was that the stratiform clouds LCH01 was looking at did not originate in cumulonimbus convection. Rather, the paper claimed that clouds are associated with extratropical systems that intrude into the tropics. In LCH01, however, SST was weighted by cloud cover since clouds would not respond to SST over the whole domain, but SST beneath the cloud. As shown by Lindzen and Christy (2020), global mean temperatures bear little relation to local temperatures.
In fact, the assertion of Hartmann and Michelsen (2002a) was grounded by correlation analyses between cloud amount for 'each grid' and cloud weighted SST averaged 'over the whole domain'. The pitfall of this grid-based analysis will be shown in contrast with the area-averaging analysis in the next section. What the paper showed was a strong negative correlation between cloudiness and SST over some subtropical regions, but no such relation near the equator, was associated with strong equatorward winds from the exterior of the tropical western Pacific (TWP; 120°E-155°W, 30°S-30°N) studied in LCH01. Lindzen et al. (2002) pointed out that the analysis by Hartmann and Michelsen (2002a) is not appropriate to test the iris effect since the way they analyzed the data cannot capture clouds popping up in different locations over the TWP and moving to the next grid. Hartmann and Michelsen (2002b), in their reply to Lindzen et al. (2002), agreed that "the specific meteorological origin of the cloud anomalies is less relevant" and admitted their result is consistent with Lindzen et al. (2002), suggesting "the deep convective core clouds do exist outside the near equatorial region and over colder SST". However, Hartmann and Michelsen (2002b) reasserted that the data do not support the iris effect due to too high sensitivity of the statistics to cloud weighted SST.
In any event, the Hartmann-Michelsen assertion was that a tropical-subtropical contrast distorts climate sensitivity. Their paper would remain doubtful for at least two reasons: (i) The iris effect would diminish as we restricted ourselves to regions closer to the equator if Michelsen (2002a, 2002b) were correct. If anything, the effect increased. (ii) Close examination of GMS5 data permitted one to identify cumulonimbus cores associated with cirrus in all parts of the domain.
Del Genio and Kovari (2002) noted that when they looked at a warmer region and a cooler region in the tropics, they found more upper-level cloud cover in the warmer region in contrast to what the iris effect would suggest. However, this simply reflected the general phenomenon that "convection preferentially occurs over the warmest surface temperatures" (Mauritsen and Stevens 2015). Thus, it may easily appear as if upper-level cloud cover increases with temperature. Unfortunately, they failed to distinguish cloud cover from cloud cover per unit cumulus tower (area with 11-μm brightness temperature less than 220 K in LCH01). The iris effect refers to the latter. As Rondanelli and Lindzen (2008) noted later, when the normalization was applied, they actually found a stronger iris effect than LCH01 had estimated. As we will explain later, the movement of convection to the warmest surface will eventually reduce the actual area of anvil cirrus for the given TWP domain.
Two other critical papers (Fu et al. (2002), and Lin et al. (2002)) argued against the LCH01's treatment of the radiative effects of the upper level cirrus. As noted in the responses (Chou et al. 2002a(Chou et al. , 2002b, both these criticisms argued that LCH01 had misrepresented the radiative properties of upper level cirrus by overemphasizing their infrared effect and understating their visible effect. Lin et al. (2002) concluded that the main impact of upper level clouds was in the visible.
The disparate findings may arise from small differences in albedo over tropical moist and dry regions: 0.27 and 0.21 for LCH01 versus 0.31 and 0.15 for Lin et al. (2002), and more importantly, from different fractional coverages of the tropical high cloud within a moist region: 0.44 for LCH01 versus 0.2 for Lin et al. (2002). Note that dry regions do not only include clear-sky albedo (0.13), but also low-level cumulus or stratocumulus that have a high albedo (0.42). LCH01 were taking Ac(260) (i.e., area of 11-μm brightness temperature < 260 K 1 ) as a surrogate for all tropical upper level cirrus. Later, Choi et al. (2005) show that Ac (260) takes just about half of the MODIS high-cloud amount determined by 1.38-mm reflectances (>1.1%) since the minimum detectable value of optical depths by the GMS5 and MODIS measurements is around 2.0 and 0.1, respectively. Therefore, much of very thin cirrus clouds were missing in Ac(260) and the high cloud fraction of 0.44 for LCH01 may not be too large (Mace et al. 2009). For the very thin cirrus, thicker cirrus with and without underlying low clouds (which total-column optical depth <~9), the infrared effect was found to dominate over the visible effect in satellite observations from CERES and MODIS (Choi and Ho 2006). Likewise, virtually all papers dealing with the impact of clouds on climate agree that for upper level cirrus, the primary radiative impact is in the infrared (e.g., Liou 1986). Thus, LCH01 focused on the clouds in the immediate vicinity of the cumulus cores which are atypically thick and do have a high albedo. LCH01 have already noted that separating thick clouds from thinner clouds did not change the magnitude of the effect. Neither did changing the reference brightness temperature. In brief, these criticisms managed to change the results obtained by LCH01 by claiming that LCH01 had followed an inappropriate procedure different from the one actually employed.
Perhaps the most misleading criticism was that due to Lin et al. (2004). As noted by Chou and Lindzen (2005), Lin et al. (2004) assumed differences from LCH01 imply the absence of negative feedbacks, but the negative feedbacks implied by their results are larger than those in LCH01. Similarly, they assumed that the discrepancies between their specification of cloud radiative properties and the observations is due to their use of the iris result that cloud cover diminishes 22% per degree. What they failed to note was that to bring their values into line with observations, they would need much more than 22%.
Several studies simply could not find direct evidence supporting the iris effect with currently available satellite data: direct energy observing CERES data , and TRMM data (Lin et al. 2006;Rapp et al. 2005). TRMM data were expected to offer precipitation efficiency in the deep convective cloud system, so as to allow exploration of the mechanism behind the iris effect. However, as Su et al. (2008) noted, cloud information derived from TRMM data may not capture the upper tropospheric thin cirrus clouds. Their suggestion was that the iris effect should also take into account radiative effects of upper tropospheric cloud water content variation related to cloud optical thickness, in addition to cirrus fraction. Likewise, Kubar and Jiang (2019) used the 2-year MISR data combined with the other satellite data to show overall thinning and top-heaviness of clouds with increasing SST over TWP, but failed to find the iris effect. As shown in Igel et al. (2014) and Choi et al. (2017), however, recent satellite observations very likely do contain the iris effect.
Though several rebuttal points were briefly reviewed here, for further information, readers may look over Lindzen et al. (2002) Su et al. (2008). The results of these comment papers are not again refuted in the literature since then. In the next section, subsequent papers viewed somewhat different, but relevant points, leading to the gradual evolution of the argument for the iris effect beyond the original argument by LCH01. Rondanelli and Lindzen (2008;2010a) used the ground-based radar at Kwajalein as well as Tropical Rainfall Measuring Mission (TRMM) satellite data to examine aspects of the iris effect. In particular, the iris effect implied that precipitation efficiency in cumulonimbus towers would increase with increasing SST, and that the area of stratiform detrainment from these towers would decrease. The Kwajalein ground radar offered good resolution and coverage. As usual, one needed to be careful about adequate time intervals. At the least, one had to avoid sampling only one part of the time evolution of the cloud system. However, they found that at any given time the radar was sensing a large number of systems at different stages of evolution, and this minimized the issue. Rondanelli and Lindzen (2008) found that convective precipitation increased 10%/K while stratiform area decreased~22%/Kconsistent with the earlier results from LCH01. TRMM data provided coverage over the whole tropics, but with poor time resolution and other problems: perhaps most notably an inability to detect low levels of precipitation. Still, the TRMM data also showed increases of convective precipitation (~7%/K) and decreases of stratiform area (~5%/K) albeit at a reduced level (though TRMM showed great regional variability with some regions showing much larger values). Thus, a variety of independent approaches (Del Genio and Kovari 2002;Horváth and Soden 2008;LCH01;Rondanelli and Lindzen 2008;Rondanelli and Lindzen 2010a;Igel et al. 2014;Igel and van den Heever 2015;Choi et al. 2017;Saint-Lu et al. 2020) have all confirmed the iris effect, though there is uncertainty as to the exact magnitude of the effect (Sherwood et al. 2020). As already noted, Del Genio and Kovari (2002) required normalization to show the effect. The reason for increased precipitation efficiency with increasing SST is, of course, not resolved by these studies.

Advancing the Iris Effect Study
The issue of normalization warrants some discussion. The normalization used by Rondanelli and Lindzen (2008) was total precipitation. LCH01 recommended cumulus mass flux. Whatever the choice is, as noted by Rondanelli and Lindzen (2008), "the normalization is a necessary condition to obtain a meaningful result since to the first order the SST distribution organizes the spatial variation of convection in the tropics. Cloud properties such as the area of detrainment, total rainfall, and cloud radiative effects can be considered extensive with respect to the amount of convection. The normalization simply takes into account the fact that the amount of convection in a region of a given temperature is not indicative of the amount of convection in a climate with the same temperature, since the global amount of convection is determined by global energy balance considerations. It remains an open question whether total precipitation, mass convective flux, or some other measure of convection is the most appropriate normalization factor with regard to climate effects." LCH01 noted that if relative humidity in the tropical boundary layer remained constant as climate changed, the cumulus mass flux was the appropriate normalization. However, small changes in relative humidity would lead to different choices, and until we have a good basis for predicting such changes, the exact choice of normalization remains uncertain. There is no guarantee that the normalized results will provide a reliable quantitative climate effect (Su et al. 2008) without relating it to radiative effects, but not normalizing the cloud variations by a measure of convection will surely provide a wrong answer (Rondanelli and Lindzen 2010a).
As an application of the iris effect, Rondanelli and Lindzen (2010b) considered its possible role in resolving the Early Faint Sun Paradox (Sagan and Mullen 1972). Briefly, Sagan and Mullen (1972) noted that according to the standard model of the sun, solar output should have been 20-30% less than today's between 2 and 3.5 billion years ago (remember that a doubling of CO 2 represents only a 2% change in radiative balance.). This should have led to an ice-covered planet, but the geological evidence indicated flowing water and possibly the complete absence of ice. The most commonly proposed solutions involved large quantities of greenhouse gases, but each of these had profound difficulties. For example, the amount of CO 2 required was at least ten times greater than permitted by geological evidence from paleosols and other proxies. Rondanelli and Lindzen (2010b) examined the potential of upper level cirrus in the tropics to resolve the problem and found that this was a possible solution. Invoking the iris effect with a magnitude of −22%/K allowed resolution without any additional greenhouse gas, while −5%/K allowed for an ice-free tropics while also allowing complete resolution with the addition of CO 2 well within what was allowed by geological data.
The assertion by Hartmann and Michelson (2002a) "the upper-level cirrus fraction in LCH10 intruded into the TWP from the subtropics have led a negative correlation with SST." was tested with the 15-years of accumulated satellite data by Choi et al. (2017). They found that the assertion was groundless. The variation in upper-level cirrus arise from the equatorial convective cloud systems in the sub-area of TWP. As to the cloud radiative forcing, the direct response appears in SST tendency (the difference between the two consecutive months) rather than SST itself. Choi et al. (2017) showed that the cloud influence in SST tendency is trivial over the TWP, but significant over the tropical eastern Pacific. This means that the upper level cirrus over the TWP was influenced by SST, but not vice versa.
As already mentioned in Section 1, one important implication of the highly variable tropical upper-level cirrus (regardless of the iris effect), is that it is no longer possible to isolate the water vapor feedback from the infrared cloud feedback associated with upper-level cirrus. It is, rather, more appropriate to simply look for the infrared (longwave) feedback. Similarly, it is reasonable to look for a shortwave feedback as well. Choi (2009, 2011) using satellite data from the ERBE (Earth Radiation Budget Experiment) and from both ERBE and CERES (Clouds and the Earth's Radiant Energy System) respectively, investigated the dependence of variations in both longwave and shortwave radiation as seen from space on SST. We note that criticisms (Chung et al. 2010;Murphy 2010;Trenberth et al. 2010;Trenberth et al. 2011) came out against Lindzen and Choi (2009). They pointed out our misinterpretation of the observed tropical radiative dependence on SST with respect to 'global climate sensitivity'. In response to this criticism, Lindzen and Choi (2011) corrected the issue and again concluded that the sensitivity was low. Subsequently, Bates (2016) examined this in greater detail, attempting to improve the simple procedure of simply averaging feedbacks concentrated in the tropics over all latitudes.
As is common in feedback studies in other fields, it was necessary to restrict oneself to fluctuations over time intervals long compared to the time scale for the feedback processes while short compared to the time for climate equilibration (Lindzen and Choi 2009). Simply regressing over the whole time record (as in Forster and Gregory 2006) can and often does lead to biases in feedback estimation (explicitly shown in Figs. 6 and 7 of Lindzen and Choi 2011). What we found was a fairly clear negative longwave feedback, an ambiguous signal for the shortwave feedback, and large differences between model responses to specified SSTs compared to what was observed.
The negative longwave feedback is compatible with the results of LCH01. Mauritsen and Stevens (2015) in introducing the iris effect into their model (by making precipitation efficiency increase due to clustering of convection in mesoscale systems) also found that the longwave feedback was essentially eliminated, but found that, in their model, shortwave feedback still provided a net positive feedback that brought climate sensitivity to the low end of the range in Intergovernmental Panel on Climate Change (IPCC) reports. The enhanced clustering of convection (and high cloud fraction decrease consistent with the iris effect) modeled in Mauritsen and Stevens (2015) has been found in several recent observational investigations (Zelinka and Hartmann 2011;Tobin et al. 2012;Choi et al. 2017;Saint-Lu et al. 2020).
A recent study by Choi et al. (2017) found the clusteringinduced observed negative relations between MODIS cirrus clouds and TRMM convective precipitations normalized by total precipitations (that is mainly due to clustering of convective towers) would give important constraints on equilibrium climate sensitivity (to 2 × CO 2 ) lower than 3 K in current climate models (CMIP5) while leaving uncertain the actual value. Thus, it should be important to extract cloud feedback effects as reliably as possible for climate sensitivity. Current climate models simulate the climate sensitivity between 1.5 and 4.5 K, and much of the uncertainty arises from cloud feedback (Sherwood et al. 2020). Note that some studies (including some from major data centers) also confirm the implausibility of higher climate sensitivities (Fyfe et al. 2013;Stott et al. 2013;Lewis and Curry, 2015;Lindzen, 2019).

Signal-to-Noise Problems
Studying feedback effects on clouds inevitably involves at least two intrinsic difficulties: measurements are always to some extent indirect (that is to say one does not in general measure cumulus mass flux or stratiform cloud area directly), and, perhaps more importantly, clouds depend on many factors other than surface temperature so that attempts to relate cloud cover to temperature have large scatters and relatively low correlations. As Rondanelli and Lindzen (2008) noted: "small values of the regression coefficients as the ones we observe do not imply a small effect but rather the presence of noise." This indicates that sufficient noise can render signal detection difficult, and, in some cases impossible (Choi et al. 2014).
This issue can be simply tested by relating the top-of-the atmosphere radiative flux anomalies (R) to the sea surface temperature anomalies (T). It is straightforward to take the regression coefficients of R onto T (hereafter, R/T) for confirmation of the R-T feedback relation. We calculated the lagged R/T based on monthly R from CERES and T from NOAA OI SST in 1°× 1°resolution for 20 years (March 2000 to February 2020). Figure 1 basically implies that observational R/Ts are very variable with the methods used for calculation of R/T. Figure 1a regressed R onto T area-averaged over the tropics (20°S-20°N) ([R]/[T]), while Fig. 1b regressed R onto T from each ith grid point within the tropics (R i /T i ).
The feedback relation between R and T is manifested at a simultaneous time, not a lagged time. Thus, the magnitude of the regression coefficient R/T should be maximized at zero lag, but weaken with lag. This has been explicitly shown in Choi et al. (2014), and in their simulations, only the peaked R/ T at zero lag securely represented a reliable feedback parameter. We do see a clear zero-lag peak of [R]/[T] is −4 W m −2 K −1 in the longwave (the red line in Fig. 1a) a negative longwave feedback consistent with many previous studies (Lindzen and Choi 2011;Choi et al. 2014;Mauritsen and Stevens 2015). Choi et al. (2014) also noted that the unclear regression coefficient R/T is also induced by numerous non-feedback noise. For example, processes other than SST caused changes in clouds. That is, autonomous cloud-induced R influences T (but not vice versa). Also, the third factors (diurnal or seasonal solar illumination, ENSO, etc.) readily change both R and T regardless of R-T feedback relation. These cases altogether make the R/T peak at lagged time (e.g., ENSO effect shown in Trenberth et al. 2011), and the zero-lag R/T far from a false feedback parameter. We reconfirm that it is hardly avoided by the simple regression in the shortwave (the blue line in Fig. 1a). Choi et al. (2014) noted that the noise dominated the attempted measurements of shortwave feedback leaving the feedback essentially undetectable.
In the grid-to-grid regression analysis (Fig. 1b), the sign of R i /T i was even the opposite to [R]/[T] in Fig. 1a. Moreover, very sharp R i /T i peaks appear at 1-month lag, and do not represent the feedback effect. This may be just due to the movement of convective clouds to the next grid (i + 1) where the surface has been warmer. R i + 1 due to the clouds should have been related to T i where the clouds popped up. Otherwise, the grid-to-grid relation readily misinterprets the feedback parameter.
For better assessment of the iris effect, some studies pointed out methodological advances are needed. Noting that part of the problem was that cloud changes could produce rapid changes in skin temperature of the ocean (which is what AVHRR is essentially measuring), Cho et al. (2012) found that using SST weighted according to the absence of clouds could reduce noise but not sufficiently to allow a reliable determination of the shortwave feedback.
Nevertheless, because this has led to more efficient terrestrial emission but less solar reflection (Trenberth and Fasullo 2009), its net radiative feedback effects continue to be questioned. Sherwood et al. (2020), in their review, inferred that the net effect may be too small to substantially affect climate sensitivity (about −0.2 W/m 2 /K with a standard deviation 0.2 W/m 2 /K). We note that this value is based on Williams and Pierrehumbert (2017) which is a grid-to-grid analysis. They have also reviewed the possibility for the net effect to be positive Chambers et al. 2002), and neutral (Zelinka and Hartmann 2011). The continual disagreements should be understood in line with Fig. 1. Studies without confirming the lagged relation between radiative fluxes and SST should be treated with extreme caution.

Suggestions for Further Study
As seen in subsequent studies, the iris effect remains a viable mechanism. To be sure, there remain such questions as the exact mechanism for precipitation efficiency (Sui et al. 2020), and the proper normalization for application to climate change. However, neither of these appears to be capable of changing the major conclusions.
Normalization (by either cumulus flux, total precipitation, or some combination) is still needed so as to look at the reaction of anvil cirrus per unit cumulus flux to SST warming. Otherwise, the iris effect would be underestimated (or veiled) because cirrus fraction is readily observed to increase with more cumulus towers concentrated on warmer areas. If possible, we recommend the normalization by total rain since that by the area of cumulus tower may be unresolved by grid-scale climate data.
Nevertheless, without the normalization, one can qualitatively find the reduction of cirrus with SST warming in various (both passive and active) satellite observations (Igel et al. 2014). Based on Choi et al. (2017), the best domain for the iris effect study would be the TWP that is used in LCH01. An area smaller than this domain could not stably capture the iris effect . The sub-areas of previous studies (Horváth and Soden 2008;Rondanelli and Lindzen 2008) would require normalization for discerning the iris effect. Therefore, the choice of domains is important, and comparison of smaller or larger tropical regions should be useful. Too small a domain may cut off total cirrus response, while too large domain invites subtropical influence. Dynamical coupling may be involved with the mechanism since several convective cloud systems embody the iris effect often by their dynamic aggregations (Wing and Emanuel 2014;Muller and Bony 2015).
It is advisable to use instrument platforms better suited to the measurements we have described. Microwave sensors on a geostationary satellite would be a significant help, but this does not seem likely for the immediate future. Hourly variation of upper-level clouds and moist areas in a geostationary satellite in relation to SST should be studied more carefully. According to LCH01 and the rebuttal papers, this geostationary imagery is important for studying the iris effect and may be a key to revealing something that cannot be seen from a sun-synchronous satellite in most previous papers.
Considerable uncertainty in the iris effect arises simply from SST per se. SST should be carefully selected especially when quantifying the iris effect. As LCH01 noted, cloudweighted SSTs beneath cumulonimbus towers are theoretically ideal, but those SSTs in sampling are heavily affected by Fig. 1 The regression coefficients of the top-of-the atmosphere radiative flux anomalies (R) onto the sea surface temperature anomalies (T), based on the area-averages ([]) over the tropics (20°S-20°N) (a), and ith grid point values within the tropics (b). The negative (positive) lag means that R precedes T (T precedes R). The red and blue line indicates the longwave (LW) and shortwave (SW), respectively. Only the zero-lag regression coefficient of −4 W m −2 K −1 in the longwave represents a reliable negative longwave feedback parameter consistent with many previous studies (Lindzen and Choi 2011;Choi et al. 2014;Mauritsen and Stevens 2015). The other regression coefficients are false feedback parameters. See the text for more details the clouds themselves. SSTs often drop due to cloud shielding effects. Using this SST may obscure the clear delineation of anvil response to SST. To obtain SST as a pre-condition of convective cloud system, at least hourly cloud data are needed to screen out bad pixels. A more detailed discussion of possible methodology may be found in Cho et al. (2012).
Further observational evidence, clarifying the mechanism, should help models to incorporate the iris effect. It is also important to identify the differences between current climate modeling of cirrus clouds and the observations. Most of modeled cirrus cloud may not be those 'detrained' from cumulonimbus towers, but entirely long-lived condensates in situ at the upper levels (Seeley et al. 2019). This means that the iris mechanism is not in the model lacking cloud resolving schemes from the very beginning (Retsch et al. 2019). Modelers would have to select cloud fractions on the top one to three vertical levels Li et al. 2019); the more cloud levels selected, the larger the cloud fraction and albedo calculated. In addition, the response of upper-level clouds to SST is dependent on the vertical resolution of the model (Ohno et al. 2019). Practically, the cirrus calculated have higher albedo than observed (Li et al. 2019). This difference certainly amplifies modeled shortwave radiative impact for any iris effect experiments. However, by appropriately implementing the iris effect, it is also true that climate models (including radiative-convective equilibrium models) may be useful to estimate the impact of the iris effect on global climate (Zhao et al. 2016). Recent studies use radiative-convective equilibrium models with the multiple cloud layers for assessment of the iris effect (Bony et al. 2016;Ohno and Satoh 2018). The issue of estimating the impact of tropical cloud feedbacks on the globe as a whole, should probably be reassessed in the light of recent work on baroclinic equilibration (baroclinic instability being the major source of meridional heat transport in the atmosphere (Jansen and Ferrari 2013;Lindzen 2020). In this connection, recent studies have shown that poleward energy transport in the atmosphere tends to prevent excessive Arctic warming (Hwang et al. 2019. Again, the iris effect is basically an infrared feedback due to the large radiative effect of high cirrus in the infrared. Adding the shortwave effect to this process amplifies uncertainty considerably. In the combinations of MODIS and CERES satellite observations, the outgoing longwave radiation clearly decreased with the increase in cirrus fraction . In this calculation, the dependence of absorbed solar radiation on cirrus fraction was much weaker. Observational calculation of shortwave cloud radiative effects is commonly by clear-sky minus total-sky absorbed solar radiation (Williams and Pierrehumbert 2017). However, that approach gives rise to overestimation of the radiative effects of cirrus, since the calculated cloud radiative effect is not just of high cirrus clouds, but of all total-column clouds including much of low-level clouds below cirrus (Choi and Ho 2006).
Radiative transfer simulations show that variable low clouds below the cirrus clouds induce large variation in cloud radiative effects, even for a given fixed top pressure of high clouds and the fixed total column cloud optical thickness (Kang et al. 2020). The iris effect with multiple cloud layers may well be investigated with a simple 3.5-box model of the sort employed in LCH01. As the iris effect is the immediate response of clouds to SST, this shortwave component should be treated very cautiously. As noted above, the uncertainties in the shortwave feedback involve contributions other than those associated with the iris effect per se.