4D flow MRI-derived variables from Marfan patients are highly heterogeneous. Our aim was to identify distinct Marfan patient subgroups based on aortic 4D flow MRI and Z-score for stratification of distinct hemodynamic profiles and clinical features by means of hierarchical cluster analysis. One hundred Marfan patients underwent baseline aortic 4D flow MRI at 3 T. Z-scores, degree of helical and vortical flow, wall shear stress, flow displacement, and peak velocity were determined in the ascending aorta. Sex, age, BMI, antihypertensive medication, and dural ectasia were recorded. Hierarchical cluster analysis was performed using 4D flow MRI variables and Z-scores as input. Cluster analysis resulted in three distinct clusters characterized by different Z-scores (mean ± SD); cluster 1: 0.4 ± 1.1 vs. cluster 2: 3.1 ± 1.1 vs. cluster 3: 3.6 ± 1.9. The three clusters delivered differences in helical and vortical flow patterns (global p = 0.003 and p < 0.001, respectively), wall shear stress (0.49 ± 0.11 vs. 0.44 ± 0.12 vs. 0.37 ± 0.09 N/m2, global p < 0.001), flow displacement (0.11 ± 0.05 vs. 0.16 ± 0.08 vs. 0.15 ± 0.07, global p = 0.006), and peak velocity (76.3 ± 9.0 vs. 60.1 ± 7.3 vs. 56.0 ± 7.8 cm/s, global p < 0.001). Patients in cluster 1 and 2 were relevantly younger than in cluster 3 (32.3 ± 13.8 vs. 32.8 ± 12.6 vs. 40.2 ± 15.0 years, all pairwise ∆p < 0.0297). Hierarchical cluster analysis based on aortic 4D flow MRI and Z-score revealed three distinct subgroups of Marfan patients, each characterized by specific hemodynamic profiles and clinical features. Follow-up of our patients is warranted to assess if 4D flow MRI- and Z-score-based stratification can predict future aortic diameter growth and ultimately improve outcomes. A combination of Z-score and 4D flow MRI-derived parameters may help identify subgroups of Marfan patients representing different stages or phenotypes of aortic disease, which require specific management strategies. Four-dimensional (4D) flow MRI-derived variables of Marfan patients are highly heterogeneous across varying Z-scores. Cluster analysis based on 4D flow MRI and Z-score revealed three distinct subgroups of Marfan patients. A combination of Z-score and 4D flow MRI-derived parameters may identify different stages of aortic disease in Marfan patients.