A GCSS Boundary-Layer Cloud Model Intercomparison Study Of The First Astex Lagrangian Experiment

Three single-column models (all with an explicit liquid water budget and compara-tively high vertical resolution) and three two-dimensional eddy-resolving models (including one with bin-resolved microphysics) are compared with observations from the first ASTEX Lagrangian experiment. This intercomparison was a part of the second GCSS boundary-layer cloud modelling workshop in August 1995.In the air column tracked during the first ASTEX Lagrangian experiment, a shallow subtropical drizzling stratocumulus-capped marine boundary layer deepens after two days into a cumulus capped boundary layer with patchy stratocumulus. The models are forced with time varying boundary conditions at the sea-surface and the capping inversion to simulate the changing environment of the air column.The models all predict the observed deepening and decoupling of the boundary layer quite well, with cumulus cloud evolution and thinning of the overlying stratocumulus. Thus these models all appear capable of predicting transitions between cloud and boundary-layer types with some skill. The models also produce realistic drizzle rates, but there are substantial quantitative differences in the cloud cover and liquid water path between models. The differences between the eddy-resolving model results are nearly as large as between the single column model results. The eddy resolving models give a more detailed picture of the boundary-layer evolution than the single-column models, but are still sensitive to the choice of microphysical and radiative parameterizations, sub-grid-scale turbulence models, and probably model resolution and dimensionality. One important example of the differences seen in these parameterizations is the absorption of solar radiation in a specified cloud layer, which varied by a factor of four between the model radiation parameterizations.     

Evaluation of the Diurnal Cycle in the Atmospheric Boundary Layer Over Land as Represented by a Variety of Single-Column Models: The Second GABLS Experiment

We present the main results from the second model intercomparison within the GEWEX (Global Energy and Water cycle EXperiment) Atmospheric Boundary Layer Study (GABLS). The target is to examine the diurnal cycle over land in today??s numerical weather prediction and climate models for operational and research purposes. The set-up of the case is based on observations taken during the Cooperative Atmosphere-Surface Exchange Study-1999 (CASES-99), which was held in Kansas, USA in the early autumn with a strong diurnal cycle with no clouds present. The models are forced with a constant geostrophic wind, prescribed surface temperature and large-scale divergence. Results from 30 different model simulations and one large-eddy simulation (LES) are analyzed and compared with observations. Even though the surface temperature is prescribed, the models give variable near-surface air temperatures. This, in turn, gives rise to differences in low-level stability affecting the turbulence and the turbulent heat fluxes. The increase in modelled upward sensible heat flux during the morning transition is typically too weak and the growth of the convective boundary layer before noon is too slow. This is related to weak modelled near-surface winds during the morning hours. The agreement between the models, the LES and observations is the best during the late afternoon. From this intercomparison study, we find that modelling the diurnal cycle is still a big challenge. For the convective part of the diurnal cycle, some of the first-order schemes perform somewhat better while the turbulent kinetic energy (TKE) schemes tend to be slightly better during nighttime conditions. Finer vertical resolution tends to improve results to some extent, but is certainly not the solution to all the deficiencies identified.     

The First Observed Cloud Echoes and Microphysical Parameter Retrievals by China's 94-GHz Cloud Radar

By using the cloud echoes first successfully observed by China's indigenous 94-GHz SKY cloud radar, the macrostructure and microphysical properties of drizzling stratocumulus clouds in Anhui Province on 8 June2013 are analyzed, and the detection capability of this cloud radar is discussed. The results are as follows.(1) The cloud radar is able to observe the time-varying macroscopic and microphysical parameters of clouds,and it can reveal the microscopic structure and small-scale changes of clouds.(2) The velocity spectral width of cloud droplets is small, but the spectral width of the cloud containing both cloud droplets and drizzle is large. When the spectral width is more than 0.4 m s-1, the radar reflectivity factor is larger(over –10 dBZ).(3) The radar's sensitivity is comparatively higher because the minimum radar reflectivity factor is about–35 dBZ in this experiment, which exceeds the threshold for detecting the linear depolarized ratio(LDR) of stratocumulus(commonly –11 to –14 dBZ; decreases with increasing turbulence).(4) After distinguishing of cloud droplets from drizzle, cloud liquid water content and particle effective radius are retrieved. The liquid water content of drizzle is lower than that of cloud droplets at the same radar reflectivity factor.     

Large-Eddy Simulation Of A Nocturnal Stratocumulus-Topped Marine Atmospheric Boundary Layer: An Uncertainty Analysis

A large-eddy simulation (LES) model has been used to study a nocturnalstratocumulus-topped marine atmospheric boundary layer. The main objectivesof our study have been first to investigate the statistical significance of LES-derived data products. Second, to test the sensitivity of our LES results with respect to the representation of subgrid-scale mixing and microphysical processes, and third to evaluate and to quantify the parametric uncertainty arising from the incomplete knowledge of the environmental parameters that are required to specify the initial and boundary conditions of a particular case study. Model simulations were compared with observations obtained in solid stratocumulus during the third flight of the first 'Lagrangian' experiment of the Atlantic Stratocumulus Transition Experiment (ASTEX). Based on these simulations the following conclusions could be drawn. Resolution(50 × 50 × 25 m³) and domain size (3.2 × 3.2 × 1.5 km³) of the LES calculations were adequate from a numerical point of view to represent the essential features of the stratocumulus-topped boundary layer. However, the ensemble runs performed in our study to investigate the statistical significance of LES-derived data products demonstrate that the area-time averaging procedure for the second-order moments produces only a low degree of statistical reliability in the model results. This illustratesthe necessity of having LES model results that are not only of adequate resolution but also of sufficiently large domain. The impact of different subgrid schemes was small, but the primary effects of drizzle were found to influence the boundary-layer structure in a climatologically significant way. The parametric uncertainty analysis revealed that the largest contribution to the variance of the LES-derived data products is due to theuncertainties in the cloud-top jump of total water mixing ratio and the net radiative forcing. The differences between the model and measurements for most of the simulated quantities were within the modelling uncertainties, but the calculated precipitation rate was found to differ significantly from that derived in the observations.     

层积云覆盖的边界层的数值模拟研究(Ⅱ): 数值模拟结果分析

着重分析了文献[1]新建立的层积云覆盖的边界层的数值模式对大涡模拟对比试验中个例模拟的结果。并将数值模拟结果与一些已有的观测和数值模拟结果进行了对比。分析结果表明,该边界层模式能较合理地模拟海区层积云覆盖的边界层的微湍流结构和云的微物理过程。     

Parametric representation of heat and moisture fluxes in cloud-topped mixed layers

Aircraft measurements were made from the NCAR Electra in stratus and stratocumulus clouds off the coast of California in June 1976. Several types of cloud conditions were observed, including (1) a broken layer less than 100 m thick, capped by an inversion at ~1000 m, (2) a broken stratocumulus layer ~300 m thick with an inversion at ~500 m, and (3) a solid stratocumulus layer ~250 m thick with an inversion at ~500 m. Although these observations indicate that a variety of cloud conditions may exist in mixed layers, simple one-dimensional mixed-layer models implicitly assume a solid cloud layer with no unsaturated region within the cloud. In order to generalize these simple models, a parametric representation of the heat and moisture fluxes is considered. In this scheme, the fluxes are parameterized in terms of the product of a cloud mass flux and the characteristic difference between the thermodynamic properties of an updraft-downdraft circulation. This representation allows for an explicit representation of the buoyancy flux when the downdraft has no liquid water.Data collected during these flights were used to calculate heat and moisture fluxes and to obtain the mean difference in the thermodynamic properties of the updrafts and downdrafts at a given level. The mass flux was calculated using updraft-downdraft differences and the fluxes. The mass fluxes obtained using various thermodynamic quantities are examined for consistency. The vertical distribution of the mass flux is determined. These results indicate that a mass flux formulation could prove to be useful in modeling applications where cloud conditions may vary between solid and broken.     

层积云覆盖的海洋边界层云详细微物理过程的数值模拟

文中建立了一个含显式分档的云微物理模式和辐射传输模式的一维 3阶湍流闭合模式 ,该模式可用于研究海洋边界层云中气溶胶和云的相互作用过程 ,同时提出了一种新的动力模式和微物理模式耦合方法 ,该方法可使动力模式中液态水相关项可以直接由微物理模式变量计算得到。作为模式的初步应用模拟了 2 0 0 1年APEX/ACE Asia在西太平洋上所观测到的一个个例。模拟结果和观测资料比较表明该模式基本上模拟出层积云覆盖的海洋边界层的基本结构     

The Sensitivity Of A Stratocumulus Transition: Model Simulations Of The Astex First Lagrangian

In this study, a one-dimensional ensemble-average model is used to simulatethe Atlantic Stratocumulus Transition Experiment firstLagrangian, where the same airmass was followed from the subtropical high pressure region en route towards the trade wind region. The airmass experiences increasing sea-surface temperature and achange from subsidence to weak ascent on its way south. Thiscauses the marine boundary layer (MBL)to grow and the cloud deck to change from a solid stratocumulus deck tomore broken stratocumulus clouds with cumulus cloudsdeveloping beneath, and reaching up into the stratocumulus clouds.A control run is analyzed and compared in detail with theobservations. Both a statistical evaluation and a more subjective evaluation are performed, where both establish confidencein the model performance. The model captures the MBL growth and the cloudliquid water, as well as the drizzle flux, is well predicted by the model.A sensitivity study was performed with the objective of examining theMBL and the cloud response to external and internal 'forces'.The results show that, if drizzle formation is not allowed,unrealistically high cloud liquid water mixing ratios are predicted. Even though the drizzle flux is very small, it is still important for the water budget of the MBL and for the boundary-layer dynamics.We also found that the sea-surface temperature increase is more important for the increasing cloud top height than the synoptic-scale divergence fields. However, the synoptic-scale subsidence is crucial during the first day, when the sea-surface temperature was constant, in keepingthe cloud top at a constant height. Drizzle evaporation below the cloud base seems to be important for below-cloud condensation. The drizzle predictions are significantly altered when the prescribed cloud droplet and/or drizzle drop numbers are altered.     

Observation of a Self-Limiting,Shear-Induced Turbulent Inversion Layer Above Marine Stratocumulus

High-resolution measurements of thermodynamic, microphysical, and turbulence properties inside a turbulent inversion layer above a marine stratocumulus cloud layer are presented. The measurements are performed with the helicopter-towed measurement payload Airborne Cloud Turbulence Observation System (ACTOS), which allows for sampling with low true air speeds and steep profiles through cloud top. Vertical profiles show that the turbulent inversion layer consists of clear air above the cloud top, with nearly linear profiles of potential temperature, horizontal wind speed, absolute humidity, and concentration of interstitial aerosol. The layer is turbulent, with an energy dissipation rate nearly the same as that in the lower cloud, suggesting that the two are actively coupled, but with significant anisotropic turbulence at the large scales within the turbulent inversion layer. The turbulent inversion layer is traversed six times and the layer thickness is observed to vary between 37 and 85 m, whereas the potential temperature and horizontal wind speed differences at the top and bottom of the layer remain essentially constant. The Richardson number therefore increases with increasing layer thickness, from approximately 0.2 to 0.7, suggesting that the layer develops to the point where shear production of turbulence is sufficiently weak to be balanced by buoyancy suppression. This picture is consistent with prior numerical simulations of the evolution of turbulence in localized stratified shear layers. It is observed that the large eddy scale is suppressed by buoyancy and is on the order of the Ozmidov scale, much less than the thickness of the turbulent inversion layer, such that direct mixing between the cloud top and the free troposphere is inhibited, and the entrainment velocity tends to decrease with increasing turbulent inversion-layer thickness. Qualitatively, the turbulent inversion layer likely grows through nibbling rather than engulfment.     

Coamps<Superscript>®</Superscript>-Les: Model Evaluation and Analysis of Second-and Third-Moment Vertical Velocity Budgets

The Naval Research Laboratory Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS^®) has been extended to perform as a large-eddy simulation (LES) model. It has been validated with a series of boundary-layer experiments spanning a range of cloud nighttime, and includes a nighttime stratocumulus case, a trade wind cumulus layer, shallow cumulus convection over land, and a mixed regime consisting of cumulus clouds under broken stratocumulus. COAMPS-LES results are in good agreement with other models for all the cases simulated. Exact numerical budgets for the vertical velocity second\((\overline{w^{'2}})\) and third moment\((\overline{w^{'3}})\) have been derived for the stratocumulus and trade wind cumulus cases. For the\(\overline{w^{'3}}\) budget in the stratocumulus, the buoyancy contribution from the updraughts and downdraughts largely cancel each other due to their similar magnitudes but opposite signs. In contrast, for the cumulus layer, the negative buoyancy contribution from the environmental downdraughts is negligible and the positive contribution from the updraughts completely dominates due to the conditional instability in the environment. As a result,\(\overline{w^{'3}}\) is significantly larger in the cumulus than in the stratocumulus layer.     

An Intercomparison of Large-Eddy Simulations of the Stable Boundary Layer

Results are presented from the first intercomparison of large-eddy simulation (LES) models for the stable boundary layer (SBL), as part of the Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study initiative. A moderately stable case is used, based on Arctic observations. All models produce successful simulations, in as much as they generate resolved turbulence and reflect many of the results from local scaling theory and observations. Simulations performed at 1-m and 2-m resolution show only small changes in the mean profiles compared to coarser resolutions. Also, sensitivity to subgrid models for individual models highlights their importance in SBL simulation at moderate resolution (6.25 m). Stability functions are derived from the LES using typical mixing lengths used in numerical weather prediction (NWP) and climate models. The functions have smaller values than those used in NWP. There is also support for the use of K-profile similarity in parametrizations. Thus, the results provide improved understanding and motivate future developments of the parametrization of the SBL.     

A note on the top-down and bottom-up gradient functions over a forested site

The dimensionless bottom-up and top-down gradient functions in the convective boundary layer (CBL) are evaluated utilizing long-term well-calibrated carbon dioxide mixing ratio and flux measurements from multiple levels of a 447-m tall tower over a forested area in northern Wisconsin, USA. The estimated bottom-up and top-down functions are qualitatively consistent with those from large-eddy simulation (LES) results and theoretical expectations. Newly fitted gradient functions are proposed based on observations for this forested site. The integrated bottom-up function over the lowest 4% of the CBL depth estimated from the tower data is about five times larger than that from LES results for a ‘with-canopy’ case, and is smaller than that from LES results for a ‘no-canopy’ case by a factor of 0.7. We discuss the uncertainty in the evaluated gradient functions due to stability, wind direction, and uncertainty in the entrainment flux and show that while all of these have a significant impact on the gradient functions, none can explain the differences between the modelled and observed functions. The effects of canopy features and atmospheric stability may need to be considered in the gradient function relations.     

Current problems in the stratocumulus-topped atmospheric boundary layer

Extended sheets of stratocumulus (Sc) in the upper part of the atmospheric boundary layer (ABL) often occur under appropriate meteorological conditions. These cloud decks are important both in climate studies and in weather forecasting. We review the current knowledge of the turbulent structure of the ABL capped by a cloud deck, in the light of recent observations and model studies. The most important physical processes determining this structure are longwave radiative cooling at cloud top, shortwave radiative wanning by absorption in the cloud, surface buoyancy flux, and wind shear in the ABL. As a result, turbulence can cause entrainment against the buoyancy jump at cloud top. In cases where only longwave radiative fluxes and surface buoyancy fluxes are important, the turbulent structure is relatively well understood. When shortwave radiative fluxes and/or wind shear are also important, the resulting turbulent structure may change considerably. A decoupling of the cloud from the sub-cloud layer or of the top of the cloud from the rest of the ABL is then regularly observed. In no cases are the details of the entrainment at cloud top understood well enough to derive a relatively simple formulation that is consistent with observations. Cloud-top entrainment instability may lead to the break-up of a cloud deck (but also to cloud deepening). The role of mesoscale circulations in determining fractional cloudiness is not yet well understood.     

A first order mixing–condensation scheme for nocturnal stratocumulus

A simple closure scheme for nocturnal stratocumulus is proposed. The scheme is formulated in conserved variables. Cloud fraction and cloud water amount are diagnosed assuming a top-hat distribution for total water. Conversion of cloud water into rain water is parameterized in terms of cloud water and the incoming rain flux. Turbulence transport in the cloud layer is accounted for by a first-order vertical diffusion scheme with a profile-type diffusivity. The length scale corresponds to the thickness of the cloud layer. The turbulent velocity scale is directly related to the long wave radiative flux divergence in the cloud. Entrainment at cloud top is implicitly treated by extending the in-cloud mixing profile slightly beyond cloud top. The excess height is derived from the buoyancy frequency at cloud top and a radiative–convective velocity scale. The scheme is capable of simulating realistic profiles of the conserved variables and cloud parameters for a case of nocturnal stratocumulus prepared on the basis of ASTEX data.     

Study of near-surface models for large-eddy simulations of a neutrally stratified atmospheric boundary layer

In large-eddy simulations (LES) of the atmospheric boundary layer (ABL), near-surface models are often used to supplement subgrid-scale (SGS) turbulent stresses when a major fraction of the energetic scales within the surface layer cannot be resolved with the temporal and spatial resolution at hand. In this study, we investigate the performance of both dynamic and non-dynamic eddy viscosity models coupled with near-surface models in simulations of a neutrally stratified ABL. Two near-surface models that are commonly used in LES of the atmospheric boundary layer are considered. Additionally, a hybrid Reynolds- averaged/LES eddy viscosity model is presented, which uses Prandtl’s mixing length model in the vicinity of the surface, and blends in with the dynamic Smagorinsky model away from the surface. Present simulations show that significant portions of the modelled turbulent stresses are generated by the near-surface models, and they play a dominant role in capturing the expected logarithmic wind profile. Visualizations of the instantaneous vorticity field reveal that flow structures in the vicinity of the surface depend on the choice of the near-surface model. Among the three near-surface models studied, the hybrid eddy viscosity model gives the closest agreement with the logarithmic wind profile in the surface layer. It is also observed that high levels of resolved turbulence stresses can be maintained with the so-called canopy stress model while producing good agreement with the logarithmic wind profile.     

A Mixed-Layer Model Of The Well-Mixed Stratocumulus-Topped Boundary Layer

Recently a range of sophisticated large-eddy simulations of thecloud-topped boundary layer have been intercompared and furthercompared with observations and single column models. Here we comparethese results with perhaps the simplest model of the cloud-toppedboundary layer, namely a mixed-layer model. Results from the model aredescribed with two aims in mind. Firstly, the good results act as areminder of the success of simple models, and, secondly, we suggestthat a simple mixed-layer model could be used as a baseline for futuremodel intercomparisons.The mixed-layer model is based on two assumptions that follow previousstudies. Firstly, the liquid-water potential temperature and the total waterspecific humidity are assumed to be constant with height in the boundarylayer. Secondly, turbulence entrains air across the inversion into the boundarylayer at a rate that is assumed to be proportional to the jump in radiative flux at the cloud top and inversely proportional to the jump in buoyancy at the inversion. The constant of proportionality is called the entrainment efficiency.Results from the model for the entrainment rate and height evolutionof the boundary layer are compared with the observations and modelsconsidered in a EUCREM intercomparison study. Thepresent mixed-layer model accurately predicts the observed heightevolution of the boundary layer, but over-estimates the entrainmentrate to a similar degree as the large-eddy simulations. We show that,if the subsidence rate is reduced to the value given by observationsrather than the value used in the EUCREM intercomparison study,then the model agrees well with observed value of the entrainment rateif the entrainment efficiency is taken to be 0.6. With this value, themodel also agrees well with a further case study byBechtold et al. An entrainment efficiency of 0.6 is a little higherthan suggested by large eddy simulations, but such simulations do notcurrently resolve the entrainment events explicitly. Hence this pointdeserves further study.     

Stratocumulus-capped mixed layers derived from a three-dimensional model

Results of a three-dimensional numerical model are analysed in a study of turbulence and entrainment within mixed layers containing stratocumulus with or without parameterized cloud-top radiative cooling. The model eliminates most of the assumptions invoked in theories of cloud-capped mixed layers, but suffers disadvantages which include poor resolution and large truncation errors in and above the capping inversion.For relatively thick mixed layers with relatively thick capping inversions, the cloud-top radiative cooling is found to be lodged mostly within the capping inversion when the cooling is confined locally to the upper 50 m or less of the cloud. It does not then contribute substantially towards increased buoyancy flux and turbulence within the well mixed layer just below.The optimal means of correlating the entrainment rate, or mixed-layer growth rate, for mixed layers of variable amounts of stratocumulus is found to be through functional dependence upon an overall jump Richardson number, utilizing as scaling velocity the standard deviation of vertical velocity existing at the top of the mixed layer (near the center of the capping inversion). This velocity is found to be a fraction of the generalized convective velocity for the mixed layer as a whole which is greater for cloud-capped mixed layers than for clear mixed layers.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

A multi-limit formulation for the equilibrium depth of a stably stratified boundary layer

Currently no expression for the equilibrium depth of the turbulent stably-stratified boundary layer is available that accounts for the combined effects of rotation, surface buoyancy flux and static stability in the free flow. Various expressions proposed to date are reviewed in the light of what is meant by the stable boundary layer. Two major definitions are thoroughly discussed. The first emphasises turbulence and specifies the boundary layer as a continuously and vigorously turbulent layer adjacent to the surface. The second specifies the boundary layer in terms of the mean velocity profile, e.g. by the proximity of the actual velocity to the geostrophic velocity. It is shown that the expressions based on the second definition are relevant to the Ekman layer and portray the depth of the turbulence in the intermediate regimes, when the effects of static stability and rotation essentially interfere. Limiting asymptotic regimes dominated by either stratification or rotation are examined using the energy considerations. As a result, a simple equation for the depth of the equilibrium stable boundary layer is developed. It is valid throughout the range of stability conditions and remains in force in the limits of a perfectly neutral layer subjected to rotation and a rotation-free boundary layer dominated by surface buoyancy flux or stable density stratification at its outer edge. Dimensionless coefficients are estimated using data from observations and large-eddy simulations. Well-known and widely used formulae proposed earlier by Zilitinkevich and by Pollard, Rhines and Thompson are shown to be characteristic of the above interference regimes, when the effects of rotation and static stability (due to either surface buoyancy flux, or stratification at the outer edge of the boundary layer) are roughly equally important.     

Aerosol-Cloud-Precipitation Interactions in a Closed-cell and Non-homogenous MBL Stratocumulus Cloud

A closed-cell marine stratocumulus case during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) aircraft field campaign is selected to examine the heterogeneities of cloud and drizzle microphysical properties and the aerosol-cloud-precipitation interactions. The spatial and vertical variabilities of cloud and drizzle microphysics are found in two different sets of flight legs: Leg-1 and Leg-2, which are parallel and perpendicular to the cloud propagation, respectively. The cloud along Leg-2 was close to adiabatic, where cloud-droplet effective radius and liquid water content linearly increase from cloud base to cloud top with less drizzle. The cloud along Leg-1 was sub-adiabatic with lower cloud-droplet number concentration and larger cloud-droplet effective, but higher drizzle droplet number concentration, larger drizzle droplet median diameter and drizzle liquid water content. The heavier drizzle frequency and intensity on Leg-1 were enhanced by the collision-coalescence processes within cloud due to strong turbulence. The sub-cloud precipitation rate on Leg-1 was significantly higher than that along Leg-2. As a result, the sub-cloud accumulation mode aerosols and CCN on Leg-1 were depleted, but the coarse model aerosols increased. This further leads to a counter-intuitive phenomenon that the CCN is less than cloud-droplet number concentration for Leg-1. The average CCN loss rates are ?3.89 \begin{document}$\mathrm{c}{\mathrm{m}}^{-3}\;{\mathrm{h}}^{-1}$\end{document} and ?0.77 \begin{document}$\mathrm{c}{\mathrm{m}}^{-3}\;{\mathrm{h}}^{-1}$\end{document} on Leg-1 and Leg-2, respectively. The cloud and drizzle heterogeneities inside the same stratocumulus can significantly alter the sub-cloud aerosols and CCN budget. Hence it should be treated with caution in the aircraft assessment of aerosol-cloud-precipitation interactions.