Conventional and organic production systems mainly differ in feeding strategies, outdoor and pasture access, and the use of antibiotic treatments. These environmental differences could lead to a genotype by environment interaction (G × E) and a requirement for including G × E in breeding decisions. The objectives of this study were to estimate variance components and heritabilities for conventional and organic production systems and investigate G × E under these 2 production systems for female fertility traits in Danish Holsteins. The analyzed traits included the interval from calving to first insemination (ICF), the interval from first to last insemination, number of inseminations per conception (NINS), and non-return rate within 56 d after the first insemination. Records of female fertility in heifers and the first 3 lactations in cows as well as grass ratio of feed at herd level were collected during the period from 2011 to 2016. The performances of a trait in heifers and cows (lactation 1 to 3) were considered as different traits. The (co)variance components and the resulting heritabilities and genetic correlations were estimated using 2 models. One was a bivariate model treating performances of a trait under organic and conventional production systems as 2 different traits using a reduced data set, and the other was a reaction norm model with random regression on the production system and the grass ratio of feed using a full data set. The full data set comprised records of 37,836 females from 112 organic herds and 513,599 females from 1,224 conventional herds, whereas the reduced data set comprised records from all these 112 organic herds and 92,696 females from 185 convention herds extracted from the full data set with grass ratio of feed lower than 0.20. All female fertility performances of the organic production system were superior to those of the conventional production system. Besides, heterogeneities in additive genetic variances and heritabilities were observed between conventional and organic production systems for all traits. Furthermore, genetic correlations between these 2 production systems ranged from 0.607 to 1.000 estimated from bivariate models and from 0.848 to 0.999 estimated from reaction norm models. Statistically significant G × E were observed for NINS in heifers, non-return rate within 56 d after the first insemination in heifers, and ICF from the bivariate model, and for ICF and NINS in cows from the reaction norm model.