Ability to visualize breast lesion vascularity and quantify the vascular heterogeneity using contrast-enhanced 3-D harmonic (HI) and subharmonic (SHI) ultrasound imaging was investigated in a clinical population. Patients (n = 134) identified with breast lesions on mammography were scanned using power Doppler imaging, contrast-enhanced 3-D HI, and 3-D SHI on a modified Logiq 9 scanner (GE Healthcare). A region of interest corresponding to ultrasound contrast agent flow was identified in 4D View (GE Medical Systems) and mapped to raw slice data to generate a map of time-intensity curves for the lesion volume. Time points corresponding to baseline, peak intensity, and washout of ultrasound contrast agent were identified and used to generate and compare vascular heterogeneity plots for malignant and benign lesions. Vascularity was observed with power Doppler imaging in 84 lesions (63 benign and 21 malignant). The 3-D HI showed flow in 8 lesions (5 benign and 3 malignant), whereas 3-D SHI visualized flow in 68 lesions (49 benign and 19 malignant). Analysis of vascular heterogeneity in the 3-D SHI volumes found benign lesions having a significant difference in vascularity between central and peripheral sections (1.71 ± 0.96 vs. 1.13 ± 0.79 dB, p < 0.001, respectively), whereas malignant lesions showed no difference (1.66 ± 1.39 vs. 1.24 ± 1.14 dB, p = 0.24), indicative of more vascular coverage. These preliminary results suggest quantitative evaluation of vascular heterogeneity in breast lesions using contrast-enhanced 3-D SHI is feasible and able to detect variations in vascularity between central and peripheral sections for benign and malignant lesions.