Single crystals of ice subjected to primary creep in torsion exhibit a softening behavior: the plastic strain rate increases with time. In a cylindrical sample, the size of the radius affects this response. The smaller the radius of the sample becomes while keeping constant the average shear stress across a section, the softer the response. The size-dependent behavior is interpreted by using a field dislocation theory, in terms of the coupled dynamics of excess screw dislocations gliding in basal planes and statistical dislocations developed through cross slip occurring in prismatic planes. The differences in the results caused by sample height effects and variations in the initial dislocation microstructure are discussed.