Abrupt reduction in shipping emission as an inadvertent geoengineering termination shock produces substantial radiative warming.

Human activities affect the Earth's climate through modifying the composition of the atmosphere, which then creates radiative forcing that drives climate change. The warming effect of anthropogenic greenhouse gases has been partially balanced by the cooling effect of anthropogenic aerosols. In 2020, fuel regulations abruptly reduced the emission of sulfur dioxide from international shipping by about 80% and created an inadvertent geoengineering termination shock with global impact. Here we estimate the regulation leads to a radiative forcing of +0.2±0.11Wm-2 averaged over the global ocean. The amount of radiative forcing could lead to a doubling (or more) of the warming rate in the 2020 s compared with the rate since 1980 with strong spatiotemporal heterogeneity. The warming effect is consistent with the recent observed strong warming in 2023 and expected to make the 2020 s anomalously warm. The forcing is equivalent in magnitude to 80% of the measured increase in planetary heat uptake since 2020. The radiative forcing also has strong hemispheric contrast, which has important implications for precipitation pattern changes. Our result suggests marine cloud brightening may be a viable geoengineering method in temporarily cooling the climate that has its unique challenges due to inherent spatiotemporal heterogeneity.

Observational evidence of strong forcing from aerosol effect on low cloud coverage.

Aerosols cool Earth's climate indirectly by increasing low cloud brightness and their coverage (Cf), constituting the aerosol indirect forcing (AIF). The forcing partially offsets the greenhouse warming and positively correlates with the climate sensitivity. However, it remains highly uncertain. Here, we show direct observational evidence for strong forcing from Cf adjustment to increased aerosols and weak forcing from cloud liquid water path adjustment. We estimate that the Cf adjustment drives between 52% and 300% of additional forcing to the Twomey effect over the ocean and a total AIF of -1.1 ± 0.8 W m-2. The Cf adjustment follows a power law as a function of background cloud droplet number concentration, Nd. It thus depends on time and location and is stronger when Nd is low. Cf only increases substantially when background clouds start to drizzle, suggesting a role for aerosol-precipitation interactions. Our findings highlight the Cf adjustment as the key process for reducing the uncertainty of AIF and thus future climate projections.

Global reduction in ship-tracks from sulfur regulations for shipping fuel.

Ship-tracks are produced by ship-emitted aerosols interacting with low clouds. Here, we apply deep learning models on satellite data to produce the first global climatology map of ship-tracks. We show that ship-tracks are at the nexus of cloud physics, maritime shipping, and fuel regulation. Our map captures major shipping lanes while missing others because of background conditions. Ship-track frequency is more than 10 times higher than a previous survey, and its interannual fluctuations reflect variations in cross-ocean trade, shipping activity, and fuel regulations. Fuel regulation can alter both detected frequency and shipping routes due to cost. The 2020 fuel regulation, together with the coronavirus disease 2019 pandemic, reduced ship-track frequency to its lowest level in recent decades across the globe and may have ushered in an era of low frequency. The regulation reduces the aerosol indirect forcing from ship emissions by 46% or between 0.02 and 0.27 W m-2 given its current estimates.

Exploring aerosols near clouds with high-spatial-resolution aircraft remote sensing during SEAC4RS.

Since aerosols are important to our climate system, we seek to observe the variability of aerosol properties within cloud systems. When applied to the satellite-borne Moderate-resolution Imaging Spectroradiometer (MODIS), the Dark Target (DT) retrieval algorithm provides global aerosol optical depth (AOD at 0.55 µm) in cloud-free scenes. Since MODIS' resolution (500 m pixels, 3 km or 10 km product) is too coarse for studying near-cloud aerosol, we ported the DT algorithm to the high-resolution (~50 m pixels) enhanced-MODIS Airborne Simulator (eMAS), which flew on the high-altitude ER-2 during the Studies of Emissions, Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC4RS) Airborne Science Campaign over the U.S. in 2013. We find that even with aggressive cloud screening, the ~0.5 km eMAS retrievals show enhanced AOD, especially within 6 km of a detected cloud. To determine the cause of the enhanced AOD, we analyze additional eMAS products (cloud retrievals and degraded-resolution AOD), co-registered Cloud Physics Lidar (CPL) profiles, MODIS aerosol retrievals, and ground-based Aerosol Robotic Network (AERONET) observations. We also define spatial metrics to indicate local cloud distributions near each retrieval, and then separate into near-cloud and far-from-cloud environments. The comparisons show that low cloud masking is robust, and unscreened thin cirrus would have only a small impact on retrieved AOD. Some of the enhancement is consistent with clear-cloud transition zone microphysics such as aerosol swelling. However, 3D radiation interaction between clouds and the surrounding clear air appears to be the primary cause of the high AOD near clouds.

Cloud Products from the Earth Polychromatic Imaging Camera (EPIC): Algorithms and Initial Evaluation.

This paper presents the physical basis of the EPIC cloud product algorithms and an initial evaluation of their performance. Since June 2015, EPIC has been providing observations of the sunlit side of the Earth with its 10 spectral channels ranging from the UV to the near-IR. A suite of algorithms has been developed to generate the standard EPIC Level 2 Cloud Products that include cloud mask, cloud effective pressure/height, and cloud optical thickness. The EPIC cloud mask adopts the threshold method and utilizes multichannel observations and ratios as tests. Cloud effective pressure/height is derived with observations from the O2 A-band (780 nm and 764 nm), and B-band (680 nm and 688 nm) pairs. The EPIC cloud optical thickness retrieval adopts a single channel approach where the 780 nm and 680 nm channels are used for retrievals over ocean and over land, respectively. Comparison with co-located cloud retrievals from geosynchronous earth orbit (GEO) and low earth orbit (LEO) satellites shows that the EPIC cloud product algorithms are performing well and are consistent with theoretical expectations. These products are publicly available at the Atmospheric Science Data Center at the NASA Langley Research Center for climate studies and for generating other geophysical products that require cloud properties as input.

Observations of Local Positive Low Cloud Feedback Patterns and Their Role in Internal Variability and Climate Sensitivity.

Modeling studies have shown that cloud feedbacks are sensitive to the spatial pattern of sea surface temperature (SST) anomalies, while cloud feedbacks themselves strongly influence the magnitude of SST anomalies. Observational counterparts to such patterned interactions are still needed. Here we show that distinct large-scale patterns of SST and low-cloud cover (LCC) emerge naturally from objective analyses of observations and demonstrate their close coupling in a positive local SST-LCC feedback loop that may be important for both internal variability and climate change. The two patterns that explain the maximum amount of covariance between SST and LCC correspond to the Interdecadal Pacific Oscillation and the Atlantic Multidecadal Oscillation, leading modes of multidecadal internal variability. Spatial patterns and time series of SST and LCC anomalies associated with both modes point to a strong positive local SST-LCC feedback. In many current climate models, our analyses suggest that SST-LCC feedback strength is too weak compared to observations. Modeled local SST-LCC feedback strength affects simulated internal variability so that stronger feedback produces more intense and more realistic patterns of internal variability. To the extent that the physics of the local positive SST-LCC feedback inferred from observed climate variability applies to future greenhouse warming, we anticipate significant amount of delayed warming because of SST-LCC feedback when anthropogenic SST warming eventually overwhelm the effects of internal variability that may mute anthropogenic warming over parts of the ocean. We postulate that many climate models may be underestimating both future warming and the magnitude of modeled internal variability because of their weak SST-LCC feedback.

Biomass smoke from southern Africa can significantly enhance the brightness of stratocumulus over the southeastern Atlantic Ocean.

Marine stratocumulus clouds cover nearly one-quarter of the ocean surface and thus play an extremely important role in determining the global radiative balance. The semipermanent marine stratocumulus deck over the southeastern Atlantic Ocean is of particular interest, because of its interactions with seasonal biomass burning aerosols that are emitted in southern Africa. Understanding the impacts of biomass burning aerosols on stratocumulus clouds and the implications for regional and global radiative balance is still very limited. Previous studies have focused on assessing the magnitude of the warming caused by solar scattering and absorption by biomass burning aerosols over stratocumulus (the direct radiative effect) or cloud adjustments to the direct radiative effect (the semidirect effect). Here, using a nested modeling approach in conjunction with observations from multiple satellites, we demonstrate that cloud condensation nuclei activated from biomass burning aerosols entrained into the stratocumulus (the microphysical effect) can play a dominant role in determining the total radiative forcing at the top of the atmosphere, compared with their direct and semidirect radiative effects. Biomass burning aerosols over the region and period with heavy loadings can cause a substantial cooling (daily mean -8.05 W m-2), primarily as a result of clouds brightening by reducing the cloud droplet size (the Twomey effect) and secondarily through modulating the diurnal cycle of cloud liquid water path and coverage (the cloud lifetime effect). Our results highlight the importance of realistically representing the interactions of stratocumulus with biomass burning aerosols in global climate models in this region.

Differences in liquid cloud droplet effective radius and number concentration estimates between MODIS Collections 5.1 and 6 over global oceans.

Differences in cloud droplet effective radius and cloud droplet number concentration (CDNC) estimates inferred from the Aqua MODIS Collections 5.1 and 6 cloud products (MYD06) are examined for warm clouds over global oceans for the year 2008. Individual pixel level retrievals for both collections are aggregated to 1° × 1° and compared globally and regionally for the three main spectral channel pairs used for MODIS cloud optical property retrievals. Comparisons between both collections are performed for cases in which all three effective radii retrievals are classified by the MODIS Cloud Product as valid. The contribution to the observed differences of several key MYD06 Collection 6 algorithm updates are also explored, with a focus on changes to the surface reflectance model, assumed solar irradiance, above cloud emission, cloud top pressure, and pixel registration. Global results show a neutral to positive (> 50 cm-3) change for C6-derived CDNC relative to C5.1 for the 1.6 µm and 2.1 µm channel retrievals, corresponding to a neutral to -2 µm difference in droplet effective radius. For 3.7 µm retrievals, CDNC results show a negative change in the tropics, with differences transitioning toward positive values with increasing latitude spanning -25 to +50 cm-3 related to a +2.5 to -1 µm transition in effective radius. Cloud optical thickness differences were small relative to effective radius, and found to not significantly impact CDNC estimates. Regionally, the magnitude and behavior of the annual CDNC cycle are compared for each effective radius retrieval. Results from this study indicate significant intercollection differences in aggregated values of effective radius due to changes to the pre-computed retrieval lookup tables for ocean scenes, changes to retrieved cloud top pressure, solar irradiance, or above cloud thermal emission, depending upon spectral channel. The observed differences between collections may have implications for existing MODIS derived climatologies and validation studies of effective radius and CDNC.

The MODIS cloud optical and microphysical products: Collection 6 updates and examples from Terra and Aqua.

The MODIS Level-2 cloud product (Earth Science Data Set names MOD06 and MYD06 for Terra and Aqua MODIS, respectively) provides pixel-level retrievals of cloud-top properties (day and night pressure, temperature, and height) and cloud optical properties (optical thickness, effective particle radius, and water path for both liquid water and ice cloud thermodynamic phases-daytime only). Collection 6 (C6) reprocessing of the product was completed in May 2014 and March 2015 for MODIS Aqua and Terra, respectively. Here we provide an overview of major C6 optical property algorithm changes relative to the previous Collection 5 (C5) product. Notable C6 optical and microphysical algorithm changes include: (i) new ice cloud optical property models and a more extensive cloud radiative transfer code lookup table (LUT) approach, (ii) improvement in the skill of the shortwave-derived cloud thermodynamic phase, (iii) separate cloud effective radius retrieval datasets for each spectral combination used in previous collections, (iv) separate retrievals for partly cloudy pixels and those associated with cloud edges, (v) failure metrics that provide diagnostic information for pixels having observations that fall outside the LUT solution space, and (vi) enhanced pixel-level retrieval uncertainty calculations. The C6 algorithm changes collectively can result in significant changes relative to C5, though the magnitude depends on the dataset and the pixel's retrieval location in the cloud parameter space. Example Level-2 granule and Level-3 gridded dataset differences between the two collections are shown. While the emphasis is on the suite of cloud optical property datasets, other MODIS cloud datasets are discussed when relevant.

Seasonally Transported Aerosol Layers over Southeast Atlantic are Closer to Underlying Clouds than Previously Reported.

From June to October, low-level clouds in the Southeast (SE) Atlantic often underlie seasonal aerosol layers transported from African continent. Previously, the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) 532 nm lidar observations have been used to estimate the relative vertical location of the above-cloud aerosols (ACA) to the underlying clouds. Here, we show new observations from NASA's Cloud-Aerosol Transport System (CATS) lidar. Two seasons of CATS 1064 nm observations reveal that the bottom of the ACA layer is much lower than previously estimated based on CALIPSO 532nm observations. For about 60% of CATS nighttime ACA scenes, the aerosol layer base is within 360 m distance to the top of the underlying cloud. Our results are important for future studies of the microphysical indirect and semi-direct effects of ACA in the SE Atlantic region.

Validation of quasi-invariant ice cloud radiative quantities with MODIS satellite-based cloud property retrievals.

Similarity relations applied to ice cloud radiance calculations are theoretically analyzed and numerically validated. If τ(1-ϖ) and τ(1-ϖg) are conserved where τ is optical thickness, ϖ the single-scattering albedo, and g the asymmetry factor, it is possible that substantially different phase functions may give rise to similar radiances in both conservative and non-conservative scattering cases, particularly in the case of large optical thicknesses. In addition to theoretical analysis, this study uses operational ice cloud optical thickness retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2 Collection 5 (C5) and Collection 6 (C6) cloud property products to verify radiative similarity relations. It is found that, if the MODIS C5 and C6 ice cloud optical thickness values are multiplied by their respective (1-ϖg) factors, the resultant products referred to as the effective optical thicknesses become similar with their ratio values around unity. Furthermore, the ratios of the C5 and C6 ice cloud effective optical thicknesses display an angular variation pattern similar to that of the corresponding ice cloud phase function ratios. The MODIS C5 and C6 values of ice cloud similarity parameter, defined as [(1-ϖ)/(1-ϖg)]1/2, also tend to be similar.

Ice cloud backscatter study and comparison with CALIPSO and MODIS satellite data.

An invariant imbedding T-matrix (II-TM) method is used to calculate the single-scattering properties of 8-column aggregate ice crystals. The II-TM based backscatter values are compared with those calculated by the improved geometric-optics method (IGOM) to refine the backscattering properties of the ice cloud radiative model used in the MODIS Collection 6 cloud optical property product. The integrated attenuated backscatter-to-cloud optical depth (IAB-ICOD) relation is derived from simulations using a CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite) lidar simulator based on a Monte Carlo radiative transfer model. By comparing the simulation results and co-located CALIPSO and MODIS (Moderate Resolution Imaging Spectroradiometer) observations, the non-uniform zonal distribution of ice clouds over ocean is characterized in terms of a mixture of smooth and rough ice particles. The percentage of the smooth particles is approximately 6% and 9% for tropical and midlatitude ice clouds, respectively.

Uncertainties in cloud phase and optical thickness retrievals from the Earth Polychromatic Imaging Camera (EPIC).

This paper presents an investigation of the expected uncertainties of a single channel cloud optical thickness (COT) retrieval technique, as well as a simple cloud temperature threshold based thermodynamic phase approach, in support of the Deep Space Climate Observatory (DSCOVR) mission. DSCOVR cloud products will be derived from Earth Polychromatic Imaging Camera (EPIC) observations in the ultraviolet and visible spectra. Since EPIC is not equipped with a spectral channel in the shortwave or mid-wave infrared that is sensitive to cloud effective radius (CER), COT will be inferred from a single visible channel with the assumption of appropriate CER values for liquid and ice phase clouds. One month of Aqua MODIS daytime granules from April 2005 is selected for investigating cloud phase sensitivity, and a subset of these granules that has similar EPIC sun-view geometry is selected for investigating COT uncertainties. EPIC COT retrievals are simulated with the same algorithm as the operational MODIS cloud products (MOD06), except using fixed phase-dependent CER values. Uncertainty estimates are derived by comparing the single channel COT retrievals with the baseline bi-spectral MODIS retrievals. Results show that a single channel COT retrieval is feasible for EPIC. For ice clouds, single channel retrieval errors are minimal (< 2%) due to the particle size insensitivity of the assumed ice crystal (i.e., severely roughened aggregate of hexagonal columns) scattering properties at visible wavelengths, while for liquid clouds the error is mostly limited to within 10%, although for thin clouds (COT < 2) the error can be higher. Potential uncertainties in EPIC cloud masking and cloud temperature retrievals are not considered in this study.

Retrieval of ice cloud properties using an optimal estimation algorithm and MODIS infrared observations. Part I: Forward model, error analysis, and information content.

An optimal estimation (OE) retrieval method is developed to infer three ice cloud properties simultaneously: optical thickness (τ), effective radius (reff ), and cloud-top height (h). This method is based on a fast radiative transfer (RT) model and infrared (IR) observations from the MODerate resolution Imaging Spectroradiometer (MODIS). This study conducts thorough error and information content analyses to understand the error propagation and performance of retrievals from various MODIS band combinations under different cloud/atmosphere states. Specifically, the algorithm takes into account four error sources: measurement uncertainty, fast RT model uncertainty, uncertainties in ancillary datasets (e.g., atmospheric state), and assumed ice crystal habit uncertainties. It is found that the ancillary and ice crystal habit error sources dominate the MODIS IR retrieval uncertainty and cannot be ignored. The information content analysis shows that, for a given ice cloud, the use of four MODIS IR observations is sufficient to retrieve the three cloud properties. However, the selection of MODIS IR bands that provide the most information and their order of importance varies with both the ice cloud properties and the ambient atmospheric and the surface states. As a result, this study suggests the inclusion of all MODIS IR bands in practice since little a priori information is available.

Positive Low Cloud and Dust Feedbacks Amplify Tropical North Atlantic Multidecadal Oscillation.

The Atlantic Multidecadal Oscillation (AMO) is characterized by a horseshoe pattern of sea surface temperature (SST) anomalies and has a wide range of climatic impacts. While the tropical arm of AMO is responsible for many of these impacts, it is either too weak or completely absent in many climate model simulations. Here we show, using both observational and model evidence, that the radiative effect of positive low cloud and dust feedbacks is strong enough to generate the tropical arm of AMO, with the low cloud feedback more dominant. The feedbacks can be understood in a consistent dynamical framework: weakened tropical trade wind speed in response to a warm middle latitude SST anomaly reduces dust loading and low cloud fraction over the tropical Atlantic, which warms the tropical North Atlantic SST. Together they contribute to appearance of the tropical arm of AMO. Most current climate models miss both the critical wind speed response and two positive feedbacks though realistic simulations of them may be essential for many climatic studies related to the AMO.

MODIS Collection 6 shortwave-derived cloud phase classification algorithm and comparisons with CALIOP.

Cloud thermodynamic phase (ice, liquid, undetermined) classification is an important first step for cloud retrievals from passive sensors such as MODIS (Moderate-Resolution Imaging Spectroradiometer). Because ice and liquid phase clouds have very different scattering and absorbing properties, an incorrect cloud phase decision can lead to substantial errors in the cloud optical and microphysical property products such as cloud optical thickness or effective particle radius. Furthermore, it is well established that ice and liquid clouds have different impacts on the Earth's energy budget and hydrological cycle, thus accurately monitoring the spatial and temporal distribution of these clouds is of continued importance. For MODIS Collection 6 (C6), the shortwave-derived cloud thermodynamic phase algorithm used by the optical and microphysical property retrievals has been completely rewritten to improve the phase discrimination skill for a variety of cloudy scenes (e.g., thin/thick clouds, over ocean/land/desert/snow/ice surface, etc). To evaluate the performance of the C6 cloud phase algorithm, extensive granule-level and global comparisons have been conducted against the heritage C5 algorithm and CALIOP. A wholesale improvement is seen for C6 compared to C5.

A framework based on 2-D Taylor expansion for quantifying the impacts of sub-pixel reflectance variance and covariance on cloud optical thickness and effective radius retrievals based on the bi-spectral method.

The bi-spectral method retrieves cloud optical thickness (τ) and cloud droplet effective radius (r e ) simultaneously from a pair of cloud reflectance observations, one in a visible or near infrared (VIS/NIR) band and the other in a shortwave-infrared (SWIR) band. A cloudy pixel is usually assumed to be horizontally homogeneous in the retrieval. Ignoring sub-pixel variations of cloud reflectances can lead to a significant bias in the retrieved τ and r e . In the literature, the retrievals of τ and r e are often assumed to be independent and considered separately when investigating the impact of sub-pixel cloud reflectance variations on the bi-spectral method. As a result, the impact on τ is contributed only by the sub-pixel variation of VIS/NIR band reflectance and the impact on r e only by the sub-pixel variation of SWIR band reflectance. In our new framework, we use the Taylor expansion of a two-variable function to understand and quantify the impacts of sub-pixel variances of VIS/NIR and SWIR cloud reflectances and their covariance on the τ and r e retrievals. This framework takes into account the fact that the retrievals are determined by both VIS/NIR and SWIR band observations in a mutually dependent way. In comparison with previous studies, it provides a more comprehensive understanding of how sub-pixel cloud reflectance variations impact the τ and r e retrievals based on the bi-spectral method. In particular, our framework provides a mathematical explanation of how the sub-pixel variation in VIS/NIR band influences the r e retrieval and why it can sometimes outweigh the influence of variations in the SWIR band and dominate the error in r e retrievals, leading to a potential contribution of positive bias to the r e retrieval. We test our framework using synthetic cloud fields from a large-eddy simulation and real observations from MODIS. The predicted results based on our framework agree very well with the numerical simulations. Our framework can be used to estimate the retrieval uncertainty from sub-pixel reflectance variations in operational satellite cloud products and to help understand the differences in τ and r e retrievals between two instruments.

Frequency and causes of failed MODIS cloud property retrievals for liquid phase clouds over global oceans.

Moderate Resolution Imaging Spectroradiometer (MODIS) retrieves cloud droplet effective radius (re ) and optical thickness (τ) by projecting observed cloud reflectances onto a precomputed look-up table (LUT). When observations fall outside of the LUT, the retrieval is considered "failed" because no combination of τ and re within the LUT can explain the observed cloud reflectances. In this study, the frequency and potential causes of failed MODIS retrievals for marine liquid phase (MLP) clouds are analyzed based on 1 year of Aqua MODIS Collection 6 products and collocated CALIOP and CloudSat observations. The retrieval based on the 0.86 µm and 2.1 µm MODIS channel combination has an overall failure rate of about 16% (10% for the 0.86 µm and 3.7 µm combination). The failure rates are lower over stratocumulus regimes and higher over the broken trade wind cumulus regimes. The leading type of failure is the "re too large" failure accounting for 60%-85% of all failed retrievals. The rest is mostly due to the "re too small" or τ retrieval failures. Enhanced retrieval failure rates are found when MLP cloud pixels are partially cloudy or have high subpixel inhomogeneity, are located at special Sun-satellite viewing geometries such as sunglint, large viewing or solar zenith angles, or cloudbow and glory angles, or are subject to cloud masking, cloud overlapping, and/or cloud phase retrieval issues. The majority (more than 84%) of failed retrievals along the CALIPSO track can be attributed to at least one or more of these potential reasons. The collocated CloudSat radar reflectivity observations reveal that the remaining failed retrievals are often precipitating. It remains an open question whether the extremely large re values observed in these clouds are the consequence of true cloud microphysics or still due to artifacts not included in this study.