Background The lumbar vertebra and paraspinal muscles play an important role in maintaining the stability of the lumbar spine. Therefore, the aim of this study was to investigate the relationship between paraspinal muscles fat infiltration and vertebral body related changes [vertebral bone quality (VBQ) score and Modic changes (MCs)] in patients with chronic low back pain (CLBP). Methods Patients with CLBP were prospectively collected in four hospitals and all patients underwent 3.0T magnetic resonance scanning. Basic clinical information was collected, including age, sex, course of disease (COD), and body mass index (BMI). MCs were divided into 3 types based on their signal intensity on T1 and T2-weighted imaging. VBQ was obtained by midsagittal T1-weighted imaging (T1WI) and calculated using the formula: SIL1-4/SICSF. The Proton density fat fraction (PDFF) values and cross-sectional area (CSA) of paraspinal muscles were measured on the fat fraction map from the iterative decomposition of water and fat with the echo asymmetry and least-squares estimation quantitation (IDEAL-IQ) sequences and in/out phase images at the central level of the L4/5 and L5/S1 discs. Results This study included 476 patients with CLBP, including 189 males and 287 females. 69% had no Modic changes and 31% had Modic changes. There was no difference in CSA and PDFF for multifidus(MF) and erector spinae (ES) at both levels between Modic type I and type II, all P values>0.05. Spearman correlation analysis showed that VBQ was weakly negatively correlated with paraspinal muscles CSA (all r values < 0.3 and all p values < 0.05), moderately positive correlation with PDFF of MF at L4/5 level (r values = 0.304, p values<0.001) and weakly positively correlated with PDFF of other muscles (all r values<0.3 and all p values<0.001). Multivariate linear regression analysis showed that age (beta = 0.141, p < 0.001), gender (beta = 4.285, p < 0.001) and VBQ (beta = 1.310, p = 0.001) were related to the total PDFF of muscles. For MCs, binary logistic regression showed that the odds ratio values of age, BMI and COD were 1.092, 1.082 and 1.004, respectively (all p values < 0.05). Conclusions PDFF of paraspinal muscles was not associated with Modic classification. In addition to age and gender, PDFF of paraspinal muscles is also affected by VBQ. Age and BMI are considered risk factors for the MCs in CLBP patients.