Abstract: A system and method for non-contrast imaging of pulmonary blood flow in a subject are described. In some aspects, the method includes acquiring, using a magnetic resonance imaging (“MRI”) system, image data from at least the subject's lungs during which little to no respiratory motion occurs in the subject, such as during a breath-hold. The method also includes assembling the image data into a plurality of time-series datasets representing temporal variations of magnetic resonance signals in a region of interest that contains all or part of the subject's lungs. The method further includes computing a statistical blood flow metric for each voxel in the region of interest, using respective time-series datasets, and generating at least one image representative of pulmonary blood flow in the subject using the computed statistical blood flow metrics.