Aims More data are needed about how genetic variation(GV)and envi-ronmental factors influence phenotypic variation within the natural populations of long-lived species with broad geographic distribu-tions.To fill this gap,we examined the correlations among envi-ronmental factors and phenotypic variation within and among 13 natural populations of Pinus tabulaeformis consisting of four demo-graphically distinct groups within the entire distributional range.Methods Using the Akaike’s information Criterion(AiC)model,we measured 12 morphological traits and constructed alternative candidate models for the relationships between each morphological trait and key climatic variables and genetic groups.We then compared the AiC weight for each candidate model to identify the best approximating model for ecogeographical variation of P.tabulaeformis.The partitioning of vari-ance was assessed subsequently by evaluating the independent vari-ables of the selected best models using partial redundancy analysis.Important Findings Significant phenotypic variation of the morphological traits was observed both within individual populations and among populations.Variation partition analyses showed that most of the phenotypic variation was co-determined by both GV and climatic factors.GV accounted for the largest proportion of reproductive trait variation,whereas local key climatic factors(i.e.actual evapotranspiration,AET)accounted for the largest proportion of phenotypic variation in the remaining investigated traits.Our results indicate that both genetic divergence and key environmental factors affect the phenotypic variation observed among populations of this species,and that reproductive and vegetative traits adaptively respond differently with respect to local environmental conditions.This partitioning of factors can inform those making predictions about phenotypic variation in response to future changes in climatic conditions(particularly those affecting AET).