Radiation damage of natural and synthetic zircon grains is evaluated by Raman spectroscopy to understand annealing and stability of fission tracks. Analyses focus on a suite of 338 Paleozoic detrital zircon grains from metamorphosed strata in the Hellenic forearc that were variably annealed by a Miocene thermal event, as well as a suite of 97 synthetic zircon grains. The Raman wavenumber shift of ν3 [SiO4] reveals that radiation damage and damage distribution in this suite mainly depends on uranium concentration. In zircon with similar uranium concentration, the Raman wavenumber shift allows for the determination of radiation damage in different crystals, which is a function of effective accumulation time. Nine detrital zircons grains with moderate radiation damage were stepwise annealed at 1000 and 1400 °C, which resulted in progressive removal of radiation damage revealed in an increase of ν3 [SiO4] peak positions. For a partly reset sample that was brought to temperatures of ~350 °C in a geologic setting (Hellenic forearc), we use the Raman measurements and uranium determination to estimate a Zircon Damage Discrimination Factor (ZRDD), which is our attempt to estimate only radiation damage in single grains by accounting for affects of the uranium atom in the Raman wavenumber. This discrimination allows for a separation of zircon fission track (ZFT) ages of single ages based on grains that have a low track retention (high damage, fully reset grain), thus refining the age determination of cooling in a rock that shows variable resetting.