Polarization relaxation was studied in Pb(Zn1/3Nb2/3)O3-PbTiO3 (PZN-PT) single crystals that show fatigue anisotropy. To excite prepoled crystals, a modest dc voltage (<1/2 of the coercive field) was applied along the poling direction. Upon removal of the voltage, the polarization decay in the time domain was measured. Experimental data were modeled with a stretched exponential function. Stretching exponent (β〈hkl〉) and characteristic time (τ〈hkl〉) constants for polarization relaxation were determined from data over four decades in the time domain at different stages of bipolar cycling. β〈hkl〉 values after 101 cycles were 0.146±0.002 and 0.247±0.0004 in the 〈001〉 and 〈111〉 orientations, respectively. The β〈111〉 constant increased up to 0.453±0.104 after 105 cycles in 〈111〉 oriented crystals that show fatigue. However, much less change is observed in β〈001〉 as a function of cycling for 〈001〉 crystals. Characteristic time constants for relaxation (τ〈hkl〉) were calculated for 〈001〉 and 〈111〉 orientations as 0.401±0.048 s and 57.46±0.10 s, respectively. These results suggest a faster polarization relaxation in 〈001〉 than in the 〈111〉 orientation of rhombohedral PZN-PT ferroelectric crystals.