Accurate prediction of genomic breeding values for Timothy was possible using genomic best linear unbiased prediction. Timothy (Phleum pratense L.) is a grass species of great importance for Finnish agricultural production systems. Genotyping-by-sequencing along with genomic prediction methods offer the possibility to develop breeding materials efficiently. In addition, knowledge about the relationships among traits may be used to increase rates of genetic gain. Still, the quality of the genotypes and the validation population may affect the accuracy of predictions. The objectives of the study were (i) to estimate variance components for yield, winter damage and digestibility traits, and (ii) to assess the accuracy of genomic predictions. Variance components were estimated using genomic residual maximum likelihood where the genomic relationship matrix was scaled using a novel approach. Genomic breeding values were estimated using genomic best linear unbiased prediction in single- and multiple-trait settings, and for different marker filtering criteria. Estimates of heritability ranged from 0.13 ± 0.03 to 0.86 ± 0.05 for yield at first cut and organic matter digestibility at second cut, respectively. Genetic correlations ranged from -0.72 ± 0.12 to 0.59 ± 0.04 between yield at first cut and winter damage, and between digestibility at first and second cuts, respectively. Accuracy of prediction was not severely affected by the quality of genotyping. Using family cross-validation and single-trait models, predictive ability ranged from 0.18 to 0.62 for winter damage and digestibility at second cut, respectively. In addition, validation using forward prediction showed that estimated genomic breeding values were moderately accurate with little dispersion. Thus, genomic prediction constitutes a valuable tool for improving Timothy in Finland.