The purpose of the current study was to develop and validate a three-dimensional (3D) free-breathing cardiac T 1 -mapping sequence using SAturation-recovery and Variable-flip-Angle (SAVA). SAVA sequentially acquires multiple electrocardiogramtriggered volumes using a multishot spoiled gradient-echo sequence. The first volume samples the equilibrium signal of the longitudinal magnetization, where a flip angle of 2 is used to reduce the time for the magnetization to return to equilibrium. The succeeding three volumes are saturation prepared with variable delays, and are acquired using a 15 flip angle to maintain the signal-to-noise ratio. A diaphragmatic navigator is used to compensate the respiratory motion. T 1 is calculated using a saturation-recovery model that accounts for the flip angle. We validated SAVA by simulations, phantom, and human subject experiments at 3 T. SAVA was compared with modified Look-Locker inversion recovery (MOLLI) and saturation-recovery single-shot acquisition (SASHA) in vivo. In phantoms, T 1 by SAVA had good agreement with the reference (R 2 = 0.99). In vivo 3D T 1 mapping by SAVA could achieve an imaging resolution of 1.25 Â 1.25 Â 8 mm 3 . Both global and septal T 1 values by SAVA (1347 ± 37 and 1332 ± 42 ms) were in between those by SASHA (1612 ± 63 and 1618 ± 51 ms) and MOLLI (1143 ± 59 and 1188 ± 65 ms). According to the standard deviation (SD) and coefficient of variation (CV), T 1 precision measured by SAVA