Iron deposition and mitochondrial dysfunction are closely associated with the genesis and progression of Parkinson's disease (PD). This study aims to extract susceptibility and oxygen extraction fraction (OEF) values of deep grey matter (DGM) to explore spatiotemporal progression patterns of brain iron-oxygen metabolism in PD.Ninety-five PD patients and forty healthy controls (HCs) were included. Quantitative susceptibility mapping (QSM) and OEF maps were computed from MRI multi-echo gradient echo data. Analysis of covariance (ANCOVA) was used to compare mean susceptibility and OEF values in DGM between early-stage PD (ESP), advanced-stage PD (ASP) patients and HCs. Then Granger causality analysis on the pseudo-time-series of MRI data was applied to assess the causal effect of early altered nuclei on iron content and oxygen extraction in other DGM nuclei.The susceptibility values in substantia nigra (SN), red nucleus, and globus pallidus (GP) significantly increased in PD patients compared with HCs, while the iron content in GP did not elevate obviously until the late stage. The mean OEF values for the caudate nucleus, putamen, and dentate nucleus were higher in ESP patients than in ASP patients or/and HCs. We also found that iron accumulation progressively expands from the midbrain to the striatum. These alterations were correlated with clinical features and improved AUC for early PD diagnosis to 0.824.Abnormal cerebral iron deposition and tissue oxygen utilization in PD measured by QSM and OEF maps could reflect pathological alterations in neurodegenerative processes and provide valuable indicators for disease identification and management.Noninvasive assessment of cerebral iron-oxygen metabolism may serve as clinical evidence of pathological changes in PD and improve the validity of diagnosis and disease monitoring.• Quantitative susceptibility mapping and oxygen extraction fraction maps indicated the cerebral pathology of abnormal iron accumulation and oxygen metabolism in Parkinson's disease. • Iron deposition is mainly in the midbrain, while altered oxygen metabolism is concentrated in the striatum and cerebellum. • The susceptibility and oxygen extraction fraction values in subcortical nuclei were associated with clinical severity.© 2023. The Author(s), under exclusive licence to European Society of Radiology.