Near-surface shear-velocity structure can be inferred from multimode dispersion data. Several methods have been developed to isolate the different modes from seismic signals observed on linear arrays of sensors. Most techniques analyze the wavefield through a frequency-wavenumber lass="inline-formula" id="inline-formula-1">lass="math mml" alt="Formula" src="R13/embed/mml-math-1.gif" /> transform, paying little attention to group-delay-time information. Moreover, classical analyses are generally restricted to fundamental-mode dispersion, limiting the resolution power at depth. We have overcome the limitations of classical lass="inline-formula" id="inline-formula-2">lass="math mml" alt="Formula" src="R13/embed/mml-math-2.gif" /> analysis by using a wavefield representation in the group-velocity/phase-velocity lass="inline-formula" id="inline-formula-3">lass="math mml" alt="Formula" src="R13/embed/mml-math-3.gif" /> domain. We have then set up a nonlinear inversion procedure, easily tractable on a common field computer, to constrain the 1D vertical profile of shear velocities. Applications to synthetic data and to a set of actual records show that lass="inline-formula" id="inline-formula-4">lass="math mml" alt="Formula" src="R13/embed/mml-math-4.gif" /> diagrams greatly help to separate dispersion information between different modes, even when they are not detectable on usual lass="inline-formula" id="inline-formula-5">lass="math mml" alt="Formula" src="R13/embed/mml-math-5.gif" /> diagrams. Tests on synthetic and actual data confirm that the inversion procedure quickly converges to the expected model.