We numerically determine subleading scaling terms in the ground-state entanglement entropy of several two-dimensional (2D) gapless systems, including a Heisenberg model with NÃ©el order, a free Dirac fermion in the Ï-flux phase, and the nearest-neighbor resonating-valence-bond wave function. For these models, we show that the entanglement entropy between cylindrical regions of length x and Lâx, extending around a torus of length L, depends on the dimensionless ratio x/L. This can be well approximated on finite-size lattices by a function ln(sin(Ïx/L)), akin to the familiar chord-length dependence in one dimension. We provide evidence, however, that the precise form of this bulk-dependent contribution is a more general function in the 2D thermodynamic limit.