The potential expression of crop productivity is a constant challenge for plant breeders in the face of oscillations in environmental variables, making selection difficult and increasing the operational and methodological costs in obtaining the ideal cultivar for a determined region. The objective of the research was to investigate the effects of genotype x environment interaction in the selection of cowpea lines with high grain yield, adaptability, and genotypic stability simultaneously, in environments of the Northwest region of Rio de Janeiro State, using mixed models (Maximum Restricted Likelihood - REML/Best Linear Unbiased Prediction - BLUP). 27 cowpea genotypes were evaluated in six environments by combining location (Cambuci-RJ and Bom Jesus do Itabapoana-RJ) and year (2016, 2017 and 2018). A randomized block design with four replications was performed. The genetic parameters were estimated via the REML/BLUP procedure, and the selection was based on the harmonic mean method of the relative performance of the genetic values (HMRPGV). A meaningful difference was observed for environment and genotype x environment interaction. The residual variance and the variance of the genotype x environment interaction formed the most significant fractions of the phenotypic variance. Lines 4 (Bico-de-ouro 1-5-24), 6 (Pingo-de-Ouro 1-5-4), 5 (Pingo-de-Ouro 1-5-26), and 9 (Pingo-de-Ouro 1-5-8), stood out as superior in terms of stability and adaptability and grain yield by HMRPGV. Commercial cultivars 13 (BRS-Tumucumaque) and 26 (BRS-Itaim) had high grain yield, adaptability, and specific stability in the evaluated edaphoclimatic conditions.