Abstract The variation and plasticity of leaf morphology plays a pivotal role in the response to environmental changes for plant individuals. Discovering the large-scale pattern of such variation can reveal plants' general adaptive strategies. We analyzed leaf morphology of three widespread woody species in the northern hemisphere using specimen data from the iDigBio and GBIF databases, to investigate the variations in the individual mean traits, in the inter- and intra-individual variability of traits, and in the allometry between traits, along climatic gradients. We found that larger and wider leaves were associated with warmer, wetter and low-sunlight habitats, while smaller but wider leaves are linked to higher wind speed, indicating the response of leaf morphology to multiple climate stresses. The inter-individual variation in leaf area was smaller in colder and windier conditions, suggesting the trait convergence among individuals under environmental filtering, while the intra-individual variation in leaf relative width was smaller in warmer habitats, indicating the similar growth optimum of leaves within one individual in more favorable conditions. Finally, the allometric exponent between leaf length (X-axis) and width (Y-axis) became greater under lower solar radiation and higher wind speed, while the squared correlation coefficient (r2) indicating phenotypic integration showed a decoupling trend under colder conditions, indicating that climate affected the variation tendency of leaf relative width during leaf enlargement. These results reveal the common patterns of leaf morphology responding to climate variation spatially and underscore the necessity to consider inter- and intra-individual variability when examining plant responses to environmental changes.