Purpose: To assess radiomics features efficacy obtained by arterial and portal MRI phase in the prediction of clinical outcomes in the colorectal liver metastases patients, evaluating recurrence, mutational status, pathological characteristic (mucinous and tumor budding) and surgical resection margin. Methods: This retrospective analysis was approved by the local Ethical Committee board, and radiological databases were used to select patients with colorectal liver metastases with pathological proof and MRI study in a pre-surgical setting after neoadjuvant chemotherapy. The cohort of patients included a training set (51 patients with 61 years of median age and 121 liver metastases) and an external validation set (30 patients with single lesion with 60 years of median age). For each segmented volume of interest on MRI by two expert radiologists, 851 radiomics features were extracted as median values using the PyRadiomics tool. Non-parametric Kruskal-Wallis test, intraclass correlation, receiver operating characteristic (ROC) analysis, linear regression modelling and pattern recognition methods (support vector machine (SVM), k-nearest neighbors (KNN), artificial neural network (NNET), and decision tree (DT)) were considered. Results: The best predictor to discriminate expansive versus infiltrative tumor growth front was wavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasis extracted on portal phase with accuracy of 82%, sensitivity of 84%, and specificity of 77%. The best predictor to discriminate tumor budding was wavelet_LLH_firstorder_10Percentile extracted on portal phase with accuracy of 92%, a sensitivity of 96%, and a specificity of 81%. The best predictor to differentiate the mucinous type of tumor was the wavelet_LLL_glcm_ClusterTendency extracted on portal phase with accuracy of 88%, a sensitivity of 38%, and a specificity of 100%. The best predictor to identify the recurrence was the wavelet_HLH_ngtdm_Complexity extracted on arterial phase with accuracy of 90%, a sensitivity of 71%, and a specificity of 95%. The best linear regression model was obtained in the identification of mucinous type considering the 13 textural significant metrics extracted by arterial phase (accuracy of 94%, sensitivity of 77% and a specificity of 99%). The best results were obtained in the identification of tumor budding with the eleven textural significant features extracted by arterial phase using a KNN (accuracy of 95%, sensitivity of 84%, and a specificity of 99%). Conclusions: Our results confirmed the capacity of radiomics to identify as biomarkers and several prognostic features that could affect the treatment choice in patients with liver metastases in order to obtain a more personalized approach.