This study analyzes the effects of bioclimate and masticatory factors on the regional variability of human cranial forms across 150 ethnic groups worldwide. Morphometric variables were generated using principal component analysis applied to 3D homologous models. Relationships between cranial form and bioclimate (temperature and precipitation) and masticatory factors (infratemporal space) were tested considering sampling bias due to past population movements during the late Pleistocene and/or early- to mid-Holocene. Cranial size correlated with thermal conditions, consistent with Bergmann's rule. The length/breadth proportion of the neurocranium aligned with Allen's rule for thermal adaptation, while no relationship with masticatory stress was found. Facial form responded to either climate or masticatory conditions, although the primary factor was unclear due to the high correlation between stresses. However, masticatory stress was identified as an equally significant factor behind facial flatness in cold regions, else than the effect of Allen's rule. High narrowness of nasal and orbital openings correlated significantly with cold temperatures and cranial size, suggesting not only functional but also allometric effect. This study demonstrated the complexity of environmental influences on cranial form diversity, nonetheless suggested reduction of selective pressure on cranial form caused by natural environmental stress due to the development of civilization.