We investigate the consequences of a nonzero bulk viscosity coefficient on the transverse momentum spectra, azimuthal momentum anisotropy, and multiplicity of charged hadrons produced in heavy ion collisions at LHC energies. The agreement between a realistic 3D hybrid simulation and the experimentally measured data considerably improves with the addition of a bulk viscosity coefficient for strongly interacting matter. This paves the way for an eventual quantitative determination of several QCD transport coefficients from the experimental heavy ion and hadron-nucleus collision programs.