The majority of large-eddy simulation (LES) studies of the planetary boundary layer (PBL) decouple the land surface and the atmosphere by adopting constant surface boundary conditions. To investigate the temporal and spatial interaction between the soil surface and the overlying atmosphere, the NOAH LSM (NCEP/OSU/Air Force/Office of Hydrology Land Surface Model) has been coupled with the NCAR LES. The simulations to be discussed will involve horizontal heterogeneity with scales up to tens of kilometers while retaining sufficient horizontal resolution (~50m) to capture the essential PBL physics. Our investigation focuses on large-scale horizontal heterogeneity in soil moisture and it's influence on wetting and drying PBLs. Large-scale heterogeneity is shown to generate organized motions that alters PBL statistics. As an example, entrainment rates are shown to be dependent on the scale of the heterogeneity.