Ultrasound-based blood flow monitoring is vital in the diagnosis and treatment of a variety of cardiovascular and neurologic conditions. Finite spatial resolution of clinical color flow (CF) systems, however, has hampered measurement of vessel cross-section areas. We propose a resolution enhancement technique that allows reliable determination of blood flow in small vessels. We leverage sparsity in the spatial distribution of the frequency spectrum of routinely collected CF data to blindly determine the point spread function (PSF) of the imaging system in a robust manner. The CF data are then deconvolved with the PSF, and the volumetric flow is computed using the resulting velocity profiles. Data were collected from phantom blood vessels with diameters between 2 and 6 mm using a clinical ultrasound system at 2 MHz insonation frequency. The proposed method yielded a flow estimation bias of 0 mL/min, standard deviation of error (SDE) of 22 mL/min, and a root-mean-square error (RMSE) of 22 mL/min over a 150 mL/min range of mean flows. Recordings were also obtained in low SNR conditions using a skull mimicking element, resulting in an estimation bias of −13 mL/min, SDE of 23 mL/min, and an RMSE of 26 mL/min. The effect of insonation frequency was also investigated by obtaining recordings at 4.3 MHz, yielding an estimation bias of −16 mL/min, SDE of 16 mL/min, and an RMSE of 22 mL/min. The results indicate that our technique can lead to clinically acceptable flow measurements across a range of vessel diameters in high and low SNR regimes.