Abstract. We explore the impacts of the diurnal cycle, free-tropospheric (FT) humidity values, and interactive surface fluxes on the cloud system evolution of non-precipitating marine stratocumuli based on a large ensemble of large-eddy simulations. Cases are separated into three categories based on their degree of decoupling and cloud liquid water path (LWPc). A new budget analysis method is proposed to analyze the evolution of LWPc under both coupled and decoupled conditions. More coupled clouds start with relatively low LWPc and cloud fraction (fc) but experience the least decrease in LWPc and fc during the daytime. More decoupled clouds undergo greater daytime reduction in LWPc and fc, especially those with higher LWPc at sunrise because they suffer from faster weakening of a net radiative cooling. During the nighttime, a positive correlation between FT humidity and LWPc emerges, consistent with higher FT humidity reducing both radiative cooling and the humidity jump, both of which reduce entrainment and increase LWPc. The time rate of change in the LWPc is more likely to be negative for higher LWPc and greater inversion base height (zi), conditions under which entrainment dominates as turbulence develops. In the morning, the rate of the LWPc reduction depends on the LWPc at sunrise, zi, and the degree of decoupling, with distinct contributions from subsidence and radiation. Under well-mixed conditions, it takes about 10 h for the surface fluxes to offset 15 % of the changes in entrainment warming and drying, assuming no changes in transfer coefficients or surface wind speed.