Reflection and transmission of electromagnetic waves at the boundaries of periodic composites (electromagnetic/optical metamaterials) depends in general on both bulk and surface waves. We investigate the interplay of these two contributions using three-dimensional full-wave numerical simulations and a recently developed non-asymptotic homogenization theory.