The paper shows that the domain walls formed in heavy QCD axion models can generate a stochastic background of gravitational waves within the reach of current and future gravitational wave detectors (e.g. LIGO-VIRGO-KAGRA, LISA). Remarkably, the signal can be disentangled from other sources due to the prediction of correlated signals at near-future neutron and proton electric dipole moments experiments.