A regional shear-wave velocity (V_S) model has been developed for the Groningen gas field in the Netherlands as the basis for seismic microzonation of an area of more than 1000 km². The V_S model, extending to a depth of almost 1 km, is an essential input to the modelling of hazard and risk due to induced earthquakes in the region. The detailed V_S profiles are constructed from a novel combination of three data sets covering different, partially overlapping depth ranges. The uppermost 50 m of the V_S profiles are obtained from a high-resolution geological model with representative V_S values assigned to the sediments. Field measurements of V_S were used to derive representative V_S values for the different types of sediments. The profiles from 50 to 120 m are obtained from inversion of surface waves recorded (as noise) during deep seismic reflection profiling of the gas reservoir. The deepest part of the profiles is obtained from sonic logging and V_P–V_S relationships based on measurements in deep boreholes. Criteria were established for the splicing of the three portions to generate continuous models over the entire depth range for use in site response calculations, for which an elastic half-space is assumed to exist below a clear stratigraphic boundary and impedance contrast encountered at about 800 m depth. In order to facilitate fully probabilistic site response analyses, a scheme for the randomisation of the V_S profiles is implemented.