Unresolved sub-grid processes, those which are too small or dissipate too quickly to be captured within a model's spatial resolution, are not adequately parameterized by conventional numerical climate models. Sub-grid heterogeneity is lost in parameterizations that quantify only the ‘bulk effect’ of sub-grid dynamics on the resolved scales. A unique solution, one unreliant on increased grid resolution, is the employment of stochastic parameterization of the sub-grid to reintroduce variability. We administer this approach in a coupled land-atmosphere model, one that combines the single-column Community Atmosphere Model (CAM-SC) and the single-point Community Land Model (CLM-SP), by incorporating a stochastic representation of sub-grid latent heat flux to force the distribution of precipitation. Sub-grid differences in surface latent heat flux arise from the mosaic of Plant Functional Types (PFT) that describe terrestrial land cover. With the introduction of a stochastic parameterization framework to affect the distribution of sub-grid PFT's, we alter the distribution of convective precipitation over regions with high PFT variability. The stochastically forced precipitation probability density functions (pdf) show lengthened tails, demonstrating the retrieval of rare events. Through model data analysis we show that the stochastic model increases both the frequency and intensity of rare events in comparison to conventional deterministic parameterization.