3.1 Flight track and meteorological conditions

ORACLES research flight tracks included in situ cloud sampling during individual ascents or descents through cloud or during a series of ascents and descents through cloud along a constant heading (sawtooth maneuvers). A case study of the fifth P-3 research flight (PRF5) flown on 6 September 2016 was used to examine aerosol and cloud properties sampled under conditions of both contact and separation between the aerosol layer and cloud tops. PRF5 was selected because it had the highest cloud profiling time among the six PRFs with at least eight cloud profiles (Table 2). Four sawtooth maneuvers (S1–S4) were flown during PRF5 (Fig. 1) along with four individual cloud profiles (P1–P4). Each sawtooth maneuver consisted of four to six individual profiles (Table 2), which were numbered sequentially (S1-1, S1-2, etc.). South-southeasterly winds (5–8 m s⁻¹) were observed at the surface and at 925 mb (Fig. 2a and b). This wind field was associated with a surface low-pressure system east of the study region centered around 17^∘ S, 13^∘ E, which resulted in advection of low clouds toward the northwest. Open and closed cells of marine stratocumulus persisted along with pockets of open cells (POCs) (Fig. 1). S1, S2, and S3 were flown along 9^∘ E in closed cells of marine stratocumulus. S4 was flown closer to the coast in a shallow boundary layer with thin closed-cell stratocumulus (Fig. 1) later in the day compared to S1–S3 (Fig. 3). Ambient temperature sampled by the aircraft sensor was 3 to 6 ^∘C higher during S2 and S3 compared to S1 because the 500 mb geopotential height and relative humidity (RH) were higher toward the north (Fig. 2b). Cloud-top height (Z_T) is identified as the highest altitude satisfying the criteria used to define cloud (N_c > 10 cm⁻³ and bulk LWC > 0.05 g m⁻³). S1, S2, and S3 had higher Z_T compared to S4 (Fig. 3) due to the advection of cold, dry continental air from the southeast and low RH (< 70 %) where S4 was flown, which resulted in cloud thinning and a shallower boundary layer (Fig. 2b and c).

Figure 4Snapshots of the boundary layer sampled below (a) S1 showing shallow cumulus and stratocumulus layers with varying bases, and (b) S4 showing stratocumulus clouds with a uniform base (NSRC/NASA Airborne Science Program).

The aircraft intermittently entered and exited cumulus clouds below the stratocumulus layer during 33 of the 71 cloud profiles flown during the IOP (Table 2), which resulted in fluctuating values of N_c and R_e, with bulk LWC < 0.05 g m⁻³. For example, during S1-3, N_c varied between 10 and 240 cm⁻³, and R_e varied between 3 and 12 µm up to 130 m below where the stratocumulus base was identified with bulk LWC > 0.05 g m⁻³. Images from a forward-facing camera on the aircraft contrast a boundary layer with multiple cloud layers (Fig. 4a; image taken at 08:53 UTC) during S1-3 and a shallow, well-mixed boundary layer capped by stratocumulus (Fig. 4b; image taken at 13:16 UTC) during S4-1. It is likely the stratocumulus layer was decoupled from the surface where S1-3 was flown because the boundary layer was deepened by the entrainment of free-tropospheric air. Subsequently, the sub-cloud layer was well-mixed with the surface and topped by shallow cumulus similar to observations by Wood (2012). The cloud base height (Z_B) for the 33 profiles was determined as the lowest altitude with N_c > 10 cm⁻³ and bulk LWC > 0.05 g m⁻³ above which a continuous cloud layer was sampled. S4 had lower Z_B (195–249 m) compared to S1 (676–691 m), S2 (534–598 m), and S3 (501–775 m) (Fig. 3).

Figure 5Five-day back trajectories from the Hybrid Single-Particle Lagrangian Integrated Trajectory model for sawtooth maneuvers flown on 6 September 2016 (a) ending at 10:00 UTC for S1–S3 at 500 $\mathrm{m}\mathrm{a}.\mathrm{m}.\mathrm{s}.\mathrm{l}.$; (b) ending at 10:00 UTC for S1–S3 at 1000 $\mathrm{m}\mathrm{a}.\mathrm{m}.\mathrm{s}.\mathrm{l}.$; and (c) ending at 13:00 UTC for S4 at 200, 500, and 2500 $\mathrm{m}\mathrm{a}.\mathrm{m}.\mathrm{s}.\mathrm{l}$.

Table 3The total (OA + ${\mathrm{SO}}_{\mathrm{4}}^{\mathrm{2}+}$ + ${\mathrm{NH}}_{\mathrm{4}}^{+}$ + ${\mathrm{NO}}_{\mathrm{3}}^{-}$) and OA M_a, PCASP N_a, and rBC and CO concentrations sampled up to 100 m below cloud base and 100 m above cloud top during four sawtooth maneuvers (S1–S4) flown on 6 September 2016. These values correspond to averages across the individual profiles flown during S1–S4. AOD was sampled during constant altitude flight legs and corresponds to the atmospheric column above the aircraft.

NA: not available.

Download Print Version | Download XLSX

3.2 Above- and below-cloud aerosol composition

For each sawtooth maneuver, the above- and below-cloud air mass source region was identified using 5 d back trajectories computed using the NOAA Hybrid Single-Particle Lagrangian Integrated Trajectory model (Stein et al., 2015) applied to the National Centers for Environmental Prediction Global Data Assimilation System model (Fig. 5). The concentrations listed in Table 3 indicate measurements up to 100 m above and below the clouds averaged across the cloud profiles for each sawtooth maneuver. The variability in above-cloud M_a and N_a for S1–S4 was driven by the above-cloud air mass source region. The above-cloud air mass sampled near S1 and S4 originated from the boundary layer from the southeast, and the above-cloud air mass sampled near S2 and S3 descended from higher altitudes over the African continent (Fig. 5b and c). The above-cloud OA M_a and N_a for S2 and S3 were over 5 times higher than the corresponding values for S1 and S4 (Table 3). The below-cloud air mass sampled during S1–S4 was advected from the boundary layer from the southeast (Fig. 5a and c). During S1 and S4, the above- and below-cloud rBC and CO concentrations were similar (Table 3) since the above-cloud air mass also originated from the southeast (Fig. 5b and c). During S2 and S3, the continental above-cloud air mass had much higher rBC and CO (over 500 cm⁻³ and 190 ppb) compared to the below-cloud air mass from the southeast (below 150 cm⁻³ and 120 ppb). Since OA, rBC, and CO are indicators of combustion, this suggests the continental above-cloud air mass had greater exposure to biomass burning products compared to the air masses from the southeast. S2 and S3 also had higher below-cloud rBC and CO compared to S1 and S4 (Table 3), which suggests the BBAs with high N_a within 100 m above clouds could be mixing into the cloud layer and polluting the boundary layer. This is also likely to be associated with the history of entrainment mixing of polluted free-tropospheric air into the boundary layer prior to these observations (Diamond et al., 2018).

3.3 Cloud profile classification

Every sawtooth maneuver was preceded by a 5–10 min constant-altitude flight leg about 100 m above the cloud layer to retrieve the above-cloud AOD using 4STAR. Average above-cloud AOD at 550 nm within 50 km of the sampling locations for S1–S4 ranged between 0.33 and 0.49, indicating a BBA layer was located at some altitude above the clouds sampled during S1–S4. During S1, above-cloud N_a < 500 cm⁻³ was sampled up to 200 m above cloud tops (Fig. 3), which indicates the BBA layer was separated from cloud tops. During S4, the level of above-cloud N_a > 500 cm⁻³ was identified over 200 m above cloud tops, indicating a similar separation. Therefore, cloud profiles flown during S1 and S4 were classified as separated profiles. During S2 and S3, the level of above-cloud N_a > 500 cm⁻³ was located within 100 m above cloud tops, and the BBA layer was likely in contact with the cloud tops. Therefore, cloud profiles flown during S2 and S3 were classified as contact profiles. In a previous study, a significantly higher threshold (PCASP N_a = 1000 cm⁻³) was used to identify the BBA layer above stratocumulus clouds off the coast of California (Mardi et al., 2018). The sensitivity of the threshold chosen in this study is examined in Appendix A, and using a threshold of 1000 cm⁻³ would have no significant impact on the results presented in this study.

Figure 6Vertical profiles of (a) N_c, (b) R_e, and (c) LWC and aLWC as a function of Z_N for the four sawtooth maneuvers. Maneuvers with contact (separation) between the biomass burning aerosol layer and cloud tops shown in blue (red).

Download

3.4 Vertical profiles of N_c, R_e, and LWC

Since Z_B and cloud thickness (H) varied between profiles, N_c, R_e, and LWC were examined as a function of normalized height above cloud base (Z_N), where Z_N = $(Z-Z_{\mathrm{B}})/(Z_{\mathrm{T}}-Z_{\mathrm{B}})$ and varied from 0 (cloud base) to 1 (cloud top). Measurements from the four sawtooth maneuvers were compared following McFarquhar et al. (2007) and divided into 10 Z_N bins, where each bin represented 10 % of the cloud layer (Fig. 6). For example, the bin with 0 < Z_N < 0.1 (represented by the midpoint, Z_N = 0.05) included data collected over the bottom 10 % of the cloud layer. For separated profiles, droplet nucleation occurred near cloud base with the median N_c increasing up to Z_N = 0.25 (S1: 132 to 179 cm⁻³; S4: 23 to 85 cm⁻³). The impact of droplet nucleation decreased above cloud base (Z_N = 0.25 to 0.75), and median N_c increased by up to 30 cm⁻³ for S1 and decreased by up to 15 cm⁻³ for S4 (Fig. 6a). Condensational growth occurred over these levels as the median R_e increased with Z_N (Fig. 6b). The median N_c decreased near cloud top (Z_N = 0.75 to 0.95) due to droplet evaporation resulting from cloud-top entrainment mixing between cloudy and non-cloudy air. Contact profiles (S2 and S3) had higher median N_c at cloud base compared to separated profiles, which decreased with height up to Z_N = 0.25 (S2: 190 to 169 cm⁻³, S3: 180 to 131 cm⁻³). The median N_c for S2 and S3 increased by up to 43 cm⁻³ over Z_N = 0.25 to 0.75 and decreased near cloud top due to droplet evaporation. S4 had the lowest N_c at cloud base because the below-cloud M_a and N_a for S4 were over a factor of 3 lower than the corresponding values for S1–S3 (Table 3).

Figure 7Vertical profiles of (a) T and (b) q_T as a function of distance from cloud top. Each line corresponds to an individual ascent through cloud during a sawtooth. The profiles flown during S2 and S3 (S1 and S4) had contact (separation) between the above-cloud biomass burning aerosol layer and cloud tops.

Download

Consistent with condensational growth and collision–coalescence, median R_e increased with Z_N from cloud base to top, from 6.0 to 6.7 µm, 4.6 to 6.9 µm, 4.9 to 8.3 µm, and 8.7 to 9.9 µm for S1–S4, respectively (Fig. 6b). S1 and S4 had higher median R_e at cloud base due to higher drizzle (droplets with diameters larger than 50 µm) concentrations (41 and 31 L⁻¹) compared to S2 and S3 (14 and 18 L⁻¹). For S4, drizzle concentration decreased from Z_N = 0.05 to 0.25, which led to the decrease in median R_e over these heights. The median LWC increased with height up to at least Z_N = 0.75 and decreased near cloud tops due to droplet evaporation (Fig. 6c). The LWC for each sawtooth maneuver was lower than the adiabatic LWC (aLWC) due to cloud-top entrainment mixing, and the ratio of LWC to aLWC was used to quantify the degree of mixing. Lower LWC / aLWC (averaged over the cloud layer) for S2 and S3 (0.37 and 0.41) compared to S1 and S4 (0.51 and 0.55) indicated that contact profiles had greater mixing between cloudy and non-cloudy air in the cloud layer, on average. The boundary layer was capped by an inversion with warmer, drier air above the clouds. During S1–S4, the temperature increased above cloud top by 10.3, 9.3, 8.9, and 1.5 ^∘C, and the total water mixing ratio decreased by 6.2, 5.4, 2.3, and 0.4 g kg⁻¹, respectively (Fig. 7). The decreases in N_c and LWC near stratocumulus tops have been attributed to cloud-top entrainment of the overlying warm and sub-saturated air (Wood, 2012). Droplet evaporation due to the entrainment mixing resulted in decreases of 14, 28, 12, and 26 % in the median N_c near cloud tops during S1–S4, respectively.

3.5 Evidence of the aerosol indirect effect

N_c and R_e were compared between sawtooth maneuvers, and the differences reported hereafter refer to 95 % confidence intervals for the difference in the variable means (based on a two-sample t test, p < 0.02). Between the contact profiles, S2 had significantly higher N_c (differences of 37 to 56 cm⁻³) compared to S3. This was despite having statistically insignificant differences in below-cloud N_a, a greater fractional decrease in median N_c near cloud top compared to S3, and greater entrainment mixing (lower LWC / aLWC). S2 had significantly higher above-cloud N_a compared to S3 and the mixing of above-cloud air with high N_a likely resulted in droplet nucleation above cloud base, where the median N_c for S2 increased from 169 to 220 cm⁻³ over Z_N = 0.25 to 0.75. Between the separated profiles, S1 had significantly higher N_c (differences of 108 to 126 cm⁻³), which could be attributed to significantly higher above-cloud N_a and greater entrainment mixing during S1 compared to S4. However, these differences could also be due to the meteorological differences at their sampling locations (lower boundary layer height, RH, and 500 mb geopotential height for S4 along with a smaller decrease in T and q_T across cloud tops) or the significantly higher below-cloud N_a for S1 compared to S4.

Contact profiles had significantly higher N_c (differences of 45 to 61 cm⁻³) and lower R_e (differences of 1.4 to 2.0 µm) compared to separated profiles. Contact profiles also had significantly higher above-cloud N_a and greater entrainment mixing in the cloud layer (lower LWC / aLWC). These microphysical changes would also impact cloud reflectance (Twomey, 1991) as seen by the significantly higher cloud optical thickness (τ) of contact profiles compared to separated profiles (differences of 2.5 to 8.2). The increase in τ and the cloud reflectance provides observational evidence of the aerosol indirect effect over the Southeast Atlantic due to contact between above-cloud BBAs and the stratocumulus clouds.

However, contact profiles also had significantly higher below-cloud N_a (differences of 145 to 190 cm⁻³), which contribute to the higher N_c relative to separated profiles. Therefore, a statistical analysis was conducted with a larger number of profiles in an attempt to attribute these differences in N_c and R_e to the vertical distance between the above-cloud BBA layer and cloud tops. Building on this case study, 71 cloud profiles flown on six flights between 6 and 25 September 2016 were examined, and the impact of above-cloud BBAs on the free-tropospheric humidity and buoyancy across cloud tops was explored. Sixty-one contact and separated profiles were further classified as low-N_a or high-N_a profiles based on the below-cloud N_a. This was done to quantify the differences in N_c and R_e between contact and separated profiles within boundary layers with similar below-cloud N_a.

Figure 8Flight tracks from PRFs 5, 7, 8, 9, 11, and 12 flown on 6, 10, 12, 14, 20, and 25 September 2016 with green segments indicating location of cloud profiles (flight tracks from PRFs 7 and 8 coincide with PRF13 and hence are not visible).

Download

Table 4The range of time, latitude, longitude, above-cloud AOD, and cloud-top height (Z_T) for cloud profiles flown during the six flights. The lowest altitude where above-cloud N_a > 500 cm⁻³ occurred during the flight (Z₅₀₀) is in the far-right column.

Download Print Version | Download XLSX

4.1 Meteorological conditions and above-cloud aerosols

Six flights (including PRF5) are included in the statistical analysis. On 10, 12, and 25 September, the P-3 took off from Walvis Bay, Namibia (23^∘ S, 14.6^∘ E), and flew northwest from 23^∘ S, 13.5^∘ E toward 10^∘ S, 0^∘ E, returning along the same track (Fig. 8). Different tracks were followed on 6, 14, and 20 September, which included meridional legs along 9, 7.5 and 9^∘ E, and 9 and 10.5^∘ E, respectively. Meteorological conditions on 10, 12, and 14 September were similar to the conditions described for the case study. South-southeasterly surface winds were associated with a surface low-pressure system over Africa. The surface wind speeds varied between 5 and 10 m s⁻¹ depending on the pressure gradient between the continental low and a surface high toward the southwest. A region of 925 mb RH < 60 % persisted along the coast due to dry-air advection from Africa. A different meteorological setup on 20 September had westerly surface winds and easterly winds at 925 mb. The aerosol plume was sampled immediately above the boundary layer (600 m) as warm surface air was overlaid by drier, polluted air from the continent. The continental surface low was located farther south on 25 September compared to other flight days with the region of low 925 mb RH to the south of the flight track. The study region had RH > 60 % with south-southeasterly surface winds and southerly 925 mb winds. The BBA layer with above-cloud N_a > 500 cm⁻³ was sampled during each flight with variability in its vertical location (Table 4). Only separated profiles were flown on 10 and 14 September (Table 2), when the BBA layer and cloud tops were separated by over 600 and 1500 m, respectively (Table 4). On 12 September, profile 1 (P1) had N_a > 500 cm⁻³ within 75 above cloud tops and was classified as a contact profile, while P2 and S1 were classified as separated profiles. On 20 September, each profile had above-cloud AOD > 0.4 and was classified as a contact profile. On 25 September, the profiles had above-cloud AOD > 0.27, and each profile (except from a sawtooth near 11^∘ S, 1^∘ E) was classified as a contact profile.

Figure 9Cloud base and top heights for contact (blue) and separated (red) profiles flown during the six PRFs.

Download

Figure 10Kernel density estimates (indicated by the width of shaded area) and boxplots showing the 25th (Q1), 50th (white point), and 75th (Q3) percentile for (a) N_c, (b) R_e, and (c) LWC as a function of Z_N for contact (blue) and separated (red) profiles.

Download

4.2 N_c, R_e, and LWC for contact and separated profiles

Since clouds sampled on different flight days had variable Z_B and Z_T (Fig. 9), vertical profiles of N_c, R_e, and LWC from the contact and separated profiles were compared as a function of Z_N. The frequency distributions of N_c, R_e, and LWC as a function of Z_N are examined in Fig. 10 using violin plots (Hintze and Nelson, 1998; Wang et al., 2020), where the width of the shaded area represents the proportion of data there. The average N_c for contact profiles was significantly higher than the average N_c for separated profiles (differences of 60 to 68 cm⁻³). During separated profiles, the median N_c had little variability up to Z_N = 0.75 (114 to 122 cm⁻³) and decreased thereafter with Z_N to 73 cm⁻³ due to droplet evaporation (Fig. 10a). During contact profiles, the median N_c decreased slightly up to Z_N = 0.25 (183 to 174 cm⁻³), increased to 214 cm⁻³ at Z_N = 0.75, and decreased near cloud top to 157 cm⁻³ due to droplet evaporation. Contact profiles had significantly lower R_e than the separated profiles (differences of 1.1 to 1.3 µm), and the median R_e increased with Z_N from 4.9 to 7.0 µm for contact and from 6.6 to 8.6 µm for separated profiles (Fig. 10b). The differences in R_e were likely due to the significantly lower drizzle concentrations for contact profiles (differences of 5 to 20 L⁻¹).

The average LWC for contact and separated profiles were within 0.01 g m⁻³, and the median LWC increased with Z_N to 0.23 g m⁻³ at Z_N = 0.85 for contact and 0.21 g m⁻³ at Z_N = 0.75 for separated profiles (Fig. 10c). Contact profiles had lower LWC / aLWC in the cloud layer (0.45) compared to separated profiles (0.57), which suggests there was greater entrainment mixing during contact profiles, on average. However, droplet evaporation near cloud top had a stronger impact on separated profiles as the median LWC decreased to 0.16 g m⁻³ for separated and 0.20 g m⁻³ for contact profiles (Fig. 10c). Separated profiles had a greater decrease in LWC / aLWC near cloud top (0.41 to 0.26) compared to contact profiles (0.38 to 0.30) and greater fractional decreases in median N_c and LWC (40 and 16 %) compared to contact profiles (25 and 9 %). The stronger impact of droplet evaporation during separated profiles contributed to the differences between N_c for contact and separated profiles.

4.3 Cloud-top evaporative cooling

Buoyancy and humidity across cloud tops were determined to explore the cloud-top entrainment mechanisms resulting in the differential impact of droplet evaporation for these profiles. Cloud-top instability is the dominant source of turbulence in stratocumulus, with evaporative cooling being a key driver of instability (Mellado, 2017). Recent studies have shown there is strong correlation between above-cloud AOD and water vapor within air masses originating from the African continent (Deaconu et al., 2019; Pistone et al., 2021). Longwave cooling by water vapor within the BBA layer leads to decreased cloud-top cooling, and cloud-top dynamics are influenced by distinct radiative contributions from water vapor and absorbing aerosols. Evaporative cooling in a mixture of dry and cloudy air near cloud top generates negatively buoyant air mixtures, which further enhances mixing and leads to an entrainment feedback called cloud top entrainment instability, or CTEI (Kuo and Schubert, 1988). Under such conditions, negative buoyancy leads to an unstable feedback, unlike the conventional association of negative buoyancy with atmospheric stability. The critical condition for cloud-top stability is given by Kuo and Schubert (1988) as

where k is the CTEI parameter, θ_e is the equivalent potential temperature, L_v is the latent heat of vaporization, and C_p is the specific heat capacity of air at constant pressure. The Δ operator represents gradients across the cloud top, defined here as the difference between θ_e (or q_T) measured 100 m above cloud top and the vertical average of θ_e (or q_T) over the top 100 m of the profile. Following Eq. (13) from Kuo and Schubert (1988), k > 0.23 indicates negative buoyancy across cloud tops. Water vapor mixing ratio measured by the chilled-mirror hygrometer was used to calculate θ_e and q_T. Since lower Δq_T was sampled during descents into cloud due to condensation on the hygrometer, k values for descents were determined to be measurement artifacts and not usable here.

Figure 11Difference between equivalent potential temperature (θ_e) and total water mixing ratio (q_T) measured within cloud and 100 m above cloud top for contact (blue) and separated (red) profiles (only ascents through cloud shown).

Download

All separated profiles (except PRF5 S1-3 and S4-1, S4-3, and S4-5) laid within the region of cloud-top instability (k > 0.23) on a Δθ_e–Δq_T plane (Fig. 11) and showed negative buoyancy across cloud tops. During PRF5 S1-3, low Δθ_e was sampled due to higher above-cloud humidity associated with the presence of N_a > 100 cm⁻³ within 50 m above cloud tops. During PRF5 S4, a weak cloud-top inversion led to positive Δθ_e and Δq_T < −2 g kg⁻¹ (Fig. 7). For the remaining separated profiles, negative buoyancy across cloud tops led to forced descent of dry free-tropospheric air into the clouds. Since the free-tropospheric air was warmer and drier than the cloudy air, droplet evaporation led to the decreases in median N_c and LWC near cloud top. The positive evaporative cooling feedback and greater Δq_T compared to contact profiles (Fig. 11) explain the stronger impact of droplet evaporation on median N_c and LWC for separated profiles. While evaporative cooling triggered the CTEI feedback, the clouds persisted, consistent with cloud-top radiative cooling or surface evaporation leading to boundary layer moistening (Lock, 2009; Mellado, 2017).

All contact profiles (except PRF13 S1-3) laid within the region of cloud-top stability and showed positive buoyancy across cloud tops. Entrainment mixing for these profiles likely occurred when the clouds penetrated the inversion. This is consistent with significantly higher average H (267 m) for contact profiles compared to separated profiles (213 m). Braun et al. (2018) found a negative correlation between H and adiabaticity (ratio of the measured and the adiabatic liquid water path), which is consistent with contact profiles having lower LWC / aLWC and higher H compared to separated profiles. In the presence of above-cloud BBAs, the above-cloud air was more humid, and the above-cloud N_a was significantly higher compared to separated profiles (differences of 768 to 831 cm⁻³). Contact profiles had greater entrainment mixing compared to separated profiles, and the median N_c increased with height over Z_N = 0.25 to 0.75. It is likely the entrainment of BBAs into clouds resulted in additional droplet nucleation over these Z_N levels. Therefore, weaker droplet evaporation near cloud top and additional droplet nucleation above cloud base in the presence of above-cloud BBAs likely contributed to the differences between N_c for contact and separated profiles.

Table 5Aerosol and cloud properties were averaged across all contact and separated profiles flown in low-N_a and-high N_a boundary layers. These averages were compared between contact and separated profiles. The values listed below represent the 95 % confidence intervals (from a two-sample t test) when the differences were statistically significant. Positive values indicate the average for contact profiles was higher, and “insignificant” denotes the differences were statistically insignificant.

Download Print Version | Download XLSX

4.4 N_c, R_e, and LWC in boundary layers with similar N_a

Contact profiles had significantly higher below-cloud N_a (differences of 93 to 115 cm⁻³) and below-cloud CO (differences of 13 to 16 ppb) in addition to higher above-cloud N_a (differences of 768 to 831 cm⁻³) compared to separated profiles. Enhanced aerosol loading within the boundary layer is consistent with BBAs immediately above cloud tops entraining into the cloud layer and polluting the boundary layer. This is consistent with higher above-cloud CO (240 ppb) sampled for contact profiles with below-cloud CO > 100 ppb compared to above-cloud CO (104 ppb) for profiles with below-cloud CO < 100 ppb. The correlations between above- and below-cloud aerosols could be partly due to the history of entrainment mixing between free-tropospheric and boundary layer air masses (Diamond et al., 2018). To investigate the contribution of below-cloud N_a relative to the impact of above-cloud BBAs on cloud properties, 28 contact and 33 separated profiles were classified into four new regimes defined as follows: contact high N_a (C-H), separated high N_a (S-H), contact low N_a (C-L), and separated low N_a (S-L), where high- and low-N_a boundary layers were separated using a threshold concentration of 350 cm⁻³. Cloud microphysical properties and above- and below-cloud N_a were compared between 20 C-H and 11 S-H profiles and between 8 C-L and 22 S-L profiles (Table 5) to compare contact and separated profiles with minor differences in below-cloud N_a.

Within low-N_a boundary layers, C-L and S-L profiles had insignificant differences in below-cloud N_a despite significantly higher above-cloud N_a for C-L profiles (differences of 592 to 669 cm⁻³), higher N_c (differences of 22.8 to 34.9 cm⁻³), and lower R_e (differences of 0.5 to 1.0 µm) compared to S-L profiles. Within high-N_a boundary layers, C-H profiles had significantly higher below-cloud N_a compared to S-H profiles (differences of 39.1 to 70.5 cm⁻³), but the differences were much smaller than those in the above-cloud N_a (differences of 738 to 884 cm⁻³). Further, the C-H profiles had significantly higher N_c (differences of 75.5 to 88.5 cm⁻³) and lower R_e (differences of 1.1 to 1.3 µm) than the S-H profiles. Previous studies have argued the changes in N_c due to the impact of BBAs are more strongly correlated with below-cloud N_a compared to above-cloud N_a (Diamond et al., 2018; Mardi et al., 2019). However, these results suggest that, although the differences in N_c were lower than the differences in above-cloud N_a, significant changes in N_c and R_e were associated with contact with above-cloud BBAs, and these changes were independent of the below-cloud aerosol loading.

Figure 12Boxplots representing vertical profiles of (a) N_c, (b) R_e, and (c) LWC as a function of Z_N for contact (blue) and separated (red) profiles within boundary layers with high N_a (> 350 cm⁻³) (darker) or low N_a (< 350 cm⁻³) (lighter). The number of 1 Hz measurements within each regime is listed within parentheses.

Download

Vertical profiles of N_c, R_e, and LWC are examined (Fig. 12) to further investigate the microphysical changes due to contact with above-cloud BBAs. Within low-N_a boundary layers, there were minor deviations in N_c with Z_N up to Z_N = 0.75 (Fig. 12a). Over the top 20 % of the cloud layer, S-L profiles had a decrease in median N_c (32 cm⁻³), with a smaller change for C-L profiles (8 cm⁻³) over the same levels. There was also a weaker decrease in water vapor mixing ratio across cloud tops for contact profiles. Thus, cloud-top entrainment of more humid air likely occurred for the C-L profiles. This is consistent with higher median R_e and LWC over Z_N = 0.75 to 0.95 for C-L profiles compared to S-L profiles despite having lower R_e and LWC closer to cloud base (Fig. 12b and c). Thus, the microphysical differences between contact and separated profiles within low-N_a boundary layers (where most separated profiles were sampled) are consistent with the processes of cloud-top entrainment and droplet evaporation.

The differences between below-cloud N_a for C-H profiles and that for S-H profiles (39.1 to 70.5 cm⁻³) were lower than the corresponding differences in N_c (75.5 to 88.4 cm⁻³). C-H profiles had significantly higher N_c and lower R_e compared to S-H profiles throughout the cloud layer (Fig. 12a and b). There was a significant increase in median N_c for C-H profiles over Z_N = 0.25 to 0.75, which was accompanied by higher median LWC for C-H profiles in the top half of the cloud layer. This is consistent with additional droplet nucleation above cloud base during C-H profiles. Additionally, there was a stronger decrease in N_c near cloud top for S-H profiles (N_c decreased by 66 cm⁻³) compared to C-H profiles (N_c decreased by 29 cm⁻³) likely due to cloud-top entrainment. It is difficult to separate the impact of changes in droplet nucleation on differences in N_c between C-H and S-H profiles from the impact of changes in droplet evaporation due to cloud-top entrainment. Therefore, it is speculated the microphysical changes within high-N_a boundary layers were likely driven by the combination of higher below-cloud N_a, potential droplet nucleation above cloud base, and weaker droplet evaporation near cloud tops in the presence of above-cloud BBAs. The sensitivity of these results to using different thresholds to locate BBAs (other than 500 cm⁻³), to define “separation” between the aerosol and cloud layers (other than 100 m), and to define a “high-N_a boundary layer” (other than 350 cm⁻³) is discussed in Appendix A but does not affect the qualitative findings.