We examine the annual, seasonal, monthly, and diurnal climate responses to the land use change (LUC) in eastern United States and Cuba during four epochs (1650, 1850, 1920, and 1992) with ensemble simulations conducted with the RegCM4 regional climate model that includes the Biosphere Atmosphere Transfer Scheme (BATS1e) surface physics package (Dickinson et al., 1993). We derived the land use (LU) data sets by harmonizing a previous reconstruction (Steyaert & Knox, 2008) with updated observations and modeled potential vegetation. The eight‐member ensembles for each epoch were driven with randomly perturbed 1990–2002 atmospheric boundary conditions derived from the National Center for Environmental Prediction global reanalysis. LUC induces statistically significant climate responses across all epochs; the largest changes occur between 1850 and 1920 with the widespread conversion of forests in the United States and forests, grassland, and woody wetlands in Cuba to agriculture. The atmospheric feedback from the aggregated grid‐cell responses attributed to physical and biophysical parameters in BATS1e alters the circulation in the lower atmosphere, thereby propagating the LUC regionally. Depending on the season and location, the altered circulation reinforces, attenuates, or has little effect on surface responses. Relative to pre‐settlement (1650), the 1992 LU produces colder mean annual air temperature (−0.09 ± 0.16 °C) and increased precipitation (0.08 ± 0.09 mm day⁻¹) over the United States, warmer (0.08 °C) and wetter (0.03 mm day⁻¹) conditions over Florida, and warming (0.32 °C) and drying (−0.03 mm day⁻¹) over Cuba, indicating that LUC has played a varying role in climate change over the region.