To enhance comprehension of the brightening effect induced by cirrus thinning, it is better to illustrate a comparative analysis between this brightening effect and corresponding cooling effect. To facilitate communication, the symbol "Δ" is employed to represent the discrepancies ("△") in relation to the REF simulation from cirrus thinning simulations (HET or SEED). In accordance with previous studies [
23,
58,
81], the cooling effect is quantified by anomalies in CRE at the top of the atmosphere (ΔCRE
TOA). Analogous to the cooling effect, the brightening effect is quantified by anomalies in CRE at Earth's surface (ΔCRE
bri). It is noteworthy that ΔCRE
TOA has longwave (Earth radiation) net flux (including both downward and upward irradiance) and shortwave (solar radiation) net flux (ΔCRE
TOAlw and ΔCRE
TOAsw) whereas ΔCRE
bri solely takes into account the component of downward solar radiation. Here, we not only demonstrate the brightening effect and cooling effect but also place significant emphasis on understanding the corresponding mechanisms. These mechanisms could yield more valuable insights for leveraging the brightening effect.
3.1. Impacts on Cloud Properties
The changes in ICs number concentration caused by cirrus thinning (i.e., both HET and SEED simulations) is analyzed firstly (
Figure 2). After artificially turning off homogeneous freezing (i.e., HET simulation), the average number concentration of newly formed ICs under cirrus conditions (
Ninuc) drastically decreases to a very low level, especially in the Southern Hemisphere where INPs are scarce. Compared to the HET simulation,
Ninuc is obviously increased in the SEED simulation due to seeding INPs. However,
Ninuc from the SEED simulation is also much lower than that from the REF simulation. In cirrus clouds, the ICs number concentration (i.e.,
Ni) is primarily influenced by the process of ice nucleation (i.e.,
Ninuc) [
69,
74]. As expected, the zonal mean
Ni from both HET and SEED simulations is significantly decreased above the –37 ℃ isotherms (i.e., cirrus clouds). All these three simulations show that the
Ni in mixed-phase clouds at mid-to-high latitudes is relatively substantial. This might be due to convective detrainment, which provides a lot of ICs. Because there is no homogeneous nucleation in mixed-phase cloud scheme, the
Ni in mixed-phase clouds does not significantly decrease in either the HET or SEED simulations. In alignment with the noticeable decrease in
Ni (mainly in cirrus clouds), the vertically integrated
Ni (i.e., column
Ni) also significantly decreases in both HET and SEED simulations. Taken overall, these cirrus thinning simulations have successfully achieved their objective of significantly reducing
Ni.
The decrease in
Ni (i.e., cirrus thinning) impacts not only the cloud water in cirrus clouds but also the cloud water in mixed-phase and liquid clouds (
Figure 3). The ice water content (IWC) in cirrus clouds from both the HET and SEED simulations shows a notable decrease due to lower
Ni. Conversely, the IWC in mixed-phase clouds shows an increase across most regions. The main reason might be that cirrus thinning reduces atmospheric stability through its impact on the radiation budget, thereby instigating increased convective activity, which brings more water to mixed-phase cloud layers. The ice water path (IWP) is significantly decreased in most regions because the decrease in IWC in cirrus clouds is stronger than the increase in IWC in mixed-phase clouds. It is worth noting that, in certain regions (e.g., middle Africa and northern Brazil), the IWP is increased because the decreases in IWC within cirrus clouds are slight (which is consistent with the slight decreasing
Ni,
Figure 2) and these decreases are even less significant than the increases in IWC within mixed-phase clouds there. The liquid water content (LWC) and liquid water path (LWP) are also impacted by the thinning of cirrus clouds, as shown in both HET and SEED simulations. However, these changes in LWC and LWP (i.e., ΔLWC and ΔLWP) are not as noticeable as the ΔIWC and ΔIWP. Despite the overall less noticeable changes in LWC and LWP, it's important to highlight that, in some low- and mid-latitude regions, there are statistically significant decreases in both LWC and LWP. Furthermore, in terms of global mean values, both LWC and LWP are also significantly decreased. One possible reason is that the larger cirrus cloud ICs (associated with cirrus thinning, not shown) fall into the lower mixed-phase and liquid cloud layers and enhances the efficiency of converting cloud droplets into precipitation [
59,
78]. Another possible reason is that cirrus thinning results in more convective activities and convective precipitation which would consume more cloud water [
42,
58,
59,
82]. The above analyses suggest that cirrus thinning also has considerable impacts on the lower mixed-phase and liquid clouds.
The changes in CRE mainly depend on the changes in cloud optical depth (COD).
Figure 4 shows the changes in COD. Here, the COD in both longwave and shortwave bands (COD
lw and COD
sw) are shown. Generally, the changes in COD
lw and COD
sw (i.e., ΔCOD
lw and ΔCOD
sw) show almost the same special pattern because the change in cloud water is the primary determinant for both COD
lw and COD
sw. The COD from cirrus clouds (iCOD) is significantly decreased in both longwave and shortwave bands (iCOD
lw and iCOD
sw) over most regions. The global mean values of iCOD
lw and iCOD
sw are decreased by more than half, especially for the HET simulation. The ΔiCOD
lw and ΔiCOD
sw pass the significance test over most regions except for middle Africa and northern Brazil. Compared to ΔiCOD, the changes in COD from mixed-phase and liquid clouds (ΔmlCOD) become complicated. This is in agreement with the complex changes in cloud water within mixed-phase and liquid clouds (
Figure 3). Both ΔmlCOD
lw and ΔmlCOD
sw show statistically significant decreases (i.e., negative values) over some low- and mid-latitude regions. In terms of global mean values, ΔmlCOD
lw is stronger than ΔiCOD
lw, and ΔmlCOD
sw is stronger than ΔiCOD
sw. In short, cirrus thinning leads to a noticeable and consistent decrease in iCOD across most regions. Additionally, it also results in a substantial reduction in mlCOD over some low- and mid-latitude regions.
3.2. Brightening Effect and Cooling Effect
In this section, we quantify the brightening effect and cooling effect of cirrus thinning using CREbri variables (e.g., ΔCREbri and ΔmlCREbri) and CRETOA variables (e.g., ΔiCRETOA and ΔmlCRETOA), respectively. A positive value of the CREbri variables indicates a brightening effect, while a negative value suggests the opposite. Similarly, a negative value of the CRETOA variables signifies a cooling effect and a positive value implies a warming effect.
Firstly, we analyze the CRE
TOA variables and CRE
bri variables solely from cirrus clouds (
Figure 5). The positive iCRE
TOAlw indicates that cirrus clouds contribute to the reduction of Earth's outgoing longwave radiation, thereby keeping our planet warm. The negative iCRE
TOAsw indicates that cirrus clouds are not conducive to the Earth's absorption of solar radiation, resulting in a cooling effect. Overall, the warming effect (i.e., positive iCRE
TOAlw) is evidently stronger than the cooling effect (i.e., the absolute value of negative iCRE
TOAsw). The globally averaged iCRE
TOA (iCRE
TOAlw + iCRE
TOAsw) from the REF simulation is 6.53 W m
−2 (net warming effect). This value falls within the potential range reported in recent studies (4.5 to 6.8 W m
−2) [
23,
32,
56,
80,
83]. The negative iCRE
bri suggests that cirrus clouds cause a dimming effect on the Earth's surface. The value of iCRE
bri is close to the value of iCRE
TOAsw, but a little weaker (less negative). This can be explained by that, in the absence of cirrus clouds, more downward solar irradiance can enter the mixed-phase and liquid cloud layers. Although the mixed-phase and liquid clouds scatter and absorb some solar radiation, most of it could reach the Earth's surface causing a brightening effect. All these radiative fluxes (i.e., iCRE
TOA, iCRE
TOAlw, iCRE
TOAsw, and iCRE
bri) show a similar spatial pattern that aligns with the COD of cirrus clouds (i.e., iCOD
lw and iCOD
sw). After cirrus clouds become thin (i.e., the HET and SEED simulations), the net warming effect and surface dimming effect from cirrus clouds also weaken. In other words, cirrus thinning leads to a cooling and brightening effect. The globally averaged ΔiCRE
TOA from the HET and SEED simulations are −3.56 ± 0.04 and −2.46 ± 0.04 W m
−2, respectively. The globally averaged ΔiCRE
bri from the HET and SEED simulations are 2.78 ± 0.03 and 2.19 ± 0.03 W m
−2, respectively. These global mean values suggest that the cirrus cloud net warming effect and surface dimming effect from the REF simulation are reduced by about half. Both the HET and SEED simulations show that the cooling effect (negative ΔiCRE
TOA) and brightening effect (positive ΔiCRE
bri) are statistically significant over most regions.
Secondly, the radiative effects of mixed-phase and liquid clouds are analyzed (
Figure 6). Similar to cirrus clouds, mixed-phase and liquid clouds also exhibit longwave warming effect (positive mlCRE
TOAlw) and shortwave cooling effect (negative mlCRE
TOAsw). Compared to cirrus clouds (iCRE
TOAlw and iCOD
lw), mlCRE
TOAlw is only increased by about half despite a roughly twenty-fold increase in mlCOD
lw. This weaker efficiency (i.e., CRE/COD) primarily stems from the relatively small temperature difference between Earth's surface and these clouds. Unlike mlCRE
TOAlw, mlCRE
TOAsw is approximately ten times stronger than iCRE
TOAsw. Therefore, the shortwave cooling effect (negative mlCRE
TOAsw) is much stronger than the longwave warming effect (positive mlCRE
TOAlw) over most regions. Mixed-phase and liquid clouds show a net cooling effect (negative mlCRE
TOA). Mixed-phase and liquid clouds also make the Earth's surface dimmer (negative mlCRE
bri). Here, mlCRE
bri (downward solar irradiance at the Earth’s surface) is a little stronger (more negative) than mlCRE
TOAsw (net radiative flux). The main reason is that a portion of mlCRE
bri (surface albedo) is reflected back into the atmosphere. Because the impact of cirrus thinning on mixed-phase and liquid clouds is complex (
Figure 3 and
Figure 4), the changes in radiative effects caused by cirrus thinning are statistically insignificant over approximately half of the Earth. These two cirrus thinning simulations (i.e., HET and SEED simulations) show that ΔmlCRE
TOAlw is generally negative (cooling effect) and ΔmlCRE
TOAsw is generally positive (warming effect). The positive ΔmlCRE
TOAsw is obviously stronger than the absolute value of ΔmlCRE
TOAlw over most low- and mid-latitude regions where solar radiation is relatively more dominant. Therefore, the ΔmlCRE
TOA (ΔmlCRE
TOAlw + ΔmlCRE
TOAsw) from cirrus thinning simulations are positive over most low- and mid-latitude regions. Over the Western Pacific Warm Pool and nearby regions, the net warming effect (i.e., positive ΔmlCRE
TOA) is considerable and statistically significant. The globally averaged ΔmlCRE
TOA from the HET and SEED simulations are 1.35 ± 0.18 and 1.25 ± 0.16 W m
−2, respectively. This warming effect caused by the alterations in mixed-phase and liquid clouds would counteract, to some extent, the cooling effect derived from the thinning of cirrus clouds alone (−3.56 ± 0.04 and −2.46 ± 0.04 W m
−2). The next paragraph will analyze this influence. The globally averaged ΔmlCRE
bri from the HET and SEED simulations are 3.05 ± 0.25 and 2.50 ± 0.21 W m
−2, respectively. The brightening effect from mixed-phase and liquid clouds is a little larger than from cirrus clouds (2.78 ± 0.03 and 2.19 ± 0.03 W m
−2). In short, after implementing the cirrus thinning approach, the cooling and dimming effects of mixed-phase and liquid clouds become weaker. This is consistent with the decrease in mlCOD caused by cirrus thinning.
Finally, the brightening effect and cooling effect caused by cirrus thinning are quantified by the changes in radiative effects of entire clouds (ice, mixed-phase, and liquid clouds;
Figure 7). In terms of solar radiation, the entire cloud has a globally averaged shortwave cooling effect (CRE
TOAsw is −56.42 W m
−2) and dimming effect (CRE
bri is −66.91 W m
−2). Cirrus thinning not only significantly decreases the shortwave cooling effect and dimming effect of cirrus clouds but also leads to substantial reductions (less negative in radiative effects) within the lower mixed-phase and liquid clouds. Hence, the entire cloud exhibits considerable reductions in both the shortwave cooling effect and the dimming effect. The globally averaged ΔCRE
TOAsw from HET and SEED simulations are 6.02 ± 0.21 and 4.85 ± 0.17 W m
−2, respectively. The globally averaged ΔCRE
bri from the HET and SEED simulations are 5.83 ± 0.26 and 4.69 ± 0.21 W m
−2, respectively. The brightening effect (i.e., positive ΔCRE
bri) is close to the shortwave warming effect (positive ΔCRE
TOAsw), and they have a similar spatial pattern. The ΔCRE
bri are predominantly significant across most low- and mid-latitude regions, primarily attributable to the intense solar radiation present in these areas. In terms of longwave radiation, the cloud has a global averaged warming effect (CRE
TOAlw is 28.27 W m
−2). Although the globally averaged CRE
TOAlw is roughly half of the absolute value of CRE
TOAsw, ΔCRE
TOAlw is generally stronger than ΔCRE
TOAsw due to the dominant contribution from cirrus clouds (i.e., ΔiCRE
TOAlw). In terms of the sum of shortwave and longwave radiation, the cloud has a globally averaged cooling effect (CRE
TOA is −28.15 W m
−2). The ΔCRE
TOA caused by cirrus thinning also shows cooling effect. The ΔCRE
TOA from the HET and SEED simulations are −2.21 ± 0.18 and −1.21 ± 0.19 W m
−2, respectively. Unlike the brightening effect, the cooling effect (i.e., negative ΔCRE
TOA) is statistically significant over most high-latitude regions because solar radiation is weak at high latitudes. In comparison to artificially halting homogeneous nucleation (i.e., the HET simulation), the globally averaged cooling effect resulting from the seeding method decreases by approximately half. Unlike the cooling effect, the relative reduction of the brightening effect is considerably less, falling by about one-fourth. In the SEED simulation, the globally averaged brightening effect is approximately four times stronger than the globally averaged cooling effect. This suggests that cirrus thinning via seeding INPs not only cools our planet but also provides much stronger brightening effect.