Changes in phenotypic traits, such as mollusc shells, are indicative of variations in selective pressure along environmental gradients. Recently, increased sea surface temperature (SST) and ocean acidification (OA) due to increased levels of carbon dioxide in the seawater have been described as selective agents that may affect the biological processes underlying shell formation in calcifying marine organisms. The benthic snail Concholepas concholepas (Muricidae) is widely distributed along the Chilean coast, and so is naturally exposed to a strong physical–chemical latitudinal gradient. In this study, based on elliptical Fourier analysis, we assess changes in shell morphology (outlines analysis) in juvenile C. concholepas collected at northern (23°S), central (33°S) and southern (39°S) locations off the Chilean coast. Shell morphology of individuals collected in northern and central regions correspond to extreme morphotypes, which is in agreement with both the observed regional differences in the shell apex outlines, the high reclassification success of individuals (discriminant function analysis) collected in these regions, and the scaling relationship in shell weight variability among regions. However, these extreme morphotypes showed similar patterns of mineralization of calcium carbonate forms (calcite and aragonite). Geographical variability in shell shape of C. concholepas described by discriminant functions was partially explained by environmental variables (pCO2, SST). This suggests the influence of corrosive waters, such as upwelling and freshwaters penetrating into the coastal ocean, upon spatial variation in shell morphology. Changes in the proportion of calcium carbonate forms precipitated by C. concholepas across their shells and its susceptibility to corrosive coastal waters are discussed.