Articles | Volume 25, issue 13
https://doi.org/10.5194/acp-25-6957-2025
https://doi.org/10.5194/acp-25-6957-2025
Research article
 | 
08 Jul 2025
Research article |  | 08 Jul 2025

Analysis of ship emission effects on clouds over the southeastern Atlantic using geostationary satellite observations

Nikos Benas, Jan Fokke Meirink, Rob Roebeling, and Martin Stengel
Abstract

This study investigates the impact of ship emissions on clouds over a shipping corridor in the southeastern Atlantic, employing geostationary-based observations, which have not been previously used in studies of this kind. Based on CLAAS-3, the 20-year (2004–2023) CLoud property dAtA set using SEVIRI (the geostationary Spinning Enhanced Visible and InfraRed Imager), the diurnal, seasonal and long-term corridor effects on clouds are examined. Results show a significant impact of ship emissions on cloud microphysics, consistent with the Twomey effect: an increase in cloud droplet number concentration (Nd) and a decrease in effective radius (re). Additionally, cloud liquid water path (W) decreases, though changes in cloud fraction are more subtle. No clear impact on cloud optical thickness is found, implying an overall minor radiative effect of the ship emissions, although methodological limitations to detect changes in the corridor cannot be excluded. Seasonal and diurnal variations of the impact are evident, influenced by regional conditions and by the cloud thinning during the day, respectively. The long-term analysis reveals a weakening of the shipping corridor effect on Nd and re, presumably following the International Maritime Organization's 2020 stricter regulations on sulfur emissions, and broader regional changes in W and cloud fraction, associated with sea surface temperature variations. Focusing on a climatically important cloud regime, and including novel aspects, namely the diurnal and full seasonal cycle analyses, this study highlights the advantages and potential of geostationary satellite-based cloud observations for studying aerosol–cloud interactions.

Share
1 Introduction

Clouds are a major regulator of the Earth's energy balance, modulating the amount of solar radiation reflected back to space and the amount of outgoing thermal radiation. They are also a critical component of the water cycle, redistributing water across the planet through advection. A large part of the global effect of clouds on the atmospheric radiation budget is attributed to low warm clouds, in particular stratocumulus clouds, which form extensive decks that cover large parts of the subtropical oceans. These clouds induce a cooling effect by reflecting solar radiation back to space and a warming effect by preventing thermal radiation from escaping. Since they have a similar temperature as the underlying surface, their longwave radiative effect is rather small, and their overall impact on the Earth's climate is a cooling effect (Wood, 2012).

Since the net radiative effect of stratocumulus clouds is dominated by the shortwave component, it depends on their albedo and consequently on their microphysical properties. These, in turn, vary depending on the availability of atmospheric aerosols that act as cloud condensation nuclei. However, the ways in which aerosols interact with clouds remain highly uncertain (Bellouin et al., 2020). The instantaneous, microphysical change that fine-mode aerosols induce in liquid clouds is an increase in the cloud droplet number concentration (Nd) and a decrease in their effective radius (re). This process, known as the Twomey or cloud albedo effect (Twomey, 1974), occurs on a timescale of minutes to a few hours (Gryspeerdt et al., 2021). It leads to an increase in cloud albedo and consequently to an increase in radiation reflected back to space. However, subsequent adjustments can take two different paths that lead to either positive (warming) or negative (cooling) radiative effects. In the first case, the smaller droplets give rise to enhanced entrainment and evaporation, leading to less liquid water in the cloud and smaller cloud fractions, thus reducing the cloud albedo (Bretherton et al., 2007). In the second case, smaller droplets lead to suppression of precipitation and longer cloud lifetime, thus increasing the cloud fraction and albedo (Albrecht, 1989; Christensen et al., 2020). Both mechanisms have been documented and studied extensively (see, e.g., the review study of Bellouin et al., 2020). They depend on the local meteorological conditions and occur on timescales of several hours (Feingold et al., 2024; Gryspeerdt et al., 2021).

In the study of aerosol interactions with clouds, ship emissions have been used extensively in the past (see, e.g., Christensen et al., 2022, and references therein). They can be considered localized sources of aerosols in otherwise undisturbed (or uniformly disturbed) environments, constituting good examples of opportunities to investigate the aerosol effects on clouds. In recent years, this research area has gained momentum for two additional reasons:

  1. In 2020, the International Maritime Organization (IMO) of the United Nations implemented new regulations to limit the use of sulfur in ship fuels (IMO, 2019). Since then, various studies have documented the consequences of these new regulations on ship emissions and their effects on clouds and radiative forcing (Gettelman et al., 2024; Yuan et al., 2024; Diamond, 2023; Yuan et al., 2022; Watson-Parris et al., 2022).

  2. As the effects of global warming become increasingly evident, there are also discussions to consider deliberate interventions in the climate system as a means to gain time in the effort to lower CO2 emissions. One possible mechanism for this, usually called “marine cloud brightening”, would be to increase the cloud albedo by injecting aerosol particles in marine clouds in order to increase their albedo and the radiation they reflect back to space (Feingold et al., 2024).

Relevant studies can be broadly categorized into campaign-based, which rely on in situ measurements from dedicated flights, and satellite-based, which use retrievals of cloud properties from satellite sensors to investigate aerosol–cloud interactions caused by ship emissions. The satellite-based approach in studying ship emission effects on clouds can be further categorized into ship track and shipping corridor studies. Individual ship tracks have been analyzed, starting already with the first satellite images in the 1960s. Until recently, however, difficulties in their identification were a limiting factor in estimating their global effects on clouds (Christensen et al., 2022). This is gradually changing, with new and more sophisticated methods increasing the number of detected ship tracks and improving the relevant global estimates of these effects (Yuan et al., 2023; Manshausen et al., 2022; Watson-Parris et al., 2022). Shipping corridor studies avoid individual ship track detection limitations by focusing on areas where ship traffic and its emissions are quasi-continuous. For example, Peters et al. (2011, 2014) examined three shipping corridors based on satellite and climate model output data. With the winds blowing mainly perpendicular to the corridors, the regions upstream and downstream of each corridor were investigated assuming clear and polluted conditions, respectively. No statistically significant effects were found, attributed to the large meteorological variability and the unknown cloud conditions in the absence of anthropogenic emissions. Diamond et al. (2020), on the other hand, found significant effects of ship emissions on cloud Nd and re, examining a corridor in the southeast (SE) Atlantic Ocean, where winds typically blow parallel to the corridor and thus confine ship emissions. Examining the same corridor, Hu et al. (2021) also found significant positive effects of ship emissions on cloud albedo and fraction under low background Nd.

The use of geostationary data to analyze ship tracks has been so far limited to ship track identification based on reflectances (level 1 data; Larson et al., 2022; Schreier et al., 2010). Other studies combined level 1 data from geostationary orbiters with cloud retrievals from polar orbiters (Goren and Rosenfeld, 2012, 2015). Cloud properties from geostationary satellites, based on recalibrated level 1 data, have not been previously used to analyze ship tracks. Instead, cloud properties from polar-orbiting sensors have been used so far. Most studies of ship emission effects on clouds have used cloud properties retrieved from the MODerate resolution Imaging Spectroradiometer (MODIS) (e.g., Yuan et al., 2023, 2022; Diamond, 2023; Manshausen et al., 2022; Gryspeerdt et al., 2021; Diamond et al., 2020). Contrary to geostationary-based data, the use of cloud properties from polar orbiters provides only two overpass times during the day, thus limiting the possibility of a detailed diurnal variation analysis.

In this study, we analyze the effects of ship emissions on clouds focusing on the shipping corridor that crosses the SE Atlantic Ocean, mainly connecting Europe with countries in southern and SE Asia and vice versa. The suitability of this corridor in terms of local conditions (winds typically parallel to the corridor and an overlying persistent stratocumulus deck) for similar analyses has already been shown (Diamond et al., 2020; Hu et al., 2021; Diamond, 2023). For the first time in such a study, we use cloud properties retrieved from a geostationary sensor, the Spinning Enhanced Visible and InfraRed Imager (SEVIRI), namely the latest, third, edition of the Cloud property dAtAset using SEVIRI (CLAAS-3; Benas et al., 2023; Meirink et al., 2022), which is provided by the Satellite Application Facility on Climate Monitoring (CM SAF) of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). The use of geostationary-based cloud retrievals from CLAAS-3 allows the shipping corridor effect on clouds to be analyzed on a range of temporal scales, including diurnal variations. Compared to polar orbiters, which typically observe a region once or twice per day, the CLAAS-3 high frequency of observations from the same sensor improves the assessment of these diurnal variations substantially. Taking advantage of the 20-year long (2004–2023) CLAAS-3 data record, we also estimate long-term changes in the corridor effect, focusing especially on the implications of the stricter regulations implemented by IMO in 2020.

The paper is structured as follows: in Sect. 2 we describe the data and the methodology used in the analysis – the CLAAS-3 cloud data; the ship density data for the identification of the SE Atlantic shipping corridor; and the approaches used to identify the corridor, estimate the effects on clouds, and propagate relevant uncertainties when calculating spatial and temporal averages. Results are presented in Sect. 3, separated based on temporal resolution: time series averages, seasonal cycles, diurnal cycles and long-term changes. A summary with main conclusions is given in Sect. 4.

2 Data and methodology

2.1 The CLAAS-3 cloud data set

CLAAS-3 provides detailed information on various cloud properties, offering high spatial and temporal resolution (3 km × 3 km at nadir, every 15 min) and covering the period 2004–present. As mentioned before, CLAAS-3 retrievals are based on observations from SEVIRI, which flies on board geostationary Meteosat Second Generation (MSG) satellites Meteosat-8, Meteosat-9, Meteosat-10 and Meteosat-11, covering the period 2004–present. For the CLAAS-3 processing, SEVIRI solar channel observations were calibrated based on the latest (collection 6.1) Aqua MODIS reflectances, following the approach described in Meirink et al. (2013). Extended to include all four MSG satellites, this methodology ensures a temporal stability of CLAAS-3 that is suitable for long-term analyses (see, e.g., the CLAAS-3 validation report, available in Meirink et al., 2022). The instantaneous (level 2) retrievals, which are based on the calibrated (level 1) reflectances, are aggregated into (level 3) daily and monthly averages at 0.05° × 0.05° and monthly diurnal averages at 0.25° × 0.25°. The areas covered include Africa, Europe, the Atlantic Ocean, parts of South America, the Indian Ocean and the Middle East. CLAAS-3 includes cloud fractional coverage (fc), cloud phase (liquid and ice), cloud top height, pressure and temperature, cloud optical thickness (τ), effective radius and water path (W), and Nd. Details on the retrieval algorithms and the validation of CLAAS-3 are given in Benas et al. (2023) and in dedicated technical documentation available in Meirink et al. (2022).

Regarding the cloud properties examined in this study, Nd is retrieved from re and τ, assuming an idealized stratiform boundary layer cloud, as described in Bennartz and Rausch (2017), while re and τ are retrieved based on a look-up table (LUT) of pre-calculated reflectances (Roebeling et al., 2006; Nakajima and King, 1990). W is also calculated from re and τ, assuming vertically homogeneous water content (e.g., Stephens, 1978). Since the 0.6 µm and the 3.9 µm channels are used for the retrieval of re and τ, it follows that re, τ, W and Nd are only available during daytime, which is defined as solar zenith angle values lower than 75°. For consistency, fc is also provided in daytime time slots (fc,day) in addition to the all-day retrievals fc.

All data used in this study come from level 3 monthly averages, either of all time slots or diurnally resolved, and refer to liquid clouds only, as aggregated from level 2. Thus, the successful exclusion of ice clouds from this analysis depends on the performance of the CLAAS-3 cloud phase retrieval algorithm. Evaluation results of CLAAS-3 cloud phase, available in Meirink et al. (2022), show good performance and very good agreement with reference data sets over the wider SE Atlantic region.

Occurrences of cloud edges and broken clouds, which are associated with retrieval biases in re and τ (and consequently W and Nd), can also affect level 3 averages. In CLAAS-3 level 3 aggregation algorithms, level 2 results of re and τ falling outside the LUT, which are usually associated with such cases, are excluded from level 3 averages. While this process does not ensure complete removal of the mentioned retrieval biases, it is expected that their impact on level 3 values will be minor. Additionally, different biases in and outside the shipping corridor are not expected.

2.2 Identification of the shipping corridor

For the identification of the shipping corridor we use shipping traffic density data available from a collaboration of the World Bank with the IMO. This data set is based on all hourly ship positions received by the Automatic Identification System (AIS) between January 2015 and February 2021 (Cerdeiro et al., 2020). The data represent the total number of AIS positions reported by ships in grid cells of dimensions 0.005° × 0.005° (approximately 500 × 500 m at the Equator).

For the present study, the original ship density values are aggregated to a coarser resolution of 0.25° using 50 × 50 grid cell blocks. The threshold for ship density classification is determined by calculating the mean plus 1 standard deviation of the ship density values within the wider SE Atlantic region shown in Fig. 1. Using this threshold, the ship density values are converted to a binary flag, with a flag value of 1 for cells with ship density values greater than the threshold and 0 for cells with values below the threshold. Finally, the flag is downscaled to the original CLAAS-3 level 3 resolution of 0.05°, in order to determine the corridor center in CLAAS-3 coordinates. This is done by selecting, in each latitudinal row, the central grid cell from the ones marked as “corridor”, i.e., having a flag value of 1.

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f01

Figure 1Shipping corridors crossing the southeastern Atlantic, derived based on data from the Automatic Identification System and scaled to the CLAAS-3 resolution, as described in the text. The dark-blue regions on the right correspond to high ship density near the African coast. The red rectangle denotes the focus region of this study.

While the threshold used is arbitrary, it is efficient in distinguishing the major shipping corridors from lower “background” ship density values, as Fig. 1 shows. Two corridors are identified in the study area (20–10° S, 10° W–10° E), which is shown in red. We focus on the denser one, which produces notable differences in the Nd and re values (Fig. 2a and b). This corridor is also further away from the coast, hence least subject to the influence of terrestrial emission sources.

2.3 Estimation of the shipping corridor effect

The methodology used for the quantification of the shipping corridor effect on cloud properties takes advantage of the good alignment of prevailing winds in the region with the corridor orientation. This alignment constrains the emissions and renders their effects on clouds more pronounced (Diamond et al., 2020).

To assess how shipping corridor emissions impact cloud properties, we first calculate the angle between the shipping corridor and the south-to-north direction. Then, for each point positioned at the corridor's center, we determine the line perpendicular to the corridor, based on this angle. We then identify the grid cells along this line and calculate their distances from the corridor center (measured in kilometers). This process yields cloud property data variations with distance from the corridor's center on both sides. By averaging these values along the corridor line, we create a distribution centered on the midpoint of the corridor. This distribution provides insights into how cloud properties change relative to their distance from the corridor center on both sides.

The underlying assumption in this approach is that the corridor effect will manifest as a deviation from an underlying smooth distribution. Motivated by the analysis of Diamond et al. (2020), in which the shipping-affected area spans 5.0° in the east–west direction, we define a distance of 250 km on both sides of the corridor center as a reasonable estimate of the affected area. To quantify the corridor effect, we apply a cubic fit to the 500 km affected range, based on data from the ranges 250–400 km away from the corridor on both sides. These ranges, and the fitted values within the affected area, represent a scenario without the shipping corridor, and they are subtracted from the actual values in order to estimate the corridor effect. When reporting an average corridor effect based on these differences, we focus on the corridor core, defined as a 150 km wide area (±75 km) around its center.

To assess the statistical significance of the corridor effect on a cloud variable, we examine the distribution of the corridor effect values within the corridor core region, defined above. If the 95 % confidence level (2σSD range, assuming a normal distribution) does not contain zero, we conclude that the corridor effect is statistically significant at the 0.05 level.

2.4 Propagation of uncertainties in monthly averages

All CLAAS-3 variables come with a pixel-based estimated uncertainty at instantaneous (level 2) retrievals. This uncertainty is propagated to level 3 daily and monthly averages. The implementation of the error propagation follows Stengel et al. (2017). For the daily averages, two cases are provided in CLAAS-3, with uncertainty correlations c equal to 0.1 and 1.0. For the monthly averages, which are computed based on the daily averages, the same two scenarios are available in CLAAS-3, using the daily mean uncertainty correlation c=0.1. For consistency, here we use the scenario where c=0.1 also for the monthly mean uncertainty. When averaging monthly mean data, we calculate the uncertainty of the averaged data (σx) based on

(1) σ x 2 = 1 N σ SD 2 + c σ i 2 ,

where N is the number of grid cells, σSD the standard deviation of the monthly mean values and σi the average uncertainty of these values (based on the c=0.1 scenario for monthly means).

For the quantification of the uncertainty of the no-ship scenarios, we repeat the cubic fit approach described above five times, varying the distance of the data ranges used in the fit from 150–300 to 350–500 km. The uncertainty of the no-ship scenario is then given by the standard deviation of these five runs. As mentioned before, the range 250–400 km is used for the quantification of the corridor effect. The corresponding uncertainty of the corridor effect is finally calculated as the combined uncertainty of the actual profile and the no-ship scenario.

3 Results

3.1 Average corridor effects

Figure 2a and b show maps of Nd and liquid re averages over the study region, calculated from the entire CLAAS-3 monthly time series (2004–2023). The shipping corridor appears as a straight line in a NW–SE direction in both Nd and re. Consistently with the Twomey effect, Nd in the corridor is higher than in the surroundings, and re is lower. This is shown even more clearly in the across-corridor average distributions of Nd (Fig. 2c) and re (Fig. 2d) and the ensuing differences from the no-ship scenarios (Fig. 2e and f, respectively). Estimation of the average corridor effect, as described in the previous section, yields an increase in Nd by 6.1 ± 1.5 cm−3 on average over the affected core region and a decrease in re by 0.21 ± 0.04 µm. The 2σSD values of the distribution of the corridor effect in the corridor core are 2.3 cm−3 and 0.08 µm, for Nd and re, respectively, showing that the effect on both is statistically significant at the 95 % confidence level. The average effect values agree well with those given by Diamond et al. (2020) for approximately the same region, even though they use a different approach and a wider region to determine their “NoShip” scenario. For Nd, they found a corridor effect of  5 cm−3 (based only on Aqua MODIS), while the effect found for re was 0.29 µm for Aqua MODIS and 0.28 µm for Terra MODIS.

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f02

Figure 2Maps of monthly time series averages of (a) Nd and (b) liquid re over the study region (with corresponding uncertainties given in Fig. S1 in the Supplement). Plots (c) and (d) show the corresponding across-corridor average distributions, with the propagated uncertainties (1σ) in light-blue shades and the no-ship scenarios in dotted curves. The standard deviation of the no-ship scenarios (1σSD) are shown in grey shades. Plots (e) and (f) show the corridor effect, calculated as the differences between the actual distribution and the no-ship scenarios, and their combined uncertainties. The average effects, calculated from the corridor core in absolute and relative terms, are also displayed. Dotted vertical lines denote the corridor center. In plots (e) and (f) the zero lines are also shown.

Download

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f03

Figure 3As in Fig. 2 but for the liquid water path (W) and the daytime cloud fraction (fc,day).

Download

Corresponding results for W and fc,day, the two cloud variables associated with possible adjustment mechanisms, are shown in Fig. 3. In the W case, the corridor perturbation is still apparent (Fig. 3c and e), with an effect of 1.6 ± 0.5 g m−2 in the core region. The 2σSD variation of the effects in the core is 0.5 g m−2, implying that the corridor effect on W is statistically significant. The mean effect on W (1.6 g m−2) is closer to the Aqua MODIS-based estimate of 1.3 g m−2 from Diamond et al. (2020) than to the estimate from Terra MODIS (0.5 g m−2). The response of W to Nd changes is complicated and depends strongly on the local meteorological and cloud conditions (see, e.g., Gryspeerdt et al., 2019, and references therein). The negative response found here is consistent with studies of similar cloud regimes (Manshausen et al., 2022; Gryspeerdt et al., 2019). The higher Nd values resulting from the ship emissions lead to an increase in the entrainment rate of dry air at the cloud top, thus decreasing W (Bretherton et al., 2007). Globally, however, results are not conclusive. A positive W adjustment is reported in Manshausen et al. (2022), while Wall et al. (2022, 2023) find decreases in W over global oceans when sulfate concentrations increase. However, these two studies are not limited to ship tracks, so their findings may not be directly comparable to the previous one. More recently, Tippett et al. (2024) showed that the W response to aerosol perturbations from ship emissions is weak on average, after correcting for biases in prior research based on tracking ship-affected air masses (Manshausen et al., 2022, 2023), related to correlations between wind and cloud properties. Several studies analyzing visible ship tracks also indicate weak W response to emitted aerosol (e.g., Toll et al., 2019).

In the case of fc,day there is a weak feature that appears centered along the corridor, with slightly higher values on the sides (between 100 and 200 km) compared to the center (Fig. 3d). Estimation of the corridor effect shows a weak increase of 0.03 ± 0.02 %, which is statistically insignificant at the 95 % confidence level (2σSD=0.08). Diamond et al. (2020) also report weak, non-significant increases in daytime cloud fraction from MODIS on both Terra (0.05 %) and Aqua (0.15 %) over roughly the same region. However, based on the shape of the across-corridor fc,day distribution (Fig. 3d) and the corresponding effect (Fig. 3f), it is possible that the corridor emissions induce a decrease in fc,day in a narrower range; in fact, the estimated uncertainty of the no-ship scenario (grey band in Fig. 3d) includes both possibilities on the sign of the corridor effect: an increase or a decrease in fc,day. This is also shown in Fig. S2d, where the no-ship scenarios based on data from different distance ranges from the corridor center are shown: results for Nd, re and W show that this method is robust for calculating the corridor effect on these variables. In the fc,day case, however, if we assume that unaffected areas start closer to the corridor center, and we apply the no-ship scenario fit accordingly, the corridor effect switches from positive to negative and leads to a decrease in fc,day by 0.22 ± 0.23 % on average. While the reason that the corridor would affect fc,day in a narrower range compared to other cloud properties is not obvious, a decrease in fc,day would be more in line with the decrease in W. Apart from the expected error in the method, another possible explanation of Fig. 3d and f is that the effect is more complicated, with a decrease closer to the corridor center and an increase on its sides. Investigation of such a scenario, however, requires modeling underlying processes beyond the scope of this study.

The same analysis was also conducted for τ, the cloud optical thickness. Based on the relationship of τ with Nd and W (see, e.g., Eq. (2) in Bennartz and Rausch, 2017), and the changes over the corridor of the two latter, a small decrease should be expected in τ. However, there is no apparent corridor effect in the across-corridor τ distributions, unlike the cloud properties examined before (see Fig. S3). This result may be due to τ variability across the corridor dominating the signal and highlights a limitation in our methodology: if the corridor does not manifest as a deviation from an otherwise smooth background, no conclusive remarks on its effect can be made. Two opposing mechanisms, namely an increase in τ due to the Twomey effect and a decrease due to the decreasing W, may also be acting here. Finally, this across-corridor τ distribution, which is a time series average, does not exclude the possibility of discernible corridor effects appearing in the seasonal and diurnal results. Thus, τ is still examined in the analyses that follow.

Apart from differences in the methodology for the shipping corridor identification and the estimation of its effects, the different spatial resolutions of MODIS and SEVIRI are also a factor that may play an important role when comparing our results with those from Diamond et al. (2020), who focused on almost the same region. A coarser resolution will lead to decreased τ and increased re values; this is the result of an increased probability of pixels containing a mixture of cloudy and clear scenes, which will lead to lower reflectances in both the visible and the shortwave infrared channel used in the retrievals. Thus, the coarser resolution of SEVIRI ( 3 km) can lead to biases with corresponding MODIS results ( 1 km resolution). However, it is beyond the setup of this study to disentangle such biases from differences originating in the methodology of the two studies and in the two sensors themselves. Additionally, the different resolution is not expected to affect the corridor effect significantly: the effect is estimated based on retrievals from the corridor and surrounding areas, which will be similarly affected by the coarser resolution of SEVIRI compared to MODIS.

3.2 Seasonal cycles

Figure 4 shows the seasonal variability in Nd, liquid re (Fig. 4a), W and fc,day (Fig. 4b) spatial averages over the study domain and the corresponding corridor effects. Average monthly values of the effect, along with their statistical significance at the 95 % confidence level, are shown in Table 1. The corridor effect on Nd, re (Fig. 4c) and W (Fig. 4d) is consistent in sign and detected throughout the year, albeit with varying strength (see Fig. S4b, d and f for respective monthly profiles).

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f04

Figure 4Seasonal variability of the spatially averaged (area-weighted over study domain as in Fig. 1) Nd and liquid re (a), W and fc,day (b), and corresponding monthly corridor effects (c, d). The corridor effect values are calculated by averaging differences between the actual and the no-ship scenario values in a 150 km wide area centered on the corridor (see also Sect. 2.3). Light-colored bands show the uncertainties (1σ) for each parameter. Dotted lines in (c) and (d) show the zero lines for each variable.

Download

Table 1Monthly averages (±2σSD) of the corridor effect on Nd, re, W and fc,day. Statistically significant effects at the 95 % confidence level are shown in bold.

Download Print Version | Download XLSX

The highest Nd values appear in May, June and July and secondarily in March and April (Fig. 4a). Lower values occur in local spring and summer, with the lowest occurring in September. Liquid re is anti-correlated with Nd, with the highest values occurring in August and September and the lowest in April–June. The seasonal cycle of Nd is very similar to corresponding results from previous studies based on MODIS (Bennartz and Rausch, 2017; Grosvenor et al., 2018), with CLAAS-3 exhibiting higher values. The peak in July and the rapid drop in August and September are prominent in all data sets. A similar Nd seasonal cycle peaking in summer is also found in Li et al. (2018) over a larger region and using both MODIS and CALIPSO data. Zeng et al. (2014) also find Nd peaks west of Namibia in April and July, based on MODIS and CALIOP data. Previous studies report that the seasonal cycle of Nd in this region is affected mainly by aerosol emissions from southern Africa (Grosvenor et al., 2018; Li et al., 2018; Bennartz and Rausch, 2017). The presence of absorbing aerosols from biomass burning activities above the cloud layer, which is typical in this region for certain months, could also have an artificial effect on the retrieved τ and re: by darkening the scene, especially in the 0.6 µm channel, smoke aerosols over clouds can lead to the retrieval of lower τ as well as lower re. However, when using the 3.9 µm channel reflectance, the effect on re is – in comparison to shorter wavelengths – relatively small (Haywood et al., 2004). Since the Nd retrieval depends weakly on τ and much more strongly on re (see, e.g., Eq. (2) in Bennartz and Rausch, 2017), the absorbing aerosol effect on Nd is also expected to be modest compared to retrievals based on smaller wavelengths (e.g., 1.6 µm). Furthermore, smoke aerosols in the region are expected to have a similar effect in both corridor-affected and unaffected parts, thus influencing the corridor effect calculations even less than area-wide average values.

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f05

Figure 5Diurnal variability of the spatially averaged (area-weighted) Nd (a) and liquid re (c) and corresponding diurnal corridor effects (b, d). Light-blue bands show the standard deviation for each parameter.

Download

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f06

Figure 6As in Fig. 5 but for W and fc. The two vertical dashed lines in the diurnal variation of fc (c) and corresponding corridor effect (d) represent the start and end time of the day when Nd, re and W retrievals are also available.

Download

The strong seasonal cycle of the (in-cloud) W follows that of fc,day (Fig. 4b), which peaks between August and November. In turn, the seasonal peak in fc,day coincides with the season of maximum static stability, which follows the seasonal cycle of sea surface temperature (SST; Klein and Hartmann, 1993) inversely. Thus, thicker clouds in September–November can explain the peak in W but also the higher re values, which are expected to grow with height. This season is also reported by Diamond et al. (2020) as giving the strongest corridor signal. In our case, the corridor effect on Nd is consistently positive and statistically significant throughout the year (meaning more cloud droplets in observations compared to a scenario without ships) and appears stronger from July to December (Fig. 4c and Table 1). Similarly, for re, there was a consistently significant negative corridor effect, meaning smaller cloud particles in reality compared to a no-ship scenario. The presence of the corridor also causes a decrease in W throughout the year (Fig. 4d and Table 1). In the case of fc,day, the corridor signal is ambiguous and appears to change sign during the year (Fig. 4d and Table 1). However, a closer examination of monthly across-corridor profiles of effects on fc,day (Fig. S4h) shows that only in austral summer months (mainly October, December, February and March) do the profiles appear as deviations roughly symmetrical across the corridor center, indicating a distinguishable corridor effect. In these cases, actual values are lower than in the no-ship scenario. Thus, while the time series average corridor effect on fc,day is ambiguous (Fig. 3f), a clearer signal of fc,day decreasing over the corridor appears in austral summer. In other months, it is difficult to draw any conclusion on the shipping corridor effect on fc,day. Therefore, a statistically significant effect in these months (i.e., May–September values in Table 1), which is not accompanied by an across-corridor profile of changes as described before, should not be interpreted as a meaningful corridor effect.

Results of across-corridor τ profiles per month are shown in Fig. S4i and corresponding profiles of corridor effect in Fig. S4j. The latter, however, are only provided for completeness: the corridor does not manifest convincingly in any month examined, so the calculated effects in Fig. S4j cannot be used to draw any conclusion.

3.3 Diurnal cycles

Figure 5 shows the diurnal variability of the Nd and re averages over the study area (Fig. 5a, c) and the corresponding average corridor effects (Fig. 5b, d). Corresponding across-corridor profiles of effects for every time slot when data are available are shown in the Fig. S5. Since Nd and re can only be retrieved during the day, it was required that for each time slot, the study region is completely covered with data from all months in the time series. In this way, time slots at the beginning and end of the day, when high solar zenith angles may lead to problematic retrievals, were omitted. For the same reason, the two extremes of the remaining available time slots (08:00 and 15:00 UTC) were also removed from the analysis. A statistical significance analysis of the diurnal cycles was not performed due to the coarser resolution of these data, which leads to very few points inside the corridor core.

Nd is nearly constant during the available time slots of the day (Fig. 5a). The effect of the corridor on Nd also appears stable during the day, ranging between 5 and 6 cm−3, with somewhat higher values in the afternoon (Fig. 5b). Liquid re decreases slightly during the day (Fig. 5c), in consistency with previous findings (Seethala et al., 2018). The corridor effect on re also gets weaker as the day progresses, stabilizing early in the afternoon (Fig. 5d).

Diamond et al. (2020) report larger corridor effects on Nd during the Terra overpass ( 10:30 local time), compared to Aqua ( 13:30 local time), contrary to this study. Their estimation of the relative change in Terra Nd is based on respective changes in W and re because of the absence of a published Terra Nd product. For re they report similar effects during the Terra and Aqua overpasses, as also found here (Fig. 5d).

Corresponding results for W show a strong decrease during the day, accompanied by a weakening of the corridor effect (Fig. 6a and b). In the afternoon, the corridor effect on W is about half compared to the morning value. The decrease in W during the day is in good agreement with results of Seethala et al. (2018): they also report peak values of W early in the morning and a decrease until the afternoon. A small decrease is also reported for re during the day. The decrease in W is associated with the cloud thinning, which is apparent in the fc diurnal plot (Fig. 6c).

The estimated corridor effect on W by Diamond et al. (2020) is lower in the morning (Terra overpass) than in the afternoon (Aqua overpass). This contrasts with our findings, which show similar decreases in the two relevant time slots (10:30 and 13:30 UTC in Fig. 6b), with a somewhat weaker effect during the Aqua overpass. Other studies of ship tracks in the northeastern Pacific, where the cloud conditions are similar, show even larger W decreases for Terra MODIS than for Aqua MODIS, with the former being almost double the latter (Christensen et al., 2009; Segrin et al., 2007). These results hinder concrete conclusions regarding the differences in the corridor effect between morning and afternoon from being drawn and highlight the value of using observations in more frequent time intervals, such as CLAAS-3, for similar analyses. However, diurnal observations should also be treated with caution, since they may suffer from diurnally varying biases related to scattering geometry (Smalley and Lebsock, 2023).

Table 2Time series averages (±2σSD) of the corridor effect on Nd, re, W and fc,day during the periods before 2020 and from 2020 onwards and their difference. Statistically significant effects at the 95 % confidence level are shown in bold.

Download Print Version | Download XLSX

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f07

Figure 7Corridor effects on Nd (a), re (b), W (c) and fc,day (d) in the periods 2004–2019 and 2020–2023, calculated as the differences between the actual average distribution and the no-ship scenarios for these periods. The light-colored bands show the associated propagated uncertainties (1σ). Dotted vertical lines denote the corridor center. The zero lines are also shown.

Download

Since fc can also be retrieved at nighttime, the full 24 h cycle variation is available, showing the cloud thinning during day and the thickening during the night (Fig. 6c). It is interesting to note that the corridor effect on fc changes sign between day and night (Fig. 6d). It is negative during the night and becomes positive at day, contrary to what would be expected from the diurnal cycles of precipitation and evaporation, which maximize at night and day, respectively (e.g., Sandu et al., 2008). A closer examination of the fc difference profiles per time slot, however, shows that the positive day time differences have a local minimum at the corridor center (Fig. S5d). This is reminiscent of the pattern found also in the seasonal analysis of fc,day, suggesting that selection of a narrower range for the corridor-affected area during the day would lead to lower (or even negative) effects on fc. However, since this indication of a narrower affected area does not appear in other variables, no concrete conclusion can be drawn regarding the effect of ship emissions on fc during the day. It can be safely concluded, however, that during the night the shipping corridor exerts a negative effect on fc, as can be shown by examining shipping corridor and no-ship scenario profiles of fc on an individual time slot basis (Fig. S6).

Results on the diurnal variation of τ over the region, across-corridor profiles and corresponding estimates of corridor effect are shown in Fig. S7. Average values of τ decrease during the day, as expected (Fig. S7a). The diurnal variation of the corridor effect shows practically no effect in the morning, with slightly positive values in the afternoon (Fig. S7b). However, these results should be interpreted with caution, since corresponding across-corridor profiles have high uncertainty and show little or no appearance of corridor-centered perturbations (Fig. S7c and d–i).

3.4 Long term changes

In order to examine the strength of the corridor effect before and after 2020, when the stricter IMO regulations on sulfur-containing emissions were implemented, we repeated the analysis described in Sect. 3.1 for two time periods: 2004–2019 and 2020–2023. While the two time periods are not equal, the longer period before 2020 is more representative of the entire CLAAS-3 time range and reduces the risk of being affected by large year-to-year variability. The across-corridor distributions of the effect on Nd, re, W and fc,day are shown in Fig. 7. The effect of the shipping corridor on Nd, re and W has weakened notably from 2020 onward, acquiring almost half the values of the 2004–2019 period. On average, before 2020 the ship emissions caused an increase in Nd by 7.00 ± 1.49 cm−3, which weakened to 2.73 ± 1.50 cm−3 from 2020 onward. Similarly for re, the average decrease by 0.25 ± 0.04 µm due to ship emissions before 2020 weakened to 0.09 ± 0.04 µm from 2020 onward. Corresponding average effects on W are 1.80 ± 0.47 g m−2 and 0.99 ± 0.52 g m−2. In the fc,day case, results are less clear. The across-corridor distribution before 2020 shows some resemblance to the full period distribution (Fig. 3f), although with weaker characteristics. From 2020 onward, no corridor effect is apparent based on the shape of the distribution.

The corridor effects on Nd, re and W appear statistically significant at the 95 % confidence level in both time periods. Similar results were obtained when the differences between the corridor effect profiles before and after 2020 were examined: significant differences were found for Nd, re and W but not for fc,day. These results are summarized in Table 2.

In the case of τ, since the cubic fit approach does not produce useful results in either of the two periods examined, we analyzed the across-corridor profiles of absolute values. As shown in Fig. S8a and b, a significant decrease in τ is apparent in the period after 2020. This is analyzed further below.

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f08

Figure 8Maps of Nd (a), liquid re (b), W (c) and fc,day (d) differences between the periods after and before 2020. Corresponding uncertainties are shown in Fig. S9.

Download

To further examine how the recent changes in ship emissions affected the cloud properties along the corridor compared to less affected areas, we calculated, for each cloud property and grid cell in the study area, the difference between the averages from the periods from 2020 onward and before 2020. Results are shown in Fig. 8.

The shipping corridor is clearly visible in the Nd and re maps. In the case of Nd it manifests as stronger decreases, compared to an overall background decrease which is stronger closer to the coast. In the case of re, it appears as larger increases in the southern part of the corridor and smaller increases in the northern part, compared to background decreases. In the maps of W and fc,day differences, the corridor is not apparent. Instead, for both variables large decreases appear over the entire area. The same is true for τ (Fig. S8c). As discussed in Sect. 3.2, SST is the main driver for the seasonal variability of fc,day, and this is also the case for a longer term. Thus, a possible explanation for this apparent absence of a corridor effect change during the examined period is that changes due to ship emissions are hidden by a wider and stronger effect of SST changes on fc,day. This scenario is examined in more detail below.

The statistical significance of the differences shown in Fig. 8 was also assessed by comparing the average differences with their combined uncertainties. In the cases of Nd and re, differences are only significant in the southeastern part of the region (Fig. S9a and b), while they cover larger parts in the W, fc,day (Fig. S9c, d), and τ (Fig. S8d) cases. In no case is the shipping corridor highlighted from adjacent areas in terms of significance of differences. Similar findings over the region are also reported in Diamond (2023). The corridor, however, is clearly visible in the Nd and re differences maps in Fig. 8, which is a strong indication that these differences originate in ship emission changes.

https://acp.copernicus.org/articles/25/6957/2025/acp-25-6957-2025-f09

Figure 9Time series of monthly average Nd (a), liquid re (c), W (e) and fc,day (g) over the shipping corridor in 2004–2023, deseasonalized using a 12-month running average. Corresponding monthly corridor effects are plotted in (b), (d), (f) and (h). The light-blue bands show the propagated uncertainty (1σ) in each case. In all plots, dotted vertical lines denote the beginning of each year. In the corridor effect plots, the zero lines are also shown.

Download

As a further step in the assessment of the long-term changes in corridor averages and effects, the time series of monthly average Nd, re, W and fc,day over the corridor and the corresponding corridor effects were also analyzed. Figure 9 shows the resulting time series of averages (left column) and evolution of the corridor effects (right column). All time series were deseasonalized and smoothed using a 12-month running average.

Results show a considerable variability in Nd and re during the study period (Fig. 9a and c). In more recent years there has been a tendency of Nd to decrease and re to increase, which may be due to natural variability, anthropogenic causes or a combination thereof. The corridor effect on these two variables has notably weakened in recent years, with values that are not found before 2020 (Fig. 9b and d). Features before 2020 may also be linked to fluctuations in shipping activity. Examining a region in the NE Pacific, Yuan et al. (2022) report dips in the ship track density in 2009–2010 (after-effect of the 2008 financial crisis) and 2014–2016, likely caused by a strong slowdown in the Chinese economy. Abrupt decreases in the corridor effect on Nd in 2010 and 2014 (Fig. 9b) may be associated with these circumstances.

It is also worth noting that fc,day decreases considerably from 2020 onward compared to the period before, with values returning closer to long-term averages in 2023 (Fig. 9g). Similar features also appear in the W time series (Fig. 9e). The corresponding corridor effects on fc,day and W, however, are less notably different from 2020 onward compared to the years before (Fig. 9h and f). As in the overall changes over the region examined in Fig. 8, the time series of τ shows very similar characteristics to the series of fc,day and W (Fig. S8b).

As mentioned before, SST is one of the main drivers of cloud variability in marine low-cloud regimes (Andersen et al., 2023), and this can largely explain these findings. Diamond (2023) also reports low cloud albedo values in 2020–2022 (which are consistent with the lower fc,day values found here) and suggests that these may be related to unusually warm SST during these years. Interestingly, Gettelman et al. (2024) suggest that a large part of the unusually high SSTs in 2023, especially in the Northern Hemisphere, can be attributed to the IMO regulations and their effect on clouds, claiming that in 2021–2023 cloud anomalies were more likely to drive SST anomalies than the other way around. While an analysis of these causation mechanisms is beyond the scope of this study, an examination of SST data from Aqua MODIS (Fig. S10) yields a correlation coefficient R=-0.67 between SST and fc,day, confirming their correlation over the study region.

4 Summary and conclusions

In this study we analyzed the effect of ship emissions on cloud properties over a busy shipping corridor that crosses the SE Atlantic. The analysis covered the 20-year period 2004–2023, when data from CLAAS-3, the cloud data set based on the geostationary SEVIRI imager, are available. Taking advantage of the CLAAS-3 temporal resolution, the corridor effect on the stratocumulus clouds of the region was quantified on diurnal and seasonal bases, while long-term changes were also examined.

Results show a clear effect of shipping emissions on the cloud properties associated with the Twomey effect, i.e., an increase in Nd and a decrease in re. Subsequent adjustments reveal a decrease in W, while the cloud fraction changes are more subtle and limited. No clear impact on τ was found, suggesting an overall minor radiative effect of the shipping emissions, although limitations in our method to detect changes in the corridor cannot be excluded. The effects vary seasonally and diurnally, depending on corresponding regional conditions in the former case and on the cloud thinning during the day in the latter.

In the long term, a weakening of the corridor effect is clearly seen in Nd and re from 2020 onward, presumably due to the stricter IMO regulations on sulfur emissions that were implemented at the beginning of that year. Changes in W and the cloud fraction are also detected over the wider region, associated mainly with corresponding SST variations.

Taking advantage of the good alignment between the SE Atlantic shipping corridor and prevailing winds, the approach used here for the quantification of ship emission effects on clouds cannot be directly implemented in other regions and shipping corridors, without prior adjustments in the methodology. The analysis, however, is valuable considering the climatic importance of the extensive SE Atlantic stratocumulus region. It also highlights the great potential of using geostationary-based cloud observations in similar studies, which has not been exploited so far.

Code availability

All code for the data analysis associated with this study is available at https://doi.org/10.5281/zenodo.14724159 (Benas, 2025).

Data availability

CLAAS-3 data were obtained from the CM SAF Web User Interface (https://doi.org/10.5676/EUM_SAF_CM/CLAAS/V003, Meirink et al., 2022). Global shipping traffic density data were obtained from the World Bank Data Catalog (https://datacatalog.worldbank.org/search/dataset/0037580/Global-Shipping-Traffic-Density, Cerdeiro et al., 2020).

Supplement

The supplement related to this article is available online at https://doi.org/10.5194/acp-25-6957-2025-supplement.

Author contributions

NB and JFM designed the study and performed the analysis. NB wrote the code for the analysis and the first draft of the manuscript. RR and MS provided input on the structure and contents of the manuscript. All authors reviewed and edited the manuscript.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

Acknowledgements

This work was performed within the EUMETSAT CM SAF framework, and all authors acknowledge the financial support of the EUMETSAT member states.

Financial support

This research has been supported by the European Organization for the Exploitation of Meteorological Satellites (funded by the EUMETSAT ground segment budget (Satellite Application Facilities)).

Review statement

This paper was edited by Graham Feingold and reviewed by Michael Diamond and two anonymous referees.

References

Albrecht, B. A.: Aerosols, Cloud Microphysics, and Fractional Cloudiness, Science, 245, 1227–1230, https://doi.org/10.1126/science.245.4923.1227, 1989. 

Andersen, H., Cermak, J., Douglas, A., Myers, T. A., Nowack, P., Stier, P., Wall, C. J., and Wilson Kemsley, S.: Sensitivities of cloud radiative effects to large-scale meteorology and aerosols from global observations, Atmos. Chem. Phys., 23, 10775–10794, https://doi.org/10.5194/acp-23-10775-2023, 2023. 

Bellouin, N., Quaas, J., Gryspeerdt, E., Kinne, S., Stier, P., Watson-Parris, D., Boucher, O., Carslaw, K. S., Christensen, M., Daniau, A.-L., Dufresne, J.-L., Feingold, G., Fiedler, S., Forster, P., Gettelman, A., Haywood, J. M., Lohmann, U., Malavelle, F., Mauritsen, T., McCoy, D. T., Myhre, G., Mülmenstädt, J., Neubauer, D., Possner, A., Rugenstein, M., Sato, Y., Schulz, M., Schwartz, S. E., Sourdeval, O., Storelvmo, T., Toll, V., Winker, D., and Stevens, B.: Bounding Global Aerosol Radiative Forcing of Climate Change, Rev. Geophys., 58, e2019RG000660, https://doi.org/10.1029/2019RG000660, 2020. 

Benas, N.: nikosbenas/Clouds_over_SEA_shipping_corridor: ACP paper code, Zenodo [code], https://doi.org/10.5281/zenodo.15222065, 2025. 

Benas, N., Solodovnik, I., Stengel, M., Hüser, I., Karlsson, K.-G., Håkansson, N., Johansson, E., Eliasson, S., Schröder, M., Hollmann, R., and Meirink, J. F.: CLAAS-3: the third edition of the CM SAF cloud data record based on SEVIRI observations, Earth Syst. Sci. Data, 15, 5153–5170, https://doi.org/10.5194/essd-15-5153-2023, 2023. 

Bennartz, R. and Rausch, J.: Global and regional estimates of warm cloud droplet number concentration based on 13 years of AQUA-MODIS observations, Atmos. Chem. Phys., 17, 9815–9836, https://doi.org/10.5194/acp-17-9815-2017, 2017. 

Bretherton, C. S., Blossey, P. N., and Uchida, J.: Cloud droplet sedimentation, entrainment ef?ciency, and subtropical stratocumulus albedo, Geophys. Res. Lett., 34, L03813, https://doi.org/10.1029/2006GL027648, 2007. 

Cerdeiro, D. A., Komaromi, A., Liu, Y., and Saeed, M.: IMF’s World Seaborne Trade monitoring system, World Bank Data Catalog [data set], https://datacatalog.worldbank.org/search/dataset/0037580/Global-Shipping-Traffic-Density (last access: 7 July 2025), 2020 

Christensen, M. W., Coakley, J. A., and Tahnk, W. R.: Morning-to-Afternoon Evolution of Marine Stratus Polluted by Underlying Ships: Implications for the Relative Lifetimes of Polluted and Unpolluted Clouds, J. Atmos. Sci., 66, 2097–2106, https://doi.org/10.1175/2009JAS2951.1, 2009. 

Christensen, M. W., Jones, W. K., and Stier, P.: Aerosols Enhance Cloud Lifetime and Brightness along the Stratus-to-Cumulus Transition, P. Natl. Acad. Sci. USA, 117, 17591–17598, https://doi.org/10.1073/pnas.1921231117, 2020. 

Christensen, M. W., Gettelman, A., Cermak, J., Dagan, G., Diamond, M., Douglas, A., Feingold, G., Glassmeier, F., Goren, T., Grosvenor, D. P., Gryspeerdt, E., Kahn, R., Li, Z., Ma, P.-L., Malavelle, F., McCoy, I. L., McCoy, D. T., McFarquhar, G., Mülmenstädt, J., Pal, S., Possner, A., Povey, A., Quaas, J., Rosenfeld, D., Schmidt, A., Schrödner, R., Sorooshian, A., Stier, P., Toll, V., Watson-Parris, D., Wood, R., Yang, M., and Yuan, T.: Opportunistic experiments to constrain aerosol effective radiative forcing, Atmos. Chem. Phys., 22, 641–674, https://doi.org/10.5194/acp-22-641-2022, 2022. 

Diamond, M. S.: Detection of large-scale cloud microphysical changes within a major shipping corridor after implementation of the International Maritime Organization 2020 fuel sulfur regulations, Atmos. Chem. Phys., 23, 8259–8269, https://doi.org/10.5194/acp-23-8259-2023, 2023. 

Diamond, M. S., Director, H. M., Eastman, R., Possner, A., and Wood, R.: Substantial Cloud Brightening from Shipping in Subtropical Low Clouds, AGU Advances, 1, e2019AV000111, https://doi.org/10.1029/2019av000111, 2020. 

Feingold, G., Ghate, V. P., Russell, L. M., Blossey, P., Cantrell, W., Christensen, M. W., Diamond, M. S., Gettelman, A., Glassmeier, F., Gryspeerdt, E., Haywood, J., Hoffmann, F., Kaul, C. M., Lebsock, M., McComiskey, A. C., McCoy, D. T., Ming, Y., Mülmenstädt, J., Possner, A., Prabhakaran, P., Quinn, P. K., Schmidt, K. S., Shaw, R. A., Singer, C. E., Sorooshian, A., Toll, V., Wan, J. S., Wood, R., Yang, F., Zhang, J., and Zheng, X.: Physical science research needed to evaluate the viability and risks of marine cloud brightening, Sci. Adv., 10, eadi8594, https://doi.org/10.1126/sciadv.adi8594, 2024. 

Gettelman, A., Christensen, M. W., Diamond, M. S., Gryspeerdt, E., Manshausen, P., Stier, P., Watson-Parris, D., Yang, M., Yoshioka, M., and Yuan, T.: Has reducing ship emissions brought forward global warming? Geophys. Res. Let., 51, e2024GL109077, https://doi.org/10.1029/2024GL109077, 2024. 

Goren, T. and Rosenfeld, D.: Satellite Observations of Ship Emission Induced Transitions from Broken to Closed Cell Marine Stratocumulus over Large Areas, J. Geophys. Res.-Atmos., 117, D17206, https://doi.org/10.1029/2012jd017981, 2012. 

Goren, T. and Rosenfeld, D.: Extensive Closed Cell Marine Stratocumulus Downwind of Europe – A Large Aerosol Cloud Mediated Radiative Effect or Forcing?, J. Geophys. Res.-Atmos., 120, 6098–6116, https://doi.org/10.1002/2015JD023176, 2015. 

Grosvenor, D. P., Sourdeval, O., Zuidema, P., Ackerman, A., Alexandrov, M. D., Bennartz, R., Boers, R., Cairns, B., Chiu, C., Christensen, M., Deneke, H., Diamond, M., Feingold, G., Fridlind, A., Hünerbein, A., Knist, C., Kollias, P., Marshak, A., McCoy, D., Merk, D., Painemal, D., Rausch, J., Rosenfeld, D., Russchenberg, H., Seifert, P., Sinclair, K., Stier, P., B. van D., Wendisch, M., Werner, F., Wood, R., Zhang, Z., and Quaas, J.: Remote sensing of cloud droplet number concentration in warm clouds: A review of the current state of knowledge and perspectives, Rev. Geophys., 56, 409–453, https://doi.org/10.1029/2017RG000593, 2018. 

Gryspeerdt, E., Goren, T., Sourdeval, O., Quaas, J., Mülmenstädt, J., Dipu, S., Unglaub, C., Gettelman, A., and Christensen, M.: Constraining the aerosol influence on cloud liquid water path, Atmos. Chem. Phys., 19, 5331–5347, https://doi.org/10.5194/acp-19-5331-2019, 2019. 

Gryspeerdt, E., Goren, T., and Smith, T. W. P.: Observing the timescales of aerosol–cloud interactions in snapshot satellite images, Atmos. Chem. Phys., 21, 6093–6109, https://doi.org/10.5194/acp-21-6093-2021, 2021. 

Haywood, J. M., Osborne, S. R., and Abel, S. J.: The effect of overlying absorbing aerosol layers on remote sensing retrievals of cloud effective radius and cloud optical depth, Q. J. Roy. Meteor. Soc., 130, 779–800, https://doi.org/10.1256/qj.03.100, 2004. 

Hu, S., Zhu, Y., Rosenfeld, D., Mao, F., Lu, X., Pan, Z., Zang, L., and Gong, W.: The dependence of ship-polluted marine cloud properties and radiative forcing on background drop concentrations, J. Geophys. Res.-Atmos., 126, e2020JD033852, https://doi.org/10.1029/2020JD033852, 2021. 

IMO: RESOLUTION MEPC.320(74) 2019 GUIDELINES FOR CONSISTENT IMPLEMENTATION OF THE 0.50 % SULPHUR LIMIT UNDER MARPOL ANNEX V, London, UK, https://wwwcdn.imo.org/localresources/en/OurWork/Environment/Documents/ResolutionMEPC.320(74).pdf (last access: 7 May 2024), 2019. 

Klein, S. A. and Hartmann, D. L.: The seasonal cycle of low stratiform clouds, J. Climate, 6, 1587–1606, 1993. 

Larson, K. M., Shand, L., Staid, A., Gray, S., Roesler, E. L., and Lyons, D.: An optical flow approach to tracking ship track behavior using GOES-R satellite imagery, IEEE J. Sel. Top. Appl. Earth Obs., 15, 6272–6282, https://doi.org/10.1109/JSTARS.2022.3193024, 2022 

Li, J., Jian, B., Huang, J., Hu, Y., Zhao, C., Kawamoto, K., Liao, S., and Wu, M.: Long-term variation of cloud droplet number concentrations from space-based Lidar, Remote Sens. Environ., 213, 144–161, doi.org/10.1016/j.rse.2018.05.011, 2018. 

Manshausen, P., Watson-Parris, D., Christensen, M. W., Jalkanen, J.-P., and Stier, P.: Invisible ship tracks show large cloud sensitivity to aerosol, Nature, 610, 101–106, https://doi.org/10.1038/s41586-022-05122-0, 2022. 

Manshausen, P., Watson-Parris, D., Christensen, M. W., Jalkanen, J.-P., and Stier, P.: Rapid saturation of cloud water adjustments to shipping emissions, Atmos. Chem. Phys., 23, 12545–12555, https://doi.org/10.5194/acp-23-12545-2023, 2023. 

Meirink, J. F., Roebeling, R. A., and Stammes, P.: Inter-calibration of polar imager solar channels using SEVIRI, Atmos. Meas. Tech., 6, 2495–2508, https://doi.org/10.5194/amt-6-2495-2013, 2013. 

Meirink, J. F., Karlsson, K.-G., Solodovnik, I., Hüser, I., Benas, N., Johansson, E., Håkansson, N., Stengel, M., Selbach, N., Schröder, M., and Hollmann, R.: CLAAS-3: CM SAF CLoud property dAtAset using SEVIRI – Edition 3, EUMETSAT [data set], https://doi.org/10.5676/EUM_SAF_CM/CLAAS/V003, 2022. 

Nakajima, T. and King, M. D.: Determination of the optical thick-ness and effective particle radius of clouds from reflected solar radiation measurements, part 1: Theory, J. Atmos. Sci., 47, 1878–1893, https://doi.org/10.1175/1520-0469(1990)047<1878:DOTOTA>2.0.CO;2, 1990. 

Peters, K., Quaas, J., and Graßl, H.: A search for largescale effects of ship emissions on clouds and radiation in satellite data, J. Geophys. Res.-Atmos., 116, D24205, https://doi.org/10.1029/2011jd016531, 2011. 

Peters, K., Quaas, J., Stier, P., and Graßl, H.: Processes Limiting the Emergence of Detectable Aerosol Indirect Effects on Tropical Warm Clouds in Global AerosolClimate Model and Satellite Data, Tellus B, 66, 24054, https://doi.org/10.3402/tellusb.v66.24054, 2014. 

Roebeling, R. A., Feijt, A. J., and Stammes, P.: Cloud property retrievals for climate monitoring: implications of differences between SEVIRI on METEOSAT-8 and AVHRR on NOAA-17, J. Geophys. Res., 111, D20210, https://doi.org/10.1029/2005JD006990, 2006. 

Sandu, I., Brenguier, J.-L., Geoffroy, O., Thouron, O., and Masson, V.: Aerosol Impacts on the Diurnal Cycle of Marine Stratocumulus, J. Atmos. Sci., 65, 2705–2718, https://doi.org/10.1175/2008jas2451.1, 2008. 

Schreier, M., Joxe, L., Eyring, V., Bovensmann, H., and Burrows, J. P.: Ship track characteristics derived from geostationary satellite observations on the west coast of southern Africa, Atmos. Res., 95, 32–39, https://doi.org/10.1016/j.atmosres.2009.08.005, 2010 

Seethala, C., Meirink, J. F., Horváth, Á., Bennartz, R., and Roebeling, R.: Evaluating the diurnal cycle of South Atlantic stratocumulus clouds as observed by MSG SEVIRI, Atmos. Chem. Phys., 18, 13283–13304, https://doi.org/10.5194/acp-18-13283-2018, 2018. 

Segrin, M. S., Coakley, J. A., and Tahnk, W. R.: MODIS Observations of Ship Tracks in Summertime Stratus off the West Coast of the United States, J. Atmos. Sci., 64, 4330–4345, https://doi.org/10.1175/2007JAS2308.1, 2007. 

Smalley, K. M. and Lebsock, M. D.: Corrections for Geostationary Cloud Liquid Water Path Using Microwave Imagery, J. Atmos. Ocean. Tech., 40, 1049–1061, https://doi.org/10.1175/JTECH-D-23-0030.1, 2023. 

Stephens, G. L.: Radiation profiles in extended water clouds: II. Parameterization schemes, J. Atmos. Sci., 35, 2123–2132, https://doi.org/10.1175/1520-0469(1978)035<2123:RPIEWC>2.0.CO;2, 1978. 

Stengel, M., Stapelberg, S., Sus, O., Schlundt, C., Poulsen, C., Thomas, G., Christensen, M., Carbajal Henken, C., Preusker, R., Fischer, J., Devasthale, A., Willén, U., Karlsson, K.-G., McGarragh, G. R., Proud, S., Povey, A. C., Grainger, R. G., Meirink, J. F., Feofilov, A., Bennartz, R., Bojanowski, J. S., and Hollmann, R.: Cloud property datasets retrieved from AVHRR, MODIS, AATSR and MERIS in the framework of the Cloud_cci project, Earth Syst. Sci. Data, 9, 881–904, https://doi.org/10.5194/essd-9-881-2017, 2017. 

Tippett, A., Gryspeerdt, E., Manshausen, P., Stier, P., and Smith, T. W. P.: Weak liquid water path response in ship tracks, Atmos. Chem. Phys., 24, 13269–13283, https://doi.org/10.5194/acp-24-13269-2024, 2024. 

Toll, V., Christensen, M., Quaas, J., and Bellouin, N.: Weak average liquid-cloud-water response to anthropogenic aerosols, Nature, 572, 51–55, https://doi.org/10.1038/s41586-019-1423-9, 2019. 

Twomey, S.: Pollution and the Planetary Albedo, Atmos. Environ., 8, 1251–1256, 1974. 

Wall, C. J., Norris, J. R., Possner, A., McCoy, D. T., McCoy, I. L., and Lutsko, N. J.: Assessing Effective Radiative Forcing from Aerosol–Cloud Interactions over the Global Ocean, P. Natl. Acad. Sci. USA, 119, 2210481119, https://doi.org/10.1073/pnas.2210481119, 2022. 

Wall, C. J., Storelvmo, T., and Possner, A.: Global observations of aerosol indirect effects from marine liquid clouds, Atmos. Chem. Phys., 23, 13125–13141, https://doi.org/10.5194/acp-23-13125-2023, 2023. 

Watson-Parris, D., Christensen, M. W., Laurenson, A., Clewley, D., Gryspeerdt, E., and Stier, P.: Shipping regulations lead to large reduction in cloud perturbations, P. Natl. Acad. Sci. USA, 119, e2206885119, https://doi.org/10.1073/pnas.2206885119, 2022. 

Wood, R.: Stratocumulus clouds, Mon. Weather Rev., 140, 2373–2423, https://doi.org/10.1175/MWR-D-11-00121.1, 2012. 

Yuan, T., Song, H., Wood, R., Wang, C., Oreopoulos, L., Platnick, S. E., von Hippel, S., Meyer, K., Light, S., and Wilcox, E.: Global reduction in ship-tracks from sulfur regulations for shipping fuel, Sci. Adv., 8, eabn7988, https://doi.org/10.1126/sciadv.abn7988, 2022. 

Yuan, T., Song, H., Wood, R., Oreopoulos, L., Platnick, S. E., Wang, C., Yu, H., Meyer, K., and Wilcox, E.: Observational evidence of strong forcing from aerosol effect on low cloud coverage, Sci. Adv., 9, eadh7716, https://doi.org/0.1126/sciadv.adh7716, 2023. 

Yuan, T., Song, H., Oreopoulos, L., Wood, R., Bian, H., Breen, K., Chin, M., Yu, H., Barahona, D., Meyer, K., and Platnick, S.: Abrupt reduction in shipping emission as an inadvertent geoengineering termination shock produces substantial radiative warming, Commun. Earth Environ., 5, 1–8, https://doi.org/10.1038/s43247-024-01442-3, 2024.  

Zeng, S., Riedi, J., Trepte, C. R., Winker, D. M., and Hu, Y.-X.: Study of global cloud droplet number concentration with A-Train satellites, Atmos. Chem. Phys., 14, 7125–7134, https://doi.org/10.5194/acp-14-7125-2014, 2014. 

Download
Short summary
This study examines how ship emissions affect clouds over a shipping corridor in the southeastern Atlantic. Using satellite data from 2004 to 2023, we find that ship emissions increase the number of cloud droplets while reducing their size and slightly decrease cloud water content. Effects on seasonal and daily patterns vary based on regional factors. The impact of emissions weakened after stricter regulations were implemented in 2020.
Share
Altmetrics
Final-revised paper
Preprint