This work addresses the formation of adjoining twin pairs (ATPs) at grain boundaries (GBs) in hexagonal close-packed (hcp) metals from the perspective of the co-nucleation (CN) of pairs of deformation twins. A continuum defect mechanics model is proposed to investigate the energetic feasibility of the CN of ATPs resulting from the dissociation of GB dislocations. Application of the model to a set of GBs in Mg reveals that CN is largely preferred over the nucleation of a single twin variant for low misorientation angle GBs, which is consistent with previous experimental analyses. The CN events are subsequently analyzed considering the GB character, and the twin system alignment (m′) and Schmid factors to obtain better insights into the conditions governing the formation of ATPs via CN. The likelihood of CN even at low m′ values reveals that CN events could be responsible for the formation of ATP pairs in these scenarios.