Rotational spectra of small-sized (H(2))(N)-OCS clusters with N = 2-7 were measured using a pulsed-jet Fourier transform microwave spectrometer. These include spectra of pure (para-H(2))(N)-OCS clusters, pure (ortho-H(2))(N)-OCS clusters, and mixed ortho-H(2) and para-H(2) containing clusters. The rotational lines of ortho-H(2) molecules containing clusters show proton spin-proton spin hyperfine structure, and the pattern evolves as the number of ortho-H(2) molecules in the cluster increases. Various isotopologues of the clusters were investigated, including those with O(13)CS, OC(33)S, OC(34)S, and O(13)C(34)S. Nuclear quadrupole hyperfine structures of rotational transitions were observed for (33)S (nuclear spin quantum number I = 3/2) containing isotopologues. The (33)S nuclear quadrupole coupling constants are compared to the corresponding constant of the OCS monomer and those of the He(N)-OCS clusters. The assignment of the number of solvating hydrogen molecules N is supported by the analyses of the proton spin-proton spin hyperfine structures of the mixed clusters, the dependence of line intensities on sample conditions (pressure and concentrations), and the agreement of the (para-H(2))(N)-OCS and (ortho-H(2))(N)-OCS rotational constants with those from a previous infrared study [J. Tang and A. R. W. McKellar, J. Chem. Phys. 121, 3087 (2004)].