Adaptation to drought involves complex interactions of traits that vary within and among species. To date, few data are available to quantify within-species variation in functional traits and they are rarely integrated into mechanistic models to improve predictions of species response to climate change. We quantified intraspecific variation in functional traits of two Hakea species growing along an aridity gradient in southeastern Australia. Measured traits were later used to parameterise the model SurEau to simulate a transplantation experiment to identify the limits of drought tolerance. Embolism resistance varied between species but not across populations. Instead, populations adjusted to drier conditions via contrasting sets of trait trade-offs that facilitated homeostasis of plant water status. The species from relatively mesic climate, Hakea dactyloides, relied on tight stomatal control whereas the species from xeric climate, Hakea leucoptera dramatically increased Huber value and leaf mass per area, while leaf area index (LAI) and epidermal conductance (g(min)) decreased. With trait variability, SurEau predicts the plasticity of LAI and g(min) buffers the impact of increasing aridity on population persistence. Knowledge of within-species variability in multiple drought tolerance traits will be crucial to accurately predict species distributional limits.