Temporal relaxation of density fluctuations in supercooled liquids near the glass transition occurs in multiple steps. Using molecular dynamics simulations for three model glass-forming liquids, we show that the short-time beta relaxation is cooperative in nature. Using finite-size scaling analysis, we extract a growing length scale associated with beta relaxation from the observed dependence of the beta relaxation time on the system size. We find, in qualitative agreement with the prediction of the inhomogeneous mode coupling theory, that the temperature dependence of this length scale is the same as that of the length scale that describes the spatial heterogeneity of local dynamics in the long-time alpha-relaxation regime.