Abstract
Aim Southern Africa lacks the palaeo‐archives needed to infer regional vegetation history. The megadiverse subtropical thicket of South Africa, locally termed the Albany subtropical thicket ( AST ), is postulated to have experienced severe reductions in range during the Pleistocene glacial periods. Here we test this glacial contraction hypothesis using spatial predictions from climate envelope models to explore the distributions of the AST vegetation subtypes at the Last Glacial Maximum ( LGM ). Location South Africa. Methods We used an ensemble approach to assess the variability of predictions due to the many methodological decisions made during the modelling process. Five sources of uncertainty were included: initial locality data sets, selected climate parameters, statistical methods, LGM global climate models and threshold criteria. Results The last three sources of uncertainty contributed significantly to the variation in projected LGM distributions. Nonetheless, estimates of gain, loss and turnover in range between present and past scenarios indicate dramatic range reduction and turnover during the LGM for all AST subtypes. Arid and valley thickets experienced the most severe reductions with an overall decline in elevation and concomitant fragmentation, whereas mesic thicket had a reduced but fairly continuous range. Main conclusions Our results suggest that the Pleistocene glacial periods saw dramatically reduced distributions of AST subtypes, supporting the glacial contraction hypothesis. Potential LGM refugia were identified and extreme bottlenecks were predicted for species from all AST subtypes during the LGM . Watersheds were also more likely to be strong barriers to gene flow during the LGM due to the retraction of thicket subtypes into basin valleys.