AbstractA series of triaxial compression tests, which has capable of measuring the volumetric strain of the sample, were conducted on columnar ice. A new testing approach of probing the experimental yield surface was performed from a single sample in order to investigate yield and hardening behaviors of the columnar ice under complex stress states. Based on the characteristic of the volumetric strain, a new method of defined the multiaxial yield strengths of the columnar ice is proposed. The experimental yield surface remains elliptical shape in the stress space of effective stress versus mean stress. The effect of temperature, loading rate and loading path in the initial yield surface and deformation properties of the columnar ice were also studied. Subsequent yield surfaces of the columnar ice have been explored by using uniaxial and hydrostatic paths. The evolution of the subsequent yield surface exhibits significant path-dependent characteristics. The multiaxial hardening law of the columnar ice was established experimentally. A phenomenological yield criterion was presented for multiaxial yield and hardening behaviors of the columnar ice. The comparisons between the theoretical and measured results indicate that this current model is capable of giving a reasonable prediction for the multiaxial yield and post-yield properties of the columnar ice subjected to different temperature, loading rate and path conditions.