<figure><img src="image_1.png" /></figure>

HAL
open science

# On the Trend in Below-Cloud Solar Irradiance in The Netherlands versus That in Aerosol Sulfate Concentration

Suzanne Crumeyrolle, Andrey Khlystov, Harry ten Brink

► To cite this version:

Suzanne Crumeyrolle, Andrey Khlystov, Harry ten Brink. On the Trend in Below-Cloud Solar Irradiance in The Netherlands versus That in Aerosol Sulfate Concentration. Atmosphere, 2023, Atmosphere, 13, 10.3390/atmos13122037. [DOI] [ HAL Id: hal-04459236 ]

Submitted on 15 Feb 2024

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# On the Trend in Below-Cloud Solar Irradiance in The Netherlands versus That in Aerosol Sulfate Concentration

Suzanne Crumeyrolle, Andrey Khlystov, Harry ten Brink

► To cite this version:

Suzanne Crumeyrolle, Andrey Khlystov, Harry ten Brink. On the Trend in Below-Cloud Solar Irradiance in The Netherlands versus That in Aerosol Sulfate Concentration. Atmosphere, 2023, Atmosphere, 13, 10.3390/atmos13122037. [https://lilloa.hal.science/hal-04459236v1](https://lilloa.hal.science/hal-04459236v1)

HAL Id: hal-04459236

Submitted on 15 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# On the Trend in Below-Cloud Solar Irradiance in The Netherlands versus That in Aerosol Sulfate Concentration

Suzanne Crumeyrolle <fnref n="1" />,*, Andrey Khlystov <fnref n="2" /> and Harry Ten Brink <fnref n="3" />

**Abstract:** The below-cloud irradiance in The Netherlands increased by over 10 Wm⁻<fnref n="2" /> in the past half century. It was hypothesized that this could be due to a decrease in the aerosol serving as cloud condensation nuclei, on which the cloud droplets form, in the following way. With unchanged macrophysics, clouds with a lower number of droplets are less reflective, in other words, they transmit more solar radiation. This hypothesis cannot be substantiated with data because of a generic absence of long-term information on cloud droplet number concentrations (CDNCs) worldwide. To assess the historic trend in CDNC, we used the Boucher and Lohmann (B&L) empirical relationship between CDNC and the mass concentration of the water-attracting hygroscopic aerosol components. The B&L parameterization was tested and validated with observations from the CHIEF cloud chamber, in which the formation of marine stratocumulus, the most frequent cloud type in The Netherlands, is simulated. This study will focus on summer periods because the irradiance governs the yearly average at this latitude. The summer trend of sulfate, the most dominant hygroscopic component of observed aerosol mass concentration, was analyzed with EBAS data from 1972 onwards. The average summer CDNCs were then calculated via the B&L parameterization, showing an upper limit of 380 cm⁻<fnref n="3" /> in the 1970s and decreasing to around 200 cm⁻<fnref n="3" /> in the 2010s. The associated increase in transmission for thin marine stratocumulus without overlying clouds would be, at most, 3.5 W m⁻<fnref n="2" />. Unobstructed stratocumuli occur only part of the time, and the change in irradiance based on the reduction in cloud droplet number is certainly small in comparison to the empirically derived trend of 10 W m⁻<fnref n="2" />.

**Keywords:** aerosol-cloud interaction; cloud chamber; stratocumulus

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Crumeyrolle, S.; Khlystov, A.; Ten Brink, H. On the Trend in Below-Cloud Solar Irradiance in The Netherlands versus That in Aerosol Sulfate Concentration. *Atmosphere* **2022**, *13*, 2037. https://doi.org/10.3390/atmos13122037

Academic Editors: Wenhao Xue, Jing Zhang, Jing Wei and Yulu Tian

Received: 28 September 2022

Accepted: 28 November 2022

Published: 4 December 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Aerosols have a significant impact on the Earth's radiation budget by scattering and absorbing solar radiation (direct effect), leading to negative (cooling) and positive (warming) radiative forcing. Moreover, interactions between aerosols and clouds (indirect effect) have a large influence on the role of aerosols on climate change. The atmospheric aerosols affect cloud formation through their role as cloud condensation nuclei (CCN). Sulfate aerosols are nowadays considered as the most dominant hygroscopic component of aerosols [1,2]. Sulfate is an important secondary aerosol formed from photochemical reactions of sulfur dioxide (SO₂) emissions in the atmosphere. Summertime peaks of sulfate concentrations are common for most areas, including The Netherlands [3], due to available solar insolation, chemical reactions facilitated in high relative humidity environments, and stagnation events e.g., [4,5]. To reduce SO₂ emissions over the world, regulatory and legislative mandates were established and total annual SO₂ emissions over Europe have decreased by 72% from 1980 to 2004 [6,7], and by 60% from 1990 to 2010 over the United States [8].

Over the past 50 years, the solar irradiance over The Netherlands significantly increased, both during clear-sky and cloudy situations [9]. The solar irradiance increase

<footnote n="1">¹ CNRS, UMR 851 LOA-Laboratoire d'Optique Atmosphérique, University Lille, 59000 Lille, France</footnote>

<footnote n="2">² Desert Research Institute DRI, Reno, NV 89512, USA</footnote>

<footnote n="3">³ Energy-Research Centre of The Netherlands ECN, 1755 ZG Petten, The Netherlands</footnote>

<footnote>* Correspondence: suzanne.crumeyrolle@univ-lille.fr</footnote>
Article

# On the Trend in Below-Cloud Solar Irradiance in The Netherlands versus That in Aerosol Sulfate Concentration

Suzanne Crumeyrolle <fnref n="1" />,*, Andrey Khlystov <fnref n="2" /> and Harry Ten Brink <fnref n="3" />

**Abstract:** The below-cloud irradiance in The Netherlands increased by over 10 Wm⁻<fnref n="2" /> in the past half century. It was hypothesized that this could be due to a decrease in the aerosol serving as cloud condensation nuclei, on which the cloud droplets form, in the following way. With unchanged macrophysics, clouds with a lower number of droplets are less reflective, in other words, they transmit more solar radiation. This hypothesis cannot be substantiated with data because of a generic absence of long-term information on cloud droplet number concentrations (CDNCs) worldwide. To assess the historic trend in CDNC, we used the Boucher and Lohmann (B&L) empirical relationship between CDNC and the mass concentration of the water-attracting hygroscopic aerosol components. The B&L parameterization was tested and validated with observations from the CHIEF cloud chamber, in which the formation of marine stratocumulus, the most frequent cloud type in The Netherlands, is simulated. This study will focus on summer periods because the irradiance governs the yearly average at this latitude. The summer trend of sulfate, the most dominant hygroscopic component of observed aerosol mass concentration, was analyzed with EBAS data from 1972 onwards. The average summer CDNCs were then calculated via the B&L parameterization, showing an upper limit of 380 cm⁻<fnref n="3" /> in the 1970s and decreasing to around 200 cm⁻<fnref n="3" /> in the 2010s. The associated increase in transmission for thin marine stratocumulus without overlying clouds would be, at most, 3.5 W m⁻<fnref n="2" />. Unobstructed stratocumuli occur only part of the time, and the change in irradiance based on the reduction in cloud droplet number is certainly small in comparison to the empirically derived trend of 10 W m⁻<fnref n="2" />.

**Keywords:** aerosol-cloud interaction; cloud chamber; stratocumulus

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Crumeyrolle, S.; Khlystov, A.; Ten Brink, H. On the Trend in Below-Cloud Solar Irradiance in The Netherlands versus That in Aerosol Sulfate Concentration. *Atmosphere* **2022**, *13*, 2037. https://doi.org/10.3390/atmos13122037

Academic Editors: Wenhao Xue, Jing Zhang, Jing Wei and Yulu Tian

Received: 28 September 2022

Accepted: 28 November 2022

Published: 4 December 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Aerosols have a significant impact on the Earth's radiation budget by scattering and absorbing solar radiation (direct effect), leading to negative (cooling) and positive (warming) radiative forcing. Moreover, interactions between aerosols and clouds (indirect effect) have a large influence on the role of aerosols on climate change. The atmospheric aerosols affect cloud formation through their role as cloud condensation nuclei (CCN). Sulfate aerosols are nowadays considered as the most dominant hygroscopic component of aerosols [1,2]. Sulfate is an important secondary aerosol formed from photochemical reactions of sulfur dioxide (SO₂) emissions in the atmosphere. Summertime peaks of sulfate concentrations are common for most areas, including The Netherlands [3], due to available solar insolation, chemical reactions facilitated in high relative humidity environments, and stagnation events e.g., [4,5]. To reduce SO₂ emissions over the world, regulatory and legislative mandates were established and total annual SO₂ emissions over Europe have decreased by 72% from 1980 to 2004 [6,7], and by 60% from 1990 to 2010 over the United States [8].

Over the past 50 years, the solar irradiance over The Netherlands significantly increased, both during clear-sky and cloudy situations [9]. The solar irradiance increase

<footnote n="1">¹ CNRS, UMR 851 LOA-Laboratoire d'Optique Atmosphérique, University Lille, 59000 Lille, France</footnote>

<footnote n="2">² Desert Research Institute DRI, Reno, NV 89512, USA</footnote>

<footnote n="3">³ Energy-Research Centre of The Netherlands ECN, 1755 ZG Petten, The Netherlands</footnote>

<footnote>* Correspondence: suzanne.crumeyrolle@univ-lille.fr</footnote>
during clear sky conditions could be attributed to a decrease in reflective aerosols within the atmosphere, such as sulfate [10,11]. Whereas, the solar irradiance decrease under cloudy conditions may have been caused by a decrease in CDNC, leading to a reduction in cloud reflectivity [12–15]. Assuming a constant liquid water path (LWP), a decrease in CDNC would be due to a lower CCN concentration. In this case, CDNC and below-cloud irradiance trends should be similar.

Over the North Sea region, stratocumulus (SC) is the most common [16] and the most frequent cloud type. Together with stratus, it covers the North Sea nearly 50% of the time [16]. SC generally forms at the top of the planetary boundary layer, and is therefore connected to aerosol observed at the surface. Moreover, previous studies, based on remote sensing, parcel modeling, and large eddy simulation (LES), have highlighted that SC clouds exhibit substantial albedo susceptibility to variations in CDNC [17,18].

Trend analyses are needed to track the evolution of solar radiation and CDNC. In a trend analysis, the annual average of the solar irradiance measured values cannot be used because of its variation, which depends on the solar angle and cloud coverage. For these reasons, cloud cover values were scaled to the yearly average of cloud cover and occurrence of clouds per solar angle [9]. The parameter obtained is defined as a "proxy" for below-cloud irradiance per year. At The Netherlands latitudes, the annual solar irradiance is determined by that in summer because of the long sunshine hours and high irradiance. This leads to a difference of a factor of ten in irradiance in summer and winter [19].

Contrary to solar irradiation, the CDNC trend cannot be evaluated from measurements because of an absence of monitoring data. This is not a local issue but a worldwide phenomenon. To assess the historic CDNC trend, de Martino et al. [20] detailed several approaches, including using the mass concentration of the hygroscopic aerosol components as a proxy. In a recent evaluation of the largest dataset in marine clouds in South California [21], the parameterization of Boucher and Lohmann (B&L) [22], based on the mass concentration of sulfate to retrieve CDNC, was found to provide the best results in comparison to the data. This was also the case in a shorter study off the Canadian coast of the North Atlantic [23].

There is very limited information on the CDNC over the North Sea region, let alone in combination with sulfate concentrations. Only some data exist on the CDNC in the 1980s [24], mostly collected in continental air masses non-representative of the background concentrations. In view of this, we derived an evaluation of the relationship between CDNC and sulfate concentration in a large cloud chamber (CHIEF) at the northwestern seacoast of The Netherlands, in which the conditions of marine stratocumulus were simulated [25].

In the present investigation, we first tested whether the B&L parameterization was appropriate for the relevant marine clouds over The Netherlands via an investigation in a large cloud chamber (CHIEF). Then, the summer trend of the sulfate concentration was reconstructed using EBAS data to infer the CDNC trend using the B&L parameterization. Finally, the CDNC trend was used to estimate the modification in solar irradiance only due to indirect effect, and will be compared to below-cloud irradiance trend. This will provide information on the magnitude of the aerosol–cloud interaction (ACI) effect on climate [3]. Before we proceed, it is important to mention that in addition to sulfate, other chemical compounds in aerosols, such as nitrate [20], might be of importance as components of the cloud nuclei. However, we showed in a recent comprehensive study [25] that, in The Netherlands, ammonium nitrate is not relevant during the summertime, and neither are organic compounds or sea salt; even though they significantly contribute to the total aerosol mass, they are not present in the size range of particles that serve as cloud nuclei.

## 2. Materials and Methods

### 2.1. Cloud Chamber

The large cloud chamber (CHIEF) used in the tests is situated at the northwestern seacoast of The Netherlands. A full description of the set-up is provided in an earlier publication [25]. Here we provide a summary of the relevant aspects. It is a 21 m³ flow-
during clear sky conditions could be attributed to a decrease in reflective aerosols within the atmosphere, such as sulfate [10,11]. Whereas, the solar irradiance decrease under cloudy conditions may have been caused by a decrease in CDNC, leading to a reduction in cloud reflectivity [12–15]. Assuming a constant liquid water path (LWP), a decrease in CDNC would be due to a lower CCN concentration. In this case, CDNC and below-cloud irradiance trends should be similar.

Over the North Sea region, stratocumulus (SC) is the most common [16] and the most frequent cloud type. Together with stratus, it covers the North Sea nearly 50% of the time [16]. SC generally forms at the top of the planetary boundary layer, and is therefore connected to aerosol observed at the surface. Moreover, previous studies, based on remote sensing, parcel modeling, and large eddy simulation (LES), have highlighted that SC clouds exhibit substantial albedo susceptibility to variations in CDNC [17,18].

Trend analyses are needed to track the evolution of solar radiation and CDNC. In a trend analysis, the annual average of the solar irradiance measured values cannot be used because of its variation, which depends on the solar angle and cloud coverage. For these reasons, cloud cover values were scaled to the yearly average of cloud cover and occurrence of clouds per solar angle [9]. The parameter obtained is defined as a "proxy" for below-cloud irradiance per year. At The Netherlands latitudes, the annual solar irradiance is determined by that in summer because of the long sunshine hours and high irradiance. This leads to a difference of a factor of ten in irradiance in summer and winter [19].

Contrary to solar irradiation, the CDNC trend cannot be evaluated from measurements because of an absence of monitoring data. This is not a local issue but a worldwide phenomenon. To assess the historic CDNC trend, de Martino et al. [20] detailed several approaches, including using the mass concentration of the hygroscopic aerosol components as a proxy. In a recent evaluation of the largest dataset in marine clouds in South California [21], the parameterization of Boucher and Lohmann (B&L) [22], based on the mass concentration of sulfate to retrieve CDNC, was found to provide the best results in comparison to the data. This was also the case in a shorter study off the Canadian coast of the North Atlantic [23].

There is very limited information on the CDNC over the North Sea region, let alone in combination with sulfate concentrations. Only some data exist on the CDNC in the 1980s [24], mostly collected in continental air masses non-representative of the background concentrations. In view of this, we derived an evaluation of the relationship between CDNC and sulfate concentration in a large cloud chamber (CHIEF) at the northwestern seacoast of The Netherlands, in which the conditions of marine stratocumulus were simulated [25].

In the present investigation, we first tested whether the B&L parameterization was appropriate for the relevant marine clouds over The Netherlands via an investigation in a large cloud chamber (CHIEF). Then, the summer trend of the sulfate concentration was reconstructed using EBAS data to infer the CDNC trend using the B&L parameterization. Finally, the CDNC trend was used to estimate the modification in solar irradiance only due to indirect effect, and will be compared to below-cloud irradiance trend. This will provide information on the magnitude of the aerosol–cloud interaction (ACI) effect on climate [3]. Before we proceed, it is important to mention that in addition to sulfate, other chemical compounds in aerosols, such as nitrate [20], might be of importance as components of the cloud nuclei. However, we showed in a recent comprehensive study [25] that, in The Netherlands, ammonium nitrate is not relevant during the summertime, and neither are organic compounds or sea salt; even though they significantly contribute to the total aerosol mass, they are not present in the size range of particles that serve as cloud nuclei.

## 2. Materials and Methods

### 2.1. *Cloud Chamber*

The large cloud chamber (CHIEF) used in the tests is situated at the northwestern seacoast of The Netherlands. A full description of the set-up is provided in an earlier publication [25]. Here we provide a summary of the relevant aspects. It is a 21 m³ flow-
through chamber in which the low supersaturation conditions of marine stratocumulus are simulated [25]. The proper setting of the chamber was derived from simultaneous CDNC measurements within the chamber and aboard a small aircraft flying in stratocumulus over the site. These results are detailed for the first time in Appendix A. The CDNC was measured using a forward scattering spectrometer probe (FSSP-100 from DMT) for droplets within the diameter range of 3 to 47 µm. Size and number concentration were calibrated according to the guidelines in the manual. The sulfate concentration was observed with a MARGA-sizer that works as follows: particles are grown to large droplets by steam condensation, and these are collected in a cyclone [26,27]. In the drain water, the hygroscopic species dissolve and are then measured online as ions. The accuracy of the MARGA system has been checked against 24-h filter measurements, e.g., [28,29], with an average deviation of less than 3%.

## 2.2. Historic Sulfate Measurements

The sulfate mass concentration has been measured since July 1972 in The Netherlands, with many of the data included in the EMEP/EBAS database [30]. A previous study [31] already conducted an analysis comparing the trend of observed sulfate mass concentration to the reflection of radiation during summer over the years 1978–1994. Since the publication of this study, the database has evolved to homogenize the mass concentration values as a function of the method used to measure it. In the early years, sulfate concentrations were measured through X-ray fluorescence (XRF), detecting sulfur, while in the period 1984–1987, a colorization reaction of the leached filter was applied [31]. Thereafter, ion chromatography was used, though data for the first 4 years of application of the latter method are not in the EBAS database. In 2009, PM₁₀ samplers were introduced within the national network. An issue with this is that these collect a larger fraction of sea salt than the original samplers. To correct this bias, the PM₁₀ samples were also analyzed for sodium, which allows the estimation of the sea salt sulfate fraction.

# 3. Results and Discussion

## 3.1. Cloud Chamber Study

Figure 1 shows the total sulfate concentrations and CDNC formed within the cloud chamber CHIEF during an intensive observation period in summer 2007. One can see that sulfate mass concentration ranges from 0.5 to 6 µg m⁻³, while the corresponding CDNC (only) spans from 180 to 400 cm⁻³. The large variability of the sulfate mass concentration is clearly not observed within the CDNC time series. Indeed, the CDNC did not significantly increase for sulfate concentrations exceeding 3 µg m⁻³. However, one can see that the tendency is similar in both figures, with two large peaks observed on the 1st and on 12 August 2007 for both parameters. On the 5 August 2007, the CDNC peak reached up to 366 cm⁻³, while the sulfate was unfortunately not measured due to technical issues. However, the sulfate concentration before and after this date is rather large (>3 µg cm⁻³), probably suggesting a peak on the 5 August.

As for the functional relationship of the CDNC and the mass concentration, we refer to Figure 2 (blue cross), in which the relationship of Leaitch et al. [23] is added. These data derive from the concentration of sulfate in the cloud water translated to the corresponding air concentration. For this translation, the liquid water content (LWC) is used, for which the average uncertainty is 20%, as we deduce from the information given in the publication. This constitutes by far the largest measurement uncertainty. The Boucher and Lohmann (B&L) parameterization (dashed lines) for this relationship is provided, which thus holds for polluted coastal aerosol at the northwestern North Atlantic coast. It should be mentioned that the asymptotic value of the CDNC at the highest concentration agrees with that measured by Martin and Johnson [24] in the North Sea area (CDNC ~ 350–400 cm⁻³), substantiating the choice for the chosen parameterization. It seems that the observed CDNC in the cloud chamber is larger than the B&L parameterization. This can be explained by the fact that CHIEF, with its residence time of 2 min, simulates the situation at the cloud
through chamber in which the low supersaturation conditions of marine stratocumulus are simulated [25]. The proper setting of the chamber was derived from simultaneous CDNC measurements within the chamber and aboard a small aircraft flying in stratocumulus over the site. These results are detailed for the first time in Appendix A. The CDNC was measured using a forward scattering spectrometer probe (FSSP-100 from DMT) for droplets within the diameter range of 3 to 47 µm. Size and number concentration were calibrated according to the guidelines in the manual. The sulfate concentration was observed with a MARGA-sizer that works as follows: particles are grown to large droplets by steam condensation, and these are collected in a cyclone [26,27]. In the drain water, the hygroscopic species dissolve and are then measured online as ions. The accuracy of the MARGA system has been checked against 24-h filter measurements, e.g., [28,29], with an average deviation of less than 3%.

## 2.2. Historic Sulfate Measurements

The sulfate mass concentration has been measured since July 1972 in The Netherlands, with many of the data included in the EMEP/EBAS database [30]. A previous study [31] already conducted an analysis comparing the trend of observed sulfate mass concentration to the reflection of radiation during summer over the years 1978–1994. Since the publication of this study, the database has evolved to homogenize the mass concentration values as a function of the method used to measure it. In the early years, sulfate concentrations were measured through X-ray fluorescence (XRF), detecting sulfur, while in the period 1984–1987, a colorization reaction of the leached filter was applied [31]. Thereafter, ion chromatography was used, though data for the first 4 years of application of the latter method are not in the EBAS database. In 2009, PM₁₀ samplers were introduced within the national network. An issue with this is that these collect a larger fraction of sea salt than the original samplers. To correct this bias, the PM₁₀ samples were also analyzed for sodium, which allows the estimation of the sea salt sulfate fraction.

# 3. Results and Discussion

## 3.1. Cloud Chamber Study

Figure 1 shows the total sulfate concentrations and CDNC formed within the cloud chamber CHIEF during an intensive observation period in summer 2007. One can see that sulfate mass concentration ranges from 0.5 to 6 µg m⁻³, while the corresponding CDNC (only) spans from 180 to 400 cm⁻³. The large variability of the sulfate mass concentration is clearly not observed within the CDNC time series. Indeed, the CDNC did not significantly increase for sulfate concentrations exceeding 3 µg m⁻³. However, one can see that the tendency is similar in both figures, with two large peaks observed on the 1st and on 12 August 2007 for both parameters. On the 5 August 2007, the CDNC peak reached up to 366 cm⁻³, while the sulfate was unfortunately not measured due to technical issues. However, the sulfate concentration before and after this date is rather large (>3 µg cm⁻³), probably suggesting a peak on the 5 August.

As for the functional relationship of the CDNC and the mass concentration, we refer to Figure 2 (blue cross), in which the relationship of Leaitch et al. [23] is added. These data derive from the concentration of sulfate in the cloud water translated to the corresponding air concentration. For this translation, the liquid water content (LWC) is used, for which the average uncertainty is 20%, as we deduce from the information given in the publication. This constitutes by far the largest measurement uncertainty. The Boucher and Lohmann (B&L) parameterization (dashed lines) for this relationship is provided, which thus holds for polluted coastal aerosol at the northwestern North Atlantic coast. It should be mentioned that the asymptotic value of the CDNC at the highest concentration agrees with that measured by Martin and Johnson [24] in the North Sea area (CDNC ~ 350–400 cm⁻³), substantiating the choice for the chosen parameterization. It seems that the observed CDNC in the cloud chamber is larger than the B&L parameterization. This can be explained by the fact that CHIEF, with its residence time of 2 min, simulates the situation at the cloud
base, where the CDNC is higher up in the cloud, as found in our cloud flight presented in Appendix A and by Leaitch et al. [23]. Hence, CDNCs measured within the chamber are indeed expected to be higher than B&L retrievals. Incidentally, the presence of a higher concentration at the cloud base is of little relevance for the reflectivity, because it concerns very small droplets that evaporate upon further upward transport.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Sulfate mass concentration and cloud droplet number concentration (CDNC) measured within the cloud chamber (CHIEF) during the indicated intensive operational period in summer 2007.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 2.** CDNC as a function of sulfate concentration measured in the CHIEF cloud chamber (blue cross) and observed during the stratocumulus flights of Leaitch et al. [23] (red squares and stars). The dashed line represents the Boucher and Lohmann [22] parameterization.</figcaption></figure>
base, where the CDNC is higher up in the cloud, as found in our cloud flight presented in Appendix A and by Leaitch et al. [23]. Hence, CDNCs measured within the chamber are indeed expected to be higher than B&L retrievals. Incidentally, the presence of a higher concentration at the cloud base is of little relevance for the reflectivity, because it concerns very small droplets that evaporate upon further upward transport.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Sulfate mass concentration and cloud droplet number concentration (CDNC) measured within the cloud chamber (CHIEF) during the indicated intensive operational period in summer 2007.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** CDNC as a function of sulfate concentration measured in the CHIEF cloud chamber (blue cross) and observed during the stratocumulus flights of Leaitch et al. [23] (red squares and stars). The dashed line represents the Boucher and Lohmann [22] parameterization.</figcaption></figure>
As the CDNC follows the B&L parameterization rather well, the B&L parameterization could be used for the CDNC trend over The Netherlands. Therefore, in the next sections, the B&L parameterization will be used to estimate the trend in CDNC based on measured sulfate mass concentrations over the past ~50 years. It is also interesting to note, in Figure 2, that the maximum CDNC values observed within the CHIEF chamber are close to 350 cm⁻³ for sulfate concentrations larger than 1 µg m⁻³. This upper limit can also be seen in Leaitch data, with a maximum observed at 380 cm⁻³.

## 3.2. Sulfate Mass Concentration and CDNC Trends

All in all, non-sea-salt (nss) sulfate mass concentrations are available for eleven sites operational at different periods, with relatively often only one location active (see supplementary material for more information). There is a good correlation of data simultaneously collected at the different locations (Figure S2), which is consistent with the previous analysis [31] for The Netherlands in the 1980s and for stations, as far apart as 200 km, in western Europe in general [32]. Median nss-sulfate mass concentration per summer from only three sites (Witteveen, Appelscha, Vredepeel) were used to reconstruct the sulfate trend presented in Figure 3a. One can clearly see that the nss-sulfate trend is clearly decreasing with time. The median summer concentration in 1972 was around 10.2 µg m⁻³, and decreased down to 1.9 µg m⁻³ in 2019. This decrease is essentially due to stakeholder decisions at the European scale, especially the gasoline restrictions [10]. Moreover, one can clearly see that the daily average variability (grey dots) is also decreasing with time due to the decrease in daily PM exceedances also due to stakeholder decisions. A similar trend is observed over the southern part of the UK. For the years 1966 to 1971, the average summer concentration of nss-sulfate in the southern UK based on data at Chilton is close to 7 µg m⁻³ [33]. Currently (2016–2021), the concentration of nss-sulfate, observed in nearby Chilbolton, is just over 2 µg m⁻³.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** (a) Time series of non-sea-salt (nss) sulfate mass concentration observed over The Netherlands during the summer season (June–August) from three sites (Witteveen, triangles; Appelscha, dots; Vredepeel, squares; see supplementary material for more details). The black signs correspond to the median over the summer season and the grey signs to the daily averages. (b) Time series of CDNC derived from the Boucher and Lohmann relationship [22] applied to the nss-sulfate concentration shown in (a).</figcaption></figure>

To estimate the summer CDNC trend, the relationship (B&L (A) and (D), shown in Figure 2) between CDNC and nss-sulfate mass concentration is used (Figure 3b). With a current-day nss-sulfate concentration of around 1.9 µg m⁻³ over The Netherlands, the CDNC should be in the order of 200 cm⁻³. Meanwhile, in the 1970s, the nss-sulfate concentration was over 10 µg m⁻³ and the CDNC was retrieved at or even over the asymptotic maximum value of 380 cm⁻³. A clear CDNC decrease, by a ratio up to 2, is observed during the last 40 years.

The trend in the measured average annual irradiance [9] is reproduced in Figure 4a. From a plateau at 76 ± 3 W m⁻² during the period 1966–1996, it increased to a level of
As the CDNC follows the B&L parameterization rather well, the B&L parameterization could be used for the CDNC trend over The Netherlands. Therefore, in the next sections, the B&L parameterization will be used to estimate the trend in CDNC based on measured sulfate mass concentrations over the past ~50 years. It is also interesting to note, in Figure 2, that the maximum CDNC values observed within the CHIEF chamber are close to 350 cm⁻³ for sulfate concentrations larger than 1 µg m⁻³. This upper limit can also be seen in Leaitch data, with a maximum observed at 380 cm⁻³.

## 3.2. Sulfate Mass Concentration and CDNC Trends

All in all, non-sea-salt (nss) sulfate mass concentrations are available for eleven sites operational at different periods, with relatively often only one location active (see supplementary material for more information). There is a good correlation of data simultaneously collected at the different locations (Figure S2), which is consistent with the previous analysis [31] for The Netherlands in the 1980s and for stations, as far apart as 200 km, in western Europe in general [32]. Median nss-sulfate mass concentration per summer from only three sites (Witteveen, Appelscha, Vredepeel) were used to reconstruct the sulfate trend presented in Figure 3a. One can clearly see that the nss-sulfate trend is clearly decreasing with time. The median summer concentration in 1972 was around 10.2 µg m⁻³, and decreased down to 1.9 µg m⁻³ in 2019. This decrease is essentially due to stakeholder decisions at the European scale, especially the gasoline restrictions [10]. Moreover, one can clearly see that the daily average variability (grey dots) is also decreasing with time due to the decrease in daily PM exceedances also due to stakeholder decisions. A similar trend is observed over the southern part of the UK. For the years 1966 to 1971, the average summer concentration of nss-sulfate in the southern UK based on data at Chilton is close to 7 µg m⁻³ [33]. Currently (2016–2021), the concentration of nss-sulfate, observed in nearby Chilbolton, is just over 2 µg m⁻³.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) Time series of non-sea-salt (nss) sulfate mass concentration observed over The Netherlands during the summer season (June–August) from three sites (Witteveen, triangles; Appelscha, dots; Vredepeel, squares; see supplementary material for more details). The black signs correspond to the median over the summer season and the grey signs to the daily averages. (b) Time series of CDNC derived from the Boucher and Lohmann relationship [22] applied to the nss-sulfate concentration shown in (a).</figcaption></figure>

To estimate the summer CDNC trend, the relationship (B&L (A) and (D), shown in Figure 2) between CDNC and nss-sulfate mass concentration is used (Figure 3b). With a current-day nss-sulfate concentration of around 1.9 µg m⁻³ over The Netherlands, the CDNC should be in the order of 200 cm⁻³. Meanwhile, in the 1970s, the nss-sulfate concentration was over 10 µg m⁻³ and the CDNC was retrieved at or even over the asymptotic maximum value of 380 cm⁻³. A clear CDNC decrease, by a ratio up to 2, is observed during the last 40 years.

The trend in the measured average annual irradiance [9] is reproduced in Figure 4a. From a plateau at 76 ± 3 W m⁻² during the period 1966–1996, it increased to a level of
approximately 87 ± 2.5 W m⁻² during the period 2000–2015. This difference (increase by 16%) is significant and far above the uncertainties. Therefore, the below-cloud solar radiation is clearly increasing over The Netherlands. As the CDNC is decreasing during the same period, the hypothesis of lower reflectivity clouds due to lower CDNC could explain this tendency. The cloud transmission would be larger, and then the total radiation below the cloud would increase as well.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** (a) Absolute cloud base proxy radiation as described by Boers et al. [9] over the indicated years averaged for five measuring sites. (b) Change in solar irradiance as deduced from the decrease in CDNC (Figure 3b).</figcaption></figure>

In order to estimate this effect, we used the relationship provided by Schwartz and Slingo [34]. They established (as demonstrated in Equation (1)), based on older studies [31,35], that the decrease in reflectivity (ΔR) is associated with CDNC decrease considering a constant liquid water path (LWP) and cloud thickness (Zc). This equation is accurate within 10% for the range of cloud-top albedo, 0.28 < R < 0.72, characteristic of the prevalent and climatically important marine stratus clouds. It should also be noted that this relationship applies in the absence of overlying higher clouds. As the cloud profiles are not observed, we assume that there are only unobstructed stratocumulus clouds over The Netherlands, which means that the calculated reflectivity evolution is an upper limit.

$$ (\Delta R)_{LWP,Zc} = 0.075 * \Delta \ln(\text{CDNC}) \quad (1) $$

Assuming a constant LWP and cloud thickness throughout the last 50 summers over The Netherlands, the estimated SC reflectivity decreases by 3.5 Wm⁻² in 2019 due to the decrease in CDNC. To add some perspective regarding the indirect forcing retrieved here, the total greenhouse gas forcing over the industrial period was about 2.5 W m⁻², with CO₂ contributing 1.5 W m⁻² [36]. The indirect forcing on cloud reflectivity is therefore not negligible. As previously stated, this value is an upper estimate, given the fact that SC clouds observed over The Netherlands are not permanently unobstructed. Moreover, this estimated value is far below the observed increase in the below-cloud irradiance (around 10 Wm⁻²). Therefore, the increase in the below-cloud irradiance can only be partly explained by the aerosol and cloud interaction in the lower layers. The rest of the increase could be due to a decrease in LWP and/or a decrease of cloud depth. Another possibility might be a systematic change in cloudiness values as a result of switching from a visual appraisal to objective probing by ceilometers [9]. This latter hypothesis could be checked via data for neighboring countries in which a similar trend in the mass concentration of sulfate occurred, but where cloudiness was very likely appreciated in a different manner over the years as compared to in The Netherlands.

## 4. Conclusions

In this study, Boucher and Lohmann parameterizations [22] were tested using cloud chamber (CHIEF) measurements of non-sea-salt (nss) sulfate mass concentration and cloud
approximately 87 ± 2.5 W m⁻² during the period 2000–2015. This difference (increase by 16%) is significant and far above the uncertainties. Therefore, the below-cloud solar radiation is clearly increasing over The Netherlands. As the CDNC is decreasing during the same period, the hypothesis of lower reflectivity clouds due to lower CDNC could explain this tendency. The cloud transmission would be larger, and then the total radiation below the cloud would increase as well.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** (a) Absolute cloud base proxy radiation as described by Boers et al. [9] over the indicated years averaged for five measuring sites. (b) Change in solar irradiance as deduced from the decrease in CDNC (Figure 3b).</figcaption></figure>

In order to estimate this effect, we used the relationship provided by Schwartz and Slingo [34]. They established (as demonstrated in Equation (1)), based on older studies [31,35], that the decrease in reflectivity (ΔR) is associated with CDNC decrease considering a constant liquid water path (LWP) and cloud thickness (Zc). This equation is accurate within 10% for the range of cloud-top albedo, 0.28 < R < 0.72, characteristic of the prevalent and climatically important marine stratus clouds. It should also be noted that this relationship applies in the absence of overlying higher clouds. As the cloud profiles are not observed, we assume that there are only unobstructed stratocumulus clouds over The Netherlands, which means that the calculated reflectivity evolution is an upper limit.

$$ (\Delta R)_{LWP,Zc} = 0.075 * \Delta \ln(\text{CDNC}) \quad (1) $$

Assuming a constant LWP and cloud thickness throughout the last 50 summers over The Netherlands, the estimated SC reflectivity decreases by 3.5 Wm⁻² in 2019 due to the decrease in CDNC. To add some perspective regarding the indirect forcing retrieved here, the total greenhouse gas forcing over the industrial period was about 2.5 W m⁻², with CO₂ contributing 1.5 W m⁻² [36]. The indirect forcing on cloud reflectivity is therefore not negligible. As previously stated, this value is an upper estimate, given the fact that SC clouds observed over The Netherlands are not permanently unobstructed. Moreover, this estimated value is far below the observed increase in the below-cloud irradiance (around 10 Wm⁻²). Therefore, the increase in the below-cloud irradiance can only be partly explained by the aerosol and cloud interaction in the lower layers. The rest of the increase could be due to a decrease in LWP and/or a decrease of cloud depth. Another possibility might be a systematic change in cloudiness values as a result of switching from a visual appraisal to objective probing by ceilometers [9]. This latter hypothesis could be checked via data for neighboring countries in which a similar trend in the mass concentration of sulfate occurred, but where cloudiness was very likely appreciated in a different manner over the years as compared to in The Netherlands.

## 4. Conclusions

In this study, Boucher and Lohmann parameterizations [22] were tested using cloud chamber (CHIEF) measurements of non-sea-salt (nss) sulfate mass concentration and cloud
droplet number concentration (CDNC). As expected, the CDNC measured in CHIEF, mimicking the cloud formation at its base, is overestimated in comparison to the B&L retrievals. Despite this overestimation, the results highlight that the B&L retrievals are reasonably comparable to our observations. The CDNC trend is then estimated from the B&L parameterizations, and the nss-sulfate mass concentration observed over The Netherlands for the past 50 years. Both the nss-sulfate and the CDNC trends are decreasing over the years. The CDNC was around $380 \text{ cm}^{-3}$ in the 1970s, and decreased down to $200 \text{ cm}^{-3}$, while the nss-sulfate mass concentration has decreased by a factor of five over the last 50 years (from 10 to $2 \text{ µg m}^{-3}$).

The CDNC evolution was then used to estimate the trend in stratocumulus cloud transmissivity. The below-cloud irradiance estimation increases by $3.5 \text{ W m}^{-2}$, at most, due to the presence of less CDNC. This increase can only explain part (35%) of the observed increase in below-cloud solar irradiance of $10 \text{ W m}^{-2}$. Therefore, there must be other parameters, such as cloud density or the appreciation of cloudiness over the years, that would explain this large modification of the solar irradiance received at the surface during cloudy days.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/atmos13122037/s1, Figure S1: Sulfate monitoring sites listed in EBAS for The Netherlands (Table S1). The grey dots correspond to the sites used to build the sulfate long term trend. The diamond corresponds to Cabauw, Table S1: Netherland sites listed in EBAS data platform and the associated instrumentation used to measure sulfate. Table S2: Sulfate concentration ($\mu\text{g m}^{-3}$) measured at EBAS-NL sites (Table S1) since 1972. Figure S2: Comparison of sulfate concentration ($\mu\text{g m}^{-3}$) measured at different EBAS-NL sites (Table S1). The dashed line corresponds to the 1:1 line. Figure S3: Vertical profiles of LWC in the vicinity of Petten. Line indicates the ideal adiabatic relation between LWC and altitude.

**Author Contributions:** Conceptualization, H.T.B., A.K. and S.C.; methodology, H.T.B. and A.K.; formal analysis, H.T.B. and A.K.; investigation, H.T.B. and S.C.; writing—original draft preparation, H.T.B. and S.C.; writing—review and editing, H.T.B. and S.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** The first author has been financially supported by the CaPPA project (Chemical and Physical Properties of the Atmosphere), which is funded by the French National Research Agency (ANR) through the PIA (Programme d'Investissement d'Avenir) under contract "ANR-11-LABX-0005-01".

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Sulfate measurements over The Netherlands are available through the EBAS database (https://ebasnilu.no, accessed on 9 November 2021).

**Acknowledgments:** This work is a contribution to the CaPPA project (Chemical and Physical Properties of the Atmosphere), funded by the ANR through the PIA (Programme d'Investissement d'Avenir) under contract ANR-11-LABX-0005-01.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Appendix A. Cloud Base and Cloud Chamber CDNC Measurements

We present here the set-up and results of an aircraft study in low clouds over the site of Petten, where the CHIEF cloud chamber is located. The aim was to assess at which setting of the chamber a similar CDNC is observed as in the cloud layer above. The campaign took place in September 1996 with specific evaluation of the flight on September 4th with a closed cloud deck. The research aircraft was a PA 31 Navajo Chieftain two engine turbo-prop. Several horizontal transects were constructed in which the altitude was kept constant. The droplet number concentration was measured with a FSSP-100 (PMS Inc., Las Cruces, NM, USA), working at a frequency of 2 Hz, similar to that installed in the cloud chamber.
droplet number concentration (CDNC). As expected, the CDNC measured in CHIEF, mimicking the cloud formation at its base, is overestimated in comparison to the B&L retrievals. Despite this overestimation, the results highlight that the B&L retrievals are reasonably comparable to our observations. The CDNC trend is then estimated from the B&L parameterizations, and the nss-sulfate mass concentration observed over The Netherlands for the past 50 years. Both the nss-sulfate and the CDNC trends are decreasing over the years. The CDNC was around $380 \text{ cm}^{-3}$ in the 1970s, and decreased down to $200 \text{ cm}^{-3}$, while the nss-sulfate mass concentration has decreased by a factor of five over the last 50 years (from 10 to $2 \text{ µg m}^{-3}$).

The CDNC evolution was then used to estimate the trend in stratocumulus cloud transmissivity. The below-cloud irradiance estimation increases by $3.5 \text{ W m}^{-2}$, at most, due to the presence of less CDNC. This increase can only explain part (35%) of the observed increase in below-cloud solar irradiance of $10 \text{ W m}^{-2}$. Therefore, there must be other parameters, such as cloud density or the appreciation of cloudiness over the years, that would explain this large modification of the solar irradiance received at the surface during cloudy days.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/atmos13122037/s1, Figure S1: Sulfate monitoring sites listed in EBAS for The Netherlands (Table S1). The grey dots correspond to the sites used to build the sulfate long term trend. The diamond corresponds to Cabauw, Table S1: Netherland sites listed in EBAS data platform and the associated instrumentation used to measure sulfate. Table S2: Sulfate concentration ($\mu\text{g m}^{-3}$) measured at EBAS-NL sites (Table S1) since 1972. Figure S2: Comparison of sulfate concentration ($\mu\text{g m}^{-3}$) measured at different EBAS-NL sites (Table S1). The dashed line corresponds to the 1:1 line. Figure S3: Vertical profiles of LWC in the vicinity of Petten. Line indicates the ideal adiabatic relation between LWC and altitude.

**Author Contributions:** Conceptualization, H.T.B., A.K. and S.C.; methodology, H.T.B. and A.K.; formal analysis, H.T.B. and A.K.; investigation, H.T.B. and S.C.; writing—original draft preparation, H.T.B. and S.C.; writing—review and editing, H.T.B. and S.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** The first author has been financially supported by the CaPPA project (Chemical and Physical Properties of the Atmosphere), which is funded by the French National Research Agency (ANR) through the PIA (Programme d'Investissement d'Avenir) under contract "ANR-11-LABX-0005-01".

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Sulfate measurements over The Netherlands are available through the EBAS database (https://ebasnilu.no, accessed on 9 November 2021).

**Acknowledgments:** This work is a contribution to the CaPPA project (Chemical and Physical Properties of the Atmosphere), funded by the ANR through the PIA (Programme d'Investissement d'Avenir) under contract ANR-11-LABX-0005-01.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Appendix A. Cloud Base and Cloud Chamber CDNC Measurements

We present here the set-up and results of an aircraft study in low clouds over the site of Petten, where the CHIEF cloud chamber is located. The aim was to assess at which setting of the chamber a similar CDNC is observed as in the cloud layer above. The campaign took place in September 1996 with specific evaluation of the flight on September 4th with a closed cloud deck. The research aircraft was a PA 31 Navajo Chieftain two engine turbo-prop. Several horizontal transects were constructed in which the altitude was kept constant. The droplet number concentration was measured with a FSSP-100 (PMS Inc., Las Cruces, NM, USA), working at a frequency of 2 Hz, similar to that installed in the cloud chamber.
The aircraft encountered a closed cloud deck that had entered The Netherlands from the north. Surface winds were from the NNE for the entire period of the measurements (6:00 to 9:00 UTC). The back trajectory indicated that the air masses were mostly coming from Arctic marine regions. The cloud layer was capped by an inversion at about 1100 m according to sounding at a nearby airport. There was no observable precipitation.

Several horizontal tracks were constructed within the cloud in which the droplet number and size distribution were determined. From these measurements, the LWC and effective diameter ($D_{eff}$) were derived. The standard deviation of both parameters within one track was rather small (see Table A1). The droplet spectrum shifts towards larger sizes at higher altitudes due to the continuing growth of the cloud droplets. The vertical profile of the LWC based on the droplet spectrum measured by the FSSP in the various horizontal tracks in the layer is shown in Figure S3. The experimental vertical profile of the LWC agrees reasonably well with the adiabatic profile of an almost ideal cloud layer.

One can see that the droplet number concentration (N) is relatively constant over the cloud height, apart from that measured at cloud base. Indeed, the CDNC value at cloud base (210 cm⁻³) is higher than that higher up in the cloud layer. This observation is similar to that of Martin and Johnson [24] within SC over the same region. The cloud base measurements (862 m) are then comparable to the onset of the cloud formation as simulated in the CHIEF chamber.

Operating the CHIEF chamber in scanning mode [37], the supersaturation is ramped up and down in a cycle of 40 min. The aerosol spectrum was monitored with an SMPS, while the droplet spectrum was assessed with the same type of FSSP as on board the measuring aircraft flying above.

At the time of the overhead flight, the aerosol number concentration was 2800 cm⁻³. The aerosol had a bimodal size distribution with an Aitken mode (Dp = 50 nm) and an accumulation mode (Dp = 150 nm). The best agreement of CDNC measured at cloud base (210 cm⁻³) and within the cloud chamber was around a supersaturation of 0.1% (Figure A1). For further details, we refer to the detailed report by Khlystov et al. [38].

<figure><img src="image_9.png" /><figcaption>**Figure A1.** CDNC measured within the CHIEF cloud chamber as a function of supersaturation.</figcaption></figure>
The aircraft encountered a closed cloud deck that had entered The Netherlands from the north. Surface winds were from the NNE for the entire period of the measurements (6:00 to 9:00 UTC). The back trajectory indicated that the air masses were mostly coming from Arctic marine regions. The cloud layer was capped by an inversion at about 1100 m according to sounding at a nearby airport. There was no observable precipitation.

Several horizontal tracks were constructed within the cloud in which the droplet number and size distribution were determined. From these measurements, the LWC and effective diameter ($D_{eff}$) were derived. The standard deviation of both parameters within one track was rather small (see Table A1). The droplet spectrum shifts towards larger sizes at higher altitudes due to the continuing growth of the cloud droplets. The vertical profile of the LWC based on the droplet spectrum measured by the FSSP in the various horizontal tracks in the layer is shown in Figure S3. The experimental vertical profile of the LWC agrees reasonably well with the adiabatic profile of an almost ideal cloud layer.

One can see that the droplet number concentration (N) is relatively constant over the cloud height, apart from that measured at cloud base. Indeed, the CDNC value at cloud base (210 cm⁻³) is higher than that higher up in the cloud layer. This observation is similar to that of Martin and Johnson [24] within SC over the same region. The cloud base measurements (862 m) are then comparable to the onset of the cloud formation as simulated in the CHIEF chamber.

Operating the CHIEF chamber in scanning mode [37], the supersaturation is ramped up and down in a cycle of 40 min. The aerosol spectrum was monitored with an SMPS, while the droplet spectrum was assessed with the same type of FSSP as on board the measuring aircraft flying above.

At the time of the overhead flight, the aerosol number concentration was 2800 cm⁻³. The aerosol had a bimodal size distribution with an Aitken mode (Dp = 50 nm) and an accumulation mode (Dp = 150 nm). The best agreement of CDNC measured at cloud base (210 cm⁻³) and within the cloud chamber was around a supersaturation of 0.1% (Figure A1). For further details, we refer to the detailed report by Khlystov et al. [38].

<figure><img src="image_8.png" /><figcaption>**Figure A1.** CDNC measured within the CHIEF cloud chamber as a function of supersaturation.</figcaption></figure>
<table><caption>**Table A1.** Average and standard deviation (STD) of droplet number concentration (N), with LWC and effective diameter (D<sub>eff</sub>) derived from the indicated horizontal tracks through a SC field in the vicinity of the CHIEF cloud chamber. Grey background is used for data measured at cloud base with higher variability in D<sub>eff</sub>.</caption><thead><tr><th colspan="2">Altitude</th><th colspan="2">N (# cm<sup>-3</sup>)</th><th colspan="2">LWC (mg m<sup>-3</sup>)</th><th colspan="2">D<sub>eff</sub> (µm)</th></tr><tr><th>Average</th><th>STD</th><th>Average</th><th>STD</th><th>Average</th><th>STD</th><th>Average</th><th>STD</th></tr></thead><tbody><tr><td>862</td><td>11</td><td>210</td><td>60</td><td>170</td><td>60</td><td>6.6</td><td>1.5</td></tr><tr><td>895</td><td>3</td><td>110</td><td>35</td><td>250</td><td>50</td><td>18.8</td><td>1.3</td></tr><tr><td>895</td><td>3</td><td>90</td><td>20</td><td>250</td><td>50</td><td>19.4</td><td>1.0</td></tr><tr><td>942</td><td>26</td><td>145</td><td>70</td><td>420</td><td>165</td><td>19.9</td><td>0.6</td></tr><tr><td>970</td><td>5</td><td>120</td><td>30</td><td>530</td><td>60</td><td>22.1</td><td>0.9</td></tr><tr><td>1020</td><td>20</td><td>115</td><td>30</td><td>730</td><td>110</td><td>24.4</td><td>0.7</td></tr><tr><td>1025</td><td>5</td><td>120</td><td>40</td><td>650</td><td>145</td><td>23.5</td><td>0.9</td></tr></tbody></table>

References

1. Boucher, O.; Randall, D.; Artaxo, P.; Bretherton, C.; Feingold, G.; Forster, P.; Kerminen, V.-M.; Kondo, Y.; Liao, H.; Lohmann, U.; et al. Climate Change 2013. In *Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change*; Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013.

2. Kim, A.-H.; Yum, S.S.; Chang, D.Y.; Park, M. Optimization of the sulfate aerosol hygroscopicity parameter in WRF-Chem. *Geosci. Model Dev.* **2021**, *14*, 259–273. [CrossRef]

3. Yoon, Y.J.; Ceburnis, D.; Cavalli, F.; Jourdan, O.; Putaud, J.P.; Facchini, M.C.; Decesari, S.; Fuzzi, S.; Sellegri, K.; Jennings, S.G.; et al. Seasonal characteristics of the physicochemical properties of North Atlantic marine atmospheric aerosols. *J. Geophys. Res. Earth Surf.* **2007**, *112*, D04206. [CrossRef]

4. Hidy, G.M.; Mueller, P.K.; Tong, E.Y. Spatial and temporal distributions of airborne sulfate in parts of the United States. *Atmos. Environ.* **1978**, *12*, 735–752. [CrossRef]

5. Tai, A.P.K.; Mickley, L.J.; Jacob, D.J. Correlations between fine particulate matter (PM2.5) and meteorological variables in the United States: Implications for the sensitivity of PM2.5 to climate change. *Atmos. Environ.* **2010**, *44*, 3976–3984. [CrossRef]

6. Engardt, M.; Simpson, D.; Schwikowski, M.; Granat, L. Deposition of sulphur and nitrogen in Europe 1900–2050. Model calculations and comparison to historical observations. *Tellus B Chem. Phys. Meteorol.* **2017**, *69*, 1328945. [CrossRef]

7. Vestreng, V.; Myhre, G.; Fagerli, H.; Reis, S.; Tarrasón, L. Twenty-five years of continuous sulphur dioxide emission reduction in Europe. *Atmos. Chem. Phys.* **2007**, *7*, 3663–3681. [CrossRef]

8. Schmalensee, R.; Stavins, R.N. The SO₂ Allowance Trading System: The Ironic History of a Grand Policy Experiment. *J. Econ. Perspect.* **2013**, *27*, 103–122. [CrossRef]

9. Boers, R.; Brandsma, T.; Siebesma, A.P. Impact of aerosols and clouds on decadal trends in all-sky solar radiation over The Netherlands (1966–2015). *Atmos. Chem. Phys.* **2017**, *17*, 8081–8100. [CrossRef]

10. Hand, J.L.; Copeland, S.A.; Day, D.E.; Dillner, A.M.; Idresand, H.; Malm, W.C.; McDade, C.E.; Moore, C.T., Jr.; Pitchford, M.L.; Schichtel, B.A.; et al. IMPROVE (Interagency Monitoring of Protected Visual Environments): Spatial and Seasonal Patterns and Temporal Variability of Haze and Its Constituents in the United States; Report V, CIRA Report; Cooperative Institute for Research in the Atmosphere: Fort Collins, CO, USA, 2011; ISSN 0737-5352-87. Available online: https://vista.cira.colostate.edu/improve/Publications/Reports/2011/2011.htm (accessed on 1 July 2012).

11. Malm, W.C. Characteristics and origins of haze in the continental United States. *Earth-Sci. Rev.* **1992**, *33*, 1–36. [CrossRef]

12. Twomey, S. *Atmospheric Aerosols*; Elsevier Scientific Publishing Company: Amsterdam, The Netherlands; Oxford, UK; New York, NY, USA, 1978.

13. Lohmann, U. Aerosol Effects on Clouds and Climate. *Space Sci. Rev.* **2006**, *125*, 129–137. [CrossRef]

14. Wood, R. Stratocumulus Clouds. *Mon. Weather Rev.* **2012**, *140*, 2373–2423. [CrossRef]

15. Berner, A.H.; Bretherton, C.S.; Wood, R. Large eddy simulation of ship tracks in the collapsed marine boundary layer: A case study from the Monterey area ship track experiment. *Atmos. Chem. Phys.* **2015**, *15*, 5851–5871. [CrossRef]

16. Duynkerke, P.G.; Driedonks, A.G.M. Turbulent Structure of a Shear-Driven Stratus-Topped Atmospheric Boundary Layer: A Comparison of Model Results with Observations. *J. Atmos. Sci.* **1988**, *45*, 2343–2351. [CrossRef]

17. Sanchez, K.J.; Russell, L.M.; Modini, R.L.; Frossard, A.A.; Ahlm, L.; Corrigan, C.E.; Roberts, G.C.; Hawkins, L.N.; Schroder, J.C.; Bertram, A.K.; et al. Meteorological and aerosol effects on marine cloud microphysical properties. *J. Geophys. Res. Atmos.* **2016**, *121*, 4142–4161. [CrossRef]

18. Platnick, S.; Twomey, S. Determining the susceptibility of cloud albedo to changes in droplet concentration with the advanced very high-resolution radiometer. *J. Appl. Meteorol.* **1994**, *33*, 334–347. [CrossRef]
<table><caption>**Table A1.** Average and standard deviation (STD) of droplet number concentration (N), with LWC and effective diameter (D<sub>eff</sub>) derived from the indicated horizontal tracks through a SC field in the vicinity of the CHIEF cloud chamber. Grey background is used for data measured at cloud base with higher variability in D<sub>eff</sub>.</caption><thead><tr><th colspan="2">Altitude</th><th colspan="2">N (# cm<sup>-3</sup>)</th><th colspan="2">LWC (mg m<sup>-3</sup>)</th><th colspan="2">D<sub>eff</sub> (μm)</th></tr><tr><th>Average</th><th>STD</th><th>Average</th><th>STD</th><th>Average</th><th>STD</th><th>Average</th><th>STD</th></tr></thead><tbody><tr><td>862</td><td>11</td><td>210</td><td>60</td><td>170</td><td>60</td><td>6.6</td><td>1.5</td></tr><tr><td>895</td><td>3</td><td>110</td><td>35</td><td>250</td><td>50</td><td>18.8</td><td>1.3</td></tr><tr><td>895</td><td>3</td><td>90</td><td>20</td><td>250</td><td>50</td><td>19.4</td><td>1.0</td></tr><tr><td>942</td><td>26</td><td>145</td><td>70</td><td>420</td><td>165</td><td>19.9</td><td>0.6</td></tr><tr><td>970</td><td>5</td><td>120</td><td>30</td><td>530</td><td>60</td><td>22.1</td><td>0.9</td></tr><tr><td>1020</td><td>20</td><td>115</td><td>30</td><td>730</td><td>110</td><td>24.4</td><td>0.7</td></tr><tr><td>1025</td><td>5</td><td>120</td><td>40</td><td>650</td><td>145</td><td>23.5</td><td>0.9</td></tr></tbody></table>

References

1. Boucher, O.; Randall, D.; Artaxo, P.; Bretherton, C.; Feingold, G.; Forster, P.; Kerminen, V.-M.; Kondo, Y.; Liao, H.; Lohmann, U.; et al. Climate Change 2013. In *Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change*; Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013.

2. Kim, A.-H.; Yum, S.S.; Chang, D.Y.; Park, M. Optimization of the sulfate aerosol hygroscopicity parameter in WRF-Chem. *Geosci. Model Dev.* **2021**, *14*, 259–273. [CrossRef]

3. Yoon, Y.J.; Ceburnis, D.; Cavalli, F.; Jourdan, O.; Putaud, J.P.; Facchini, M.C.; Decesari, S.; Fuzzi, S.; Sellegri, K.; Jennings, S.G.; et al. Seasonal characteristics of the physicochemical properties of North Atlantic marine atmospheric aerosols. *J. Geophys. Res. Earth Surf.* **2007**, *112*, D04206. [CrossRef]

4. Hidy, G.M.; Mueller, P.K.; Tong, E.Y. Spatial and temporal distributions of airborne sulfate in parts of the United States. *Atmos. Environ.* **1978**, *12*, 735–752. [CrossRef]

5. Tai, A.P.K.; Mickley, L.J.; Jacob, D.J. Correlations between fine particulate matter (PM2.5) and meteorological variables in the United States: Implications for the sensitivity of PM2.5 to climate change. *Atmos. Environ.* **2010**, *44*, 3976–3984. [CrossRef]

6. Engardt, M.; Simpson, D.; Schwikowski, M.; Granat, L. Deposition of sulphur and nitrogen in Europe 1900–2050. Model calculations and comparison to historical observations. *Tellus B Chem. Phys. Meteorol.* **2017**, *69*, 1328945. [CrossRef]

7. Vestreng, V.; Myhre, G.; Fagerli, H.; Reis, S.; Tarrasón, L. Twenty-five years of continuous sulphur dioxide emission reduction in Europe. *Atmos. Chem. Phys.* **2007**, *7*, 3663–3681. [CrossRef]

8. Schmalensee, R.; Stavins, R.N. The SO₂ Allowance Trading System: The Ironic History of a Grand Policy Experiment. *J. Econ. Perspect.* **2013**, *27*, 103–122. [CrossRef]

9. Boers, R.; Brandsma, T.; Siebesma, A.P. Impact of aerosols and clouds on decadal trends in all-sky solar radiation over The Netherlands (1966–2015). *Atmos. Chem. Phys.* **2017**, *17*, 8081–8100. [CrossRef]

10. Hand, J.L.; Copeland, S.A.; Day, D.E.; Dillner, A.M.; Idresand, H.; Malm, W.C.; McDade, C.E.; Moore, C.T., Jr.; Pitchford, M.L.; Schichtel, B.A.; et al. IMPROVE (Interagency Monitoring of Protected Visual Environments): Spatial and Seasonal Patterns and Temporal Variability of Haze and Its Constituents in the United States; Report V, CIRA Report; Cooperative Institute for Research in the Atmosphere: Fort Collins, CO, USA, 2011; ISSN 0737-5352-87. Available online: https://vista.cira.colostate.edu/improve/Publications/Reports/2011/2011.htm (accessed on 1 July 2012).

11. Malm, W.C. Characteristics and origins of haze in the continental United States. *Earth-Sci. Rev.* **1992**, *33*, 1–36. [CrossRef]

12. Twomey, S. *Atmospheric Aerosols*; Elsevier Scientific Publishing Company: Amsterdam, The Netherlands; Oxford, UK; New York, NY, USA, 1978.

13. Lohmann, U. Aerosol Effects on Clouds and Climate. *Space Sci. Rev.* **2006**, *125*, 129–137. [CrossRef]

14. Wood, R. Stratocumulus Clouds. *Mon. Weather Rev.* **2012**, *140*, 2373–2423. [CrossRef]

15. Berner, A.H.; Bretherton, C.S.; Wood, R. Large eddy simulation of ship tracks in the collapsed marine boundary layer: A case study from the Monterey area ship track experiment. *Atmos. Chem. Phys.* **2015**, *15*, 5851–5871. [CrossRef]

16. Duynkerke, P.G.; Driedonks, A.G.M. Turbulent Structure of a Shear-Driven Stratus-Topped Atmospheric Boundary Layer: A Comparison of Model Results with Observations. *J. Atmos. Sci.* **1988**, *45*, 2343–2351. [CrossRef]

17. Sanchez, K.J.; Russell, L.M.; Modini, R.L.; Frossard, A.A.; Ahlm, L.; Corrigan, C.E.; Roberts, G.C.; Hawkins, L.N.; Schroder, J.C.; Bertram, A.K.; et al. Meteorological and aerosol effects on marine cloud microphysical properties. *J. Geophys. Res. Atmos.* **2016**, *121*, 4142–4161. [CrossRef]

18. Platnick, S.; Twomey, S. Determining the susceptibility of cloud albedo to changes in droplet concentration with the advanced very high-resolution radiometer. *J. Appl. Meteorol.* **1994**, *33*, 334–347. [CrossRef]
19. Velds, C.A. *Solar Irradiance in The Netherlands (in Dutch)*; Koninklijk Nederlands Meteorologisch Instituut: De Bilt, The Netherlands, 1992; ISBN 90-5210-140-X. Available online: https://cdn.knmi.nl/knmi/pdf/bibliotheek/knmipubDIV/Zonnestraling_in_Nederland.pdf (accessed on 1 September 2022).

20. De Martino, G.; van Ulft, B.; ten Brink, H.; Schaap, M.; van Meijgaard, E.; Boers, R. *An Aerosol-Cloud Module for Inclusion in the KNMI Regional Climate Model RACMO2; Scientific Report, WR 2008-05; KNMI: De Bilt, The Netherlands, 2008.*

21. MacDonald, A.B.; Mardi, A.H.; Dadashazar, H.; Aghdam, M.A.; Crosbie, E.; Jonsson, H.H.; Flagan, R.C.; Seinfeld, J.H.; Sorooshian, A. On the relationship between cloud water composition and cloud droplet number concentration. *Atmos. Chem. Phys.* **2020**, *20*, 7645–7665. [CrossRef] [PubMed]

22. Boucher, O.; Lohmann, U. The sulfate-CCN-cloud albedo effect. *Tellus B Chem. Phys. Meteorol.* **1995**, *47*, 281–300. [CrossRef]

23. Leaitch, W.R.; Banic, C.M.; Isaac, G.A.; Couture, M.D.; Liu, P.S.K.; Gultepe, I.; Li, S.-M.; Kleinman, L.; Daum, P.H.; MacPherson, J.I. Physical and chemical observations in marine stratus during the 1993 North Atlantic Regional Experiment: Factors controlling cloud droplet number concentrations. *J. Geophys. Res. Earth Surf.* **1996**, *101*, 29123–29135. [CrossRef]

24. Martin, G.M.; Johnson, D.W.; Spice, A. The Measurement and Parameterization of Effective Radius of Droplets in Warm Stratocumulus Clouds. *J. Atmos. Sci.* **1994**, *51*, 1823–1842. [CrossRef]

25. Crumeyrolle, S.; Mensah, A.; Khlystov, A.; Kos, G.; Brink, H.T. On the importance of nitrate for the droplet concentration in stratocumulus in the North-Sea region. *Atmos. Environ.* **2021**, *252*, 118278. [CrossRef]

26. ten Brink, H.; Otjes, R.; Jongejan, P.; Slanina, S. An instrument for semi-continuous monitoring of the size-distribution of nitrate, ammonium, sulphate and chloride in aerosol. *Atmos. Environ.* **2007**, *41*, 2768–2779. [CrossRef]

27. Brink, H.T.; Henzing, B.; Otjes, R.; Weijers, E. Visibility in The Netherlands during New Year's fireworks: The role of soot and salty aerosol products. *Atmos. Environ.* **2018**, *173*, 289–294. [CrossRef]

28. Rumsey, I.C.; Cowen, K.A.; Walker, J.T.; Kelly, T.J.; Hanft, E.A.; Mishoe, K.; Rogers, C.; Proost, R.; Beachley, G.M.; Lear, G.; et al. An assessment of the performance of the Monitor for AeRosols and GAses in ambient air (MARGA): A semi-continuous method for soluble compounds. *Atmos. Chem. Phys.* **2014**, *14*, 5639–5658. [CrossRef]

29. Weijers, E.P.; Kos, G.P.A.; Blom, M.J.; Otjes, R.P.; Schaap, M.; van der Swaluw, E. *Measurements of Secondary Inorganic Aerosols in The Netherlands*; ECN-report ECN-E-12-003; National Institute for Public Health and the Environment: Utrecht, The Netherlands, 2012.

30. Tsyro, S.; Gauss, M.; Hjellbrekke, A.-G.; Aas, W. PM10, PM2.5 and Individual Aerosol Components, Supplementary Material to EMEP Status Report 1/2019; The Norwegian Meteorological Institute: Oslo, Norway, 2019; Available online: www.emep.int (accessed on 1 September 2022).

31. Twomey, S. Pollution and the planetary albedo. *Atmos. Environ.* **1974**, *8*, 1251–1256. [CrossRef]

32. Arends, B.; Baard, J.; Brink, H.T. Trends in summer sulphate in Europe. *Atmos. Environ.* **1997**, *31*, 4063–4072. [CrossRef]

33. Schaap, M.; van Loon, M.; Brink, H.M.T.; Dentener, F.J.; Builtjes, P.J.H. Secondary inorganic aerosol simulations for Europe with special attention to nitrate. *Atmos. Chem. Phys.* **2004**, *4*, 857–874. [CrossRef]

34. Salmon, L.; Atkins, D.; Fisher, E.; Healy, C.; Law, D. Retrospective trend analysis of the content of U.K. air particulate material 1957–1974. *Sci. Total Environ.* **1978**, *9*, 161–199. [CrossRef]

35. Schwartz, S.E.; Slingo, A. Enhanced shortwave cloud radiative forcing due to anthropogenic aerosols. In *Clouds, Chemistry and Climate—Proceedings of NATO Advanced Research Workshop*; Crutzen, P., Ramanathan, V., Eds.; Springer: Heidelberg, Germany, 1996; pp. 191–236. Available online: http://www.homepages.ed.ac.uk/shs/Climatechange/Data%20sources/Schwarz%20and%20Slingo.pdf (accessed on 1 September 2022).

36. Twomey, S. The Influence of Pollution on the Shortwave Albedo of Clouds. *J. Atmos. Sci.* **1977**, *34*, 1149–1152. [CrossRef]

37. Shine, K.P.; Derwent, R.G.; Wuebbles, D.J.; Morcrette, J.-J. Radiative Forcing of Climate. In *Climate Change: The IPCC Scientific Assessment*; Houghton, J.T., Jenkins, G.J., Ephraums, J.J., Eds.; Cambridge University Press: New York, NY, USA, 1990; pp. 41–68.

38. Khlystov, A.Y.; Kos, G.P.; Even, A.; ten Brink, H. *Micro-Physical Properties of Stratocumulus Clouds during the CLARA Campaign, Chapter 3 of the Final Report*. 1999. Available online: https://www.researchgate.net/publication/234922811 (accessed on 1 September 2022).
19. Velds, C.A. *Solar Irradiance in The Netherlands (in Dutch)*; Koninklijk Nederlands Meteorologisch Instituut: De Bilt, The Netherlands, 1992; ISBN 90-5210-140-X. Available online: https://cdn.knmi.nl/knmi/pdf/bibliotheek/knmipubDIV/Zonnestraling_in_Nederland.pdf (accessed on 1 September 2022).

20. De Martino, G.; van Ulft, B.; ten Brink, H.; Schaap, M.; van Meijgaard, E.; Boers, R. *An Aerosol-Cloud Module for Inclusion in the KNMI Regional Climate Model RACMO2; Scientific Report, WR 2008-05; KNMI: De Bilt, The Netherlands, 2008.*

21. MacDonald, A.B.; Mardi, A.H.; Dadashazar, H.; Aghdam, M.A.; Crosbie, E.; Jonsson, H.H.; Flagan, R.C.; Seinfeld, J.H.; Sorooshian, A. On the relationship between cloud water composition and cloud droplet number concentration. *Atmos. Chem. Phys.* **2020**, *20*, 7645–7665. [CrossRef] [PubMed]

22. Boucher, O.; Lohmann, U. The sulfate-CCN-cloud albedo effect. *Tellus B Chem. Phys. Meteorol.* **1995**, *47*, 281–300. [CrossRef]

23. Leaitch, W.R.; Banic, C.M.; Isaac, G.A.; Couture, M.D.; Liu, P.S.K.; Gultepe, I.; Li, S.-M.; Kleinman, L.; Daum, P.H.; MacPherson, J.I. Physical and chemical observations in marine stratus during the 1993 North Atlantic Regional Experiment: Factors controlling cloud droplet number concentrations. *J. Geophys. Res. Earth Surf.* **1996**, *101*, 29123–29135. [CrossRef]

24. Martin, G.M.; Johnson, D.W.; Spice, A. The Measurement and Parameterization of Effective Radius of Droplets in Warm Stratocumulus Clouds. *J. Atmos. Sci.* **1994**, *51*, 1823–1842. [CrossRef]

25. Crumeyrolle, S.; Mensah, A.; Khlystov, A.; Kos, G.; Brink, H.T. On the importance of nitrate for the droplet concentration in stratocumulus in the North-Sea region. *Atmos. Environ.* **2021**, *252*, 118278. [CrossRef]

26. ten Brink, H.; Otjes, R.; Jongejan, P.; Slanina, S. An instrument for semi-continuous monitoring of the size-distribution of nitrate, ammonium, sulphate and chloride in aerosol. *Atmos. Environ.* **2007**, *41*, 2768–2779. [CrossRef]

27. Brink, H.T.; Henzing, B.; Otjes, R.; Weijers, E. Visibility in The Netherlands during New Year's fireworks: The role of soot and salty aerosol products. *Atmos. Environ.* **2018**, *173*, 289–294. [CrossRef]

28. Rumsey, I.C.; Cowen, K.A.; Walker, J.T.; Kelly, T.J.; Hanft, E.A.; Mishoe, K.; Rogers, C.; Proost, R.; Beachley, G.M.; Lear, G.; et al. An assessment of the performance of the Monitor for AeRosols and GAses in ambient air (MARGA): A semi-continuous method for soluble compounds. *Atmos. Chem. Phys.* **2014**, *14*, 5639–5658. [CrossRef]

29. Weijers, E.P.; Kos, G.P.A.; Blom, M.J.; Otjes, R.P.; Schaap, M.; van der Swaluw, E. *Measurements of Secondary Inorganic Aerosols in The Netherlands*; ECN-report ECN-E-12-003; National Institute for Public Health and the Environment: Utrecht, The Netherlands, 2012.

30. Tsyro, S.; Gauss, M.; Hjellbrekke, A.-G.; Aas, W. PM10, PM2.5 and Individual Aerosol Components, Supplementary Material to EMEP Status Report 1/2019; The Norwegian Meteorological Institute: Oslo, Norway, 2019; Available online: www.emep.int (accessed on 1 September 2022).

31. Twomey, S. Pollution and the planetary albedo. *Atmos. Environ.* **1974**, *8*, 1251–1256. [CrossRef]

32. Arends, B.; Baard, J.; Brink, H.T. Trends in summer sulphate in Europe. *Atmos. Environ.* **1997**, *31*, 4063–4072. [CrossRef]

33. Schaap, M.; van Loon, M.; Brink, H.M.T.; Dentener, F.J.; Builtjes, P.J.H. Secondary inorganic aerosol simulations for Europe with special attention to nitrate. *Atmos. Chem. Phys.* **2004**, *4*, 857–874. [CrossRef]

34. Salmon, L.; Atkins, D.; Fisher, E.; Healy, C.; Law, D. Retrospective trend analysis of the content of U.K. air particulate material 1957–1974. *Sci. Total Environ.* **1978**, *9*, 161–199. [CrossRef]

35. Schwartz, S.E.; Slingo, A. Enhanced shortwave cloud radiative forcing due to anthropogenic aerosols. In *Clouds, Chemistry and Climate—Proceedings of NATO Advanced Research Workshop*; Crutzen, P., Ramanathan, V., Eds.; Springer: Heidelberg, Germany, 1996; pp. 191–236. Available online: http://www.homepages.ed.ac.uk/shs/Climatechange/Data%20sources/Schwarz%20and%20Slingo.pdf (accessed on 1 September 2022).

36. Twomey, S. The Influence of Pollution on the Shortwave Albedo of Clouds. *J. Atmos. Sci.* **1977**, *34*, 1149–1152. [CrossRef]

37. Shine, K.P.; Derwent, R.G.; Wuebbles, D.J.; Morcrette, J.-J. Radiative Forcing of Climate. In *Climate Change: The IPCC Scientific Assessment*; Houghton, J.T., Jenkins, G.J., Ephraums, J.J., Eds.; Cambridge University Press: New York, NY, USA, 1990; pp. 41–68.

38. Khlystov, A.Y.; Kos, G.P.; Even, A.; ten Brink, H. *Micro-Physical Properties of Stratocumulus Clouds during the CLARA Campaign, Chapter 3 of the Final Report*. 1999. Available online: https://www.researchgate.net/publication/234922811 (accessed on 1 September 2022).