The objectives of this study were to investigate changes in genetic parameters for milk yield (MY) and heat tolerance of the crossbred Thai Holstein Friesian population under different heat stress levels over time, and to investigate the threshold point of heat stress manifestation on milk production. Genetic parameters were estimated using single-step genomic REML (ssGREML) and traditional REML models. Data included 58,965 test-day MY records from 1999 to 2008 (old data) and 105,485 test-day MY records from 2009 to 2018 (recent data) from the first parity of 24,520 cows. The pedigree included 55,168 animals, of which 882 animals had genotypes. Variance components were estimated with the REMLF90 program using a repeatability model with random regressions on a function of temperature-humidity index (THI) for additive genetic and permanent environmental effects. Fixed effects included farm-calving season combination, breed group-months in milk combination, and age at first calving. Random effects included additive genetic (intercept and slope) effects, permanent environmental (intercept and slope) effects, and herd-month-year of test. The phenotypic mean for MY was 13.33 ± 4.39 kg/d in the old data, and 14.48 ± 4.40 kg/d in the recent data. Estimates over different THI levels for the intercept additive genetic variance using old data ranged from 2.61 to 2.77 and from 5.02 to 5.38 using recent data with the REML method. In ssGREML analyses (performed with recent data only) the estimates for the intercept additive genetic variance ranged from 4.71 to 5.05. Estimates for the slope additive genetic variance were close to zero in all cases, with the largest values (0.024-0.030) at the most extreme THI value (80). Using REML, the covariance between the intercept and the slope additive genetic effects (THI from 72 to 80) ranged from -0.001 to 0.019 with old data and from 0.027 to 0.060 with recent data. The same covariance ranged from 0.026 to 0.057 in ssGREML analyses. The covariance between the intercept and the slope permanent environmental effects ranged from -0.42 to -0.67 for all data and THI levels. Across THI levels, the genetic correlation between MY and heat tolerance varied from -0.06 to 0.13 with old data, from 0.16 to 0.30 with recent data in REML analyses, and from 0.15 to 0.30 in ssGREML analyses, suggesting that in the current population the top animals for MY are more resistant to heat stress. This was expected, because of the introduction of Bos indicus genes in the last years. Heritability estimates for MY ranged from 0.19 to 0.21 (old data) and from 0.33 to 0.40 (recent data) for REML analyses. Heritability estimates for MY using ssGREML ranged from 0.31 to 0.38. A decline in MY was found when the animals' breed composition had more than 97.3% of Holstein genetics, and it was greatest at THI 80. The heritability and genetic correlations observed in this study show that selection for MY is possible without a negative correlated response for heat tolerance. Although the inclusion of genomic information is expected to increase the accuracy of selection, more genotypes must be collected for successful application. Future research should address other production and fitness traits within the Thai Holstein population.