Near-field radiative heat transfer (NFRHT) between rough surfaces, due to its widespread presence in engineering practice of near-field energy utilization, requires indepth studies, especially from the perspective of physical mechanism. In this paper, an effective multilayer model is built to approach the NFRHT between random rough surfaces of silicon carbide (SiC). Using the effective medium theory (EMT), the effective dielectric function of each layer is obtained, which forms a gradient distribution of dielectric function along the depth of the medium. The influence of the effective dielectric function on surface phonon polaritons (SPhPs) is analyzed, showing that the effective layers with small filling fraction of SiC feature lower SPhP resonance frequencies than SiC bulk. The coupling of SPhPs from the gradient distribution of dielectric function produces new surface modes that dominates the NFRHT. Investigation on the effect of root mean square height (RMS height, {sigma}) reveals that the peaks of local density of states (LDOS) and spectral heat flux are red-shifted as {sigma} increases, while the spectral heat flux below the peak frequency gets larger. This can be attributed to the coupling of SPhPs inside the rough layer. We also found the total net heat flux between rough surfaces separated by an average distance exceeds that between smooth plates and increases with increasing {sigma}, which offer a new way to enhance NFRHT. This work provides a reference for the simulation and understanding of the NFRHT between rough surfaces.