The results show that whilst elastic strains dominate lattice distortion in comparison to GNDs at the loaded state i.e. at the peak load applied, these strains are largely recovered upon unloading and the contribution from GND development subsequently dominates the broadening seen. Two initial crystallographic configurations were considered. In the example in which the crystal c-axis was oriented parallel to the loading direction, the <c+a> pyramidal systems contributed most to slip with <a> basal slip system contributing to a lesser extent. However, in the example where the c-axis was oriented perpendicular to the loading direction, the <c+a> pyramidal and <a> prismatic systems were the more significant contributors to slip and the <a> basal contributing to a lesser extent. However, basal, prismatic and c+a pyramidal slip systems were found to be active in both examples and this was attributed to significant lattice rotation driven by locally high stresses which enabled otherwise badly oriented slip systems to become favourable for slip. Finally, increases in GND density were seen upon unloading for c-axis orientation parallel with loading. This was attributed to the influence of <c+a> pyramidal slip on reverse plasticity leading to diffuse GND density distributions and significant resulting lattice strains compared to that from <a> prism slip loading