This paper focuses on the variability in entrainment rate in individual cumulus clouds using the entrainment rate estimated on the scale of 5m in 186 shallow cumulus clouds from eight aircraft flights, using in situ observations from the RACORO field campaign (the routine atmospheric radiation measurement aerial facility clouds with low optical water depths optical radiative observations) over the atmospheric radiation measurement Southern Great Plains site, USA. The result shows that the mean entrainment rate of all the 186 clouds systematically decreases from the cloud edge to the cloud center. Further analysis of the fluctuation of entrainment rate shows that the probability density function of entrainment rate in each flight can be fitted by the lognormal, gamma, or Weibull distributions virtually equally well, with the Weibull distribution being the best. The parameter “standard deviation” in the lognormal distribution is weakly negatively correlated, and the other parameters in the three distributions are positively correlated with relative humidity in the entrained dry air and dilution effect, respectively. Entrainment rate is negatively correlated with droplet concentration, droplet size, and liquid water content, but positively correlated with relative dispersion. The effect of entrainment rate on the spectral shape of cloud droplet size distribution is examined and linked to the systems theory on the cloud droplet size distribution.