The analysis of data from an airborne synthetic aperture radar (SAR) campaign in the percolation zone of Greenland revealed an interferometric coherence undulation behavior with respect to vertical wavenumber, which cannot be explained with existing models. We propose a model extension that accounts for scattering from distinct layers below the surface. Simulations show that the periodicity of the coherence undulation is mainly driven by the vertical distance between dominant subsurface layers, while the amplitude of the undulation is determined by the ratio between scattering from distinct layers and scattering from the firn volume. We use the model to interpret quad-pol SAR data at X-, C-, S-, Land nd P-band. The inferred layer depths match layer detections in ground based radar data and in situ measurements. We conclude that in the percolation zone scattering from subsurface layers has to be taken into account to correctly interpret SAR data and demonstrate the potential to retrieve geophysical information about the vertical subsurface structure.