We present spatially extended gas kinematics at parsec-scale resolution for the nuclear regions of four nearby disk galaxies, and model them as rotation of a gas disk in the joint potential of the stellar bulge and a putative central black hole. The targets were selected from a larger set of long-slit spectra obtained with the Hubble Space Telescope as part of the Survey of Nearby Nuclei with STIS (SUNNS). They represents the 4 galaxies (of 24) that display symmetric gas velocity curves consistent with a rotating disk. We derive the stellar mass distribution from the STIS acquisition images adopting the stellar mass-to-light ratio normalized so as to match ground-based velocity dispersion measurements over a large aperture. Subsequently, we constrain the mass of a putative black hole by matching the gas rotation curve, following two distinct approaches. In the most general case we explore all the possible disk orientations, alternatively we constrain the gas disk orientation from the dust-lane morphology at similar radii. In the latter case the kinematic data indicate the presence of a central black hole for three of the four objects, with masses of 10^7 - 10^8 solar masses, representing up to 0.025 % of the host bulge mass. For one object (NGC2787) the kinematic data alone provide clear evidence for the presence of a central black hole even without external constraints on the disk orientation. These results illustrate directly the need to determine black-hole masses by differing methods for a large number of objects, demonstrate that the variance in black hole/bulge mass is much larger than previously claimed, and reinforce the recent finding that the black-hole mass is tightly correlated with the bulge stellar velocity dispersion.

26 pages, 11 Postscript figures, accepted for publication on ApJ