Abstract Objective 24-hour ambulatory blood pressure measurements (ABPM) is the best out-of-office BP measurement used to refine cardiovascular (CV) risk stratification. Currently, only summary ABPM indexes are used in the clinical practice while the temporal characteristics of the raw BP recordings have not been explored. Therefore, in this study we employed an unsupervised machine learning (ML) algorithm suitable for multivariate time-series analysis to identify distinct BP and heart rate (HR) profiles associated with the incidence of adverse CV events. Design and Method In 1391 community-dwelling individuals (mean age, 46.9 years; 50.97% women), we acquired 24-hour ABPM and clinical data and collected CV events (n=399) on average 22 years later. We employed multivariate dynamic time warping and k-medoids algorithm on the raw recordings of systolic and diastolic BP and HR. To assess the clinical significance of the derived clusters, we compared the clinical characteristics and the incidence of adverse events. We validated the trained model using 24-hour ABPM obtained from external population cohort (n=1137). Results The silhouette score and Dunn index showed that the optimal number of clusters was 4. Across all clusters we observed significant difference in age and sex distribution. Cluster 4 had the worst CV risk factors profiling, with higher office BP and anti-hypertensive medication intake compared to the other clusters. We observed differences in the incidence of adverse events between clusters (Figure), with cluster 1 representing the lowest risk group and cluster 4 the highest. After adjusting for traditional CV risk factors, including office systolic BP, cluster 1 had the lowest risk (HR:0.92; 95%CI:0.84-0.99; P=0.036) and cluster 4 the highest (HR:1.13; 95%CI:1.03-1.24; P=0.006) as compared to the average population risk. The analysis of the external validation cohort revealed similar ABPM patterns and clinical characteristics. Conclusion Employing an unsupervised ML model on the raw ABPM recordings we identified clinically meaningful clusters associated with an increasing cumulative incidence of CV events. This ML analysis paves the way towards utilization of temporal BP and HR information and subsequently optimization of CV risk stratification.