Abstract We present the results of a radial velocity survey of young stellar objects (YSOs) in early stages of evolution in the core of the L1688 molecular cloud. New and archival spectra obtained with four high-resolution infrared spectrographs were analyzed using Markov chain Monte Carlo techniques that simultaneously fit for the radial velocity, T eff, v sin i, and veiling by comparison with synthetic spectra. The radial velocity distribution for 32 objects, most with Class I or flat-spectrum spectral energy distributions, is marginally Gaussian, with a higher dispersion relative to optical surveys at the 2σ level. When comparing the results from both proper-motion and radial velocity surveys in L1688, there is a trend for the 1D dispersions to be higher for samples of Class I/flat-spectrum YSOs that reside in the cloud core compared to Class II/III dominated samples, which are located in the lower extinction periphery. In addition, there is a velocity gradient along the major axis of the cloud core that appears more pronounced than that derived from optically visible objects at the cloud edges. If these higher dispersions for Class I/flat-spectrum objects are confirmed by future surveys, this could imply a supervirial state for the less evolved objects in the cloud core and be a signature of the initial collapse and rebound of the cluster as suggested by recent simulations of cluster evolution.