Geoengineering: Impact of Marine Cloud Brightening Control on the Extreme Temperature Change over East Asia

We investigated the effect of artificial marine cloud brightening on extreme temperatures over East Asia. We used simulation data from five global climate models which have conducted the GeoMIP G4cdnc experiment. G4cdnc was designed to simulate an increase in the cloud droplet number concentration of the global marine lower clouds by 50% under the greenhouse gas forcing of the RCP4.5 scenario. G4cdnc decreased the net radiative forcing in the top of the atmosphere more over the ocean, alleviating the rise in mean temperature under RCP4.5 forcing. For extreme temperatures, G4cdnc reduced both the monthly minimum of daily minimum temperature (TNn) and monthly maximum of daily maximum temperature (TXx). The response of TNn was higher than that of TXx, especially in the winter, over the Sea of Okhotsk and the interior of the continent. This spatial heterogeneity and seasonality of the response were associated with sea ice–albedo and snow–albedo feedbacks. We also calculated the efficacy of warming mitigation as a measure of the relative effect of geoengineering. The efficacy for TXx was higher than that for TNn, opposite to the absolute effect. After the termination of geoengineering, both TNn and TXx tended to rapidly revert to their trend under the RCP4.5 forcing.


Introduction
Since the industrial revolution, atmospheric concentrations of greenhouse gases (GHGs) have risen owing to anthropogenic activities. In particular, the rise in concentration of carbon dioxide has accelerated, recently exceeding 400 ppm. This has resulted in global warming, which is causing significant increases in the duration, frequency, and/or intensity of various types of extreme weather and climate events over East Asia and across the world. Many efforts have been devoted to solving the problems caused by global warming. For example, the Paris Agreement was signed by the United Nations Framework Convention on Climate Change (UNFCCC) in 2015 and aimed at limiting the increase in the global mean temperature to well below 2 • C above pre-industrial levels, and pursuing efforts to limit the increase to 1.5 • C. Each nation participating in the Agreement thus set an individual goal for the reduction of greenhouse gases and consented to meet it. However, it has been suggested that the Agreement is not sufficient to alleviate the damage caused by climate change [1]. Therefore, the geoengineering field is emerging; it is aimed at mitigating the adverse effects of global warming and has received significant attention. Geoengineering includes a broad No [29] To investigate the effects of the G4cdnc experiment for extreme climate events, we focused our analyses on monthly minimum and maximum extreme temperatures instead of annual values. Analysis using monthly extremes enables us to consider monthly variations (i.e., seasonality) and identify extremes more precisely. To this end, we calculated TNn, a monthly minimum of daily minimum temperatures, and TXx, a monthly maximum of daily maximum temperatures, as defined by the Expert Team of Climate Change Detection and Indices (ETCCDI) [30]. These extreme indices have been frequently used in previous research, both for observed weather [31] and model simulations ( [32][33][34][35]. In the GeoMIP G4cdnc experiment under the RCP4.5 scenario, four models simulated daily minimum (TASmin) and daily maximum (TASmax) temperatures with regard to extreme climatic events. The BNU-ESM model, however, did not simulate these variables. Moreover, the MIROC-ESM model simulated only the increase in CDNC, without the termination of the experiment. Therefore, we excluded the BNU-ESM model from the analysis of extreme temperatures and the MIROC-ESM model from the analysis after the termination of the experiment.
As the spatial resolutions of simulation data differ between the models, as indicated in Table 1, all outputs by G4cdnc and RCP4.5 models were re-gridded to 144 × 73 (2.5 • longitude × 2.5 • latitude) for the convenience of comparison and calculation. The interpolation was carried out using the Climate Data Operators (CDO) package, in which we chose a first-order conservative remapping algorithm [36]. To analyze the changes over East Asia, we selected an analysis domain of 20-60 • N, 90-160 • E.

Methods
To investigate the effect of the GeoMIP G4cdnc geoengineering experiment over East Asia, we selected three periods for analysis: base period: 2006-2019, progress period: 2060-2069, and termination period: 2081-2090. For the base period, we used RCP4.5 data, and for the progress and termination periods, both G4cdnc and RCP4.5 data were used. Although many studies for geoengineering have adopted a 40-year period for analysis [33,35,37], we adopted only the last 10 years from the progress and termination periods to sufficiently consider the saturated conditions of forcing.
To obtain better insights into the changes in extreme temperature, we checked the changes in mean temperature through the calculation of probability density functions (PDFs) of monthly temperatures. To create PDFs of monthly temperature, the method reported by Curry et al. and Ji et al. was employed [33,35]. At every grid point in each model, we first calculated the standardized monthly anomalies of monthly mean surface temperatures in the base period under RCP4. 5; i.e., τ Base   Table 1 and for the multi-model  ensemble; equal weights were used for each model. Then, we conducted a two-sample non-parametric Kolmogorov-Smirnov test to determine if the distribution under G4cdnc differs significantly from that under the RCP4.5 scenario in the progress period [38]. For extreme temperatures, we took the results of the difference between G4cdnc and its base experiment of RCP4.5 to investigate the effect of marine cloud brightening over East Asia. The change in the extreme temperature was investigated in terms of both temporal and spatial aspects, which were displayed using a timeseries and a field of the difference, respectively. For temporal changes, the mean of the difference deduced from the multi-model ensemble over both the progress and termination periods are printed on the upper right of the time series. Extreme temperatures in summer (JJA) and winter (DJF) were calculated to examine the effects of seasonality.
In the perspectives of establishing policies on local climate change, relative effects of geoengineering to warming magnitude under RCP4.5 may be more important than absolute effects. To consider this relative effect of G4cdnc, we adopted a simple quantitative metric, efficacy used in [33,39]. To do so, we first computed the root mean square (RMS) differences between experiments at each grid point as the following equation: where X is the extreme index of interest, Exp the perturbed experiment (G4cdnc or RCP4.5 in progress period), Base the base experiment (RCP4.5 in base period), and a sum over the grid cells (dA) of a specified area A. Then, we calculated the efficacy e(X) of G4cdnc defined by the following equation: In the case of e(X) → 1, the temperature under G4cdnc approaches the temperature in the RCP4.5 base period so that the G4cdnc experiment is the most effective in offsetting the warming under the RCP4.5 scenario. On the other hand, e(X) → 0 indicates that the G4cdnd experiment has little effect in reducing the warming. Values of e(X) can be negative if RMS differences for G4cdnc exceed those for RCP4.5. Efficacy was calculated for mean temperatures (TAS, TASmax, and TASmin) as well as for the two extreme temperature indices (TNn and TXx) to check if there are significant differences between relative effects on extreme temperatures and mean temperatures. We also analyzed the efficacy not only for the progress period but also for the period after the termination of geoengineering in order to investigate how long the relative effects can last.

Radiation Budget at TOA
As the G4cdnc experiment simulates the increase in the reflection of solar radiation by marine clouds, we first examined the changes in the radiation budget at the top of the atmosphere (TOA), expecting an increase in upward shortwave radiation. Figure 1 shows the global distributions of the differences in radiative fluxes at the TOA between G4cdnc and the RCP4.5 experiments during the progress period. During this period, the upward shortwave radiation at the TOA increased by 2.95 Wm −2 in the global mean multi-model ensemble average. Although there was a decrease in certain areas, especially in the western equatorial Pacific, the reflection of solar radiation was increased overall around the globe with a large increase over the eastern parts of the Pacific and Atlantic Oceans (Figure 1a).
Atmosphere 2020, 11, 1345 5 of 18 reflective shortwave radiation at the TOA resulted in a decrease in outgoing longwave radiation (OLR) at the TOA, indicating a surface cooling not only over the global ocean but over the global land ( Figure 1b). A significant decrease in OLR over tropical rainforest regions including the Amazon, equatorial Africa, and the Indo-Pakistan region is noteworthy. In the global mean and multi-model ensemble average, the reduction of OLR (−2.27 Wm −2 ) largely compensates the energy loss by the shortwave reflection. Overall, the changes in the upward shortwave radiation and OLR tend to compensate each other. Nevertheless, the net radiation at the TOA (net shortwave radiation minus OLR) decreased largely over the oceans, except for the ITCZ area and an increased OLR and Arctic sea-ice region, and increased over most of the land areas, except for the Tibetan Plateau where the upward shortwave radiation significantly increased, unlike the other land areas (Figure 1c). The global mean multimodel ensemble average of the net radiation at the TOA is −0.68 Wm −2 with a downward positive sign during the progress period. This amount of change in the net radiation corresponds to a 15% reduction of the radiative forcing (4.5 Wm −2 ) induced by anthropogenic activities under the RCP4.5 scenario.
Other than the sea-ice/snow albedo feedback effect, the heterogeneity in regional responses to G4cdnc can also be explained by the changes in cloudiness, soil moisture, and other climate variables, but further investigation into the subject is beyond the scope of this paper.
The differences in the TOA radiation fluxes are replicated in Figure 2 and focused on East Asia, the analysis domain for this study. The marine cloud brightening increased the reflection of shortwave radiation at the TOA over most of the domain but caused a significant decrease over the Gobi Desert area (Figure 2a). There was a considerable increase over the Sea of Okhotsk, which can be related to a positive sea ice-albedo feedback. The increase over Southern China and the Tibetan Plateau can also be related to a positive snow-albedo feedback (figures are not shown). Meanwhile, OLR decreased over most of the domain area even in the area with a decreased shortwave reflection Except for some areas of the inter-tropical convergence zone (ITCZ), the overall increase in reflective shortwave radiation at the TOA resulted in a decrease in outgoing longwave radiation (OLR) at the TOA, indicating a surface cooling not only over the global ocean but over the global land ( Figure 1b). A significant decrease in OLR over tropical rainforest regions including the Amazon, equatorial Africa, and the Indo-Pakistan region is noteworthy. In the global mean and multi-model ensemble average, the reduction of OLR (−2.27 Wm −2 ) largely compensates the energy loss by the shortwave reflection.
Overall, the changes in the upward shortwave radiation and OLR tend to compensate each other. Nevertheless, the net radiation at the TOA (net shortwave radiation minus OLR) decreased largely over the oceans, except for the ITCZ area and an increased OLR and Arctic sea-ice region, and increased over most of the land areas, except for the Tibetan Plateau where the upward shortwave radiation significantly increased, unlike the other land areas (Figure 1c). The global mean multi-model ensemble average of the net radiation at the TOA is −0.68 Wm −2 with a downward positive sign during the progress period. This amount of change in the net radiation corresponds to a 15% reduction of the radiative forcing (4.5 Wm −2 ) induced by anthropogenic activities under the RCP4.5 scenario.
Other than the sea-ice/snow albedo feedback effect, the heterogeneity in regional responses to G4cdnc can also be explained by the changes in cloudiness, soil moisture, and other climate variables, but further investigation into the subject is beyond the scope of this paper.
The differences in the TOA radiation fluxes are replicated in Figure 2 and focused on East Asia, the analysis domain for this study. The marine cloud brightening increased the reflection of shortwave radiation at the TOA over most of the domain but caused a significant decrease over the Gobi Desert area (Figure 2a). There was a considerable increase over the Sea of Okhotsk, which can be related to a positive sea ice-albedo feedback. The increase over Southern China and the Tibetan Plateau can also be related to a positive snow-albedo feedback (figures are not shown). Meanwhile, OLR decreased over most of the domain area even in the area with a decreased shortwave reflection at the TOA (Figure 2b). This reduction of OLR induced by the marine cloud brightening resulted in an increase in downward net radiation at the TOA except for the Sea of Okhotsk and Tibetan Plateau areas (Figure 2c). at the TOA (Figure 2b). This reduction of OLR induced by the marine cloud brightening resulted in an increase in downward net radiation at the TOA except for the Sea of Okhotsk and Tibetan Plateau areas ( Figure 2c). Figure 1, except for over East Asia.

Mean Temperature
To check whether the global warming effect under the RCP4.5 scenario has been alleviated over East Asia through the radiative forcing changes induced by the marine cloud brightening, probability distribution functions (PDFs) for surface air mean temperature are shown in Figure 3 for each of the five models as well as their ensemble mean, using standarized monthly mean anomalies throughout the given periods. In the figure, the PDFs from the progress period of the G4cdnc and RCP4.5 results are compared with the baseline (the base period RCP4.5) results.

Mean Temperature
To check whether the global warming effect under the RCP4.5 scenario has been alleviated over East Asia through the radiative forcing changes induced by the marine cloud brightening, probability distribution functions (PDFs) for surface air mean temperature are shown in Figure 3 for each of the five models as well as their ensemble mean, using standarized monthly mean anomalies throughout the given periods. In the figure, the PDFs from the progress period of the G4cdnc and RCP4.5 results are compared with the baseline (the base period RCP4.5) results.
In Figure 3, the PDFs of the multi-model ensemble mean exhibit a modest increase in mean temperature between G4cdnc and the baseline, indicating warming despite the implementation of the G4cdnc experiment. The PDFs in the progress period, however, show that the distribution mean under G4cdnc is lower than that under RCP4.5 (G4cdnc: 1.206, RCP4.5: 3.469). The difference in PDFs between the G4cdnc and RCP4.5 experiments is statistically significant at the 95% level by the Kolmogorov-Smirnov test. Hence, this feature implies that the G4cdnc experiment significantly alleviates the additional rise in mean temperature caused by the RCP4.5 scenario over East Asia. The mean of the PDFs from all the models under the G4cdnc experiment is lower than of RCP4.5; the CSIRO-Mk3L-1-2 model resulted in the lowest mean of −0.326, and the MIRCO-ESM model resulted in the greatest mean of 1.198.
This cooling effect of mean temperature is consistent with the results of previous studies investigating the effects of other geoengineering experiments. For example, a study conducted by Curry et al. calculated the PDFs, the same as our method but for the globe, and focused on the difference in mean temperature between the GeoMIP G1 and CMIP5 abrupt4xCO 2 experiments [33]. In that study, PDFs obtained from the G1 experiment are located further left compared to the PDFs obtained from the abrupt4xCO 2 experiment, indicating a very high cooling effect caused by the G1 experiment. Ji et al. also calculated PDFs using the same algorithm for the globe and reported the difference between GeoMIP G4 and CMIP5 RCP4.5 [35]. Similarly, in that study, there was a discernable difference between the PDFs from the G4 experiment and those from RCP4.5, representing an alleviation of the increase in mean temperature caused by RCP4.5.
Atmosphere 2020, 11, x FOR PEER REVIEW 7 of 19 Figure 3. Probability distributions of standardized monthly mean anomalies for three experiments: G4cdnc (blue), RCP4.5 (red), and base period RCP4.5 (black). Probability density functions (PDFs) are over East Asia during the progress period (year 2060 to 2069). Anomalies are calculated with respect to the mean and standard deviation of the base period for the RCP4.5 scenario for each experiment. "False" shown in the upper right corner indicates that the difference in the distribution between G4cdnc and RCP4.5 is significant at the 95% level by the Kolmogorov-Smirnov test and "True" indicates vice versa.
In Figure 3, the PDFs of the multi-model ensemble mean exhibit a modest increase in mean temperature between G4cdnc and the baseline, indicating warming despite the implementation of the G4cdnc experiment. The PDFs in the progress period, however, show that the distribution mean under G4cdnc is lower than that under RCP4.5 (G4cdnc: 1.206, RCP4.5: 3.469). The difference in PDFs between the G4cdnc and RCP4.5 experiments is statistically significant at the 95% level by the Kolmogorov-Smirnov test. Hence, this feature implies that the G4cdnc experiment significantly alleviates the additional rise in mean temperature caused by the RCP4.5 scenario over East Asia. The mean of the PDFs from all the models under the G4cdnc experiment is lower than of RCP4.5; the CSIRO-Mk3L-1-2 model resulted in the lowest mean of −0.326, and the MIRCO-ESM model resulted in the greatest mean of 1.198.
This cooling effect of mean temperature is consistent with the results of previous studies investigating the effects of other geoengineering experiments. For example, a study conducted by Curry et al. calculated the PDFs, the same as our method but for the globe, and focused on the difference in mean temperature between the GeoMIP G1 and CMIP5 abrupt4xCO2 experiments [33]. In that study, PDFs obtained from the G1 experiment are located further left compared to the PDFs obtained from the abrupt4xCO2 experiment, indicating a very high cooling effect caused by the G1 experiment. Ji et al. also calculated PDFs using the same algorithm for the globe and reported the difference between GeoMIP G4 and CMIP5 RCP4.5 [35]. Similarly, in that study, there was a . Anomalies are calculated with respect to the mean and standard deviation of the base period for the RCP4.5 scenario for each experiment. "False" shown in the upper right corner indicates that the difference in the distribution between G4cdnc and RCP4.5 is significant at the 95% level by the Kolmogorov-Smirnov test and "True" indicates vice versa. Figure 3 also shows that the variance of the multi-model ensemble mean PDF obtained from the G4cdnc experiment is lower than that from the RCP4.5 experiment (G4cdnc: 1.655, RCP4.5: 2.748), which is consistent with a previous study result for the G1 experiment compared with the CMIP5 abrupt4xCO 2 experiment [33]. A greater variance in PDFs equates to a larger inter-annual variability in the temperature and probably to the more frequent occurrence of extreme temperatures. Assuming the latter case, extreme temperatures under the G4cdnc experiment may occur less frequently than under the RCP4.5 scenario. Most of the models resulted in a lower variance in PDFs under the G4cdnc experiment than under RCP4.5, with the exception of the CSIRO-Mk3L-1-2 model, presenting the opposite. Meanwhile, the skewness of the multi-model ensemble's PDFs under G4cdnc also decreases slightly compared with that under RCP4.5 (G4cdnc: 0.380, RCP4.5: 0.457).

Temporal and Spatial Response in Extreme Temperatures
To investigate the temporal changes in extreme temperatures over East Asia under the G4cdnc experiment, the differences between G4cdnc and RCP4.5 for the TNn and TXx indices have been opposite. Meanwhile, the skewness of the multi-model ensemble's PDFs under G4cdnc also decreases slightly compared with that under RCP4.5 (G4cdnc: 0.380, RCP4.5: 0.457).

Temporal and Spatial Response in Extreme Temperatures
To investigate the temporal changes in extreme temperatures over East Asia under the G4cdnc experiment, the differences between G4cdnc and RCP4.5 for the TNn and TXx indices have been analyzed and their time series are shown in Figure 4. Hereafter, ∆ denotes the multi-model ensemble mean of a difference between G4cdnc and RCP4.5 averaged over the progress period (2060-2069). . As a result, G4cdnc lowers not only the mean temperature but also extreme temperatures over East Asia by more or less than one degree. Among the models, CanESM2 produced the least change in the extreme temperature indices, which is consistent with the result for mean temperature presented in Section 3.2.
It is notable that the ∆TNn is larger than the ∆TXx, suggesting an implication for extreme cold.  . As a result, G4cdnc lowers not only the mean temperature but also extreme temperatures over East Asia by more or less than one degree. Among the models, CanESM2 produced the least change in the extreme temperature indices, which is consistent with the result for mean temperature presented in Section 3.2.
It is notable that the ∆TNn is larger than the ∆TXx, suggesting an implication for extreme cold.  [32,34,35,40,41]. Figure 5 shows the geographical distribution of the difference in extreme temperatures between the G4cdnc and RCP4.5 experiments over East Asia during the progress period. The figures demonstrate negative values of ∆TNn and ∆TXx over the entire area of East Asia. However, there appears to be evident regional heterogeneity in the response of extreme temperatures to the G4cdnc simulation.
Notably, the ∆TNn in Figure 5a has a considerably large negative value over the Sea of Okhotsk, lower than −4 • C, whereas ∆TXx has a minimum value of −2 • C over the same region. Several studies noted the larger temperature cooling effect over high-latitude oceanic regions including the Sea of Okhotsk due to SRM geoengineering [35,42,43]. This large temperature response over high-latitude oceanic regions appears to be in line with the well-known global warming effect, polar amplification [44][45][46]. Ji et al. showed that both solar dimming (G1) and stratospheric aerosol (G4) geoengineering resulted in a large ∆TNn over high-latitude oceanic regions in both hemispheres, which was not seen in TXx [35]. It has been found that this high-latitude cooling by geoengineering experiments, such as cloud brightening [42] and solar dimming [43], could restore the sea ice extent to the state of their control experiments.  [32,34,35,40,41]. Figure 5 shows the geographical distribution of the difference in extreme temperatures between the G4cdnc and RCP4.5 experiments over East Asia during the progress period. The figures demonstrate negative values of ∆TNn and ∆TXx over the entire area of East Asia. However, there appears to be evident regional heterogeneity in the response of extreme temperatures to the G4cdnc simulation. Notably, the ∆TNn in Figure 5a has a considerably large negative value over the Sea of Okhotsk, lower than −4 °C, whereas ∆TXx has a minimum value of −2 °C over the same region. Several studies noted the larger temperature cooling effect over high-latitude oceanic regions including the Sea of Okhotsk due to SRM geoengineering [35,42,43]. This large temperature response over high-latitude oceanic regions appears to be in line with the well-known global warming effect, polar amplification [44][45][46]. Ji et al. showed that both solar dimming (G1) and stratospheric aerosol (G4) geoengineering resulted in a large ∆TNn over high-latitude oceanic regions in both hemispheres, which was not seen in TXx [35]. It has been found that this high-latitude cooling by geoengineering experiments, such as cloud brightening [42] and solar dimming [43], could restore the sea ice extent to the state of their control experiments.
Although ∆TNn and ∆TXx resulted in cooling all over the area, a land-sea difference of the results is distinct in Figure 5. The G4cdnc experiment has a larger cooling effect of ∆TNn and ∆TXx over the interior of the continent compared to the ocean, with the exception of the Sea of Okhotsk. This land-sea difference in response can also be seen in previous studies that have investigated other geoengineering experiments such as the injection of stratospheric aerosols or marine sea salt [35,47,48]. Alterskjaer et al. conducted a model experiment with sea salt injection and showed that the cooling effect in mean temperature is larger over the interior of the continent than in the oceans, except for high-latitude northern hemisphere regions [47]. Volodin et al. and Ji et al. showed a larger cooling effect of stratospheric aerosol injection over tropical and mid-latitude land areas than in the corresponding oceanic areas [35,48].

Seasonality of Extreme Temperature Response
This section examines the effect of seasonality on the changes in the extreme temperature indices, TNn and TXx. To investigate the seasonality, temporal and spatial changes were analyzed and are shown in Figures 6 and 7, respectively, for East Asia in boreal summer and winter during the progress period. Although ∆TNn and ∆TXx resulted in cooling all over the area, a land-sea difference of the results is distinct in Figure 5. The G4cdnc experiment has a larger cooling effect of ∆TNn and ∆TXx over the interior of the continent compared to the ocean, with the exception of the Sea of Okhotsk. This land-sea difference in response can also be seen in previous studies that have investigated other geoengineering experiments such as the injection of stratospheric aerosols or marine sea salt [35,47,48]. Alterskjaer et al. conducted a model experiment with sea salt injection and showed that the cooling effect in mean temperature is larger over the interior of the continent than in the oceans, except for high-latitude northern hemisphere regions [47]. Volodin et al. and Ji et al. showed a larger cooling effect of stratospheric aerosol injection over tropical and mid-latitude land areas than in the corresponding oceanic areas [35,48].

Seasonality of Extreme Temperature Response
This section examines the effect of seasonality on the changes in the extreme temperature indices, TNn and TXx. To investigate the seasonality, temporal and spatial changes were analyzed and are shown in Figures 6 and 7, respectively, for East Asia in boreal summer and winter during the progress period. Figure 6 presents the time series for the entire simulation period, but now only the last 10-year averages of the progress period from 2060 to 2069 are used. In both seasons, G4cdnc resulted in cooling of the extreme temperatures. In winter, the multi-model ensemble means of ∆TNn = −1.425 • C and ∆TXx = −1.072 • C are larger than the annual response shown in Figure 4. The inter-model difference is 1.771 • C for TNn, ranging from −2.222 • C to −0.451 • C, and 0.990 • C for TXx, ranging from −1.555 • C to −0.565 • C. These inter-model differences are larger than the annual range. Unlike the winter case, G4cdnc has a larger cooling effect on TXx than on TNn in summer, having multi-model means of ∆TNn = −0.994 • C and ∆TXx = −1.119 • C. The inter-model difference is lower than the annual range: 0.362 • C for TNn, ranging from −1.356 • C to −0.994 • C, and 0.444 for TXx, ranging from −1.458 • C to −1.014 • C. In short, G4cdnc has a greater cooling effect on TNn in winter and on TXx in summer. The inter-model difference for TNn is larger than the ensemble mean annual range during winter, indicating a large uncertainty in simulating the extreme cold. Atmosphere 2020, 11, x FOR PEER REVIEW 10 of 19    Figure 5, the cooling effect caused by G4cdnc is dominant over the whole of East Asia in both winter and summer. The particularly large cooling effect over the Sea of Okhotsk shown in the annual case for TNn (Figure 5a) does not appear in summer but in winter (Figure 7a,b). An increase in sea ice lowers the air temperature, particularly the minimum temperature [35] and, according to this effect, the G4cdnc experiment demonstrated the largest cooling over the Sea of Okhotsk in winter (Figure 7a). Sea ice in the Sea of Okhotsk varies with season, reaching its maximum in March and minimum in September [49]. Many CMIP5 models projected that the extent of Arctic sea ice will decline under a GHG forcing scenario of RCP4.5 in both summer and winter seasons in the future [50]. In the meantime, a study by Moore et al. demonstrated that the G1 scenario restores the extent of Arctic sea ice from conditions under the high GHG scenario to a preindustrial state [43]. Because of the removal of snow and sea ice in summer, the TNn response to the G4cdnc experiment is weaker in summer over the East Asian continent and the Sea of Okhotsk compared to winter (Figure 7b). The sea ice over the Sea of Okhotsk also causes a larger TXx cooling compared to the rest of the ocean, but a moderate response compared to the TNn cooling in winter (Figure 7c,d).   Figure 4. The inter-model difference is 1.771 °C for TNn, ranging from −2.222 °C to −0.451 °C, and 0.990 °C for TXx, ranging from −1.555 °C to −0.565 °C. These inter-model differences are larger than the annual range. Unlike the winter case, G4cdnc has a larger cooling effect on TXx than on TNn in summer, having multimodel means of ∆TNn = −0.994 °C and ∆TXx = −1.119 °C. The inter-model difference is lower than the annual range: 0.362 °C for TNn, ranging from −1.356 °C to −0.994 °C, and 0.444 for TXx, ranging from −1.458 °C to −1.014 °C. In short, G4cdnc has a greater cooling effect on TNn in winter and on TXx in summer. The inter-model difference for TNn is larger than the ensemble mean annual range during winter, indicating a large uncertainty in simulating the extreme cold. Figure 7 presents the difference in TNn and TXx between the G4cdnc and RCP4.5 experiments over East Asia in boreal winter and summer during the progress period (2060-2069). The same as the annual case shown in Figure 5, the cooling effect caused by G4cdnc is dominant over the whole of East Asia in both winter and summer. The particularly large cooling effect over the Sea of Okhotsk shown in the annual case for TNn (Figure 5a) does not appear in summer but in winter (Figure 7a,b). An increase in sea ice lowers the air temperature, particularly the minimum temperature [35] and, according to this effect, the G4cdnc experiment demonstrated the largest cooling over the Sea of Okhotsk in winter (Figure 7a). Sea ice in the Sea of Okhotsk varies with season, reaching its maximum in March and minimum in September [49]. Many CMIP5 models projected that the extent of Arctic sea ice will decline under a GHG forcing scenario of RCP4.5 in both summer and winter seasons in the future [50]. In the meantime, a study by Moore et al. demonstrated that the G1 scenario restores the extent of Arctic sea ice from conditions under the high GHG scenario to a preindustrial state [43]. Because of the removal of snow and sea ice in summer, the TNn response to the G4cdnc experiment is weaker in summer over the East Asian continent and the Sea of Okhotsk compared to winter The land-sea difference revealed in the annual response of the extreme temperature indices in Figure 5 can also be found in the seasonal responses of TNn, especially in winter. Compared to the oceanic response, the G4cdnc experiment results in a larger cooling effect over the East Asian continent, except for Southern China in winter and summer (Figure 7a,c). This continental region with a larger cooling effect is frequently covered by snow during winter, suggesting snow-albedo feedback. A previous study showed that G3 and G4 geoengineering would restore the snow cover that is reduced by GHG forcing, which decreases the surface-sensible heat flux, raises heat capacity, and thus reduced the sensitivity of temperature to radiative forcing changes [35,51]. Figure 8 shows the difference in snow cover and cloudiness between the G4cdnc and RCP4.5 experiments over East Asia in winter during the progress period, but it has been drawn using only two model outputs (CanESM2 and MIROC-ESM) currently available from the data archive. The pattern of snow cover change with a decrease approximately north of 40 • N and an increase in the south is consistent with that of the upward shortwave radiation change at the TOA shown in Figure 2a. This result illustrates positive feedback between snow and albedo. The net incoming solar radiation increases over the northern area and it would warm the atmosphere. Then, the resultant decrease in vertical temperature gradient near the surface can cause a decrease in latent and sensible heat fluxes from the surface, followed by a decrease in cloudiness ( Figure 8b) and snowfall, daytime warming, and nighttime cooling of the surface after all. This feedback process can possibly be related to the greater cooling of TNn (Figure 7a) over the northern area with the decreased shortwave flux and the smaller cooling of TXx over the area (Figure 7c). Such positive snow-albedo feedback can be applied to the southern area with the increased upward shortwave radiation with the opposite signs as well.
heat fluxes from the surface, followed by a decrease in cloudiness ( Figure 8b) and snowfall, daytime warming, and nighttime cooling of the surface after all. This feedback process can possibly be related to the greater cooling of TNn (Figure 7a) over the northern area with the decreased shortwave flux and the smaller cooling of TXx over the area (Figure 7c). Such positive snow-albedo feedback can be applied to the southern area with the increased upward shortwave radiation with the opposite signs as well.

Relative Effect of Extreme Temperature
We calculated annual and monthly efficacies based on root mean square (RMS) differences to investigate the relative effect on extreme temperatures of the G4cdnc experiment over East Asia. The annual efficacies of extreme temperatures (TNn and TXx) and mean temperatures (TAS, TASmin,

Relative Effect of Extreme Temperature
We calculated annual and monthly efficacies based on root mean square (RMS) differences to investigate the relative effect on extreme temperatures of the G4cdnc experiment over East Asia. The annual efficacies of extreme temperatures (TNn and TXx) and mean temperatures (TAS, TASmin, and TASmax) are shown in Figure 9, and seasonal variations in efficacy are displayed using monthly values in Figure 10. In this section, we refer to the efficacy in terms of percentages.
Atmosphere 2020, 11, x FOR PEER REVIEW 13 of 19 and TASmax) are shown in Figure 9, and seasonal variations in efficacy are displayed using monthly values in Figure 10. In this section, we refer to the efficacy in terms of percentages.    The annual efficacies of G4cdnc for decreasing the warming of TNn and TXx under the RCP4.5 conditions in the multi-model ensemble mean over East Asia are 56.6% and 63.4%, respectively. In other words, the efficacy of TXx is larger than that of TNn by 7%, which implies that marine cloud brightening is a relatively effective method of alleviating the warming of TXx compared to TNn under GHG forcing. This asymmetric relative effect on extreme temperatures is the opposite to that represented by the absolute effect on extreme temperatures (∆TNn ≥ ∆TXx in Section 3.3). However, Curry et al., for G1 over the globe, showed that the difference between the efficacies of TNn and TXx is not significant (TNn: 87% and TXx: 89%) [33]. Therefore, the relative effect of marine cloud brightening over East Asia has a different asymmetry from that reported by the G1 experiment. Meanwhile, the annual efficacy for mean temperature over the progress period by the G4cdnc simulation is 60.9% for TAS, 61% for TASmin, and 60.4% for TASmax, which implies that G4cndc is similarly effective in alleviating the warming of the mean temperature. In addition, the efficacies for all mean temperatures are larger than that for TNn but lower than that for TXx.
The monthly efficacy was investigated for extreme temperatures during the progress period ( Figure 10). The efficacy of TXx was larger than that of TNn during all the months, except from The annual efficacies of G4cdnc for decreasing the warming of TNn and TXx under the RCP4.5 conditions in the multi-model ensemble mean over East Asia are 56.6% and 63.4%, respectively. In other words, the efficacy of TXx is larger than that of TNn by 7%, which implies that marine cloud brightening is a relatively effective method of alleviating the warming of TXx compared to TNn under GHG forcing. This asymmetric relative effect on extreme temperatures is the opposite to that represented by the absolute effect on extreme temperatures (∆TNn ≥ ∆TXx in Section 3.3). However, Curry et al., for G1 over the globe, showed that the difference between the efficacies of TNn and TXx is not significant (TNn: 87% and TXx: 89%) [33]. Therefore, the relative effect of marine cloud brightening over East Asia has a different asymmetry from that reported by the G1 experiment. Meanwhile, the annual efficacy for mean temperature over the progress period by the G4cdnc simulation is 60.9% for TAS, 61% for TASmin, and 60.4% for TASmax, which implies that G4cndc is similarly effective in alleviating the warming of the mean temperature. In addition, the efficacies for all mean temperatures are larger than that for TNn but lower than that for TXx.
The monthly efficacy was investigated for extreme temperatures during the progress period ( Figure 10). The efficacy of TXx was larger than that of TNn during all the months, except from August to October. The relative effect of TXx reached its maxima of 75% in December to March and its minima of 55% in April and September to November. The efficacy, however, rebounded to 65% from May to July, despite this being during the minima period. The relative effect of TXx shown by the G4cdnc experiment demonstrated a seasonal response; it reached maxima during winter and minima during spring and fall. On the other hand, the relative effect of TNn reached its maxima of 65% in February, March, August, and December and its minima of 50% in the other months. Hence, the relative effect of TNn revealed its double maximum peaks in February and August and its minimum in May, showing relatively low efficacies in spring and fall. In short, the relative effect of extreme temperatures produced by the G4cdnc experiment demonstrates seasonal dependency over East Asia, and the range between the maximum and minimum efficacies reaches up to 30%.

Durability of Effects after the Termination of G4cdnc
The last 20-year simulation period (2071 to 2090) highlights the consequences of terminating marine cloud brightening after the restoration of CDNC, and the termination period has been specified as the latter 10 years (2081-2090) in this study. Here, we analyze the changes in TNn and TXx indices over East Asia during the termination period, whereas previous studies have researched changes in mean temperatures [4,47,52] after impeding other geoengineering experiments over the globe. The multi-model ensemble mean of differences in extreme temperatures between the G4cdnc and RCP4.5 experiments over East Asia during the termination period shown in Figure 4 are ∆TNn = −0.321 • C for TNn with an inter-model range of −0.366 • C to −0.255 • C and ∆TXx = −0.278 • C for TXx with a range of −0.417 • C to −0.013 • C. Although the 10-year averaged values of TNn and TXx are still lower than the values under the RCP4.5 experiment, the year-to-year values asymptote to the RCP4.5 level till 2090, the last simulation year of the experiment. This termination effect on extreme temperatures was examined in a previous study that showed a 10th percentile global temperature increase of 0.61 • C and 90th percentile increase of 0.59 • C after shutting down injections of sea salt into the tropical marine atmospheric boundary layer [34]. In the MIROC-ESM model, which did not simulate the termination effect, the values of ∆TNn and ∆TXx remained at around −1 • C during the termination period. Figure 11 shows the difference pattern of extreme temperatures between the G4cdnc and RCP4.5 experiments over East Asia averaged during the termination period. Both ∆TNn and ∆TXx increased and approached the RCP4.5 results compared to the differences obtained from the progress period. Other studies also presented an increase in extreme temperature over the globe after terminating the forcing of their geoengineering methods [34,47]. It is noteworthy that while the cloud brightening control was implemented in the G4cdnc experiment, the extreme temperatures remained cooler than the RCP4.5 experiment results over the whole East Asian region, but after the termination of the control showing restoration of the experimental forcing CDNC, they increased and, in some regions, showed higher temperature values than the RCP4.5 results. These regionally different responses can cause a political problem between the countries taking control and other countries not participating in the control experiment. Finally, Figure 12 shows the monthly efficacy of the G4cdnc experiment for extreme temperatures over East Asia during the termination period. As the G4cndc is terminated, the relative effects on both TNn and TXx also decreased significantly. Moreover, the relative effect on TNn during the termination period is larger than TXx in all months except for February, March, and November. For TNn, the efficacy reached its maxima of 15% on average in March and November and its minima of 5% in the rest of the months. TXx reached its maxima of 20% on average in May to August and December and minima of 5% in January to April. In short, the relative effect on extreme temperature during the termination period also shows a seasonal variation but with much weaker seasonality than that during the progress period. Finally, Figure 12 shows the monthly efficacy of the G4cdnc experiment for extreme temperatures over East Asia during the termination period. As the G4cndc is terminated, the relative effects on both TNn and TXx also decreased significantly. Moreover, the relative effect on TNn during the termination period is larger than TXx in all months except for February, March, and November. For TNn, the efficacy reached its maxima of 15% on average in March and November and its minima of 5% in the rest of the months. TXx reached its maxima of 20% on average in May to August and December and minima of 5% in January to April. In short, the relative effect on extreme temperature during the termination period also shows a seasonal variation but with much weaker seasonality than that during the progress period. the termination period is larger than TXx in all months except for February, March, and November. For TNn, the efficacy reached its maxima of 15% on average in March and November and its minima of 5% in the rest of the months. TXx reached its maxima of 20% on average in May to August and December and minima of 5% in January to April. In short, the relative effect on extreme temperature during the termination period also shows a seasonal variation but with much weaker seasonality than that during the progress period.

Summary and Conclusions
We have investigated the marine cloud brightening effects on extreme temperature indices, TNn and TXx, over East Asia by the GeoMIP Phase 2 G4cdnc experiment under the RCP4.5 warming scenario. For this study, we used five climate model datasets that participated in the G4cdnc experiment and analyzed their multi-model ensemble mean and the individual model results.
Through a brief analysis of the radiation fluxes at the TOA, our results confirmed that the marine cloud brightening could increase the reflection of shortwave radiation at the TOA and decrease the outgoing longwave radiation, resulting in a net radiation loss over ocean and a gain over land, largely

Summary and Conclusions
We have investigated the marine cloud brightening effects on extreme temperature indices, TNn and TXx, over East Asia by the GeoMIP Phase 2 G4cdnc experiment under the RCP4.5 warming scenario. For this study, we used five climate model datasets that participated in the G4cdnc experiment and analyzed their multi-model ensemble mean and the individual model results.
Through a brief analysis of the radiation fluxes at the TOA, our results confirmed that the marine cloud brightening could increase the reflection of shortwave radiation at the TOA and decrease the outgoing longwave radiation, resulting in a net radiation loss over ocean and a gain over land, largely compensating the loss in East Asia. Our results also confirmed the sea-ice albedo feedback over ocean, especially over the Sea of Okhotsk and the snow albedo feedback could be a crucial factor determining the radiation response at the TOA in East Asia.
The G4cdnc experiment shows the capability of counteracting the warming of extreme temperatures as well as the mean temperature under the RCP4.5 scenario, though the amount was insufficient to remove the warming to the level of the base period (2006-2019). The efficacy of the cloud control for mean temperature reaches approximately 60%, which is a middle of the efficacies for extreme temperatures: 56% for TNn and 63% for TXx annually. The efficacies for the extreme temperatures revealed a significant seasonality, showing the highest efficacy in winter and relatively low efficacies in spring and fall during the progress period (2060-2069).
Although the cloud control was implemented over ocean only, the changes in the extreme temperatures were larger over land than over ocean in East Asia. The considerable cooling of TNn over the Sea of Okhotsk in winter implies the restoration of the extent of sea ice by the G4cdnc experiment. These findings need to be highlighted since the extent of sea ice has a great influence on the wintertime temperature in mid-latitude land areas of the Northern Hemisphere. It is also probable that the significant cooling of TNn in winter can be related to the increase in wintertime cold events.
In the context that extreme cold events as well as hot spells have occurred more frequently under global warming, how these findings can be used to act on severe weather events is an important issue to be answered prior to any implementation of the cloud control experiment in reality. Thus, more in-depth analysis on the extreme values not only for temperature but also other variables is required on a regional scale. The authors expect this study can be a trigger for these further studies.