Résumé : The physical properties of eight ice cores recovered from along a flowline at Glacier de tranfleuron. Switzerland, have led to the identification of three distinctive internal zones. We use variations in the bulk ionic chemistry of these zones to approximate their relative liquid-water concentrations and ice viscosities. Results suggest that relative bulk water concentration and ice softness vary by over an order of magnitude between the zones. Implications of this variability for predictions of the glacier's response to climate change are evaluated by incorporating these relative softnesses into a multi-layered (two-dimensional) model of ice flow. Model output is compared with that from an identical model constrained with a spatially uniform ice viscosity under advance and retreat modelling scenarios. The former scenario is used to tune viscosity by growing a glacier to its present long-section geometry, resulting in best-fit ice hardness values of 1.2 a 1bar 3 for the englacial ice in the multi-layered model and 7.0a 1bar 3 for all of the ice in the single-layered model. Both result in close approximations to the current long profile, yielding rms deviations between measured and modelled ice thicknesses that are < 5 m. In contrast, a single-layered model constrained with a hardness of 1.2a -1bar 3 overestimates the current measured long-section area by 31%, having a rms ice-thickness error of 15.0 m. Under the retreat modelling scenario, which gauges the response of the glacier to an imposed 75 m rise in equilibrium-line altitude (ELA), the multi-layered model predicts a long-section area reduction of 78%. This contrasts with a reduction of 6.4% for the single-layered model (hardness = 7.0 a 1bar 3) and 85% for the single-layered model (hardness = 1.2 a 1bar 3). © International Glaciological Society.