Aim The environmental preferences of species are an important facet of their response to changing conditions, and these have long been thought to exhibit phylogenetic conservatism. However, these bioclimatic envelopes have not previously been imputed from climate records at the date and location of occurrence, and the strength of their phylogenetic signal has not been studied at a broad scale. Here, we combine records from global climate reconstructions with contemporaneous plant occurrences for all available terrestrial plant species and test for phylogenetic niche conservatism in plant climatic traits. Location Global. Time period 1901-2018. Major taxa studied Terrestrial plants. Methods We used >100 million plant records from the Global Biodiversity Information Facility (GBIF) to produce distributions of bioclimatic envelopes for >200,000 species, using a range of climate variables. We matched species observations to historical climate reconstructions from the European Centre for Medium-Range Weather Forecasting (ECMWF) and compared this with WorldClim climate averages. We tested for phylogenetic signal in a supertree of plants using Pagel's λ. Finally, to investigate how well bioclimatic envelopes could be inferred for poorly known and rare species, we performed cross-validation by removing occurrence records for some common species to test how accurately their bioclimatic envelopes were estimated. Results We found extremely strong phylogenetic signals (λ > 0.9 in some cases) for climate variables from both climate datasets, including temperature, soil temperature, solar radiation and precipitation. We were also able to impute missing bioclimatic envelopes for artificially removed species, having a correlation with observed data of .7. Main conclusions We reconstructed plant climatic tolerances for >200,000 plant species historically recorded on GBIF using a technique that could be applied to any comparable biodiversity dataset. Although global information on most species is sparse, we explored methods for bias correction and data imputation, with positive results for both.