Fully recrystallized commercial Zirconium plates were subjected to uniaxial tension. Tests were conducted at different temperatures (123 K - 623 K) and along two plate directions. Both directions were nominally unfavorable for deformation twinning. The effect of the working temperature on crystallographic texture and in-grain misorientation development was insignificant. However, systematic variation in work hardening and in the area fraction and morphology of deformation twins was observed with temperature. At all temperatures, twinning was associated with significant near boundary mesoscopic shear, suggesting a possible linkage with twin nucleation. A binary tree based model of the polycrystal, which explicitly accounts for grain boundary accommodation and implements the phenomenological extended Voce hardening law, was implemented. This model could capture the measured stress-strain response and twin volume fractions accurately. Interestingly, slip and twin system hardness evolution permitted multiplicative decomposition into temperature-dependent, and accumulated strain-dependent parts. Furthermore, under conditions of relatively limited deformation twinning, the work hardening of the slip and twin systems followed two phenomenological laws proposed in the literature for non-twinning single-phase face centered cubic materials. © 2016 Acta Materialia Inc.