We propose an innovative workflow based on the complementary use of Rayleigh waves alongside standard P-wave refraction tomography, which better depicts the shallow part of the near-surface P-wave velocity model. Our surface-wave processing sequence led to an S-wave near-surface velocity field that can be used as a constraint for P-wave tomography and can improve P-wave statics determination. Rayleigh waves are processed in three steps. The first step consists of an accurate frequency-dependent traveltime measurement for each selected source-receiver pair in which the phase difference between two adjacent traces is used to derive the phase velocity. Then, a frequency-dependent surface-wave velocity tomography is performed from the picked traveltimes. Finally, after surface-wave tomography, the frequency-dependent phase velocity volume output by the tomography is inverted to deliver an S-wave near-surface velocity model. This model is used to constrain the first-arrival P-wave tomography. To illustrate our method, we use a 3D narrow-azimuth land data set, acquired along a river valley. Strong lateral velocity variations exist in the shallow part, with slow velocities around the unconsolidated sediments of the riverbed and faster velocities in the consolidated sediments of the surrounding hills. A combined first-arrival tomography using the S-wave velocity model, the initial unconstrained refracted P-wave velocity model, and the original first arrivals is used to obtain a more accurate near-surface P-wave velocity model. This new approach led to a constrained P-wave velocity model from which primary statics are derived and then applied, leading to an improved image with better focusing and continuity of thin layers in the shallowest part.