Sea ice thickness is an essential climate variable. Current L-Band sea ice thickness retrieval methods do not account for sea ice surface roughness that is hypothesised to be not relevant to the process. This study attempts to validate this hypothesis that has not been tested yet. To test this hypothesis, we created a physical model of sea ice roughness based on geometrical optics and merged it into the L-band emissivity model of sea ice that is similar to the one used in the operational sea ice thickness retrieval algorithm. The facet description of sea ice surface used in geometrical optics is derived from 2-D surface elevation measurements. Subsequently the new model was tested with TB measurements performed during the SMOSice 2014 field campaign. Our simulation results corroborate the hypothesis that sea ice surface roughness has a marginal impact on near-nadir TB (used in the current operational retrieval). We demonstrate that the probability distribution function of surface slopes can be approximated with a parametric function whose single parameter can be used to characterise the degree of roughness. Facet azimuth orientation is isotropic at scales greater than 4.3 km. The simulation results indicate that surface roughness is a minor factor in modelling the sea ice brightness temperature. The change in TB is most pronounced at incidence angles greater than 40° and can reach up to 8 K for vertical polarisation at 60°. Therefore current and future L-band missions (SMOS, SMAP, CIMR, SMOS-HR) measuring at such angles can be affected. Comparison of the brightness temperature simulations with the SMOSice 2014 radiometer data does not yield definite results.