AbstractThe U.S. Department of Energy Atmospheric Radiation Measurement (ARM) program's Raman lidar at the ARM Southern Great Plains site in north central Oklahoma has collected water vapor mixing ratio (q) profile data more than 90% of the time since October 2004. Three hundred (300) cases were identified where the convective boundary layer was quasi‐stationary and well mixed for a 2 h period, and q mean, variance, third‐order moment, and skewness profiles were derived from the 10 s, 75 m resolution data. These cases span the entire calendar year, and demonstrate that the q variance profiles at the mixed layer (ML) top changes seasonally and is strongly related to the gradient of q across the interfacial layer. The q variance at the top of the ML shows only weak correlations (r < 0.3) with sensible heat flux, Deardorff convective velocity scale, and turbulence kinetic energy measured at the surface. The median q skewness profile is most negative at 0.85 zi, zero at approximately zi, and positive above zi, where zi is the depth of the convective ML. The spread in the q skewness profiles is smallest between 0.95 zi and zi. The q skewness at altitudes between 0.6 zi and 1.2 zi is correlated with the magnitude of the q variance at zi, with increasingly negative values of skewness observed lower down in the ML as the variance at zi increases, suggesting that in cases with larger variance at zi there is deeper penetration of the warm, dry free tropospheric air into the ML.