Core and margin in warm convective clouds . Part I : core types and evolution 1 during a cloud ' s lifetime 2 1

Abstract. The properties of a warm convective cloud are determined by the competition
between the growth and dissipation processes occurring within it. One way to
observe and follow this competition is by partitioning the cloud to core and
margin regions. Here we look at three core definitions, namely positive vertical
velocity (Wcore), supersaturation (RHcore), and positive buoyancy (Bcore),
and follow their evolution throughout the lifetime of warm convective
clouds. Using single cloud and cloud field simulations with bin-microphysics
schemes, we show that the different core types tend to be subsets of one
another in the following order: Bcore⊆RHcore⊆Wcore. This property is seen for several different thermodynamic profile
initializations and is generally maintained during the growing and mature
stages of a cloud's lifetime. This finding is in line with previous works
and theoretical predictions showing that cumulus clouds may be dominated by
negative buoyancy at certain stages of their lifetime. The RHcore–Wcore
pair is most interchangeable, especially during the growing stages of the
cloud. For all three definitions, the core–shell model of a core (positive values)
at the center of the cloud surrounded by a shell (negative values) at the
cloud periphery applies to over 80 % of a typical cloud's lifetime. The
core–shell model is less appropriate in larger clouds with multiple cores
displaced from the cloud center. Larger clouds may also exhibit buoyancy
cores centered near the cloud edge. During dissipation the cores show less
overlap, reduce in size, and may migrate from the cloud center.



Introduction
Clouds are important players in the climate system (Trenberth et al., 2009) and currently constitute one of the largest uncertainties in climate and climate change research (IPCC, 2013).One of the reasons for this large uncertainty is the complexity created by opposing processes that occur at the same time but in different locations within a cloud.Although a cloud is generally considered to be a single entity, physically, it can be partitioned to two main regions: (i) a core region, where mainly cloud growth processes occur (i.e., condensation -accumulation of cloud mass), and (ii) a margin region, where cloud suppression processes occur (i.e., evaporation -loss of cloud mass).Changes in thermodynamic or microphysical (aerosol) conditions impact the processes in both regions (sometimes in different ways) and thus the resultant total cloud properties (Dagan et al., 2015).To better understand cloud properties and their evolution in time, it is necessary to understand the interplay between physical processes within the core and margin regions (and the way they are affected by perturbations in the environmental conditions).
Considering convective clouds, there are several objective measures that have been used in previous works for separating a cloud's core from its margins (this will be referred to as physical cores hereafter).In deep convective cloud simulations the core is usually defined by the updrafts' magnitude using a certain threshold, usually W > 1 m s −1 (Khairoutdinov et al., 2009;Kumar et al., 2015;Lebo and Seinfeld, 2011;Morrison, 2012).Studies on warm cumulus clouds have defined the clouds' core as parts with positive buoyancy and positive updrafts (Dawe and Austin, 2012;de Roode et al., 2012;Heus and Jonker, 2008;Siebesma and Cuijpers, 1995) Published by Copernicus Publications on behalf of the European Geosciences Union.
or solely regions with positively buoyancy (Heus and Seifert, 2013;Seigel, 2014).More recently, cloud partition to regions of supersaturation, and sub-saturation has been used to define the cloud core in single cloud simulations (Dagan et al., 2015).
For simplicity, we focus on warm convective clouds (only contain liquid water), avoiding the additional complexity and uncertainties associated with mixed-phase and ice-phase microphysics.The common assumption when partitioning a convective cloud to its physical core and margin is that the cloud core is at its geometrical center and the peripheral regions (i.e., edges) are the margin.Previous observational (Heus et al., 2009a;Rodts et al., 2003;Wang et al., 2009) and numerical (Heus and Jonker, 2008;Jonker et al., 2008;Seigel, 2014) works have studied the gradients of cloud thermodynamic properties from the cloud center to its edge and suggest that a cloud is best described by a core-shell model.This model assumes a core with positive vertical velocity and buoyancy, surrounded by a shell with negative vertical velocity and buoyancy.The shell is the region where mixing between cloudy and environmental air parcels occurs, leading to evaporative cooling followed by a decrease in buoyancy and then a decrease in vertical velocity.The cloud shell serves as a buffer between the core and the environment, and its extent is affected by, among others, environmental humidity, aerosol concentrations, and the magnitude and radius of the updraft creating the cloud (Dawe and Austin, 2011;Hannah, 2017;Seigel, 2014).
Based on previous findings, here we explore the partition of clouds to the core and margin using three different objective core definitions where the cloud core threshold is set to be a positive value (of buoyancy, vertical velocity, or supersaturation).Cloud buoyancy (B) can be approximated by the following formula: where θ o represents the reference-state potential temperature, q v is the water vapor mixing ratio, and q l is the liquid water content.The ( ) stands for the deviation from the reference state per height (Wang et al., 2009).Buoyancy is a measure for the vertical acceleration, and its integral is the convective potential energy.Latent heat release during moistadiabatic ascent fuels positive buoyancy and clouds' growth, while evaporation and subsequent cooling drives cloud decay (Betts, 1973;de Roode, 2007).The prevalence of negatively buoyancy parcels at the cloud edges due to mixing and evaporation is a well-known phenomenon (Morrison, 2017).Mixing diagrams have been used to assess this effect (de Roode, 2007;Paluch, 1979;Taylor and Baker, 1991) and are at the root of convective parameterization schemes (Emanuel, 1991;Gregory and Rowntree, 1990;Kain and Fritsch, 1990) and parameterizations of entrainment and detrainment in cumulus clouds (de Rooy and Siebesma, 2008;Derbyshire et al., 2011).
Neglecting cases of air flow near obstacles or air mass fronts, buoyancy is the main source for vertical momentum in the cloud.In its simplest form, the vertical velocity (w) in the cloud can be approximated by the convective available potential energy (CAPE) of the vertical column up to that height (Rennó and Ingersoll, 1996;Williams and Stanfill, 2002;Yano et al., 2005): (2) Here we define CAPE to be the vertical integral of buoyancy from the lowest level of positive buoyancy (h 0 ; initiation of vertical velocity) to an arbitrary top height (h).Usually, the CAPE serves as a theoretical upper limit, and the vertical velocity is smaller due to multiple effects (de Roode et al., 2012), most importantly the perturbation pressure gradient force (which opposes the air motion) and mixing with the environment (entrainment or detrainment; de Roode et al., 2012;Morrison, 2016a;Peters, 2016).Recent studies have shown that entrainment effects on vertical velocity are of the second order, and the forces acting on a rising thermal show a balance between buoyancy and the perturbation pressure gradient (Hernandez-Deckers and Sherwood, 2016;Romps and Charn, 2015), the latter acting as a drag force on the updrafts.Nevertheless, initial updraft and environmental conditions play a crucial role in determining the magnitude of mixing effects on buoyancy and thus also the vertical velocity profile in the cloud (Morrison, 2016a(Morrison, , b, 2017)).The supersaturation (S; where S = 1 is 100 % relative humidity) core definition (S − 1 > 0 or RH > 100 %) partitions the cloud core and margin to areas of condensation and evaporation.Since we consider convective clouds, the only driver of supersaturation during cloud growth is upward vertical motion of air.Neglecting mixing with the environment, S and w can be linked as follows: where Q 1 and Q 2 are thermodynamic factors (Rogers and Yau, 1989).The thermodynamic factors are nearly insensitive to pressure for temperature above 0 • C, and both weakly decrease (less than 15 % net change) with temperature increase between 0 and 30 • C (Pinsky et al., 2013).The first term on the right-hand side is related to the change in the supersaturation due to adiabatic cooling or heating of the moist air (due to vertical motion).The second term is related to the change in the supersaturation due to condensation of vapor or evaporation of drops.Hence, the supersaturation in a rising parcel depends on the magnitude of the updraft and on the condensation rate of vapor to drops (a sink term).The latter is proportional to the concentration of aerosols in the cloud (Reutter et al., 2009;Seiki and Nakajima, 2014), which serve as cloud condensation nuclei (CCN) for cloud droplets.In Part 2 of this work (Heiblum et al., 2019), we demonstrate some of the insights gained by investigating differences between the different cores properties and their time evolution when changing the aerosol loading.The purpose of this part of the work (Part 1) is to compare and understand the differences between the three basic definitions of cloud core (i.e., W core , RH core , and B core ) throughout a convective cloud's lifetime, using both theoretical arguments and numerical simulations.Here, all simulated clouds are analyzed.It should be noted that the binmicrophysical schemes used here calculate saturation explicitly by solving the diffusion growth equation, enabling supersaturation and sub-saturation values in cloudy pixels.This is in contrast to many other works that used bulk-microphysical schemes which rely on saturation adjustment to 100 % within the cloud (Khain et al., 2015).This difference may produce significant differences in the evolution of clouds and their cores.Specifically, we aim to answer questions such as the following: -Which core type is largest?Which is smallest?
-How do the cores change during the lifetime of a cloud?-Can different core types be used interchangeably without much effect on analysis results?
-Are the cores centered at the clouds' geometrical center, as expected from the core-shell model?
It should be noted that previous works tracking clouds throughout their lifetime (e.g., Dawe and Austin, 2012;Heiblum et al., 2016a;Heus et al., 2009b) have reported multi-pulse core growth in cumulus clouds, where multiple buoyancy cores may initiate successively near the cloud base and fuel the cloud.However, these findings did not directly track the cores and were based mainly on the largest, most long-lived clouds.The differences between the cores' evolution in time shed new light on the competition of processes within a cloud in time and space.Moreover, such an understanding can serve as a guideline for all studies that perform the partition to cloud core and margin and assist in determining the relevance of a given partition.

Single cloud model
For single cloud simulations we use the Tel Aviv University axisymmetric, non-hydrostatic, warm convective single cloud model (TAU-CM).It includes a detailed (explicit) treatment of warm cloud microphysical processes solved by the multi-moment bin method (Feingold et al., 1988(Feingold et al., , 1991;;Tzivion et al., 1989Tzivion et al., , 1994)).The warm microphysical processes included in the model are nucleation, diffusion (i.e., condensation and evaporation), collision-coalescence, breakup, and sedimentation (for a more detailed description, see Reisin et al., 1996).Convection was initiated using a thermal perturbation near the surface.A time step of 1 s is chosen for dynamical computations, and 0.5 s is chosen for the microphysical computations (e.g., condensation-evaporation). The total simulation time is 80 min.There are no radiation processes in the model.The domain size is 5 km×6 km, with an isotropic 50 m resolution.The model is initialized using a Hawaiian thermodynamic profile, based on the 91285 PHTO Hilo radiosonde at 00Z on 21 August 2007.A typical oceanic size distribution of aerosols is chosen (Altaratz et al., 2008;Jaenicke, 1988), with a total concentration of 500 cm −3 .This concentration produced clouds that are non-precipitating to weakly precipitating.In Part 2 additional aerosol concentrations are considered, including ones which produce heavy precipitation.

Cloud field model
Warm cumulus cloud fields are simulated using the System for Atmospheric Modeling (SAM) model (version 6.10.3; for details see the following web page: http://rossby.msrc.sunysb.edu/~marat/SAM.html,last access: 10 October 2018; Khairoutdinov and Randall, 2003).SAM is a non-hydrostatic, anelastic model.Cyclic horizontal boundary conditions are used together with damping of gravity waves and maintaining temperature and moisture gradients at the model top.An explicit spectral bin-microphysics (SBM) scheme (Khain et al., 2004) is used.The scheme solves the same warm microphysical processes as in the TAU-CM single cloud model and uses an identical aerosol size distribution and concentration (i.e., 500 cm −3 ) for the droplet activation process.
We use the BOMEX case study as our benchmark for shallow warm cumulus fields.This case simulates a tradewind cumulus (TCu) cloud field based on observations made near Barbados during June 1969 (Holland and Rasmusson, 1973).This case study has a well-established initialization setup (sounding, surface fluxes, and surface roughness) and large-scale forcing setup (Siebesma et al., 2003).It has been thoroughly tested in many previous studies (Grabowski and Jarecka, 2015;Heus et al., 2009b;Jiang and Feingold, 2006;Xue and Feingold, 2006).To check the robustness of the cloud field results, two additional case studies are simulated: (1) the same Hawaiian profile used to initiate the single cloud model and (2) a continental shallow cumulus convection case study (named CASS), based on long-term observations taken at the ARM Southern Great Plains (SGP) site (Zhang et al., 2017).
The soundings, large-scale forcing, and surface properties used to initialize the model are detailed in previous works (Heiblum et al., 2016a;Siebesma et al., 2003;Zhang et al., 2017) to 100 m in the horizontal direction and 40 m in the vertical direction for all simulations.For CASS, above a height of 5 km, the vertical grid size gradually increases to 1 km.The time step for computation is 1 s for all simulations, with a total runtime of 8 h for BOMEX and Hawaii and 12 h for CASS.The initial temperature perturbations (randomly chosen within ±0.1 • C) are applied near the surface during the first time step.

Physical and geometrical core definitions
A cloudy pixel is defined here as a grid box with liquid water amount that exceeds 0.01 g kg −1 .The physical core of the cloud is defined using three different definitions: (1) RH core , all grid boxes for which the relative humidity (RH) exceeds 100 % and condensation occurs; (2) B core , buoyancy (see definition in Eq. 1) above zero, where the buoyancy is determined in each time step by comparing each cloudy pixel with the mean thermodynamic conditions for all non-cloudy pixels per vertical height; and (3) W core , vertical velocity above zero.These definitions apply for both the single cloud and cloud field model simulations used here.We note that setting the core thresholds to positive values (> 0) may increase the amount of non-convective pixels which are classified as part of a physical core, especially for the W core .Indeed, taking higher thresholds for the W core (e.g., W > 0.2 m s −1 ) decreases the W core extent in the cloud and reduces the variance of W core fractions between different clouds in a cloud field (as seen in Fig. 4).Nevertheless, any threshold taken is subjective in nature, while the positive vertical velocity definition is based on the process and is objective.
The centroid (i.e., mean location in each of the axes) and center of gravity (i.e., cloud center of mass) are used here to represent the geometrical location of the total cloud (i.e., cloud geometrical core) and its specific physical cores.The distances between the total cloud and its cores (D norm ), as presented here, are normalized to the cloud size to reflect the relative distance between the two centroids or centers of gravity (COGs), where D norm = 0 indicates coincident physical and geometrical cores and D norm = 1 indicates a core located at the cloud boundary.In case more than one core exists in a cloud, D norm is calculated for each of the cores, and then a mass-weighted (for each core) mean D norm is taken to represent the entire cloud.The single cloud simulations rely on an axisymmetric model, and thus all centroids are horizontally located on the center axis while vertical deviations are permitted.For this model the distance is normalized by half the cloud's thickness.For the cloud field simulations both horizontal and vertical deviations are possible; therefore distances are normalized by the maximum distance from the centroid or COG to a pixel at the cloud's edge.

Center of gravity vs. mass (CvM) phase space
Recent studies (Heiblum et al., 2016a, b) suggested the center of gravity vs. mass (CvM) phase space as a useful approach for reducing the high dimensionally and studying results of large statistics of clouds during different stages of their lifetimes (such as seen in cloud fields).In this space, the COG height and mass of each cloud in the field at each output time step (taken here to be 1 min) are collected and projected in the CvM phase space.This enables a compact view of all clouds in the simulation during all stages of their lifetimes, with the main disadvantage being the loss of gridsize resolution information on in-cloud dynamical processes.Although the scatter of clouds in the CvM is sensitive to the microphysical and thermodynamic settings of the cloud field, it was shown that the different subspaces in the CvM space correspond to different cloud processes and stages (Heiblum et al., 2016a, b).The lifetime of a cloud can be described by a trajectory on this phase space.
A schematic illustration of the CvM space is shown in Fig. 1.Most clouds are confined between the adiabat (curved, dashed line) and the inversion layer base (horizontal dashed line).The adiabat curve corresponds to the theoretical evolution of a moist adiabat 1-D cloud column in the CvM space.The large majority of clouds form within the growing branch (yellow shade) at the bottom left part of the space, adjacent to the adiabat.Clouds then follow the growing trajectory (grow in both COG and mass) to some maximal values.The growing branch deviates from the adiabat at large masses, depending on the degree of sub-adiabaticity of the cloud field (i.e., the degree of mixing between the cloud and its surrounding environment), which depends on its thermodynamic profile.After or during the growth stage of clouds, they may undergo the following processes: they (i) dissipate via a quasi-reverse trajectory adjacent to the growing one, (ii) dissipate via a gradual dissipation trajectory (magenta shade), (iii) shed small mass cloud fragments (red shades), and (iv) in the case of precipitating clouds, they can shed cloud fragments in the sub-cloudy layer (grey shade).The former two processes form continuous trajectories in the CvM space, while the latter two processes create disconnected subspaces.

Cloud tracking
To follow the evolution of individual clouds within a cloud field, we use an automated 3-D cloud tracking algorithm (see Heiblum et al., 2016a, for details).It enables tracking of continuous cloud entities (CCEs) from formation to dissipation, even if interactions between clouds (splitting or merging) occur during that lifetime.A CCE initiates as a new cloud forming in the field and is tracked under the condition that it retains the majority (> 50 %) of its mass during an interaction event with another cloud.Thus, a CCE can terminate due to either cloud dissipation or cloud interactions.

Theoretical estimations for different core sizes
Here we propose simple physical considerations to evaluate the differences in cloud partition to core and margin using different definitions.The arguments rely on key findings from previous works (see Sect. 1), with the aim of gaining intuitive understanding of the potential differences between the core types.It is convenient to separate the analysis to an adiabatic case and then add another layer of complexity and consider the effects of mixing of cloudy and noncloudy air.In this theoretical derivation, saturation adjustment to RH = 100 % is assumed for both cases, while in the other models used in this study, transient supersaturated and sub-saturated cloudy parcels are treated (more realistic).

Adiabatic case -no mixing
Considering moist-adiabatic ascent, the excess vapor above saturation is instantaneously converted to liquid (saturation adjustment).Thus, the adiabatic cloud is saturated (S = 1) throughout its vertical profile, and only W core and B core differences can be considered.It is assumed that the adiabatic convective cloud is initiated by positive buoyancy initiating from the sub-cloudy layer.As long as the cloud is growing it should have positive CAPE and will experience positive w throughout the column even if the local buoyancy at a spe-cific height is negative.Eventually the cloud must decelerate due to negative buoyancy and reach a top height, where CAPE = 0 and w = 0. Hence, for the adiabatic column case, B core is always a proper subset of W core (i.e., B core ⊂ W core ).These effects are commonly seen in warm convective cloud fields where permanent vertical layers of negative buoyancy (but with updrafts) within clouds typically exist at the bottom and top regions of the cloudy layer (Betts, 1973;de Roode and Bretherton, 2003;Garstang and Betts, 1974;Grant and Lock, 2004;Heus et al., 2009b;Neggers et al., 2007).

Cloud parcel entrainment model
A mixing model between a saturated (cloudy) parcel and a dry (environment) parcel is used to illustrate the effects of mixing on the different core types.The details of these theoretical calculations are shown in Appendix A. The initial cloudy parcel is assumed to be saturated (part of RH core ), have positive vertical velocity (part of W core ), and experience either positive or negative buoyancy (part of B core or B margin ), as is seen for the adiabatic column case.Additionally, mixing is assumed to be isobaric and in a steady environment where the average temperature of the environment per a given height does not change.The resultant mixed parcel will have lower humidity content and lower LWC, as compared to the initial cloudy parcel, and a new temperature.In nearly all cases (beside in an extremely humid environment) the mixed parcel will be sub-saturated and evaporation of LWC will occur.Evaporation ceases when equilibrium is reached due to air saturation (S = 1) or due to complete evaporation of the droplets (which means S < 1, and the mixed parcel is no longer cloudy, since it has no liquid water content).
In addition to mixing between cloudy (core or margin) and non-cloudy parcels, mixing between core and margin parcels (within the cloud) also occurs.This mixing process can be considered to be "entrainment-like" with respect to the cloud core.Considering the changes in the W core and RH core , there is no fundamental difference in the treatment of mixing of cloudy and non-cloudy parcels or mixing between core and margin (because the margins and the environment are typically sub-saturated and experience negative vertical velocity).However, for the changes in the B core after mixing, a fundamental difference exists between mixing with the reference temperature or humidity state (in the case of mixing with the environment) and mixing given a reference temperature or humidity state (in mixing between B core and B margin ).Thus, it is interesting to check the effects of mixing between B core and B margin parcels on the total extent of the B core with respect to the other two core types.The details of this second case are shown in Appendix B.

Effects of non-cloudy entrainment on buoyancy
When mixed with non-cloudy air, the change in buoyancy of the initial cloudy parcel (which is a part of W core and RH core and either B core or B margin ) happens due to both mixing and evaporation processes.The theoretical calculations show that for all relevant temperatures (∼ 0 to 30 • C, representing warm Cu), the change in the parcel's buoyancy due to evaporation alone will always be negative (see Appendix A).This is because the negative effect of the temperature decrease outweighs the positive effects of the humidity increase and water loading decrease.Nevertheless, the total change in the buoyancy (due to both mixing and evaporation) depends on the initial temperature, relative humidity, and liquid water content of the cloudy and non-cloudy parcels.
In Fig. A1 in Appendix A a wide range of non-cloudy environmental parcels, each with their own thermodynamic conditions, are mixed with a saturated cloud parcel with either positive or negative buoyancy.The main conclusions regarding the effects of such mixing on the buoyancy are as follows: 1. To a first order, the initial buoyancy values are temperature dependent, where a cloudy parcel that is warmer (colder) by more than ∼ 0.2 • C than the environment will be positively (negatively) buoyant for common values of cloudy layer environmental relative humidity (RH > 80 %).
2. Parcels that are initially part of B core may only lower their buoyancy due to entrainment; the change results in either to positive or negative buoyancy values, depending on the environmental conditions.
3. The lower the environmental RH, the larger the probability for parcel transition from B core to B margin after entrainment.
4. Parcels that are initially part of B margin can either increase or decrease their buoyancy value but never become positively buoyant.The former case (buoyancy decrease) is expected to be more prevalent, since it occurs for the smaller range of temperature differences with the environment.
In summary, entrainment is expected to always have a net negative effect on B core extent and B margin values, while evaporation feedbacks serve to maintain RH core in the cloud.Thus, we can predict that B core should be a subset of RH core (i.e., B core ⊆ RH core ).

Effects of core and margin mixing on buoyancy
We consider the case of mixing between the B core and B margin , meaning positively buoyant and negatively buoyant cloud parcels.For simplicity, we assume that both parcels are saturated (S = 1; both included in the RH core ).As seen above, such conditions exist in both the adiabatic case and in the case where an adiabatic cloud has undergone some entrainment with the environment.The buoyancy differences between the saturated parcels are mainly due to temperature differences but also due to the increasing saturation vapor pressure with increasing temperature (see Appendix B for details).In Fig. B1 in Appendix B is it shown that the resultant mixed parcel's buoyancy can be either positive or negative, depending on the magnitude of temperature difference of each parcel (core or margin) from that of the environment.However, in all cases the mixed parcel is supersaturated.This result can be generalized: given two parcels with equal RH but different temperature, the RH of the mixed parcel is always equal to or higher than the initial value.Hence, B core can either increase or decrease in extent, while the RH core can only increase due to mixing between saturated B core and B margin parcels.This again strengthens the assumption that B core should be a subset of RH core .
We note that an alternative option for mixing between the core and margin parcels exists here, in which either or both of the parcels are subsaturated so that the mixed parcel is subsaturated as well.In this case evaporation will also occur.As seen in Appendix A, this should further reduce the buoyancy value of the mixed parcel (while increasing the RH).

Effects of entrainment on vertical velocity
The vertical velocity equation dictates that buoyancy is the main production term (de Roode et al., 2012;Romps and Charn, 2015) and is balanced by perturbation pressure gradients and mixing (on grid and sub-grid scales).Thus, all changes of magnitude (and sign) in vertical velocity should lag the changes in buoyancy.This is the basis of convective overshooting and cumulus formation in the transition layer (see Sect. 3.1).It is interesting to assess the magnitude of this effect by quantifying the expected time lag between buoyancy and vertical velocity changes.The calculations in Appendix A indicate negative buoyancy values reaching −0.1 m s −2 due to entrainment.However, measurements from within clouds show that the temperature deficiency of cloudy parcels with respect to the environment is generally restricted to less than 1 • C for cumulus clouds (Burnet and Brenguier, 2010;Malkus, 1957;Sinkevich and Lawson, 2005;Wei et al., 1998), and thus the negative buoyancy should be no larger than −0.05 m s −2 .This value is closer to current and previous simulations and also observations that show negative buoyancy values within clouds to be confined between −0.001 and −0.01 m s −2 (Ackerman, 1956;de Roode et al., 2012).
Given an initial vertical velocity of ∼ 0.5 m s −1 , the deceleration due to buoyancy (and reversal to negative vertical velocity) should occur within a typical time range of 1-10 min.These timescales are much longer than the typical timescales of evaporation (that eliminates the B core ), which range between 1-10 s (Lehmann et al., 2009).Moreover, the fact that a drag force typically balances the buoyancy acceleration (Romps and Charn, 2015) can also contribute to a time lag between effects on buoyancy and subsequent effects on vertical velocity.Therefore, the switching of the sign for vertical velocity should occur with substantial delay compared to the reduction of buoyancy, and B core should be a subset of W core (i.e., B core ⊆ W core ) during the growing and mature stages of a cloud's lifetime.

The relation between supersaturation and vertical velocity cores
Here we revisit the terms in Eq. ( 3) to explore an intuitive, first-order understanding of the relation between the vertical velocity core and the supersaturation core.A rising parcel initially has no liquid water content, with its only source of supersaturation being the updraft w, and thus initially the RH core should always be a subset of W core .In general, since the sink term dq l dt becomes a source only when S < 1 (the condition for evaporation), the only way for a convective cloud to produce supersaturation (i.e., S > 1) is by updrafts during all stages of its lifetime.Once supersaturation is achieved, the sink term becomes positive dq l dt > 0 and balances the updraft source term so that supersaturation either increases or decreases.At any stage, if downdrafts replace the updrafts within a supersaturated parcel, the consequent change in supersaturation becomes strictly negative (i.e., dS dt < 0).This negative feedback limits the possibility of finding supersaturated cloudy parcels with downdrafts.Hence, we can expect the RH core to be smaller than W core during the majority of a cloud's lifetime.

Results -single cloud simulation
The differences between the three types of core definitions are examined during the lifetime of a single cloud (Fig. 2), based on the Hawaiian profile.The cloud's total lifetime is 36 min (between t = 7 and t = 43 min of simulation).Each panel in Fig. 2 presents vertical cross-sections of the three cores (magenta -W core , green -RH core , and yellow -B core ) at four points in time (with 10 min intervals).The cloud has an initial cloud base at 850 m and grows to a maximal top height of 2050 m.The condensation rates (red shades) increase toward the cloud center, and the evaporation rates (blue shades) increase toward the cloud edges.Evaporation at the cloud top results in a large eddy below it that contributes to mixing and evaporation at the lateral boundaries of the cloud.Thus, a positive feedback is initiated which leads to cooling, negative buoyancy, and downdrafts.The dissipation of the cloud is accompanied with a rising cloud base and lowering of the cloud top.
During the growing stage (t = 10, 20 min), when substantial condensation still occurs within the cloud, all of the cores seem to be self-contained within one another, with B core be-ing the smallest and W core being the largest.During the final dissipation stages, when the cloud shows only evaporation (t = 40), W core and RH core disappear, while there is still small B core near the cloud top.Further analysis (see Part 2) shows that the entire dissipating cloud is colder and more humid than the environment, but downdrafts from the cloud top (see arrows in Fig. 2) promote heating, and by that increase the buoyancy in dissipating cloudy pixels, sometimes reaching positive values.These buoyant pockets will be discussed further in Part 2. The results indicate that the three types of physical cores of the cloud are not located around the cloud's geometrical core along the whole cloud lifetime.During cloud growth (i.e., increase in mass and size) the three types of cores surround the cloud's center, while during late dissipation the B core is offset from the cloud center.
For a more complete view of the evolution of the three core types in the single cloud case, time series of core fractions are shown in Fig. 3. Figure 3a and b show the core liquid mass (core mass / total mass − f mass ) and volume (core volume / total volume -f vol ) fractions out of the cloud's totals.The results are similar for both measures, except for the fact that core mass fractions are larger than core volume fractions.This is due to significantly higher LWC per pixel in the cores compared to the margins, which skews the core mass fraction to higher values.Core mass fractions during the main cloud growing stage (between t = 7 and t = 27 min simulation time) are around 0.7-0.85,and core volume fractions are around 0.5-0.7.The time series show that as opposed to the W core and RH core fractions, which decrease monotonically with time, B core shows a slight increase during stages of cloud growth.In addition, for most of the cloud's lifetime, the B core fractions are the smallest and the W core fractions are the largest, except for the final stage of the clouds dissipation, where downdrafts from the cloud top creates pockets of positive buoyancy.These pockets are located at the cloud's peripheral regions rather than near the cloud's geometrical center as is typically expected for the cloud's core.In the cloud's center (the geometrical core) the B core is the first one to terminate (at t = 32 min) compared to both W core and RH core , which decay together (at 36 min).
For describing the locations of the physical cores, we examine the normalized distances (D norm ) between the cloud's centroid and the cores' centroids.The evolution of these distances is shown in Fig. 3c.At cloud initiation (t = 7 min), when the cloud is very small, all of the cores' centroids coincide with the total cloud centroid location.The B core (and RH core to a much lesser degree) centroid then deviates from the cloud centroid to a normalized distance of 0.27 (t = 8 min).As cloud growth proceeds, B core grows and its centroid coincides with the cloud's centroid.All of the cores' centroids are located near the cloud centroid during the majority of the growing and mature stages of the cloud, showing normalized distances < 0.1.During dissipation (t > 27 min), the cores' centroid locations start to move away from the cloud's geometrical core; this is followed by a reduction in  Finally, in Fig. 3d the fraction of pixels of each core contained within another core is shown.It can be seen that for the majority of cloud lifetime (up to t = 33 min) B core is a subset (pixel fraction of 1) of RH core , and the latter is a subset of W core .As expected, the other three permutations of pixel fractions (e.g., W core in B core ) show much lower val-ues.The cloudy regions that are not included within B core but are included within the two other cores are exclusively at the cloud's boundaries (see Fig. 2).The same pattern is seen for cloudy regions that are included within W core but not in RH core .During the dissipation stage of the cloud its core subset property (i.e B core ⊆ RH core ⊆ W core ) breaks down.Similar temporal evolutions to those shown here are seen for the other simulated clouds (with various aerosol concentrations) in Part 2 of this work.
5 Results -cloud field simulations

Partition to different core types
To test the robustness of the observed behaviors seen for a single cloud, it is necessary to check whether they also apply to large statistics of clouds in a cloud field.The BOMEX simulation is taken for the analyses here.We discard the first 3 h of cloud field data, during which the field spins up and its mean properties are unstable.In Fig. 4 the volume (f vol ) and mass (f mass ) fractions of the three core types are compared for all clouds (at all output times -every 1 min) in the CvM space.As seen in Fig. 1, the location of specific clouds in the CvM space indicates their stage in evolution.Most clouds are confined to the region between the adiabat and the inversion layer base except for small precipitating (lower-left region) and dissipating clouds (upper-left region).The color shades of the clouds indicate whether a cloud is all core (red -core fraction 1), all margin (blue -core fraction 0), or equally divided into the core and margin (white -core fraction 0.5).The size of each point in the scatter is proportional to the cloud's mean horizontal cross-sectional area.A general increase in mean cloud area with increase in mean cloud LWP is seen (i.e., synchronous growth in the horizontal and vertical axis).
As seen for the single cloud, the core mass fractions tend to be larger than core volume fractions for all core types.This is due to the fact that LWC values in the cloud core regions are higher than in margin regions, so a cloud might be core dominated in terms of mass while being margin dominated in terms of volume.Focusing on the differences between core types, the color patterns in the CvM space imply that B core definition yields the lowest core fractions (for both mass and volume), followed by RH core with higher values and W core with the highest values.The absence of the B core is especially noticeable for small clouds in their initial growth stages after formation (COG ∼ 550 m and LWP < 1 g m −2 ).Those same clouds show the highest core fractions for the other two core definitions.This large difference can be explained by the existence of the transition layer (as discussed in Sect.3) near the lifting condensation level (LCL) in warm convective cloud fields, which is the approximated height of a convective cloud base (Craven et al., 2002;Meerkötter and Bugliaro, 2009).Within this layer parcels rising from the sub-cloudy layer are generally colder than parcels subsiding from the cloudy layer.Thus, this transition layer clearly marks the lower edge of the buoyancy core, as most convective clouds are initially negatively buoyant.
Generally, the growing cloud branch (i.e., the CvM region closest to the adiabat) shows the highest core fractions.The RH core and W core fractions decrease with cloud growth (increase in mass and COG height), while the B core initially increases, shows the highest fraction values around the middle region of the growing branch, and then decreases for the largest clouds.The transition from the growing branch to the dissipation branch is manifested by a transition from coredominated to margin-dominated clouds (i.e., transition from red to blue shades).Mixed within the margin-dominated dissipating cloud branch, a scatter of W core -dominated small clouds can be seen as well.These represent cloud fragments which shed large clouds during their growing stages with positive vertical velocity.They are sometimes RH core dominated as well but are strictly negatively buoyant.The few precipitating cloud fragments seen for this simulation (cloud scatter located below the adiabat) tend to be margin dominated, especially for the RH core .
The percentages in the panel legends (Fig. 4) indicate the fractions of clouds (out of the scatter) which are core dominated with respect to volume or mass.Only ∼ 2 % of clouds are dominated by B core in terms of cloud volume, but more than 45 % of the clouds have the majority of their mass within the B core region.These numbers increase considerably for the RH core (W core ), where 44 % (80 %) of the clouds are core dominated with respect to cloud volume and 85 % (87 %) of the clouds are core dominated with respect to cloud mass.Thus, the B core can be considered to be taking up a small portion of a typical cloud mass and volume, while the W core generally occupies most of the cloud.We note that some of the largest clouds in the field (indicated by large scatter points) show higher (lower) B core (RH core and W core ) volume fractions in comparison with smaller clouds located adjacent to them in the CvM phase space.Further analysis shows that these clouds are also precipitating to the surface.The increase in B core fractions in precipitating clouds is discussed in Part 2 of this work.

Subset properties of cores
From Fig. 4 it is clear that W core tends to be the largest and B core tends to be the smallest.To what degree, however, are the cores subsets of one another as those seen for the single cloud simulation?In Fig. 5 the pixel fraction (f pixel ) of each core type within another core type is shown for all clouds in the CvM space.A f pixel of 1 (bright colors) indicates that the pixels of the specific core in question (labeled in each panel title) are a subset of the other core (also labeled in the panel title), and a f pixel of 0 (dark colors) indicates no intersection between the two cores in the cloud.It is seen that B core tends to be a subset of both other cores, with the f pixel being around 0.75-1 for most of the growing branch area and large mass dissipating clouds which still have some positive buoyancy.The pixel fractions are higher for B core inside W core compared with B core inside RH core , but both show a decrease with an increase in growing branch cloud mass, meaning that the chance for finding a proper subset B core decreases in large clouds.
The CvM space of RH core inside W core shows an even stronger relation between these two core types.For almost all growing branch clouds, the RH core is a subset of W core (i.e., RH core ⊆ W core ).The pixel fractions tend to decrease gradually with loss of cloud mass in the dissipation branch.However, some small dissipating clouds show that f pixel = 1.These clouds also experience high core volume fractions (f vol ∼ 1), as indicated by the scatter point sizes in Fig. 5.The other three permutations of the f pixel (W core inside B core , W core inside RH core , and RH core inside B core ) give an indication of cores' sizes and of which cloud types show no overlap between different cores.As stated above, growing (dissipation) clouds show higher (lower) overlap between the different core types.The W core is almost twice as large as the B core and 30 %-40 % larger than the RH core along most of the growing branch.
To give an objective measure of the degree to which different core types can be used interchangeably, we define an interchangeable fraction (f int ), which is the multiplication of the two pixel fractions of a core pair (e.g., f pixel B in RH • f pixel RH in B ).In Fig. 6 the f int is shown for all clouds and the three core pairs.It can be seen that only a small percentage (< 5 %) of clouds can be considered to have fully interchangeable core types with f int > 0.75.The RH core -W core pair shows the highest degree of interchangeability (83 % and 54 % of clouds with f int > 0.25 and f int > 0.5, respectively), showing high f int for clouds at formation and growing stages and sometimes also late dissipation.The B core -W core pair shows the lowest degree of interchangeability (46 % and 6 % of clouds with f int > 0.25 and f int > 0.5, respectively), with mature growing clouds showing the highest f int values.The B core -RH core pair shows similar results, but with slightly higher f int values on average.

Revisiting the core-shell model
Here we test how well the core-shell model can be applied to the three types of cores in different clouds seen in a warm cumulus cloud field.We test both the location of the cores with respect to the cloud center and horizontal profiles of the three types of core parameters within the cloud.In Fig. 7 the normalized distances between the total cloud centroid and each specific physical core centroid locations (i.e., D norm,Centroid ) are evaluated.Since clouds are not always axisymmetric, we also test the distances between total cloud COG and core COG (D norm,COG ), since the COG gives a better representation for where cloud and core mass is concentrated.We take D norm < 0.2 as a threshold for cores located near the centroid or COG and D norm > 0.8 as a threshold for cores located at the cloud edges.For all core types, the large majority of clouds' cores are centered near the clouds' centroid or COG.Only less than 1 % of the clouds' cores reside at the cloud edges, mostly seen for small dissipating clouds.Distances between the cloud COG and core COG yield smaller values than for distances between centroids, implying that the mass is not equally distributed within the clouds, and hence the centroid may be "missing" the true cloud center in terms of mass distribution.Along the growing branch the clouds and physical cores tend to be centered in close proximity, while during cloud dissipation the cores tend to increase in distance from the cloud's center.This type of evolution is most prominent for the W core , which shows a clear gradient of transition from small (dark colors) to large (bright colors) distances.Focusing on D norm,COG < 0.2, the B core shows a lower chance of being in proximity to the cloud COG (76 %) than the other core types (83 %).This may be due to a larger prevalence of cloud-edge B core pixels during dissipation (see Sect. 4 here and Sect.4.2 in Part 2).Compared to the other clouds, the W core shows a slightly larger probability of being located at the cloud edge in small dissipating clouds.This can be due to the fact that during cloud dissipation, complex patterns of updrafts and downdrafts within the cloud can create scenarios where the W core is comprised of very weak updrafts and located anywhere in the cloud.Further analysis shows that most clouds with D norm,COG > 0.2 values can be attributed to the relatively larger-sized clouds which typically contain multiple cores within them (Fig. 8).For the B core , RH core , and W core , 68 %, 79 %, and 81 % of the cloud scatter analyzed (which contains a core) has a single core, respectively.Thus, most clouds have a single core.Moreover, it is more probable to find multiple buoyancy cores in a cloud than vertical velocity cores.This is surprising given our choice of "weak" W core thresholds (i.e., positive values) and indicates that vertical velocity patterns are relatively well behaved in cumulus clouds, at least for the LES scales chosen here.For clouds with a single core, the growing branch cloud COG and core COG are co-located at the same point.A gradual transition to larger distances is seen as the clouds dissipate to lower mean LWP values.In total, over 80 % of single core clouds have D norm,COG < 0.2 for all core types.For clouds with multiple cores, about 50 % of clouds show large distances (D norm,COG > 0.2), with little difference between growing branch and dissipating branch clouds.This is to be expected since a large cloud with multiple cores should have a COG somewhere between those cores, explaining the larger normalized distances.
The core-shell model assumes that the highest values (of a core parameter in question) are located at the center of the cloud (Heus and Jonker, 2008).Is this indeed the case in clouds?In Fig. 9 we observe the likelihood and shape of pre-defined categories of horizontal profiles for core pa-rameters.Profiles are taken along the horizontal plane of the cloud's COG, with distances normalized to cloud maximum horizontal size so that different cloud sizes can be averaged together.Only clouds with at least three pixels in the horizontal plane are taken.Profile categories include the following: (i) core-shell (CS) profiles, which have a positive, maximum value near the COG at D norm,COG < 0.2; (ii) displaced core-shell profiles (DCS), which have a positive, maximum value somewhere between the COG and periphery at 0.2 < D norm,COG < 0.8; (iii) periphery core (PC) profiles, which have a positive, maximum value at the cloud periphery at D norm,COG > 0.8; and (iv) no-core (NC) profiles, which are comprised of only negative values.We take only clouds with a single core (or no core), since clouds with multiple cores show more complex profiles that represent a superposition of several single core profiles.The data are further divided to growing and dissipating stages of clouds by checking if a cloud grew in mass compared to the previous time step.
For all core types, there are more single core (and no-core) growing clouds (∼ 55 %-57 %) than dissipating clouds.Generally, it can be seen that the CS category profile is the most prevalent in clouds with single cores, ranging from a maximum of 66 % of growing cloud W core profiles to a minimum of 26 % of dissipating cloud B core profiles.As seen in Figs. 4, 7, and 8, growing clouds show a relatively higher percentage of the CS and DCS categories, while dissipation clouds show relatively higher percentages of PC-and NC-category profiles.The W core and RH core profiles show similar behavior, with decreasing prevalence from the CS category to the NC category (CS > DCS > PC > NC) for growing clouds.For dissipating clouds, the partition is similar, but with the PC category being the least prevalent (CS > DCS > NC > PC).The main difference in the partition to categories in B core profiles is the increasing prevalence and dominance of the NC category, as seen in previous analyses.For example, NC profiles are almost non-existent in growing clouds for the W core and RH core definitions (<= 1 %) but second-most prevalent using the B core definition (28 %).Out of the three core types, the W core shows the highest probability for matching CS and PC categories, the RH core for the DCS category, and the B core for the NC category.
On average, the CS category profiles show a monotonic decrease in value from positive to negative values.For growing clouds, vertical velocity may stay positive throughout the horizontal profile, not necessarily showing a downdraft "shell".The DCS category profiles show positive values from the COG to more than half the cloud size (or the entire cloud size for growing cloud vertical velocity) and have a maximum at D norm ∼ 0.2, indicating that they are only marginally displaced from the cloud COG.This may indicate that most DCS profiles can actually be considered to be CS profiles, but for clouds with significant asymmetry to the core maximum, they seem to be displaced.Merging of CS and DCS categories comprises 70 %-90 % of all clouds.For both CS and DCS categories, the transition from core to margin (i.e., positive to negative values) occurs at shorter D norm for dissipat-ing clouds.This transition D norm value also decreases, gradually moving from the W core to the RH core and then to the B core , indicating smaller core sizes for the latter core types.The NC profiles show little variance with distance from the COG.
All core (and margin) types show decreasing values moving from growing to dissipating clouds.A decrease in values is also seen when comparing the maximum of the CS, DCS, and PC mean profiles, respectively.An exception is seen for the buoyancy PC category, which shows a slightly higher buoyancy peak value for dissipating clouds.Compared to the RH core and W core PC category profiles, which show positive values throughout the cloud (with little change) for smallerthan-average clouds, the B core PC category shows a transition from margin at the COG to core at the periphery for largerthan-average clouds.This transition to positive buoyancy is even more pronounced (i.e reaches higher values) for dissipating multi-core clouds (not shown here) that tend to be significantly larger.This may indicate that a non-convective process is at play in creating the B core at the cloud periphery (see Part 2).

Consistency of the cloud partition to core types
The results for cloud fields are summarized in Fig. 10, which presents the evolution of core fractions of continuous cloud entities (CCEs; see Sect.2.5 for details) from formation to dissipation.Only CCEs that undergo a complete life cycle are averaged here.These CCEs fulfill the following four conditions: they (i) form near the LCL, (ii) live for at least 10 min, (ii) reach maximum cloud mean LWP values above 10 g m −2 , and (iv) terminate with a mass value below 10 g m −2 .As a test of generality, we performed this analysis for Hawaiian and CASS warm cumulus cloud field simulations in addition to the BOMEX one.For each simulation, hundreds of CCEs are collected (see panel titles), and their core volume fractions are averaged according to their normalized lifetimes (τ ).
Consistent results are seen for all three simulations.Clouds initiate with a W core fraction of ∼ 1, RH core fraction of ∼ 0.8, and B core fraction of ∼ 0-0.15.The former two core types' volume fraction decreases monotonically with lifetime, while the latter core type's volume fraction increases up to 0.15-0.35at τ ∼ 0.3 and then monotonically decreases for increasing τ .The continental (CASS) simulation consistently shows lower buoyancy volume fractions than the oceanic simulations.This can be attributed to lower RH in the CASS cloudy layer (60 %-80 %) compared with the oceanic simulations (85 %-95 %).The lower RH increases entrainment and reduces buoyancy.The fact that clouds end their life cycle with non-zero volume fractions may indicate that some of the CCEs terminate not because of full dissipation but rather because of significant splitting or merging events.
Normalized distances (D norm ) between the CCE core and cloud are also shown in Fig. 10c, f, and i.Both distances between the core centroid (solid lines) and COG (dashed lines) to total cloud centroid and COG are shown.As seen in Fig. 7, D norm,COG shows smaller values, indicating that the COG indicates the "true" cloud center better compared with the centroid.Distances tend to monotonically increase for RH core and W core with the CCE lifetime for all simulations.The gradient of increase is larger at the later stages of the CCE lifetime.Initially the W core is closer to the geometrical core, but at later stages of the CCE lifetime (typically τ > 0.8) this switches and RH core remains the closest.As seen above, for the first (second) half of the CCE lifetime, the B core D norm decreases (increases), starting at normalized distances around 0.2 for all simulations.The physical cores' COG stay closer to the cloud COG (D norm < 0.2) for the majority of their lifetimes for the three cases.Taking again the value D norm,COG = 0.2 as a threshold for physical cores centered near the cloud COG, BOMEX, Hawaii, and CASS simulation CCEs' W core values all cross this threshold at τ = 0.9.Thus, the core-shell geometrical model is true for about 90 % of a typical cloud's lifetime.
The analysis of core subset properties (Fig. 10c, f, i) shows that the assumption B core ⊆ RH core ⊆ W core is true for the initial formation stages of a cloud.Although the corresponding pixel fractions decrease slightly during the lifetime of the CCE, they remain above 0.9 (e.g., B core is 90 % contained within RH core ).A sharp decrease in pixel fractions is seen for τ > 0.8 (τ > 0.5 for the CASS simulation), as the overlaps between the different cores are reduced during dissipation stages of the cloud.For all simulations, the highest pixel fraction values are seen for the B core inside the W core pair, followed by RH core inside the W core pair and B core inside the RH core pair, showing lower values.In addition, it can be seen that the variance of the average pixel fraction (per τ ) increases with an increase in τ .This is due to the fact the all CCEs initiate with almost identical characteristics but may terminate in very different ways.In Part 2 of this work we show that this variance is highly influenced from precipitation, which contributes to more significant interactions between clouds (Heiblum et al., 2016b).

Summary
In this paper we study the partition of warm convective clouds to core and margin according to three different definitions: (i) positive vertical velocity (W core ), (ii) relative humidity supersaturation (RH core ), and (iii) positive buoyancy (B core ), with emphasis on the differences between those definitions.Using theoretical considerations of both an adiabatic cloud column and a simple two-parcel mixing model (see Appendices A and B), we support our simulated results, as we show that the B core is expected to be the smallest of the three.This finding is in line with previous works that showed that negative buoyancy is prevalent in cumulus clouds for a wide range of thermodynamic conditions (de Roode, 2007;Paluch, 1979;Taylor and Baker, 1991).This is due to the fact that entrainment into the core (i.e., mixing with non-cloudy environment or mixing with the margin regions of the cloud) may result in sub-saturation, followed by evaporation that always has a negative net effect on buoyancy.The same process has an opposing effect on the relative humidity of the mixed parcel and acts to reach saturation.Entrainment (or mixing) also acts to decrease vertical velocity but at slower manner compared to the timescales of changes in the buoyancy and relative humidity.In addition, the supersaturation equation (Eq. 3) predicts that it is unlikely to maintain supersaturation in a cloudy volume with negative vertical velocity.Hence, W core can be expected to be the largest of the three cores.
Using numerical simulations of both a single cloud and cloud fields of warm cumulus clouds, we show that during most stages of clouds' lifetime, W core is indeed the largest of the three and B core is the smallest.Only 2 % of clouds are dominated (in volume fraction) by the B core , while 44 % and 83 % of clouds are dominated by RH core and W core , respectively.The warm convective cloud fields simulated here typically have a transition layer near the lifting condensation level (LCL).Thus, the lower parts of the clouds are negatively buoyant or even lack B core at formation.After cloud formation, internal growth processes (i.e., condensation and latent heat release) increase the B core until dissipation processes become dominant and the B core decreases quickly due to entrainment.In contrast, clouds are initially dominated by the W core and RH core (fractions close to 1).The fractions of these cores then decrease monotonically with cloud lifetime.During dissipation stages, the clouds are mostly margin dominated such that most of the small mass dissipation cloud fragments are entirely coreless.However, several small mass dissipating cloud fragments which shed large cloud entities (with large COG height) may be core dominated, especially when using the vertical velocity core definition.The same is observed for small precipitating cloud fragments which reside below the convective cloud base.We note that the results here are similar for both volume and mass core fractions out of the cloud's total fraction, with the core mass fractions being larger due to a skewed distribution of cloud LWC which favors the core regions.Moreover, we show that these results are consistent for various levels of aerosol concentrations (will be seen in Part 2) and different thermodynamic profiles used to initialize the models.
In addition to the differences in their sizes, the three cores tend to be subsets of one another in the following order: B core ⊆ RH core ⊆ W core .This property is most valid for a cloud at its initial stages and breaks down gradually during a cloud's lifetime.The decrease in overlap between different core types during dissipation implies that minor local effects enable core existence rather than cloud convection.Only during growth and mature stages can the three core definitions be used interchangeably with the least amount of difference in core sizes.Generally, the RH core -W core pair are most interchangeable, while the B core -W core pair are the least interchangeable.
With respect to cloud morphology, the majority of clouds are composed from single cores (for all core types), are located near the cloud centroid and COG, and fit the intuitive core-shell model of decreasing core parameter values from the cloud center to periphery.This is especially true during cloud growth, as during dissipation the cores may decouple from the geometrical core and often comprise just a few isolated pixels at the cloud's edges.In terms of cloud lifetime, the core-shell model applies to at least 80 % of a typical cumulus cloud lifetime.We note that using the COG as a measure for the cloud and core geometrical centers yields smaller cloud-core distances than their centroids.Thus, the COG better represents the cloud physical center.Out of the three core types, the W core (B core ) shows the highest (lowest) chance of being a single core in the cloud.This is despite choosing a low W core threshold of W > 0. Relatively large clouds tend to have multiple cores so that the mean (mass weighted) core COG location is displaced from the cloud COG.The B core COG shows the highest chance to be located away from the cloud COG.In some cases of larger clouds, the buoyancy horizontal profile may look exactly opposite to the core-shell one (i.e., maximum at the periphery and minimum at the center).This may be due to downdraft-induced heating at the clouds' edge that promotes positive buoyancy (see more in Part 2).In Part 2 of this work we use the insights gained here on clouds' partition to its core and margin to understand aerosol effects on warm convective clouds.
For a wide temperature range between 200 300 (K), dB evap is always negative.This result is not trivial because evaporation both decreases the T and increases the q v , which have opposite effects.The total change in buoyancy is taken as the sum of dB evap and dB mix .
Figure A1 presents a phase space of possible changes in cloudy pixel buoyancy due to mixing with outside air, for various thermodynamic conditions, and a mixing fraction of 0.5.The initial cloudy parcel is chosen to be saturated (S = 1) and includes an LWC of 1 g kg −1 .The pressure is assumed to be 850 mbar, and the temperature is assumed to be 15 • C.However, we note that the conclusions here apply to all atmospherically relevant values of pressure, temperature, supersaturation (values of RH > 100 %), and the LWC in warm clouds.The x axis in Fig. A1 spans a range of non-cloudy environmental relative humidity values (60 % < RH < 100 %), and the y axis spans a temperature difference range between the cloud and the environment parcels (−3 • < dT < 3 • ).The initial (B i ) and final (B f ; after entrainment) buoyancy values and the differences between them can be either positive or negative.The regions of B i > 0 (B i < 0) in fact illustrate the effects of entrainment on B core (B margin ) parcels.

Figure 1 .
Figure 1.A schematic representation of a cloud field center-ofgravity height (y axis) vs. mass (x axis) phase space (CvM in short).The majority of clouds are confined to the region between the adiabatic approximation (curved dashed line) and the inversion layer base height (horizontal dashed line).The yellow, magenta, red, and grey shaded regions represent cloud growth, gradual dissipation, cloud fragments which shed large clouds, and cloud fragments which shed precipitating clouds, respectively.The black arrows represent continuous trajectories of cloud growth and dissipation.The hatched arrows represent two possible discontinuous trajectories of cloud dissipation where clouds shed segments.

Figure 2 .
Figure 2. Four vertical cross-sections (at t = 8, 20, 30, and 40 min) during the single cloud simulation.y axis represents height (m), and x axis represents the distance from the axis (m).The black, magenta, green, and yellow lines represent the cloud, W core , RH core , and B core , respectively.The black arrows represent the wind, the background represents the condensation (red) and evaporation rate (blue; g kg −1 s −1 ), and the black asterisks indicate the vertical location of the cloud centroid.Note that in some cases the lines indicating core boundaries overlap (mainly seen for RH and W cores).

Figure 3 .
Figure 3. Temporal evolution of selected core properties, including (a) the fraction of the cores' mass from the total cloud mass (f mass ), (b) the fraction of the cores' volume from the total cloud volume (f vol ), (c) the normalized distance between the cloud centroid and core centroid (D norm ), and (d) the fraction of cores' pixels contained within another core (f pixel ), including all six permutations.See panel legends for descriptions of line colors.

Figure 4 .
Figure 4. CvM phase-space diagrams of B core (a, d), RH core (b, e), and W core (c, f) fractions for all clouds between 3 and 8 h in the BOMEX simulation.Both volume fractions (f vol ; a, b, c) and mass fractions (f mass ; d, e, f) are shown.The red (blue) colors indicate a core fraction above (below) 0.5.The size of each point in the scatter is proportional to the cloud mean area, where the smallest (largest) point corresponds to an area of 0.01 (2.36) km 2 .The percentage of clouds that are core dominated (f vol f mass > 0.5) is included in panel legends.For a general description of CvM space characteristics, the reader is referred to Sect.2.4.

Figure 5 .
Figure5.CvM phase-space diagrams of pixel fractions (f pixel ) of each of the three cores within another core, including six different permutations (as indicated in the panel titles).Bright colors indicate high-pixel fractions (large overlap between two core types), while dark colors indicate low-pixel fractions (little overlap between two core types).Only clouds with a non-zero core fraction (for the core in question) are considered (e.g., for the B core in RH core in panel a, only clouds that contain at least one pixel with positive buoyancy are considered).Scatter point size is proportional to the minimum f vol of the core pairs in question.

Figure 6 .
Figure 6.CvM phase-space diagrams of degree of interchangeability (f int ) for each of the core pairs (as indicated in the panel titles).Bright colors indicate high values (cores can be interchanged with little effect), while dark colors indicate small values (no overlap between cores).Only clouds with a core by at least one definition are considered.Scatter point size is proportional to the minimum f vol of the core pairs in question.Panel legends include percentage of points (out of the scatter) with f int above a certain threshold.

Figure 7 .
Figure 7. CvM phase-space diagrams of distances between core centroid and cloud centroid (D norm,centroid ; a, b, c), and distances between core COG and cloud COG (D norm,COG ; d, e, f) location, for the three different physical core types.The distances are normalized by the maximum distance between the cloud centroid or COG and the cloud perimeter.Bright (dark) colors indicates large (small) distances.Legends include percentage of points (out of the scatter) with D norm below a certain threshold.As seen in Fig. 5, only clouds which contain a core fraction above zero (for the core in question) are considered.

Figure 8 .
Figure 8. Same as Fig. 7 but for only distances between core COG and cloud COG (D norm,COG ).Scatter data are partitioned to clouds with a single core (a, b, c) and multiple cores (d, e, f).The size of each point in the scatter is proportional to the cloud mean area.

Figure 9 .
Figure 9. Mean horizontal profiles of core parameters from the cloud COG to cloud edge, for clouds with single cores and no cores.Data are divided into growing clouds (a, b, c) and dissipating clouds (d, e, f), where the horizontal distances are normalized by the maximum distance to cloud edge.Parameters include buoyancy (a, d), diffusion rate (b, e; taken as a proxy for the supersaturation core), and vertical velocity (c, f).The data are divided to profiles that match core-shell (CS), displaced core-shell (DCS), peripheral core (PC), or no-core (NC) categories, as indicated by the different line colors.The percentages of cloud number (N) and cloud mass (M) attributed to each category are shown in the panel legends.We note that comparing the number percentages with mass percentages for each category gives an indication for the relative sizes of the clouds (e.g., higher N% than M% indicates smaller clouds).

Figure 10 .
Figure 10.Normalized time series of CCE averaged core fractions for the BOMEX (a, b, c), Hawaii (d, e, f), and CASS (g, h, i) simulations.Both core volume fractions (f vol ; a, d, g), normalized distances between cloud and core (D norm ; b, e, h), and pixel fractions of one core within another (f pixel ; c, f, i) are considered.Normalized distance between both COG locations (solid lines) and centroid locations (dotted lines) is shown.Line colors indicated different core types (see legends), while corresponding shaded color regions indicate the standard deviation.Normalized time enables averaging together CCEs with different lifetimes, from formation to dissipation.The number of CCEs averaged together for each simulation is included in the left-column panel titles.

Figure B1 .
Figure B1.Phase space presenting the resultant buoyancy (a) and relative humidity (RH; b) when mixing B core and B margin parcels with equal RH but different temperatures.A mixing fraction of 0.5 is chosen.Both parcels are initially saturated (RH = 100 %) and have an LWC of 0.5 g kg −1 .The environment has a temperature of 15 • C and pressure of 850 mbar.The x(y) axis spans the range of temperature differences between the B core (B margin ) parcel and the environment.