Abstract
Clouds affected by solar eclipses could influence the reflection of sunlight back into space and might change local precipitation patterns. Satellite cloud retrievals have so far not taken into account the lunar shadow, hindering a reliable spaceborne assessment of the eclipse-induced cloud evolution. Here we use satellite cloud measurements during three solar eclipses between 2005 and 2016 that have been corrected for the partial lunar shadow together with large-eddy simulations to analyze the eclipse-induced cloud evolution. Our corrected data reveal that, over cooling land surfaces, shallow cumulus clouds start to disappear at very small solar obscurations (~15%). Our simulations explain that the cloud response was delayed and was initiated at even smaller solar obscurations. We demonstrate that neglecting the disappearance of clouds during a solar eclipse could lead to a considerable overestimation of the eclipse-related reduction of net incoming solar radiation. These findings should spur cloud model simulations of the direct consequences of sunlight-intercepting geoengineering proposals, for which our results serve as a unique benchmark.
Similar content being viewed by others
Introduction
Blocking part of the solar radiation incident on the Earth’s (lower) atmosphere and surface is one of the proposed strategies to counteract the current and future global temperature rise, which may be inevitable if climate change mitigation efforts prove to be insufficient1,2,3,4. This type of (solar) geoengineering is based on placing sun shields or reflecting particles in space between the Earth and the Sun5,6,7,8,9,10, or on the injection of aerosols into the stratosphere11,12. General circulation models (GCMs) suggest that an insolation reduction of 3.5–5.0% can largely undo the global temperature rise and intensified hydrological cycle associated with a quadrupled pre-industrial CO2 concentration13,14,15,16. However, those GCMs also predict latitudinal variations in temperature response, and an extra reduction of precipitation in the tropics. Moreover, although clouds play a vital role in the Earth’s radiation balance17, the impact of solar dimming on clouds is still poorly understood18,19,20,21. GCMs modeling the response to (extraterrestrial) dimming of sunlight are based on idealized scenarios22. They inherently suffer from uncertainties23, mainly focus on the long-term impact, and highly parameterize short-term and small-scale processes such as cloud formation20,23. Twice a year on average, for a few hours the opportunity arises to take measurements of the Earth experiencing gradual insolation reductions from 0 to nearly 100%, during the partial phase of a solar eclipse. Although the time scales involved with solar geoengineering will most likely not be equivalent to those of solar eclipses, these measurements can help to better understand (and test models that predict) the immediate cloud response to the deployment of sunlight intercepting material.
Ground-based meteorological observations during solar eclipses have primarily focused on fast drops in air temperature, winds and turbulence, and on changing (photo)chemistry24,25. Weather publications contain anecdotal descriptions of dissipating low-level cumulus clouds right before totality, while mid- and high-level clouds survived26,27. The timing of the eclipse-induced effect on cumulus clouds is difficult to quantify with ground-based measurements due to the chaotic nature of cumulus cloud evolution and the missing observation of the non-eclipse state, and atmospheric model studies of eclipses did not yet analyze the cumulus cloud sensitivity28,29,30,31. Additionally, the precise locations of the affected clouds can be hard to predict and are sometimes inaccessible for ground-based observers. Earth observation satellites in geostationary orbit can continuously monitor clouds in large geographical areas32,33,34 and hints of dissipating cumulus clouds have been observed by comparing satellite images before and after a total solar eclipse, without estimating the non-eclipse state26,35. However, during solar eclipses, satellite images show spatio-temporally varying darkening35,36 and satellite retrievals of cloud cover and cloud optical thickness (COT) are biased. This bias is caused by not taking into account the insolation reduction in the calculation of the top-of-atmosphere (TOA) reflectance from which cloud properties are derived. Hence, it has remained unknown how fast clouds are modified by various solar obscuration fractions.
Here, we present geostationary satellite measurements of clouds during three solar eclipses between 2005 and 2016 that have been corrected for the insolation reduction. Our corrected measurements reveal that shallow cumulus clouds start to dissipate at a solar obscuration of ~15% over cooling land surfaces, which would have been hidden in the partial lunar shadow without insolation reduction correction. Using large-eddy simulations we explain the timing of the cloud dissipation and demonstrate that the rising air parcels in the atmospheric boundary layer are already affected by the solar eclipse at even smaller obscurations. We calculate that neglecting the dissipating cloud behavior in the simulations would result in an overestimation of 20 W m−2 of the eclipse-related reduction of net incoming shortwave radiation at TOA. Finally, we discuss that the high cloud sensitivity to rather small insolation reductions should spur cloud model simulations of the short-term impact of sunlight-intercepting geoengineering concepts.
Results
Satellite observations
Figure 1 shows the uncorrected and corrected TOA visible (VIS) reflectance over East Africa and part of the Indian Ocean, at four subsequent hours during the annular solar eclipse on 3 October 2005, obtained by the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) instrument (see Methods). In the first hour, cumulus clouds were present over land. They are better visible in the corrected high resolution images of the COT zoomed in on the specific area over land at 3∘–7∘ latitude and 27∘–31∘ longitude (Fig. 1c), which we call the study area. The derived COT in the cloudy pixels was ≲5 and SEVIRI estimated a cloud top altitude of ~2 km (see Supplementary Fig. 1), indicating that they were shallow and located in the planetary boundary layer. Similar large-scale daytime shallow cumulus cloud fields over land can be found in central Africa and the Amazonian rainforest throughout the year, and in northeast America and Siberia during the boreal summer37,38. When the obscuration increased in the second hour, the cloud cover in the study area diminished, although remnants of the cloud pattern of the previous hour are still recognizable. In the third hour, the insolation increased again, but the shallow cumulus clouds stayed away while the clouds that survived had grown in size and COT. It was only during the final stage of the eclipse that shallow cumulus clouds returned throughout the study area. Over ocean, the cumulus clouds did not disappear. In Supplementary Figs. 2 and 3, we provide two more examples of vanishing shallow cumulus clouds over land during solar eclipses.
The disappearance of shallow cumulus clouds only occurred on the day of the solar eclipse. Figure 2a shows the time series of the cloud cover in the study area, compared to that during 11 comparable days which we selected based on a similar type of cloud pattern before the eclipse started (see Methods). Also depicted in Fig. 2 is the obscuration fraction on the eclipse day. The time series shows that the increasing cloud cover in the morning already halted at low obscuration fractions (~0.15), happening at around 09:30 UTC which was 30 minutes after the start of the eclipse. Secondly, there was a ~50 min time lag with respect to the instant of maximum obscuration at 10:52 UTC before the clouds started to return. During the increase in cloud cover between 12:00 and 12:30 UTC, the mean COT in the cloudy pixels decreased (see Fig. 2b), which can be attributed to the contribution of the newly formed shallow clouds. This decrease is absent in the time series of the comparable days.
Because shallow cumulus clouds in the boundary layer are generated by rising thermals originating from air close to the surface39, we collected land surface temperature (LST) measurements from space by SEVIRI, derived from infrared radiation emitted by the shallow land surface layer (see Methods). Figure 2c shows the spatially averaged LST in the study area on the eclipse day and the comparable days. The maximum LST on the comparable days was two hours delayed with respect to local noon at 10:04 UTC, which could possibly be explained by the smaller heat flux into the ground due to the warmed subsurface layer in the afternoon40. On the eclipse day, the LST drops instantly with the obscuration fraction, due to the direct response of the shallow land surface layer temperature to net radiation forcing40. We estimate a maximum LST drop of 5.8 K induced by the eclipse at 11.00 UTC (see Methods). Comparable fast drops in satellite LST measurements have been found in a study over Europe during the total solar eclipse of 20 March 2015 by ref. 41 who showed dependencies of the drop magnitude on the eclipse duration and time of the day (earlier eclipses gave larger drops), vegetation, surface height and distance to the coast. We did not detect a time lag in the LST minimum with respect to maximum solar obscuration: a time lag of ~1.5 min as reported in literature41 is not resolved in our LST data at 15 minutes intervals. Hence, the measured time lag of ~50 min before clouds return in Fig. 2a cannot be explained by a time lag in the LST. Over ocean, we found no sea surface temperature drop when the eclipse passed (see Supplementary Fig. 4), due to the large heat capacity of water and the efficient heat transport from the sea surface to deeper water layers through turbulent mixing40.
Large-eddy simulations
In order to explain the land-cloud interaction in the study area, we simulated the evolution of shallow cumulus clouds during a solar eclipse with the Dutch Atmospheric Large-Eddy Simulation (DALES) model42. Figures 2d and e show the time series of the simulated cloud cover and mean COT, respectively, as would be measured from space (see Methods), and Supplementary Fig. 5 contains snapshots of the spatially resolved cloud fields. We present the results for the solar eclipse case, using the measured LST as input, and a reference case without eclipse-induced LST drop (see Methods). Indeed, in the solar eclipse case, our simulations show a substantial decrease in cloud cover with respect to the reference case. The cloud cover already differed ~15 to ~20 min after the start of the eclipse when the obscuration was still smaller than 10%. As in the observations (Fig. 2a), there is a time lag in the instant of minimum cloud cover with respect to mid-eclipse, after which the cloud cover rapidly increases. We note that the simulated cloud cover right after the eclipse is even larger than in the reference case, while the simulated COT is lower. A larger cloud cover after the eclipse compared to the hypothetical non-eclipse scenario is difficult to prove with our observations, due to the large variability of the cloud cover on the comparable days in the afternoon (cf. Fig. 2a). During the rapid increase in cloud cover, the simulated COT also increases, which is not in agreement with the observations, but can be attributed to the absence of the deeper convective clouds in the south part of the study area due to the horizontally averaged input settings in our simulations (see Methods).
The disappearing clouds during a solar eclipse can be explained by the drop in sensible (thermal) and latent (moisture) heat fluxes from the surface to the lowest atmosphere layer (Fig. 2f), as a result of the dropping LST (see Eqs. (1) and (2)). Those heat fluxes drive the buoyancy flux (Fig. 3b, f) of relatively warm and moist air parcels from the surface, through the well-mixed boundary layer, up to the level of minimum buoyancy flux (LMBF) where the parcels are capped by a temperature inversion (see Supplementary Fig. 6). The rising parcels are cooled through adiabatic expansion which increases the parcel relatively humidity (RH) up to 100% at the lifting condensation level (LCL) where shallow cumulus clouds are formed that can extend to higher altitudes (Fig. 3a, e). During a solar eclipse, this process is suppressed, which is clear from the overall drop in buoyancy flux. The drop in surface and buoyancy fluxes is consistent with the diminished boundary layer turbulence during solar eclipses found in other studies31,43,44. It should be noted that the LMBF is still higher in our model than the LCL, but with smaller upward parcel velocities (Fig. 3c) which are controlled by the surface buoyancy flux, fewer parcels reach the LCL.
The ~15–20 min delay of the simulated vanishing cloud cover with respect to the dropping LST can be related to the ~16–24 min travel time around 09:00 UTC of the fastest rising parcels from the lowest atmospheric layer to the cloud top as shown in Fig. 3d, after which the first individual clouds could fully disappear. We note that the vertical updraft velocity of the lowest layer responded within 5 min to differences in surface fluxes. The travel time depends on the cloud top height and vertical updraft velocities (see Fig. 3c, g and Methods). Thus, the cloud response was initiated when the parcels affected by the eclipse started rising, at even smaller obscuration fractions than at which this response could be observed. The ~18 min time lag of the simulated cloud cover minimum with respect to mid-eclipse can be related to the relatively long travel time around 11:10 UTC of ~13–19 min to LCL at which the newly formed clouds started influencing the cloud cover. Furthermore, we note that the LCL and cloud base remained slightly lower after the eclipse compared to the reference case, due to the continued colder air near the surface (see Supplementary Fig. 7 and Eq. (4)). In Fig. 4 we provide a conceptual model in which we summarize the most important processes responsible for the shallow cumulus cloud behavior during a solar eclipse.
The disappearance of the shallow cumuli during the solar eclipse has a notable feedback on the solar radiative fluxes. This can be understood from Fig. 5, which shows the simulated reflected and net incoming shortwave (SW) radiative flux at TOA, for both the reference and solar eclipse cases. As a result of the solar eclipse the net incoming SW flux started to decrease at 09:00 UTC. However, this change, which has solar dimming as its main cause, is also affected by the decrease of SW radiation reflected back into space due to the clearance of the sky. Indeed, the latter indirect effect causes an opposing increase in the net incoming SW flux at TOA. Neglecting the solar eclipse-induced cloud disappearance in our simulations (as illustrated by dotted line in Fig. 5), resulted in an overestimation of 20 W m−2 of the eclipse-related reduction of net incoming SW flux at TOA at 11:22 UTC. We note that this error would further increase with longer time lags of the cloud return with respect to mid-eclipse, such as found in the satellite observations (Fig. 2a), because then more sunlight illuminates the cloud-free scenes.
Discussion
The observed response of shallow cumulus clouds to a solar eclipse at already ~15% obscuration, initiated at even smaller obscurations due to the parcel travel time, reveals the potential direct consequence of deploying sunlight intercepting material in the stratosphere or in space. We note that the duration of the cloud response is expected to depend on the speed and magnitude of the local obscuration variations, as the altered difference between the near-surface air and surface temperature, which causes the response, may possibly restore after a certain period. Diminished shallow cumulus clouds would partly oppose the objective of solar geoengineering which is to decrease the amount of net incoming solar radiation, and could prevent the growth into deeper convective and possibly precipitating clouds45. While solar geoengineering proposals aim to reduce the solar radiation reaching the (lower) atmosphere and surface globally by only a few percent (depending on the required compensation), the use of non-uniform reductions to achieve this goal could increase the locally experienced variations in obscuration10. Space-based examples are the deployment of solar reflectors in Earth orbit8,9 or in orbit around the 1st Lagrange point10, offering daily and seasonally varying shading, respectively. Injected stratospheric aerosols can also exhibit spatio-temporally varying patterns due to seasonally changing global stratospheric circulation, depending on the injection location46 and strategy employed (whether constant or step-wise)47,48. Consequently, aerosol optical depths of 0.4–0.6 could be attained locally46,47,48, causing up to ~45% of the local direct sunlight to be scattered or absorbed by the aerosols before it reaches the lower atmosphere. Our results should spur model simulations investigating the response of shallow cumulus clouds to those geoengineering concepts, particularly for scenes over land where the surface temperature can adjust quickly. Additionally, our measurements provide an opportunity to validate these models, enhancing their reliability in predicting cloud behavior under natural conditions and in a world influenced by solar geoengineering.
Methods
Cloud measurements
The primary data set used in this research consists of measurements from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) on board the Meteosat Second Generation (MSG) series of geostationary satellites operated by EUMETSAT. We used data from Meteosat-8 and Meteosat-10, for the 2005–2006 and 2016 cases, respectively. SEVIRI measures TOA radiances over the full Earth disk centered at 0∘ latitude and 0∘ longitude every 15 minutes in twelve channels across the visible and near-infrared part of the spectrum49. Eleven channels have a narrow bandwidth at low spatial resolution (3 × 3 km2 at sub-satellite point) and one channel has a broad bandwidth (0.6–0.9 μm) at high spatial resolution (1 × 1 km2), the latter being referred to as the HRV channel. Shortwave channel reflectances were obtained from the observed radiances and the calculated solar irradiance, and were calibrated with MODIS following ref. 50. For the longwave channels, the operational calibration from EUMETSAT was used.
For the surface temperature analyses at low spatial resolution (see Surface temperature measurements), cloud masks were calculated using the EUMETSAT Nowcasting and Very Short Range Forecasting Satellite Application Facility (NWC SAF) v2021 cloud algorithm51. The NWC SAF software also provided the cloud top height at low spatial resolution. For the cloud analyses at high spatial resolution, we used the low-resolution NWC SAF cloud mask as a basis and improved the spatial resolution of the cloud mask using the HRV TOA reflectances by comparing with a HRV TOA reflectance climatology, which was generated for each pixel with HRV TOA reflectance measurements in a 16-day period centered at the day of the eclipse, following the method of ref. 52. The TOA reflectances of the 0.6 and 1.6 μm channels were paired to simultaneously retrieve the COT and effective droplet radius34,53, which were downscaled to high spatial resolution using the HRV channel as described in more detail in refs. 54,55. The COT of the pixels that were not masked as cloudy was set equal to zero. For the study area (3–7∘ N latitude, 27–31∘ E longitude on 3 October 2005) and between 06:00 and 09:00 UTC, pixels with a TOA reflectance value larger than 1.0, originating from the sunglint in the rivers, were removed from the cloud product.
Solar eclipse correction
We corrected the TOA reflectance for its reduction during a solar eclipse, that was due to the ignorance of the reduced solar irradiance in its calculation, through a division by 1 − fo, with fo the solar obscuration fraction. In previous work56, we validated this approach with the TROPOMI satellite instrument, which allowed for accurate monitoring of aerosols in the partial lunar shadow up to fo = 0.92. Recently, Wen et al.57 applied a similar type of correction to images of the Earth Polychromatic Imaging Camera (EPIC) instrument on the Deep Space Climate Observatory (DSCOVR), to quantify the measured TOA reflectance error without eclipse correction of the sunlit side of the Earth during the annular solar eclipse on 21 June 2020. The value of fo is different for every pixel as it depends on measurement time, surface height, latitude and longitude, which are provided with the SEVIRI data. The so-called Besselian elements describing the temporal variation of the Moon shadow’s geometry were taken from ref. 58. The value of fo also depends on wavelength through the wavelength dependent limb darkening of the solar disk59. For the corrections of the TOA reflectances in the HRV channel and 0.6 μm, 0.8 μm and 1.6 μm channels, we used the limb darkening coefficients of Pierce et al.60 and Pierce et al.61 at 0.7 μm and the central wavelengths 0.635 μm, 0.81 μm and 1.64 μm, respectively. We refer to Trees et al.56 for more details about the solar eclipse correction of the TOA reflectance. Applying the cloud algorithms to the corrected TOA reflectance yielded the corrected cloud mask and COT.
Surface temperature measurements
For the measurements of land surface temperature (LST), we used the LST product of the Land Surface Analysis Satellite Application Facility (LSA SAF)62 derived from the 10.8 μm and 12 μm SEVIRI channels63 with an uncertainty of 1 to 2 K64 which is expected to be stable within the time scales of solar eclipses41. We computed the horizontal average LST in the study area in every 15 minutes time step, after replacing the cloudy pixels by interpolated nearest neighbor values using the corrected SEVIRI NWC SAF cloud mask (see Cloud algorithm and Solar eclipse correction). The maximum LST drop due to the eclipse was estimated with respect to the average LST of comparable days (see Selection of comparable days). For the measurements of the sea surface temperature (SST), we used the hourly SST product of the Ocean and Sea Ice Satellite Application Facility (OSI SAF)65 derived from the 10.8 μm and 12 μm SEVIRI channels with an uncertainty well below 1 K66.
Selection of comparable days
The selection of days without solar eclipse that are comparable to our study case on 3 October 2005 in East Africa was done by first selecting 100 days from September and October in 2004, 2005 and 2006 with the smallest differences with respect to our study case in the sums of the predicted sensible and latent heat fluxes by ERA567 (which does not take into account the solar eclipse effect) in the complete diurnal cycle. Subsequently, we refined the selection by removing the days for which the mean of the absolute differences in high spatial resolution cloud cover with respect to our study case between 06:00 and 09:00 UTC was larger than the threshold of 0.1. In this way, we obtained 11 days for which we find the cloud cover comparable in the morning before the lunar shadow reached the study area: 2004-09-10, 2004-10-24, 2005-09-03, 2005-09-19, 2005-09-25, 2005-09-29, 2005-10-04, 2006-09-01, 2006-09-04, 2006-09-08, 2006-10-06. We visually inspected the HRV TOA reflectance and COT for those comparable days and indeed found between 06:00 and 09:00 UTC similar types of shallow cumulus clouds throughout the scene, with low COTs (≲5) as shown in Fig. 2b.
Cloud simulations
The cloud simulations were performed with the Dutch Atmospheric Large-Eddy Simulation (DALES) model42. The setup was a horizontally cyclic domain of 50 by 50 km2 with horizontal cell sizes of 100 by 100 m2 and 209 vertical layers between 0 and 14 km altitude. The vertical extent of each layer was stretched by a factor 1.01 with respect to the layer just below, and it was 20 m in the lowest layer. The ground surface was assumed flat. The atmosphere variables in the domain were initialized with horizontally homogeneous vertical profiles of the liquid potential temperature, total water specific humidity, and horizontal wind speed and direction. Those profiles were the horizontal mean profiles at 02:00 UTC in the study area (3∘-7∘N latitude, 27∘-31∘E longitude on 3 October 2005), taken from ERA567. The pressure profile was determined with the use of the thermodynamic profiles, the gas law and mean hydrostatic balance. After 02:00 UTC, the simulation freely propagated those variables, until 15:00 UTC when the simulations ended. The time and altitude dependent mean horizontal advective tendencies of heat and moisture, in addition to the geostrophic winds, were all diagnosed from ERA5 and prescribed in the DALES runs. That is, we neglect the horizontal in- and outflow of eclipse-induced atmospheric changes at the boundaries of the domain, e.g. due to short-term disturbances in the horizontal pressure gradient68. At the top boundary during the simulation, we imposed large-scale subsidence, also taken from ERA5.
The surface fluxes of heat and moisture were computed according to the vertical difference of potential temperature θ and the water vapor specific humidity qv between the ground surface and the lowest atmospheric model layer (indicated by the subscript ‘bot’), respectively. With the appropriate conversion of potential temperature to temperature the sensible and latent fluxes can be expressed as
where ρ is the air density in kg m−3, cp = 1004 J kg−1 K−1 the heat capacity of dry air, g = 9.81 m s−2 the gravitational acceleration, Lv = 2.5 ⋅ 106 J kg−1 the latent heat for vaporization of water, and qsat the temperature dependent saturation specific humidity. The aerodynamic resistance coefficients for heat and moisture, ra,T and ra,q, respectively, depend on the atmospheric stability near the surface following Monin-Obukhov similarity theory42. The dependencies of the ra values on the surface roughness length and the actual soil water content, which for land conditions will typically be much lower than the saturated value qsat, were taken into account by multiplying them with constant correction factors. This calibration step was made to obtain sensible and heat fluxes that were consistent with the ECMWF model for the reference case. As a consequence, the surface heat fluxes were not parameterized as functions of net incoming radiation, to avoid unnecessary complexity and uncertainties. Instead, the surface temperature Tsfc was prescribed using the SEVIRI LST measurements in the study area (see Satellite measurements), and shifted to the ERA5 skin temperature at 02:00 UTC through a positive offset of 2.0274 K applied to the complete LST time series, for consistency with the initial atmospheric profiles. In the simulation without eclipse and until 09:00 UTC, Tsfc was identical to that in the simulation with eclipse, but after 09:00 UTC it was the average of the LST time series of the comparable days shifted to the ERA5 skin temperature at 02:00 UTC.
The simulated cloud cover in the study area was computed as follows. First, we regridded the DALES output liquid water specific humidity ql in kg kg−1 to a grid with cell sizes of 1 × 1 km2, for a fair comparison with the SEVIRI observations. Secondly, we computed the COT by evaluating the following integral from the surface (0) to TOA (zTOA)69:
where ρair is the air density in kg m−3, ρliq = 1000 kg m−3 is the density of liquid water, and reff is the droplet effective radius which we assumed to be constant and equal to 10−5 m. Thirdly, the columns (x, y) with a COT larger than 1 were flagged as a cloudy column. The cloud cover was computed as the number of cloudy columns divided by the total number of columns in the domain. For the mean COT time series, we computed the horizontal average COT considering the cloudy columns only.
The lifting condensation level (LCL) in m was computed from the DALES output of air temperature Tbot in K and relative humidity RHbot in % of the bottom atmospheric layer (at 10 meter altitude), using the suggested formula by70:
The level of minimum buoyancy flux (LMBF) was computed as the altitude at which the buoyancy flux \(\overline{{w}^{{\prime} }{\theta }_{v}^{{\prime} }}\) was minimum, where θv is the virtual potential temperature which depends on specific humidity qv in kg kg−1, temperature T in K, pressure p in Pa, reference pressure p0 = 105 Pa, cp, the gas constant for dry air Rd = 287.0 J kg−1 K−1, and the gas constant for water vapor Rv = 461.5 J kg−1 K−171:
and \(\overline{{w}^{{\prime} }{\theta }_{{{{{{{{\rm{v}}}}}}}}}^{{\prime} }}\) was defined as:
in which Nx and Ny are the number of grid cells in the horizontal x- and y directions, and the overbar indicates a horizontal mean value. The mean updraft virtual potential temperature, \({\overline{\theta }}_{{{{{{{{\rm{v,up}}}}}}}}}(z)\), was computed as the mean of θv(z) of the updrafts only (i.e., the grid points at a certain altitude z for which wij > 0). The velocities of the fast rising parcels were defined by the p-percentile velocities wp% (with p = 1, 3 and 5), which were the vertical velocities at a certain altitude z for which a percentage p of the w-distribution contained grid points with wij > wp% (see ref. 72). The travel time Δt of the fast rising parcels from the surface to a certain altitude z*, using vertical velocity wp%, was computed by numerically evaluating
with z* = zLCL for the travel time to LCL and z* = zCT for the travel time to the cloud top. The cloud top was defined as the highest altitude for which ql > 0.
DALES uses the Rapid Radiative Transfer Model for Global climate model applications (RRTMG) radiation scheme73,74 to compute the shortwave (SW) and longwave (LW) radiative fluxes through the atmosphere, emitted and reflected by the surface, and emerging from TOA during the simulation. In the solar eclipse case, the SW radiation incident at TOA was multiplied by (1 − fo), where fo was the horizontal mean obscuration fraction in the study area taken at 0.635 μm (see Solar eclipse correction). We used vertical profiles of the ozone mass mixing ratios in the study area from ERA567. The surface albedo was the horizontal mean white-sky albedo in the study area measured by MODIS75.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The data required to replicate the timeseries presented in this article, as well as the input and boundary conditions for our model simulations, are available in the public repository Zenodo (https://doi.org/10.5281/zenodo.10371414).
Code availability
The Dutch Atmospheric Large-Eddy Simulation (DALES) software used for this research is publicly available under the terms of the GNU GPL version 3 on https://github.com/dalesteam/dales.
References
Keith, D. W. Geoengineering. Nature 409, 420 (2001).
Lenton, T. M. & Vaughan, N. E. The radiative forcing potential of different climate geoengineering options. Atmos. Chem. Phys. 9, 5539–5561 (2009).
Shepherd, J. Geoengineering the climate: science, governance and uncertainty. Project Report https://eprints.soton.ac.uk/156647/ (2009).
Kosugi, T. Role of sunshades in space as a climate control option. Acta Astronautica 67, 241–253 (2010).
Early, J. T. Space-based solar shield to offset greenhouse effect. J. Brit. Interplanet. Soc. 42, 567–569 (1989).
Angel, R. Feasibility of cooling the Earth with a cloud of small spacecraft near the inner Lagrange point (L1). Proc. Natl Acad. Sci. 103, 17184–17189 (2006).
Fuglesang, C. & de Herreros Miciano, M. G. Realistic sunshade system at L1 for global temperature control. Acta Astronautica 186, 269–279 (2021).
Mautner, M. A Space-Based Solar Screen Against Climatic Warming. Bull. Am. Astron. Soc. 21, 993 (1989).
Pearson, J., Oldson, J. & Levin, E. Earth rings for planetary environment control. Acta Astron. 58, 44–57 (2006).
Sánchez, J.-P. & McInnes, C. R. Optimal Sunshade Configurations for Space-Based Geoengineering near the Sun-Earth L1 Point. PLoS One 10, e0136648 (2015).
Crutzen, P. J. Albedo Enhancement by Stratospheric Sulfur Injections: A Contribution to Resolve a Policy Dilemma? Clim. Change 77, 211–220 (2006).
Rasch, P. J. et al. An overview of geoengineering of climate using stratospheric sulphate aerosols. Philos. Trans. R. Soc. Lond. Ser. A 366, 4007–4037 (2008).
Lunt, D. J., Ridgwell, A., Valdes, P. J. & Seale, A. "Sunshade World”: A fully coupled GCM evaluation of the climatic impacts of geoengineering. Geophys. Res. Lett. 35, L12710 (2008).
Kravitz, B. et al. Climate model response from the Geoengineering Model Intercomparison Project (GeoMIP). J. Geophys. Res. 118, 8320–8332 (2013).
Bal, P. K., Pathak, R., Mishra, S. K. & Sahany, S. Effects of global warming and solar geoengineering on precipitation seasonality. Environ. Res. Lett. 14, 034011 (2019).
Irvine, P. J., Kravitz, B., Lawrence, M. G. & Muri, H. An overview of the Earth system science of solar geoengineering. WIREs Clim. Change 7, 815–833 (2016).
Ramanathan, V. et al. Cloud-Radiative Forcing and Climate: Results from the Earth Radiation Budget Experiment. Science 243, 57–63 (1989).
Schmidt, H. et al. Solar irradiance reduction to counteract radiative forcing from a quadrupling of CO2: climate responses simulated by four earth system models. Earth Syst. Dyn. 3, 63–78 (2012).
Russotto, R. D. & Ackerman, T. P. Changes in clouds and thermodynamics under solar geoengineering and implications for required solar reduction. Atmos. Chem. Phys. 18, 11905–11925 (2018).
Kravitz, B. et al. Comparing different generations of idealized solar geoengineering simulations in the Geoengineering Model Intercomparison Project (GeoMIP). Atmos. Chem. Phys. 21, 4231–4247 (2021).
Virgin, J. G. & Fletcher, C. G. On the Linearity of External Forcing Response in Solar Geoengineering Experiments. Geophys. Res. Lett. 49, e00200 (2022).
Kravitz, B. et al. The Geoengineering Model Intercomparison Project (GeoMIP). Atmos. Sci. Lett. 12, 162–167 (2011).
Irvine, P. J. et al. Key factors governing uncertainty in the response to sunshade geoengineering from a comparison of the GeoMIP ensemble and a perturbed parameter ensemble. J. Geophys. Res. Atmos. 119, 7946–7962 (2014).
Aplin, K. L., Scott, C. J. & Gray, S. L. Atmospheric changes from solar eclipses. Philos. Trans. R. Soc. Lond. Ser. A 374, 20150217 (2016).
Harrison, R. G. & Hanna, E. The solar eclipse: a natural meteorological experiment. Philos. Trans. R. Soc. Lond. Ser. A 374, 20150225 (2016).
Anderson, J. Meteorological changes during a solar eclipse. Weather 54, 207–215 (1999).
Hanna, E. Meteorological effects of the solar eclipse of 11 August 1999. Weather 55, 430–446 (2000).
Founda, D. et al. The effect of the total solar eclipse of 29 March 2006 on meteorological variables in Greece. Atmos. Chem. Phys. 7, 5543–5553 (2007).
Montornès, A., Codina, B., Zack, J. W. & Sola, Y. Implementation of bessel’s method for solar eclipses prediction in the wrf-arw model. Atmos. Chem. Phys. 16, 5949–5967 (2016).
Clark, P. A. Numerical simulations of the impact of the 20 March 2015 eclipse on UK weather. Philos. Trans. R. Soc. Lond. Ser. A 374, 20150218 (2016).
Buban, M. S., Lee, T. R., Dumas, E. J., Baker, C. B. & Heuer, M. Observations and Numerical Simulation of the Effects of the 21 August 2017 North American Total Solar Eclipse on Surface Conditions and Atmospheric Boundary-Layer Evolution. Boundary Layer Meteorol. 171, 257–270 (2019).
Rossow, W. B. Measuring Cloud Properties from Space: A Review. J. Clim. 2, 201–213 (1989).
Stengel, M. S. et al. CLAAS: the CM SAF cloud property data set using SEVIRI. Atmos. Chem. Phys. 14, 4297–4311 (2014).
Benas, N. et al. The MSG-SEVIRI-based cloud property data record CLAAS-2. Earth Syst. Sci. Data 9, 415–434 (2017).
Gerth, J. J. Shining light on sky cover during a total solar eclipse. J. Appl. Remote Sens. 12, 020501 (2018).
Peñaloza-Murillo, M. A. & Pasachoff, J. M. Cloudiness and Solar Radiation During the Longest Total Solar Eclipse of the 21st Century at Tianhuangping (Zhejiang), China. J. Geophys. Res. Atmos. 123, 13,443–13,461 (2018).
Rabin, R. M. & Martin, D. W. Satellite observations of shallow cumulus coverage over the central United States: An exploration of land use impact on cloud cover. J. Geophys. Res. 101, 7149–7155 (1996).
Heiblum, R. H., Koren, I. & Feingold, G. On the link between Amazonian forest properties and shallow cumulus cloud fields. Atmos. Chem. Phys. 14, 6063–6074 (2014).
Siebesma, A. P.Shallow Cumulus Convection, 441–486 (Springer Netherlands, Dordrecht, 1998). https://doi.org/10.1007/978-94-011-5058-3_19.
Stull, R. B.An introduction to boundary layer meteorology (Kluwer Academic Publishers, Dordrecht, Boston and London, 1988).
Good, E. Satellite observations of surface temperature during the March 2015 total solar eclipse. Philos. Trans. R. Soc. Lond. Ser. A 374, 20150219 (2016).
Heus, T. et al. Formulation of the dutch atmospheric large-eddy simulation (dales) and overview of its applications. Geosci. Model Dev. 3, 415–444 (2010).
Eaton, F. D. et al. Solar Eclipse Effects Observed in the Planetary Boundary Layer Over a Desert. Boundary Layer Meteorol. 83, 331–346 (1997).
Mauder, M., Desjardins, R. L., Oncley, S. P. & MacPherson, I. Atmospheric Response to a Partial Solar Eclipse over a Cotton Field in Central California. J. Appl. Meteorol. Climatol. 46, 1792–1803 (2007).
Grabowski, W. W. et al. Daytime convective development over land: A model intercomparison based on LBA observations. Q. J. R. Meteorol. Soc. 132, 317–344 (2006).
Tilmes, S. et al. Sensitivity of Aerosol Distribution and Climate Response to Stratospheric SO2 Injection Locations. J. Geophys. Res. Atmos. 122, 12,591–12,615 (2017).
Visioni, D. et al. Seasonal Injection Strategies for Stratospheric Aerosol Geoengineering. Geophys. Res. Lett. 46, 7790–7799 (2019).
Lee, W. R., MacMartin, D. G., Visioni, D. & Kravitz, B. High Latitude Stratospheric Aerosol Geoengineering Can Be More Effective if Injection Is Limited to Spring. Geophys. Res. Lett. 48, e92696 (2021).
Schmetz, J. et al. An Introduction to Meteosat Second Generation (MSG). Bull. Am. Meteorol. Soc. 83, 977–992 (2002).
Meirink, J. F., Roebeling, R. A. & Stammes, P. Inter-calibration of polar imager solar channels using SEVIRI. Atmos. Measur. Tech. 6, 2495–2508 (2013).
Kerdraon, G. & Fontaine, E. Algorithm theoretical basis document for the cloud product processors of the nwc/geo (geo-cma-v5.1 (nwc-009), geo-ct-v4.1 (nwc-016), geo-ctth-v4.1 (nwc-017) and geo-cmic-v2.1 (nwc-021). Technical Report NWC/CDOP3/GEO/MFL/SCI/ATBD/Cloud, Issue 1, Rev 0.1, Météo-France / Centre d’études en Météorologie Satellitaire (2021).
Bley, S. & Deneke, H. A threshold-based cloud mask for the high-resolution visible channel of Meteosat Second Generation SEVIRI. Atmos. Measur. Tech. 6, 2713–2723 (2013).
Nakajima, T. & King, M. D. Determination of the optical thickness and effective particle radius of clouds from reflected solar radiation measurements. I - Theory. J. Atmos. Sci. 47, 1878–1893 (1990).
Deneke, H. et al. Increasing the spatial resolution of cloud property retrievals from Meteosat SEVIRI by use of its high-resolution visible channel: Implementation and examples. Atmos. Measur. Tech. 14, 5107–5126 (2021).
Deneke, H. M. & Roebeling, R. A. Downscaling of Meteosat SEVIRI 0.6 and 0.8 μm channel radiances utilizing the high-resolution visible channel. Atmos. Chem. Phys. 10, 9761–9772 (2010).
Trees, V., Wang, P. & Stammes, P. Restoring the top-of-atmosphere reflectance during solar eclipses: a proof of concept with the UV absorbing aerosol index measured by TROPOMI. Atmos. Chem. Phys. 21, 8593–8614 (2021).
Wen, G., Marshak, A., Herman, J. & Wu, D. Reduction of spectral radiance reflectance during the annular solar eclipse of 21 june 2020 observed by epic. Front. Remote Sens. 3, https://www.frontiersin.org/articles/10.3389/frsen.2022.777314 (2022).
Espenak, F. & Meeus, J. Five millennium canon of solar eclipses : -1999 to +3000 (2000 BCE to 3000 CE) (NASA Technical Publication, NASA Center for AeroSpace Information, TP-2006-214141, Hannover, Maryland, USA, 2006).
Koepke, P., Reuder, J. & Schween, J. Spectral variation of the solar radiation during an eclipse. Meteorologische Zeitschrift 10, 179–186 (2001).
Pierce, A. K. & Slaughter, C. D. Solar limb darkening. I: lambda lambda (3033 - 7297). Solar Phys. 51, 25–41 (1977).
Pierce, A. K., Slaughter, C. D. & Weinberger, D. Solar limb darkening in the interval 7404 - 24018 A˚, II. Solar Phys. 52, 179–189 (1977).
Trigo, I. et al. Algorithm theoretical basis document for land surface temperature (lst). Technical Report SAF/LAND/IM/ATBD_MLST/1.2, Issue 2, Land surface analysis satellite application facility (LSA-SAF, 2017).
Wan, Z. & Dozier, J. A generalized split-window algorithm for retrieving land-surface temperature from space. IEEE Trans. Geosci. Remote Sens. 34, 892–905 (1996).
Freitas, S. C., Trigo, I. F., Bioucas-Dias, J. M. & Gottsche, F.-M. Quantifying the Uncertainty of Land Surface Temperature Retrievals From SEVIRI/Meteosat. IEEE Trans. Geosci. Remote Sens. 48, 523–534 (2010).
Saux-Picart, S. Algorithm theoretical basis document for msg/seviri sea surface temperature data record, version 1.3. Technical Report SAF/OSI/CDOP2/MF/SCI/MA/256, Ocean and Sea Ice Satellite Application Facility (OSI SAF, 2018).
Saux Picart, S. et al. A Sea Surface Temperature data record (2004-2012) from Meteosat Second Generation satellites. Remote Sens. Environ. 240, 111687 (2020).
Hersbach, H. et al. The era5 global reanalysis. Q. J. R. Meteorol. Soc. 146, 1999–2049 (2020).
Gross, P. & Hense, A. Effects of a Total Solar Eclipse on the Mesoscale Atmospheric Circulation over Europe - A Model Experiment. Meteorol. Atmos. Phys. 71, 229–242 (1999).
Stephens, G. L. Radiation Profiles in Extended Water Clouds. II: Parameterization Schemes. J. Atmos. Sci. 35, 2123–2132 (1978).
Lawrence, M. G. The Relationship between Relative Humidity and the Dewpoint Temperature in Moist Air: A Simple Conversion and Applications. Bull. Am. Meteorol. Soc. 86, 225–233 (2005).
Clouds and Climate: Climate Science’s Greatest Challenge (Cambridge University Press, 2020).
Siebesma, A. P., Soares, P. M. M. & Teixeira, J. A Combined Eddy-Diffusivity Mass-Flux Approach for the Convective Boundary Layer. J. Atmos. Sci. 64, 1230 (2007).
Iacono, M. J. et al. Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models. J. Geophys. Res. Atmos. 113, D13103 (2008).
Blossey, P. N. et al. Marine low cloud sensitivity to an idealized climate change: The CGILS LES intercomparison. J. Adv. Model. Earth Syst. 5, 234–258 (2013).
Schaaf, C. & Wang, Z. Mcd43c3 modis/terra+aqua brdf/albedo albedo daily l3 global 0.05deg cmg v006 data set. NASA EOSDIS Land Processes DAAC, accessed 2023-06-24 https://doi.org/10.5067/MODIS/MCD43C3.006 (2015).
Acknowledgements
The authors thank Richard Bintanja (KNMI) for the discussions during manuscript preparation and David Donovan (KNMI) for reading the manuscript before submission. V.T. acknowledges funding by the User Support Programme Space Research (GO), project number ALWGO.2018.016, of the Dutch Research Council (NWO). S.d.R. acknowledges funding by the Refreeze the Arctic Foundation for research on geoengineering techniques at TU Delft.
Author information
Authors and Affiliations
Contributions
V.T. wrote the manuscript with contributions from S.d.R., J.W. and J.F.M. V.T. computed the solar obscuration fractions. S.d.R. did the DALES computations. J.W. and J.F.M. prepared the SEVIRI observations and J.W. converted the observations to high spatial resolution. V.T., S.d.R., J.W., J.F.M., P.W., P.S. and A.P.S. reviewed the manuscript and were involved in the analysis and selection of the presented results.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Communications Earth & Environment thanks Alexander Marshak and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Sylvia Sullivan and Aliénor Lavergne. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Trees, V.J.H., de Roode, S.R., Wiltink, J.I. et al. Clouds dissipate quickly during solar eclipses as the land surface cools. Commun Earth Environ 5, 71 (2024). https://doi.org/10.1038/s43247-024-01213-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s43247-024-01213-0
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.