The effect of genotype by environment interaction is one of the limiting factor for successfully identifying superior genotypes in plant breeding programs. The study aimed to identify celery genotypes with high and stable yields across different environments by assessing genotype-environment (GE) interactions, stability indices, and associations between agronomic parameters. Twenty celery genotypes were sown in randomized complete block design field trials at three locations in Punjab, India, viz., Ludhiana, Amritsar, and Gurdaspur during 2018–19, 2019–20, and 2020–21. Phenotypic data for thirteen agronomic parameters were collected for GE interactions study. The investigation mainly focused on the significance of GGE biplots in explaining GE interactions based on seed yield-related traits. Results revealed significant differences in genotypes, locations, years, and interactions in the studied traits. The WAASB (weighted average of absolute sScores) stability index estimation was highlighted to identify superior genotypes with maximum seed yield and high stability. The investigation specified the suitability of MTSI (multi-trait stability index), which considers multiple traits through genotype-ideotype distance, for assessing stability. Seed yield showed a positive and significant correlation with primary branches and umbels per plant that would be effective in simultaneous selection for multiple traits. Based on MTSI, three genotypes, PAU1, PAU9, and PAU10, were identified as stable and high-yielding across all three tested locations. The investigation suggests that efforts should focus on identifying location-specific genotypes using multi-location and multi-year data for varietal releases rather than emphasizing the overall performance of genotypes.