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(2). 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.