Scholarly article on topic '4D flow cardiovascular magnetic resonance consensus statement'

4D flow cardiovascular magnetic resonance consensus statement Academic research paper on "Medical engineering"

Share paper
OECD Field of science

Academic research paper on topic "4D flow cardiovascular magnetic resonance consensus statement"

Dyverfeldt et al. Journal of Cardiovascular Magnetic Resonance (2015) 17:72 DOI 10.1186/s12968-015-0174-5

Journal of Cardiovascular Magnetic Resonance


Open Access

4D flow cardiovascular magnetic resonance consensus statement


Petter Dyverfeldt1'2*1, Malenka Bissel!31, Alex J. Barker4, Ann F. Bolger1,2,5, Carl-Johan Carlhäll1,2,6, Tino Ebbers1,2, Christopher J. Francios7, Alex Frydrychowicz8, Julia Geiger9, Daniel Giese10, Michael D. Hope11, Philip J. Kilner12, Sebastian Kozerke13, Saul Myerson3, Stefan Neubauer3, Oliver Wieben7,14 and Michael Markl4,15


Pulsatile blood flow through the cavities of the heart and great vessels is time-varying and multidirectional. Access to all regions, phases and directions of cardiovascular flows has formerly been limited. Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has enabled more comprehensive access to such flows, with typical spatial resolution of 1.5x1.5x1.5 - 3x3x3 mm3, typical temporal resolution of 30-40 ms, and acquisition times in the order of 5 to 25 min. This consensus paper is the work of physicists, physicians and biomedical engineers, active in the development and implementation of 4D Flow CMR, who have repeatedly met to share experience and ideas. The paper aims to assist understanding of acquisition and analysis methods, and their potential clinical applications with a focus on the heart and greater vessels. We describe that 4D Flow CMR can be clinically advantageous because placement of a single acquisition volume is straightforward and enables flow through any plane across it to be calculated retrospectively and with good accuracy. We also specify research and development goals that have yet to be satisfactorily achieved. Derived flow parameters, generally needing further development or validation for clinical use, include measurements of wall shear stress, pressure difference, turbulent kinetic energy, and intracardiac flow components. The dependence of measurement accuracy on acquisition parameters is considered, as are the uses of different visualization strategies for appropriate representation of time-varying multidirectional flow fields. Finally, we offer suggestions for more consistent, user-friendly implementation of 4D Flow CMR acquisition and data handling with a view to multicenter studies and more widespread adoption of the approach in routine clinical investigations.

Keywords: 4D Flow CMR, 4D Flow MRI, Phase-contrast magnetic resonance imaging, MR flow imaging, Hemodynamics, Flow visualization, Flow quantification, Recommendations, Clinical, Cardiovascular


Pulsatile blood flow through the cavities of the heart and great vessels is multidirectional and multidimensional. However, access to all the directions, regions and phases of such flows has been limited with cardiovascular magnetic resonance (CMR) as well as other imaging modalities. Four-dimensional (4D) flow CMR has been developed to attain more comprehensive access to blood flow through the heart and large vessels

* Correspondence: tEqualcontributors

1Division of Cardiovascular Medicine, Department of Medicaland Health Sciences, Linköping University, Linköping, Sweden 2Centerfor Medicallmage Science and Visualization, Linköping University, Linköping, Sweden

Fulllist of author information is available at the end of the article

Bio Med Central

[1-4]. This unique technique enables a wide variety of options for visualization and quantification of flow, ranging from basic aspects such as flow volume and peak velocity to more advanced features such as the estimation of hemodynamic effects at the vessel wall and myocardium, as well as visualization of flow pathways in the heart and great vessels.

"4D Flow CMR" refers to phase-contrast CMR with flow-encoding in all three spatial directions that is resolved relative to all three dimensions of space and to the dimension of time along the cardiac cycle (3D + time = 4D). For concise description and unification of nomenclature we recommend the use of the term "4D Flow CMR" or "4D Flow MRI" as will be used throughout the document. For methodological clarification, we recommend that a full-

© 2015 Dyverfeldt et al. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http:// applies to the data made available in this article, unless otherwise stated.

length description such as "three-dimensional (3D) cine (time-resolved) phase-contrast CMR with three-directional velocity-encoding" is included in the methods section of reports that employ this technique.

In 2012, leading teams in 4D Flow CMR started meeting on a regular basis to discuss the state of the technique and directions for future work. This group consists of medical physicists, physicians, and biomedical engineers who have been active in cardiovascular 4D Flow CMR research. As the field has grown, the list of invitees was extended by personal acquaintance and availability to join regular meetings. Soon, the need to find a uniform basis for terminology, potential clinical applications, and various technical aspects regarding data acquisition, data processing, visualization, and quantification became apparent. Therefore, we started a consensus initiative during the 2nd biannual 4D Flow Workshop held in Oxford, UK, in September 2013 to summarize our discussions.

This consensus statement aims to assist the understanding of the acquisition, analysis, and possible clinical applications of 4D Flow CMR in the heart and great vessels (aorta, pulmonary arteries), including all the steps involved in a 4D Flow CMR study. We will also discuss research and development goals that have yet to be achieved, in order to address current limitations and ensure data reliability and validity.

The consensus statement is based on published data and the shared experience of the 4D Flow CMR consensus group. The authors understand the manuscript as a state-of-the-art summary on the acquisition or analysis of 4D Flow CMR that creates a basis for future multi-vendor and multicenter research, can serve as a reference to established research, and provide guidance to researchers and clinicians new to the field.

Clinical and scientific significance

Flow assessment has long been used in the evaluation of cardiovascular disease. In recent decades, largely due to the advent of multidimensional flow imaging and computational fluid dynamics (CFD), the importance of improving our understanding of physiological and pathophysiological blood flow conditions is increasingly acknowledged [5, 6].

The most common clinical tool for cardiovascular flow assessment is Doppler echocardiography, which can measure the blood flow velocity component in the direction of the ultrasound beam or provide a 2D visualization of one-directional blood flow velocities using the color Doppler mode. Doppler echocardiography is often used to assess peak and mean velocities in the aorta and pulmonary artery for calculation of peak and mean pressure drops, known as gradients, via the simplified Bernoulli equation [7, 8]. This approach, however, is limited by the constraints associated with echo-Doppler imaging, which

include variable velocity assessment (due to beam alignment), limited acoustic window, and operator expertise [9—11]. Further, the calculation of mean velocities and net flow is often based on assumptions regarding the underlying flow profile and vessel cross-sectional area which may results in inaccurate flow quantification in the presence of complex flow and/or vessel geometry.

The most common CMR flow imaging technique is 2D cine phase contrast (PC) CMR with velocity-encoding in a single direction (2D cine PC-CMR) [12-19]. The single velocity-encoding direction is typically selected perpendicular to the 2D plane, which enables calculation of the volume of flow passing through the plane. 2D cine PC-CMR is arguably the gold-standard for flow volume quantification. The formerly widely used, but invasive thermodilution technique for flow quantification, is subject to inaccuracies due to underlying assumptions. Unlike Doppler echocardiography or 2D cine PC-CMR, 4D Flow CMR acquisition includes measurements representing all directions and spatial regions of flow within the boundaries of the volume imaged. Although methodologically different, computational fluid dynamics (CFD) is comparable in terms of multidirectional, volumetric representation [20-25]. Flow fields can potentially be calculated by CFD to high spatial and temporal resolution [26-28]. However, CFD requires accurate definition of geometrical and physiological boundary conditions and the ability of the computed flow fields to represent reality depends on the precision of the boundary conditions and the validity of underlying assumptions. For these reasons, CFD is currently not used in clinical decision-making. The relatively direct, voxel by voxel measurements of velocity provided by 4D Flow CMR can be complementary to the higher resolution velocity fields computed by CFD.

Clinical utility

CMR-based flow volume quantification is routinely used at many institutions to estimate shunt flows, regurgitant flows, collateral flows, etc. [29, 30]. These diagnostic tests are primarily based on 2D cine PC-CMR. A large number of studies across different institutions and MR-systems have demonstrated that 4D Flow CMR permits flow volume quantification that is comparable to 2D cine PC-CMR [31-39] and has good scan-rescan repeatability [32, 40, 41]. A recent study that assessed Qp/Qs ratios in intracardiac shunts reported that flow volumes are underestimated compared to 2D cine PC-CMR but that the Qp/Qs ratios were not different [42]. 2D cine PC-CMR has been hampered by artifacts such as background phase offsets which can lead to errors inflow volume measurements. These issues are shared by 4D Flow CMR and proper measures to compensate for them are necessary.

However, flow volume quantification with 4D Flow CMR has several advantages when compared to 2D cine PC-CMR. 4D Flow CMR permits investigation of the internal consistency of the data by employing the 'conservation of mass' principle (e.g. Qp/Qs within the same dataset). This important feature lends itself well to standardization of data quality assurance and will be discussed later. Investigators have used this feature and demonstrated that flow volume measurements with 4D Flow CMR have good internal consistency [32, 38, 41-49].

Another advantage of 4D Flow CMR is the retrospective placement of analysis planes at any location within the acquisition volume. While standard 2D cine PC techniques can easily be applied during a single breath hold, 4D Flow CMR on the other hand, offers the ability to retrospectively calculate blood flow through any planes of interest across the 3D volume. Despite longer scan times, 4D Flow CMR allows easy scan prescription (positioning of a single 3D volume) compared to the need to predetermine and accurately locate all relevant planes of 2D acquisitions. This may be especially advantageous in cases where multiple 2D cine PC-CMR scans would be needed [36]. In these situations, 4D Flow CMR may even be faster than prescribing and scanning a series of 2D cine breath-held PC CMR acquisitions, enabling a reduced period of anesthesia for younger children or of scan time in decompensated patients. Further, the option of valve tracking may improve assessment of flow through heart valves [43]. Compared to 2D cine PC-CMR, 4D Flow CMR measures velocity in all spatial directions and has superior spatial coverage and may therefore also be better at capturing the peak velocity of a stenotic jet [37]. However, one recent study suggested that the peak flow rate was lower for 4D compared to 2D flow [35]. These findings may in part be explained by relatively low temporal resolution (50-55 ms). Similarly, another recent study with 46 ms temporal resolution obtained smaller net flow volumes with 4D compared to 2D flow CMR [42]. Larger, preferably multicenter, studies with optimized protocols would be helpful to establish the comparability between 2D and 4D Flow CMR, as well as the spatial and temporal resolution needed for different applications.

In addition to the flexible retrospective quantification of conventional flow parameters, 4D Flow CMR allows for the visualization of multidirectional flow features and alterations of these associated with cardiovascular disease [50-53]. Previously reported results include the application of 4D Flow CMR for the analysis of blood flow in the ventricles [54-63] and atria [64-67] of the heart, heart valves [3, 43, 68, 69], aorta [41, 69-82], main pulmonary vessels [83-86], carotid arteries [87-90], large intracranial arteries and veins [91-98], arterial and portal venous systems of

the liver [46, 85, 99-101], peripheral arteries [102] and renal arteries [103, 104]. The intuitive flow visualizations that 4D Flow CMR offers have already found utility in several clinical studies. For example, time-resolved visualizations of blood flow have been used clinically to identify flow directionality and areas of flow acceleration in visceral abdominal blood flow [105-107]. In addition, a number of studies have shown that visualization of aortic blood flow can be helpful to quickly identify regions with high velocity flow close to the vessel wall that may indicate altered fluid mechanical effects on the vessel wall [69, 74, 108-110]. Finally, there are promising applications in complex congenital heart disease [39, 85, 111, 112]. While these examples are promising and illustrate the potential of 4D Flow analysis to better understand complex hemodynamic patterns, the clinical utility needs further evaluation in larger prospective and multi-center trials. For a more detailed overview of recent 4D Flow CMR developments and its use for 3D flow visualization and quantification throughout the human circulatory systems the reader is referred to a number of recently published review articles [113-118].

Research utility

4D Flow CMR has made it possible to investigate in-vivo cardiovascular flow fields more comprehensively than was previously possible. Multidisciplinary research teams are using the technique to 1) address gaps in the understanding of cardiovascular physiology and pathophysi-ology, 2) better understand the impact of hemodynamics on the heart and vasculature, 3) delineate further to what degree alterations of flow predispose to or result from cardiovascular disease processes such as remodeling, and 4) assess the degree to which physiological flow and pressure profiles have been restored following inter-ventional or surgical procedures. Thus, by affording visualization and quantification of flow parameters ranging from conventional parameters such as flow volume and regurgitant fraction to more advanced parameters such as flow energetics and shear stress, there are several applications where 4D Flow CMR has significant potential for advancing our knowledge and assessment of the cardiovascular system. For example, 4D Flow CMR has been used to demonstrate separation of blood that transits heart chambers according to compartmental origin and fate, retrograde flow embolization pathways from the descending aorta to the brain, and associations of valve outflow jet patterns with aortopathy [56, 74, 82, 109, 119-124]. In addition, the technique can be used to derive new physiologic and pathophysiologic hemodynamic parameters such as such as wall shear stress [125-127], pressure difference [103, 128-131], pulse wave velocity [132, 133], turbulent kinetic energy [134-137], and others [57, 58, 122, 138-141] for more differentiated

characterization of cardiovascular pathophysiology beyond simple measures of flow. The majority of these in-vivo hemodynamic measures cannot be assessed non-invasively with any other imaging technique. In all areas, further studies are required to assess the clinical impact of these measurements.

Consensus recommendations

This section provides recommendations for the use of 4D Flow CMR for 3D flow visualization and quantification of flow volume, retrograde flow and peak velocity in the heart and large vessels (aorta, pulmonary arteries). We focus on basic flow visualization and standard CMR parameters which can most easily been incorporated into routine clinical application. More advanced parameters such as pressure difference mapping, wall shear stress and turbulent kinetic energy are not addressed in this section but will be discussed in the future work section below. Table 1 lists analysis parameters and visualizations that are recommended for different clinical indications. The recommendations for quantification are based on the literature described in the "Clinical utility" section above. We emphasize that the amount of supporting literature is smaller for flow visualization

compared to quantification and thus the recommendations are primarily based on consensus discussions.

Several steps are required to assess blood flow with 4D Flow CMR including proper patient preparation, choice of acquisition parameters, and data conditioning through pre-processing, and data analysis. We suggest a structured workflow for data acquisition and processing as shown in Fig. 1.

Patient preparation

4D Flow CMR requires a reliable ECG trace with detectable R-wave to ensure consistency between RR-intervals. Standard ECG positioning applies. For aortic flow assessment it is important that the surface coils are positioned high enough to also fully encompass the aortic arch, which can be quite high in some aortic pathologies. 4D Flow CMR scans can be relatively long and it is useful to inform the patient about this prior to starting the scan in order to minimize discomfort.

4D flow CMR data acquisition PC-MR signal and use of contrast agents

4D Flow CMR employs spoiled gradient echo sequences with short TR for rapid imaging. As such, the

Table 1 Recommended 4D Flow CMR analysis for different clinical indications acquisition. For a comprehensive overview of 4D Flow CMR quantification and references please see recently published review articles [113-118]

- all aspects below can be derived from a single visualization methodology including additional




Heart valve disease (stenosis, regurgitation)

Shunts and collateralvessels (Ventricular-septal defect, atrial-septaldefect, fistulae)

Complex congenitalheart disease (e.g. single ventricle physiology, Fontan circulation, Fallot's tetralogy),

Aortic disease (aneurysm, coarctation, dissection)

Flow volume

• Regurgitant flow volumes & fraction Peak velocity

• Estimated pressure gradients with modified Bernoulli equation

Flow volume

• Shunt flow volume

• Qp/Qs Flow volume

• Regurgitant flow volumes & fraction

• Flow distribution (e.g. left vs right pulmonary artery, relative SVC/IVC flow)

• Collateralflow volume Peak velocity

Flow volume

• Regurgitant flow volumes & fraction

• Relative flows in true & false lumen Peak velocity

• Identification of regurgitant and stenotic jets using streamlines and pathlines

• Peak velocity location by systolic streamlines or maximum intensity projections of speed images

• Outflow patterns using streamlines

• Time course of flow curve

• Identification of shunt flow and flow directionality using pathlines

• Flow directionality using pathlines

• Shunt flow using pathlines

• Flow connectivity and distribution using pathlines

• Peak velocity location by systolic streamlines or maximum intensity projections of speed images

• Identification of flow in false lumen and potentialentry/ exit sites

• Identification of highly disrupted flow patterns (likely to reduce forward flow) in tortuous aortic conditions

aThe amount of supporting literature is smaller for visualization compared to quantification

Fig. 1 Recommended workflow for clinicalapplication of 4D Flow CMR with the main components of 1) patient preparation, 2) data acquisition in the magnet, 3) data reconstruction, 4) pre-processing of the reconstructed data, and 5) data analysis

signal magnitude from the blood is weighted inversely with the T1 relaxation time. This allows for the generation of PC angiograms without the need for an external contrast agent [107, 142]. Although 4D Flow CMR does not require any contrast agents, it is often used as part of a comprehensive CMR study that does requires the use of T1 shortening gadolinium-based contrast agents, for example for perfusion MRI or late gadolinium enhancement imaging. In such cases, acquiring the 4D Flow CMR data after the study that requires contrast administration takes advantage of the enhanced signal-to-noise ratio (SNR) and thus velocity to noise ratio (VNR) as well as contrast between blood and surrounding tissue [107, 143, 144]. However, contrast agents that wash out during the 4D Flow scan can result in time-varying blood T1 times and the effects of this variability on PC-CMR velocity data is not fully known. As a result, the effect of contrast agents depend on the type (extravascular vs. intravascular) and SNR, VNR can vary depending on the timing of the contrast agent administration.

Scan parameters

Table 2 provides an overview of the most important scan and reconstruction parameters and lists recommended values for each of them. The recommendations are primarily based on the experience of the authors of this consensus document and the literature that has investigated 4D Flow-based flow volume quantification [31-38, 40-48]. We propose this list of parameter settings as a baseline 4D Flow CMR protocol against which alterative protocols can be compared. Specialized applications such as measurements in a paediatric population or analysis of advanced flow parameters may require optimized parameter choices. Table 2 also indicates the potential advantages of 'ideal' parameters, the factors that may prevent their achievement and a value of the parameter that has been shown to be practicable and sufficient for most flow quantification in 4D Flow CMR covering the heart and/or large vessels (aorta, pulmonary artery) in healthy adults. The main technical goal of a 4D Flow CMR acquisition is accuracy. This is largely determined by sufficient spatial and temporal resolution, and adequate SNR (hence VNR). The accuracy

Table 2 4D Flow CMR scan parameters



Limiting factor

Consensus value

Acquisition Parameters Field of view Max

SNR, coverage


Velocity encoding timing (beat- vs. TR-interleaved)

k-space segmentation factor

Temporalresolutionc ECG synchronizationd

Respiratory motion compensationf

Maximum, at least Accuracy

5-6 voxels across vessel diameter of interest13, isotropic resolution.



100 % acceptance, motion correction

Partialk-space coverage Fullk-space coverage in phase- and slice-encoding directions

Flip Angleg

Parallellmaging k-t undersamplingh Venc

Ernst angle: a = acos(e-TR/T1)

No parallelimaging

No k-t under sampling

Maximum expected velocity, multiple vencs

Postprocessing Parameters Maxwellcorrection Yes Eddy current correction Yes

Phase unwrapping

Gradient non-linearity Yes correction

Avoid inter-cycle variability

Accuracy (temporal resolution)


Cover entire ECG cycle, avoid sequence interruption

Scan time, reduction of breathing artifacts

SNR, resolution


VNR, avoid aliasing

Accuracy Accuracy



Scan time, system imperfections

Scan time, SNR

Temporal resolution

Scan time Scan time

Reconstruction complexity

Reconstruction complexity, robustness, breathing artefacts (ghosting and blurring)

Scan time

Contrast vs. SNR

Scan time Scan time Scan time

Different methods and their validity and robustness

Different methods and their Yes validity and robustness

Availability If available

Cover region of interest

<2.5x2.5x2.5 mm3 for aorta or pulmonary artery

<3.0x3.0x3.0 mm3 for whole heart and greater vessels


<40 ms

If available: retrospective Else: Prospectivee

If available: Leading or trailing MR navigator on liver/diaphragm interface, 6 mm window size, typically resulting in 50 % acceptance rate.

Otherwise: Bellows with 50 % acceptance rate.

If available: Ellipticalk-space Otherwise: Half scan 75 % x 75 % (y x z)

Ernst angle

R = 2-3 (depends on #channels in coil array) If available: R = 4-5

Single venc, 10 % higher than maximum expected velocity

Yes Yes

aAlways indicate the effectively acquired resolution in combination with the interpolated resolution

bStudies have demonstrated that 5-6 voxels across the vessel diameter is sufficient for flow volume quantification [165]

cAlways indicate the effectively acquired resolution. If a temporal interpolation is performed, also indicate the interpolated temporal resolution along with the interpolation method used

dSo called self-gating techniques have been evaluated and may become an alternative to the ECG [32]

eFor prospective gating, analyses that involve integration over the whole cardiac cycle needs to be accompanied with a description of how the incomplete temporal coverage was handled

fDifferent types of respiratory navigators exist; variants include approaches that allow less motion in the central parts of k-space. Always describe the method that has been used and indicate the mean navigator efficiency in percent as well as the navigator acceptance window in mm. For fix window sizes and no k-space reordering, 6 mm navigator window is recommended, and this typically results in 50 % navigator efficiency

gThe SNR is strongly dependent on the in-flow effect, therefore the flip angle can be and is often chosen higher than the Ernst angle. When using contrast agents, the Ernst angle further increases (due to lower T1

hk-t undersampling factor 4-5 in combination with conventional parallel imaging factor 2-3 is not recommended

is also affected by artefacts. The practical achievability of the optimum parameter values is mostly hampered by total scan time, the methods availability and validity as well as system imperfections.

4D Flow CMR requires the user to define an upper velocity limit, termed the velocity encoding range (venc), similar to 2D cine PC-CMR techniques. Venc is defined as the (positive or negative) velocity that gives a phase

shift of n radians. Since phase is a cyclic entity, phase shifts greater than n radians result in velocity aliasing, which are visible as phase wraps in flow images. Higher venc results in lower VNR. We recommend choosing a venc slightly greater than the maximum velocity expected in the territory of interest. In stenotic and regurgitant flows, a multi-venc approach can be useful.

The total scan time available for adding 4D Flow CMR to a routine clinical CMR exam is often the most important limiting factor. If the total 4D Flow scan time is limited, e.g. not more than 5-8 min, the following tradeoffs may be useful:

1. Acquire free-breathing 4D Flow CMR without respiratory gating to increase scan efficiency (studies have demonstrated reasonably accurate flow volume quantification without compensation for respiratory motion) [32, 34].

2. Reduce temporal resolution by increasing the k-space segmentation factor to 3. This decreases the temporal resolution from approximately 40 ms to 60 ms and may result in reduced accuracy of peak velocity and flow volume quantification.

3. Reduce spatial resolution and SNR by acquiring 65 % x 65 % of ky and kz phase encoding lines

Employing these parameter adjustments can result in a substantial reduction of scan time. However, these changes will result in decreased spatiotemporal resolution and SNR and increased artifacts, which negatively impact flow quantification and visualization accuracy. Deviations from a validated standard protocol should be followed up by additional quality control.

In order to achieve comparability between different studies and to facilitate reproducibility of previously published work, a crucial requirement is the inclusion of all major scan and post-processing parameters in published reports. We recommend listing all scan parameters included in Table 2, and we encourage authors to specify the employed flow-encoding scheme, such as symmetric, asymmetric, or Hadamard 4-point encoding, 5-point-encoding, multipoint encoding etc. [137, 145-147]. The total scan time should be listed as the total scan time including respiratory gating efficiency or as the total scan time excluding navigator efficiency in combination with the respiratory gating efficiency.

Data pre-processing

4D Flow CMR data processing usually involves the use of automated or semi-automated corrections of known artefacts and often requires calculation of a geometric representation of the underlying 3D cardiac or vascular geometry through segmentation. Several sources of error can compromise 4D Flow CMR analysis and need to be

addressed prior to flow quantification and visualization. Similar to 2D cine PC-CMR, the major sources of errors include eddy current effects [148], concomitant gradient field effects (Maxell terms) [149], gradient field non-linearity [150, 151], and phase wraps resulting in velocity aliasing [152, 153]. Correction strategies have been presented in the literature and should be applied and evaluated to ensure accurate flow quantification and visualization [148-151, 154]. Investigators have also explored various types of image enhancement methods (noise filtering, divergence free corrections, etc.) to improve data quality. The use of such methods should be clearly reported in manuscripts, as they can also affect data quality negatively. Details and recommendations for the most common types of data processing are provided below. We emphasize that optimal approaches for data processing, especially corrections for background phase offsets, may vary between MR systems, sequences, protocols and applications.

Background phase offsets, concomitant gradient fields

Concomitant gradient fields, also referred to as Maxwell fields, lead to spatially varying background phase offsets in any type of PC-CMR acquisition. Correction factors for the concomitant gradient field correction can be directly derived from the gradient waveforms used for the data acquisition [149]. This correction scheme is implemented on MR systems as part of the standard PC-CMR image reconstruction engine.

Background phase offsets, Eddy currents The switching of time-varying magnetic field gradients result in changes in magnetic flux which in turn induce eddy currents in the conducting parts of the scanner system. These eddy currents alter the strengths and durations of the desired gradients and thus result in spatially and temporally varying phase offsets in any type of PC-CMR [155, 156]. Modern MR scanners have pre-emphasis systems that adjust the gradient waveforms by incorporating predictions of eddy currents effects. However, not all eddy current effects can be compensated for and there currently is no definite solution to remove all eddy current induced background phase offsets. We recommend the approach of fitting polynomials through the phase of tissue known to be static [148]. It should be noted that the order of the polynomial and the approach to detect static tissue may be vendor, sequence, and application specific. Assessment of heart-phase dependent differences is recommended.

Phase wraps, velocity aliasing

Blood flow velocities that exceed the velocity sensitivity (venc) value result in velocity aliasing, or phase wraps. We recommend that the venc is set higher than the maximum expected velocity. However, such a venc setting can cause

insufficient VNR in interesting flow regions with low velocity. Also, it is not always possible to predict the maximum velocity. We therefore recommend the use of a phase-unwrapping algorithm. The phase-unwrapping algorithm should be robust and not risk introducing additional errors. Identification of abrupt phase shifts in the temporal domain is a commonly used approach [153]. It should be noted that the visual perception and optimal phase-unwrapping strategies are different for different flow-encoding schemes.

Phase-Contrast Magnetic Resonance Angiography (PC-MRA)

4D Flow CMR data can be used to derive time-averaged 3D phase-contrast MR angiography (PC-MRA) based on the combination of velocity and magnitude data [142, 143, 157, 158]. The 3D PC-MRA can be used to guide anatomic orientation for flow visualization and regional flow quantification.

Data analysis Flow visualization

We generally recommend users of 4D Flow CMR to engage in visualizations and learn to interact with the data. Multiple options for the visualization of volumetric, time-resolved velocity vector fields on a 2D screen exist and none is entirely representative of the rich underlying data. It is a matter of choosing the proper visualization approach or combination of approaches that best address a particular question.

Visualization techniques commonly used with 4D Flow CMR include vector maps, streamlines and pathlines as well as maximum intensity projections, isosurfaces and volume renderings (see Figs. 2 and 3) [50, 51, 54, 159]. The choice of one technique over another, and the choice of a color map, depends on the application in question, the display medium and the time available for processing, among other factors. Visualization of 4D Flow data can be time-consuming and often benefits from informed user interaction. Interactive user-guided visualizations are valuable for generation of flow-based hypotheses. Efforts toward more standardized automated flow visualization approaches could be helpful in certain applications and would minimize operator-dependent variation, although users should understand the principles, strengths and limitations of different techniques. Further descriptions of the various approaches and their applications can be found in recent review articles [113-118,159, 160].

Relatively thorough interrogation of a flow field can be achieved by interactively browsing through the slices of a 4D dataset using vector map displays and/or color coding according to speed (magnitude of the velocity vector) (Fig. 2a) or by creating a maximum 'intensity' projection (MIP) for speed (Fig. 2b). Instantaneous streamlines can be used to represent the directions of flow throughout

a cavity or vessel at a given cardiac phase. Comparable instantaneous visualizations are achieved by 3D velocity vector maps, which display the magnitude and direction of blood velocity of each voxel. Instantaneous streamlines are traces through a 3D velocity field, parallel to the velocity vectors at all spatial points along their length. They represent a specific temporal phase (see Fig. 2c and Fig. 3c-e). An instantaneous streamline does not, in a pulsatile flow field, correspond with the path traveled by any given blood cell [161]. In contrast, path-lines (Fig. 2d and Fig. 3a-b), follow the paths of virtual massless particles. Streamlines and pathlines of a pulsatile flow field differ from each other, and each should be interpreted accordingly. To explore the path of blood through space and time, pathlines are likely to be more telling, whereas instantaneous streamlines would be more suitable for depictions of instantaneous flow features. Inclusion of the adjective "instantaneous" (instantaneous streamlines) helps to avoid confusion.

Many analysis parameters, such as flow speed, vorticity and turbulent kinetic energy, are scalar fields that can be visualized using MIP images or isosurface and volume rendering techniques. Isosurfaces and volume renders can be combined with vector graphs, streamlines or pathlines to create visualizations of multiple parameters. 4D Flow CMR visualizations may also be fused with other types of MR images, such as contrast-enhanced MR angiography and balanced steady-state free-precession (bSSFP) cine images to display anatomy. Such combinations can provide additional integration of cardiovascular morphology and function (see example in Fig. 3).

Flow quantification

Flow volume and retrograde flow quantification in 4D Flow CMR is similar to that used in conventional 2D cine PC-CMR. However, a few important differences exist. As mentioned above and illustrated in Fig. 4, the volumetric coverage of 4D Flow CMR offers retrospective positioning of planes for flow volume measurements at any location within the acquired data volume [37, 47, 49, 162-164]. The use of a 3D or 4D PC-MR angiogram derived from the 4D Flow CMR data is recommended for anatomical orientation and identification of cross-sectional analysis planes for flow quantification. This may be combined with streamlines or pathlines visualizations for further guidance of plane positioning. Segmentation of the lumen can be done with similar approaches as in 2D cine PC-CMR, but due to smaller inflow effects, the contrast between blood and surrounding tissue is inferior in 4D Flow CMR. Another option is to perform 3D or 4D segmentation during data processing and use this segmentation as geometrical boundaries in the flow volume calculation. This approach can be advantageous for quantification of peak velocity in

Fig. 2 Examples of 4D Flow CMR visualization techniques. Allexamples are based on data acquired in the aorta of a healthy volunteer. In these examples, flow visualization is overlaid onto a segmentation of the aorta. a An oblique slice that transects the aorta has been color-coded by flow speed and combined with a graph of velocity vectors which here displays the speed and direction of blood velocity in black arrows at a coarser grid than the acquired voxels. This type of visualization provides a quick overview of velocity fields. b A maximum intensity projection (MIP) image of flow speed permits identification of areas of elevated velocity and the point of peak velocity while displaying the peak velocities of the whole volume projected onto this single slice image. c Streamlines are instantaneously tangent to the velocity vector field and are usefulto visualize 3D velocity fields at discrete time points. Here, the peak systolic velocity field is shown. d Pathlines are the trajectories that massless fluid particles would follow through the dynamic velocity field. Pathlines are suitable for studies of the path of pulsatile blood flow over time. This example shows pathlines emitted from a plane in the ascending aorta at the onset of systole and traced to early systole (left), peak systole (middle) and late systole (right). Allfigures have been color-coded based on flow speed using the same color-window settings according to the scale shown in (b) and (d). In a, c and d, the visualizations have been combined with a PC-MRA isosurface which has been derived from the 4D Flow CMR data

an entire vessel segment rather than relying on 2D analysis planes which do not coincide with the location of the maximum systolic velocity. Studies with 2D PC-CMR have shown that at least 5-6 voxels across the vessel lumen are needed for accurate flow volume quantification [165].

Quality control

Quality control is important for every clinical and research study. The versatility of 4D Flow CMR allows

several approaches to be used, that can be included in imaging and post-processing without excessive additional effort.

Screening of 4D Flow CMR source images can reveal phase wraps, background phase offsets (by using narrow color-window), fold-over, and other image artifacts. Further, with its volumetric coverage, 4D Flow CMR offers several opportunities for control of the internal data consistency. For flow volume quantification, the conservation-of-mass

Pathlines emitted from the mitral valve at approx. time of peak A-wave (left panel) and traced to early systole (right panel)

Fig. 3 Examples of 4D Flow CMR visualization techniques, demonstrated on intracardiac flow data acquired in a healthy volunteer. In these examples, flow visualization is overlaid onto a 2D bSSFP acquisition in a three-chamber view. a Pathlines are the trajectories that massless fluid particles would follow through the dynamic velocity field and are suitable for studies of the path of pulsatile blood flow over time. Here, the transit of blood through the left ventricle (LV) is shown by pathlines emitted from the mitral valve at the time point of peak A-wave and traced to the time point of early systole systole. The timing of the ECG (TECG) is included for reference. b-d Streamlines are instantaneously tangent to the velocity vector field and are usefulto visualize 3D velocity fields at discrete time points. Here, streamlines generated in a long-axis plane show parts of the intracardiac velocity field at the time points of b peak early filling (E-wave), c peak late filling (A-wave), and d peak systole

principle can be employed to assess pulmonary vs. systemic flow volume ratios and flow volume in vs. out of the left ventricle [32, 38, 41-49]. The conservation-of-mass principle can also be used for quality control of pathlines analysis as the number of pathlines that enter and leave a specified region of interest should be the same (e.g. cardiac ventricles) [57]. Another complementary approach is to screen data for streamlines or pathlines that abruptly change direction or slowly drift out of the lumen, which can be indicative of phase wraps or uncompensated background phase offsets, respectively. Similarly, the presence of uncompensated background phase offsets can be suspected if pathlines emitted from the chest or back move in a non-random fashion.

The following approaches are recommended for general data quality control:

1) Visual inspection of source images

2) Quantitative quality control that targets the parameter of interest. For example, the conservation-of-mass principle is an excellent option when assessing the quality of flow volume

quantification. We emphasize that requirements on the data depend on the analysis approach; sufficient data quality for accurate estimation of parameter A (e.g. peak velocity) does not necessarily imply accurate estimation of parameter B (e.g. flow volume). When the quantitative quality control method matches the analysis parameter of interest, this may be used as the first-in-line quality control step.

3) When the quantitative quality control signals poor data quality, we recommend performing additional visual inspection of the source images, as well as inspection of pathlines emitted from static tissue such as the chest or back.

Controversies and recommendations for future work

The imaging sequences as well as data processing and analysis methods described in the recommendations section above constitute the current state of the technique as it is available at a large number of institutions. However, the field of 4D Flow CMR is rapidly evolving with improvements in imaging acquisition methods as

Fig. 4 Illustration of retrospective flow quantification. For retrospective quantification of flow parameters based on 2D analysis, planes can be positioned at any anatomic location. In this example, an isosurface of 3D PC-MRA data derived from the 4D Flow CMR data (gray shaded) has been used to guide positioning analysis planes throughout the thoracic aorta. For each analysis plane, the vesselcontours are segmented for all cardiac time frames to calculate flow volume, peak velocity and retrograde fraction

well as data processing and analysis techniques. These development efforts increase the diversity of 4D Flow CMR. We encourage this trend, but we also see a need for improved conformity across sites and companies that develop and use 4D Flow CMR methods. There are also several limitations of current and emerging methods, some of which are not fully understood, and it is important to acknowledge the limitations and develop our understanding of them, so that improvements can be made where possible. This section outlines some of the more advanced techniques and areas for development in improving data quality and simplifying wide-scale clinical applicability.

Advanced analysis parameters

Beyond the basic flow parameters discussed in the recommendation section above, a variety of more advanced analysis parameters are currently used in research settings, including wall shear stress, pressure difference mapping, kinetic energy, turbulent kinetic energy, energy

dissipation, differential flow analysis, flow angles, flow displacement, and pulse wave velocity [40, 57, 59, 60, 103, 122, 123, 125-139, 166-168]. The quantitative information that these parameters provide can distinguish normal from abnormal blood flow, as well as differentiating types of abnormal hemodynamics. However, the application of these parameters can be complex. The effects of underlying assumptions, the impact of the quality of 4D Flow CMR data on the parameter, and the physiological meaning of the parameter as estimated with 4D Flow CMR should be taken into consideration. While a complete list of analysis parameters is beyond the scope of this consensus document, Table 3 contains some of the most common parameters and describes what they are, potential applications, the controversies and unmet needs.

Understanding the limits of the technique

Further work is needed to understand the accuracy and precision of existing and new 4D Flow CMR methods,

Table 3 Commonly used advanced analysis parameters

Target parameter



Requirements and uncertainties

WallShear Stress (WSS) [93, 125, 126, 193]

Pulse Wave Velocity (PWV) [132, 133, 168]

Turbulent Kinetic Energy (TKE) [134, 136, 137]

Relative Pressure Fields [128, 129, 131]

Viscous shear forces of flowing blood acting tangentially to the vesselwall

Propagation speed of systolic pressure pulse in the arterial system

Energy content of turbulent flow and direction-independent measure of intensity of turbulent velocity fluctuations

Relative blood pressure field

Volume and Kinetic Energy of Ventricular Flow Components or Compartments [56, 57, 59]

Separation of blood that transits heart chambers according to compartmentalorigin and fate

Indicator for impact of flow alterations on endothelial celland extracellular matrix function and risk for vesselwall remodelling

Marker of arterialstiffness and predictive of cardiovascular disease.

Estimate of turbulence-related loss of energy or pressure. Indicator of impact of turbulent flow on blood constituents or vessel wall.

Noninvasive estimation of pressure differences

Indicator of ventricular dysfunction. Risk stratification and optimization and individualization of treatment heart failure

Dependent on spatialresolution. Relationship to actualWSS values are unclear [127, 173]. Limited longitudinal data that demonstrates its predictive value for risk stratification.

Requires high temporalresolution. Sensitive to artifacts.

The effect of intravoxelmean velocity variations affects the estimation of low TKE values. Is based on information from signal magnitude data from each individualflow-encoding segment, which are usually not obtained in standard reconstructions.

Pressure field calculations based on MR velocity data do not take turbulence effects into account and do therefore not reflect turbulence-related pressure losses that occur in stenotic flows. Computation of pressure fields is associated with severalpitfalls and a best strategy has not been established.

Pathlines used to map the transit of blood through the chambers accumulate errors that are inversely related to the quality of velocity data. Mixing effects are unknown.

including sequences, reconstruction methods and analysis parameters. The assessment of spatiotemporal fidelity and noise propagation of image acquisition, reconstruction and analysis methods is of key importance. Besides localization in space and time, any bias or noise-related uncertainty requires careful consideration, as it not only depends on the CMR experiment but also on MR system settings and tuning as shown in a recent 2D PC-CMR study [155].

Spatial and temporal resolution

The acquisition of 4D Flow CMR data is, in a certain sense, complete. All dimensions and directions of the cyclically changing flow field are covered, albeit with spatial and temporal resolution that does not resolve all features of the flow. If partial k-space acquisitions are used, the method used for reconstruction, e.g. zero-filling or Margosian/homodyne reconstruction, should be reported [169]. The method needs to be chosen with respect to its impact on the phase of the MR signal. While spatial resolution is typically quoted as the ratio of field-of-view to acquisition matrix, it needs to be emphasized that the effective spatial resolution can be less. Likewise, the ability to resolve temporal features of flow may not be appropriately captured by quoting the number of acquired heart phases and any methods for temporal interpolation or view sharing should be reported. Utilizing the concept of spatiotemporal point-spread function (PSFxt) or transfer function is recommended for detailed investigations of a method's ability to portray

information [170-173]. Choices of spatial and temporal resolution need to be made according to the degree of spatial localization and temporal bandwidth required to sufficiently describe, depict, and measure the flow feature of interest. We emphasize that resolution is driven by application, and recommendations for measuring parameters such as flow volume may not be sufficient for quantities such as wall shear stress or pulse wave velocity. Careful choices and investigations are required if quantities are derived from the measured velocity vector fields including spatial and temporal velocity derivatives as required for assessing wall shear stresses, relative pressure fields or pulse wave velocities, for example. Many parameters are directly affected by the choice of temporal and spatial resolution and therefore the impact of spatial and temporal resolution on the accuracy and precision of a given parameter should be considered. It is recommended to assess if a different resolution would produce a different result. As a way of avoiding direction-dependent estimates, the acquisition of iso-tropic voxels is recommended. If this is not possible, the effect of voxel anisotropy should be investigated.

Mean flow and small-scale variation in velocity

The time-resolved velocity fields measured with 4D Flow CMR are mean velocity fields and should be viewed as such. Spatial averaging occurs over the spatial extent of the voxel, and each measured cardiac phase (time frame) represents flow fields effectively averaged (phase-averaged) over multiple cardiac cycles extending over several

minutes. The spatiotemporal resolution and effective averaging over multiple cardiac cycles limits the size of the flow features that can be characterized with velocity mapping techniques. However, the measured mean velocity field is accurate and corresponds very well to the actual mean velocity field [174-176]. In disturbed and turbulent flows, a fluctuating velocity field is superimposed on the mean velocity field. These small-scale velocity fluctuations are thus not resolved by 2D or 4D Flow CMR velocity mapping. In fact, resolving all scales of velocity is not a realistic goal for 4D Flow CMR velocity mapping, as this would require <0.1 mm spatial resolution and <1 ms real-time temporal resolution. However, this aspect of flow can be addressed by a complementary 4D Flow CMR technique referred to as intravoxel velocity standard deviation (IVSD) mapping, or turbulence mapping. This technique, which can be viewed as a flow-analogue to diffusion-weighted imaging, is based on an MR signal model that describes the relationship between the amplitude (not phase) of the PC-CMR signal and the range of velocities that are present in a voxel. The IVSD mapping technique permits the estimation of the intensity of turbulent velocity fluctuations and turbulent kinetic energy in stenotic flows [135, 174, 176, 177]. Its application in flows with only minor fluctuations may be hampered by the fact that laminar flow effects such as shear also give rise to intravoxel velocity variations that contribute to the measured IVSD. However, this effect appears to be small compared to intravoxel velocity variations caused by unstable fluctuations [60, 135].

Noise propagation and confidence

Noise remains a limitation of the technique. An important parameter with respect to noise is the venc parameter that determines the velocity sensitivity of a 4D Flow CMR acquisition. The VNR is inversely proportional to the venc. Consequently, for a given venc, the estimation of low blood flow velocities < < venc is less reliable compared to flow velocities closer to venc. This can particularly be a limiting factor for multi-purpose flow analysis (e.g. quantification of both high flow velocities in a stenotic aorta and low flow velocities in a cardiac shunt in the same patient). New sequences are under development that permit the use of two or more venc's [137,147,178]. In addition, other strategies can be employed to maximize SNR and VNR. This includes optimizing the experimental setup, including main magnetic field strength and receive-coil instrumentation, as well as protocol modifications.

Further work is required to understand the impact of noise, and we recommend the method of pseudo replicas [179, 180] to study and assess SNR and VNR dependencies and noise propagation. Accordingly, different noise realizations of same statistics are added to the

original MR raw data and image reconstruction or parameter calculation is repeated to provide confidence intervals of velocity values. In a similar fashion, postprocessing strategies including the impact of region-of-interest analysis can be tested and referenced.

Systematic errors

Systematic errors causing unwanted bias of the measured velocity field are typically related to gradient induced eddy-currents, concomitant gradient fields and gradient non-linearity [148-151, 154]. While the latter two sources of error are corrected/calibrated with sufficient accuracy by clinical MR systems, eddy-currents depend on a range of parameters including the pre-emphasis settings of an individual MR system, gradient performance, orientation of the image volume and temperature of the gradient mount. Accordingly, prediction of the bias is often impossible and correction methods need to be applied retrospectively during image post-processing.

It is recommended to carefully study potential bias in a static gel phantom under identical experimental conditions, including navigators. This includes analysis of the spatial order of eddy-current induced background phase errors for each acquired heart phase [181]. If phantom calibration is applied routinely to subtract potential background phase offsets, fitted functions should be used to avoid compromising SNR/VNR of the original data upon subtraction. If background phase errors are corrected for by fitting polynomials through the phase of tissue known to be static, the fit error needs to be weighted against the degrees of freedom of the fit function to avoid over- or underfitting. Assessment of heart-phase dependent differences is recommended.


In-vivo comparison against current gold-standard methods is lacking for many areas of 4D Flow CMR, often due to the lack of such a gold standard for in-vivo assessment. The entire chain of data acquisition, reconstruction and image processing should, if possible, be evaluated for accuracy and precision. It is helpful to compare with existing techniques, where they exist. However, 4D Flow CMR may provide more accurate quantification and so potentially become the new gold-standard, or it may be the only technique capable of assessing certain parameters. For areas where an in-vivo gold-standard is lacking, controlled steady and pulsatile flow phantom experiments with accurate reference quantification can be used to assess accuracy. In view of the range of commercial and custom-built phantoms available, it should be feasible to validate applications by simulating flow rates and pulsatility (e.g. Reynolds and Womersley number), cycle-to-cycle variation and presence of sufficient static tissue for background correction. Reference methods (Particle Tracking or Image Velocimetry,

Laser Doppler Anemometry) can also be used to establish baseline data in-vitro. Numerical phantoms providing idealized model data are important for the study of certain aspects of data reconstruction and processing in flow fields that are fully known. Evaluation of precision should include testing sensitivity especially to spatiotemporal resolution and SNR, which can be done in-vivo, in-vitro, or using simulations.

Status of implementation and standardization Sequence

In addition to the sequence settings described in Table 2, several options exist for non-Cartesian 4D Flow CMR [182-186]. Different acquisition strategies (Cartesian, spiral, radial, EPI, bSSFP, etc.) have different strengths and weaknesses and thus the optimal acquisition strategy depends on the targeted application and analysis parameter. Moreover, in addition to standard parallel imaging, more advanced acceleration techniques have shown promising results and there exists many options for reduction of 4D Flow CMR scan times [39, 49, 187-192]. Reduced scan times are particularly relevant to applications in smaller vessels where higher resolution is needed.

At the time of writing, none of the major MR systems manufacturers (GE, Philips, Siemens) routinely provide 4D Flow CMR sequences or packages to researchers or clinical users. On Philips scanners, however, the necessary sequence exists and users can set up a 4D Flow CMR protocol ('exam card'), similar to the consensus protocol, on a standard commercial system without any software modifications. Siemens offer a 'work-in-progress' package to selected users. Due to the lack of widely available commercial 4D Flow CMR sequences, a large number of studies and applications are still based on 4D Flow pulse sequences that individual research groups have developed in-house and shared with collaborators worldwide. 4D Flow CMR pulse sequences implemented by research groups exist in a variety of flavors (different k-space trajectories, acceleration methods, etc.) and for all major MR platforms (GE, Philips, Siemens). The lack of standardization across MR platforms (even for the same vendor) and data output formats, as well as the absence of commercial 4D Flow CMR sequences and protocols are limiting factors that hinder introduction of the technique to the clinical environment.

Software for pre-processing, visualization and flow quantification

Pre-processing, visualization and flow quantification is being performed using in-house developed tools, early-stage commercial packages, or manufacturer prototypes. The field would benefit from greater standardization in data analysis methods, workflows, and data output formats,

which in turn affect the use of analysis tools. Wide clinical utility would benefit from the availability of user-friendly tools that are integrated in MR scanner consoles and workstations, as well as PACS systems. This would ideally include the following capabilities: 1) retrospective flow quantification on the scanner console and/or workstations and/or PACS system, 2) analysis and representation of clinically relevant parameters such as flow waveforms and cardiac output in DICOM format 3) 4D Flow visualizations on MR scanner console and/or workstations and/or PACS systems, and 4) animations in DICOM format to store and display in PACS. We encourage vendors and third-party developers to consider implementing these key features as a basis for more routine clinical use of 4D Flow CMR.


Relatively easy scan prescription and retrospective placement of analysis planes makes 4D Flow CMR a potentially advantageous tool in the clinical setting, particularly if several regions and directions of flow merit investigation. Conventional flow parameters can be obtained at any location in the data volume where the employed parameter settings provide sufficient accuracy. At the same time, 4D Flow CMR visualizations offer more versatile and comprehensive depictions of flow fields than any other in-vivo imaging technique. Further, advanced 4D Flow CMR analysis parameters are currently used in the research setting but require testing for clinical utility. Widespread clinical usage would be facilitated by further integration into the standard MR environment. Multicenter studies are necessary to establish the repeatability of various aspects of the technique across centers.

Competing interests

The authors declare that they have no competing interests. Authors' contributions

Allauthors reviewed and provided input before and during the drafting of the manuscript. PD, MB, PJK, and MM drafted and edited the manuscript. All authors read and approved the finalmanuscript.


Petter Dyverfeldt acknowledges funding from the Swedish Research Council, the Medical Research Council of Southeast Sweden, and Linkoping University. Malenka Bissellacknowledges funding from the British Heart Foundation Centre of Research Excellence and the Oxford NIHR BiomedicalResearch Centre.

Alex Barker acknowledges funding from NIH K25HL119608.

Carl-Johan Carlhallacknowledges funding from the Swedish Heart and Lung

Foundation, and Linkoping University.

Tino Ebbers acknowledges funding from the Swedish Research Counciland the European Research Council(HEART4FLOW, 310612). MichaelHope acknowledges funding from the RadiologicalSociety of North America (RSNA) Research Scholar Grant.

Philip Kilner acknowledges funding from the NIHR Cardiovascular Biomedical Research Unit at the RoyalBrompton and Harefield NHS Foundation Trust and ImperialCollege London.

Saul Myerson acknowledges funding from the British Heart Foundation Centre of Research Excellence and the Oxford NIHR Biomedical Research Centre.

Stefan Neubauer acknowledges funding from the British Heart Foundation Centre of Research Excellence and the Oxford NIHR Biomedical Research Centre.

Oliver Wieben acknowledges funding from the National Institutes of Health (NIH) National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) R01 R01DK096169.

Michael Markl acknowledges funding from National Institute of Health (NIH) National Heart, Lung, and Blood Institute (NHLBI) grant R01HL115828.

Author details

1 Division of Cardiovascular Medicine, Department of Medical and Health Sciences, Linköping University, Linköping, Sweden. 2Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden. 3Division of Cardiovascular Medicine, Radcliffe Department of Medicine, University of Oxford Centre for Clinical Magnetic Resonance Research, Oxford, UK. 4Department of Radiology, Northwestern University, Chicago, USA. 5Department of Medicine, University of California San Francisco, San Francisco, CA, United States. 6Department of Clinical Physiology, Department of Medical and Health Sciences, Linköping University, Linköping, Sweden. 7Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA. 8Klinik für Radiologie und Nuklearmedizin, University Hospital Schleswig-Holstein, Campus Lübeck, Lübeck, Germany. 9Department of Radiology, University Children's Hospital Zurich, Zurich, Switzerland. 10Department of Radiology, University Hospital of Cologne, Cologne, Germany. "Department of Radiology, University of California San Francisco, San Francisco, CA, United States. 12NIHR Cardiovascular Biomedical Research Unit, Royal Brompton and Harefield NHS Foundation Trust, National Heart and Lung Institute, Imperial College, London, UK. 13Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland. 14Department of Medical Physics, University of Wisconsin, Madison, Wisconsin, USA. 15Department of Biomedical Engineering, Northwestern University, Chicago, IL, USA.

Received: 24 March 2015 Accepted: 17 July 2015 Published online: 10 August 2015


1. Firmin DN, Gatehouse PD, Konrad JP, Yang GZ, Kilner PJ, Longmore DB. Rapid 7-dimensional imaging of pulsatile flow. Computers in Cardiology IEEE computer Society, London. 1993;14:353-6.

2. Wigstrom L, Sjoqvist L, Wranne B. Temporally resolved 3D phase-contrast imaging. Magn Reson Med. 1996;36:800-3.

3. Kozerke S, Hasenkam JM, Pedersen EM, Boesiger P. Visualization of flow patterns distal to aortic valve prostheses in humans using a fast approach for cine 3D velocity mapping. J Magn Reson Imaging. 2001;13:690-8.

4. Markl M, Chan FP, Alley MT, Wedding KL, Draney MT, Elkins CJ, et al. Time-resolved three-dimensional phase-contrast MRI. J Magn Reson Imaging. 2003;17:499-506.

5. Richter Y, Edelman ER. Cardiology Is Flow. Circulation. 2006;113:2679-82.

6. Carlhall CJ, Bolger A. Passing strange: flow in the failing ventricle. Circ Heart Fail. 2010;3:326-31.

7. Vahanian A, Alfieri O, Andreotti F, Antunes MJ, Baron-Esquivias G, Baumgartner H, et al. Guidelines on the management of valvular heart disease (version 2012) The Joint Task Force on the Management of Valvular Heart Disease of the European Society of Cardiology (ESC) and the European Association for Cardio-Thoracic Surgery (EACTS). Eur Heart J. 2012;33:2451-96.

8. Nishimura RA, Otto CM, Bonow RO, Carabello BA, Erwin JP, Guyton RA, et al. 2014 AHA/ACC guideline for the management of patients with valvular heart disease: a report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines. J Am Coll Cardiol. 2014;63:e57-e185.

9. Fyrenius A, Wigström L, Bolger AF, Ebbers T, Öhman KP, Karlsson M, et al. Pitfalls in Doppler evaluation of diastolic function: insights from 3-dimensional magnetic resonance imaging. J Am Soc Echocardiogr. 1999;12:817-26.

10. Biner S, Rafique A, Rafii F, Tolstrup K, Noorani O, Shiota T, et al. Reproducibility of proximal isovelocity surface area, vena contracta, and

regurgitant jet area for assessment of mitral regurgitation severity. JACC Cardiovasc Imaging. 2010;3:235-43.

11. Bach DS. Echo/Doppler evaluation of hemodynamics after aortic valve replacement: principles of interrogation and evaluation of high gradients. JACC Cardiovasc Imaging. 2010;3:296-304.

12. Pelc NJ, Herfkens RJ, Shimakawa A, Enzmann DR. Phase contrast cine magnetic resonance imaging. Magn Reson Q. 1991;7:229-54.

13. Nayler GL, Firmin DN, Longmore DB. Blood flow imaging by cine magnetic resonance. J Comput Assist Tomogr. 1986;10:715-22.

14. Burt CT. NMR measurements and flow. J Nucl Med. 1982;23:1044-5.

15. Moran PR. A flow velocity zeugmatographic interlace for NMR imaging in humans. Magn Reson Imaging. 1982;1:197-203.

16. Bryant DJ, Payne JA, Firmin DN, Longmore DB. Measurement of flow with NMR imaging using a gradient pulse and phase difference technique.

J Comput Assist Tomogr. 1984;8:588-93.

17. Firmin DN, Nayler GL, Klipstein RH, Underwood SR, Rees RS, Longmore DB. In vivo validation of MR velocity imaging. J Comput Assist Tomogr. 1987;11:751-6.

18. Underwood SR, Firmin DN, Rees RS, Longmore DB. Magnetic resonance velocity mapping. Clin Phys Physiol Meas. 1990;11:37-43.

19. Atkinson DJ, Edelman RR. Cineangiography of the heart in a single breath hold with a segmented turboFLASH sequence. Radiology. 1991;178:357-60.

20. Taylor CA, Draney MT, Ku JP, Parker D, Steele BN, Wang K, et al. Predictive medicine: computational techniques in therapeutic decision-making. Comput Aided Surg. 1999;4:231-47.

21. Steinman DA. Image-based computational fluid dynamics modeling in realistic arterial geometries. Ann Biomed Eng. 2002;30:483-97.

22. Yoganathan AP, Chandran K, Sotiropoulos F. Flow in prosthetic heart valves: state-of-the-art and future directions. Ann Biomed Eng. 2005;33:1689-94.

23. Taylor CA, Steinman DA. Image-based modeling of blood flow and vessel wall dynamics: applications, methods and future directions. Ann Biomed Eng. 2010;38:1188-203.

24. Khalafvand SS, Ng EY, Zhong L. CFD simulation of flow through heart: a perspective review. Comput Methods Biomech Biomed Engin. 2011;14:113-32.

25. Chan BT, Lim E, Chee KH, Abu Osman NA. Review on CFD simulation in heart with dilated cardiomyopathy and myocardial infarction. Comput Biol Med. 2013;43:377-85.

26. Jin S, Oshinski J, Giddens DP. Effects of wall motion and compliance on flow patterns in the ascending aorta. J Biomech Eng. 2003;125:347-54.

27. Arzani A, Dyverfeldt P, Ebbers T, Shadden SC. In vivo validation of numerical prediction for turbulence intensity in an aortic coarctation. Ann Biomed Eng. 2012;40:860-70.

28. Lantz J, Ebbers T, Engvall J, Karlsson M. Numerical and experimental assessment of turbulent kinetic energy in an aortic coarctation. J Biomech. 2013;46:1851-8.

29. Srichai MB, Lim RP, Wong S, Lee VS. Cardiovascular applications of phase-contrast MRI. Am J Roentgenol. 2009;192:662-75.

30. Chai P, Mohiaddin R. How we perform cardiovascular magnetic resonance flow assessment using phase-contrast velocity mapping. J Cardiovasc Magn Reson. 2005;7:705-16.

31. Brix L, Ringgaard S, Rasmusson A, Sorensen TS, Kim WY. Three dimensional three component whole heart cardiovascular magnetic resonance velocity mapping: comparison of flow measurements from 3D and 2D acquisitions. J Cardiovasc Magn Reson. 2009;11:3.

32. Uribe S, Beerbaum P, Sorensen TS, Rasmusson A, Razavi R, SchaeffterT. Four-dimensional (4D) flow of the whole heart and great vessels using real-time respiratory self-gating. Magn Reson Med. 2009;62:984-92.

33. Hope MD, Meadows AK, Hope TA, Ordovas KG, Saloner D, Reddy GP, et al. Clinical evaluation of aortic coarctation with 4D flow MR imaging. J Magn Reson Imaging. 2010;31:711-8.

34. Nordmeyer S, Riesenkampff E, Crelier G, Khasheei A, Schnackenburg B, Berger F, et al. Flow-sensitive four-dimensional cine magnetic resonance imaging for offline blood flow quantification in multiple vessels: a validation study. J Magn Reson Imaging. 2010;32:677-83.

35. Carlsson M, Toger J, Kanski M, Bloch KM, Stahlberg F, Heiberg E, et al. Quantification and visualization of cardiovascular 4D velocity mapping accelerated with parallel imaging or kt BLAST: head to head comparison and validation at 1.5 T and 3 T. J Cardiovasc Magn Reson. 2011;13:55.

36. Valverde I, Nordmeyer S, Uribe S, Greil G, Berger F, Kuehne T, et al. Systemic-to-pulmonary collateral flow in patients with palliated

univentricular heart physiology: measurement using cardiovascular magnetic resonance 4D velocity acquisition. J Cardiovasc Magn Reson. 2012;14:25.

37. Nordmeyer S, Riesenkampff E, Messroghli D, Kropf S, Nordmeyer J, Berger F, et al. Four-dimensional velocity-encoded magnetic resonance imaging improves blood flow quantification in patients with complex accelerated flow. J Magn Reson Imaging. 2013;37:208-16.

38. Hsiao A, Tariq U, Alley MT, Lustig M, Vasanawala SS. Inlet and outlet valve flow and regurgitant volume may be directly and reliably quantified with accelerated, volumetric phase-contrast MRI. J Magn Reson Imaging. 2014;41(2):376-85.

39. Giese D, Wong J, Greil GF, Buehrer M, Schaeffter T, Kozerke S. Towards highly accelerated Cartesian time-resolved 3D flow cardiovascular magnetic resonance in the clinical setting. J Cardiovasc Magn Reson. 2014;16:42. doi:10.1186/1532-429X-16-42.

40. Markl M, Wallis W, Harloff A. Reproducibility of flow and wall shear stress analysis using flow-sensitive four-dimensional MRI. J Magn Reson Imaging. 2011;33:988-94.

41. Wentland AL, Grist TM, Wieben O. Repeatability and internal consistency of abdominal 2D and 4D phase contrast MR flow measurements. Acad Radiol. 2013;20:699-704.

42. Hanneman K, Sivagnanam M, Nguyen ET, Wald R, Greiser A, Crean AM, et al. Magnetic Resonance Assessment of Pulmonary (Qp) to Systemic (Qs) Flows Using 4D Phase-contrast Imaging: Pilot Study Comparison with Standard Through-plane 2D Phase-contrast Imaging. Acad Radiol. 2014;21:1002-8.

43. Westenberg JJ, Roes SD, Ajmone Marsan N, Binnendijk NM, Doornbos J, Bax JJ, et al. Mitral valve and tricuspid valve blood flow: accurate quantification with 3D velocity-encoded MR imaging with retrospective valve tracking. Radiology. 2008;249:792-800.

44. Roes SD, Hammer S, van der Geest RJ, Marsan NA, Bax JJ, Lamb HJ, et al. Flow assessment through four heart valves simultaneously using 3-dimensional 3-directional velocity-encoded magnetic resonance imaging with retrospective valve tracking in healthy volunteers and patients with valvular regurgitation. Invest Radiol. 2009;44:669-75.

45. Nett EJ, Johnson KM, Frydrychowicz A, Del Rio AM, Schrauben E, Francois CJ, et al. Four-dimensional phase contrast MRI with accelerated dual velocity encoding. J Magn Reson Imaging. 2012;35(6):1462-71.

46. Roldan-Alzate A, Frydrychowicz A, Niespodzany E, Landgraf BR, Johnson KM, Wieben O, et al. In vivo validation of 4D flow MRI for assessing the hemodynamics of portal hypertension. J Magn Reson Imaging. 2013;37:1100-8.

47. Frydrychowicz A, Wieben O, Niespodzany E, Reeder SB, Johnson KM, Francois CJ. Quantification of Thoracic Blood Flow Using Volumetric Magnetic Resonance Imaging With Radial Velocity Encoding: In Vivo Validation. Invest Radiol. 2013;48:819-25.

48. Barker AJ, Roldan-Alzate A, Entezari P, Shah SJ, Chesler NC, Wieben O, et al. Four-dimensional flow assessment of pulmonary artery flow and wall shear stress in adult pulmonary arterial hypertension: Results from two institutions. Magn Reson Med. 2014;73(5):1904-13.

49. Tariq U, Hsiao A, Alley M, Zhang T, Lustig M, Vasanawala SS. Venous and arterial flow quantification are equally accurate and precise with parallel imaging compressed sensing 4D phase contrast MRI. J Magn Reson Imaging. 2013;37:1419-26.

50. Napel S, Lee DH, Frayne R, Rutt BK. Visualizing three-dimensional flow with simulated streamlines and three-dimensional phase-contrast MR imaging.

J Magn Reson Imaging. 1992;2:143-53.

51. Buonocore MH. Visualizing blood flow patterns using streamlines, arrows, and particle paths. Magn Reson Med. 1998;40:210-26.

52. Markl M, Draney MT, Hope MD, Levin JM, Chan FP, Alley MT, et al. Time-Resolved 3-Dimensional Velocity Mapping in the Thoracic Aorta: Visualization of 3-Directional Blood Flow Patterns in Healthy Volunteers and Patients. J Comput Assist Tomogr. 2004;28:459-68.

53. Bogren HG, Mohiaddin RH, Kilner PJ, Jimenez-Borreguero LJ, Yang GZ, Firmin DN. Blood flow patterns in the thoracic aorta studied with three-directional MR velocity mapping: the effects of age and coronary artery disease. J Magn Reson Imaging. 1997;7:784-93.

54. Wigstrom L, Ebbers T, Fyrenius A, Karlsson M, Engvall J, Wranne B, et al. Particle trace visualization of intracardiac flow using time-resolved 3D phase contrast MRI. Magn Reson Med. 1999;41:793-9.

55. Kilner PJ, Yang GZ, Wilkes AJ, Mohiaddin RH, Firmin DN, Yacoub MH. Asymmetric redirection of flow through the heart. Nature. 2000;404:759-61.

56. Bolger AF, Heiberg E, Karlsson M, Wigstrom L, Engvall J, Sigfridsson A, et al. Transit of blood flow through the human left ventricle mapped by cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2007;9:741-7.

57. Eriksson J, Carlhall CJ, Dyverfeldt P, Engvall J, Bolger AF, Ebbers T. Semiautomatic quantification of 4D left ventricular blood flow. J Cardiovasc Magn Reson. 2010;12:9.

58. Toger J, Carlsson M, Soderlind G, Arheden H, Heiberg E. Volume Tracking: A new method for quantitative assessment and visualization of intracardiac blood flow from three-dimensional, time-resolved, three-component magnetic resonance velocity mapping. BMC Med Imaging. 2011;11:10.

59. Eriksson J, Bolger AF, Ebbers T, Carlhäll C-J. Four-dimensional blood flow-specific markers of LV dysfunction in dilated cardiomyopathy. European Heart Journal-Cardiovascular Imaging. 2013;14:417-24.

60. Zajac J, Eriksson J, Dyverfeldt P, Bolger AF, Ebbers T, Carlhäll CJ. Turbulent kinetic energy in normal and myopathic left ventricles. J Magn Reson Imaging. 2014;41(4):1021-9.

61. Calkoen EE, Roest AA, Kroft LJ, van der Geest RJ, Jongbloed MR, van den Boogaard PJ, et al. Characterization and improved quantification of left ventricular inflow using streamline visualization with 4DFlow MRI in healthy controls and patients after atrioventricular septal defect correction. J Magn Reson Imaging. 2014;41(6):1512-20.

62. Elbaz MS, Calkoen EE, Westenberg JJ, Lelieveldt BP, Roest AA, van der Geest RJ. Vortex flow during early and late left ventricular filling in normal subjects: quantitative characterization using retrospectively-gated 4D flow cardiovascular magnetic resonance and three-dimensional vortex core analysis. J Cardiovasc Magn Reson. 2014;16:78.

63. Eriksson J, Bolger AF, Carlhäll C-J, Ebbers T. Spatial heterogeneity of four-dimensional relative pressure fields in the human left ventricle. Magnetic Resonance in Medicine. 2015. doi: 10.1002/mrm.25539.

64. Fyrenius A, Wigstrom L, Ebbers T, Karlsson M, Engvall J, Bolger AF. Three dimensional flow in the human left atrium. Heart. 2001;86:448-55.

65. Fluckiger JU, Goldberger JJ, Lee DC, Ng J, Lee R, Goyal A, et al. Left atrial flow velocity distribution and flow coherence using four-dimensional FLOW MRI: a pilot study investigating the impact of age and Pre- and Postintervention atrial fibrillation on atrial hemodynamics. J Magn Reson Imaging. 2013;38:580-7.

66. Arvidsson PM, Töger J, Heiberg E, Carlsson M, Arheden H. Quantification of left and right atrial kinetic energy using four-dimensional intracardiac magnetic resonance imaging flow measurements. J Appl Physiol. 2013;114:1472-81.

67. Dyverfeldt P, Kvitting JPE, Carlhäll CJ, Boano G, Sigfridsson A, Hermansson U, et al. Hemodynamic aspects of mitral regurgitation assessed by generalized phase-contrast MRI. J Magn Reson Imaging. 2011;33:582-8.

68. Kvitting JP, Dyverfeldt P, Sigfridsson A, Franzen S, Wigstrom L, Bolger AF, et al. In vitro assessment of flow patterns and turbulence intensity in prosthetic heart valves using generalized phase-contrast MRI. J Magn Reson Imaging. 2010;31:1075-80.

69. Hope MD, Hope TA, Crook SE, Ordovas KG, Urbania TH, Alley MT, et al. 4D flow CMR in assessment of valve-related ascending aortic disease. JACC Cardiovasc Imaging. 2011;4:781-7.

70. Bogren HG, Buonocore MH. 4D magnetic resonance velocity mapping of blood flow patterns in the aorta in young vs. elderly normal subjects. J Magn Reson Imaging. 1999;10:861-9.

71. Bogren HG, Buonocore MH, Valente RJ. Four-dimensional magnetic resonance velocity mapping of blood flow patterns in the aorta in patients with atherosclerotic coronary artery disease compared to age-matched normal subjects. J Magn Reson Imaging. 2004;19:417-27.

72. Kvitting JP, Ebbers T, Wigstrom L, Engvall J, Olin CL, Bolger AF. Flow patterns in the aortic root and the aorta studied with time-resolved, 3-dimensional, phase-contrast magnetic resonance imaging: implications for aortic valve-sparing surgery. J Thorac Cardiovasc Surg. 2004;127:1602-7.

73. Markl M, Draney MT, Miller DC, Levin JM, Williamson EE, Pelc NJ, et al. Time-resolved three-dimensional magnetic resonance velocity mapping of aortic flow in healthy volunteers and patients after valve-sparing aortic root replacement.

J Thorac Cardiovasc Surg. 2005;130:456-63.

74. Hope MD, Hope TA, Meadows AK, Ordovas KG, Urbania TH, Alley MT, et al. Bicuspid aortic valve: four-dimensional MR evaluation of ascending aortic systolic flow patterns. Radiology. 2010;255:53-61.

75. Frydrychowicz A, Markl M, Hirtler D, Harloff A, Schlensak C, Geiger J, et al. Aortic Hemodynamics in Patients With and Without Repair of Aortic Coarctation: In Vivo Analysis by4D Flow-Sensitive Magnetic Resonance Imaging. Invest Radiol. 2011;46:317-25.

76. Burk J, Blanke P, Stankovic Z, Barker A, Russe M, Geiger J, et al. Evaluation of 3D blood flow patterns and wall shear stress in the normal and dilated thoracic aorta using flow-sensitive 4D CMR. J Cardiovasc Magn Reson. 2012;14:84.

77. Geiger J, Arnold R, Herzer L, Hirtler D, Stankovic Z, Russe M, et al. Aortic wall shear stress in Marfan syndrome. Magn Reson Med. 2012;70(4):1137-44.

78. von Knobelsdorff-Brenkenhoff F, Trauzeddel RF, Barker AJ, Gruettner H, Markl M, Schulz-Menger J. Blood flow characteristics in the ascending aorta after aortic valve replacement-a pilot study using 4D-flow MRI. Int J Cardiol. 2014;170:426-33.

Dyverfeldt P, Hope MD, Tseng EE, Saloner D. Magnetic resonance measurement of turbulent kinetic energy for the estimation of irreversible pressure loss in aortic stenosis. JACC Cardiovasc Imaging. 2013;6:64-71. Frydrychowicz A, Berger A, Munoz Del Rio A, Russe M, Bock J, Harloff A, et al. Interdependencies of aortic arch secondary flow patterns, geometry, and age analysed by 4-dimensional phase contrast magnetic resonance imaging at 3 Tesla. Eur Radiol. 2012;22:1122-30.

81. Geiger J, Markl M, Herzer L, Hirtler D, Loeffelbein F, Stiller B, et al. Aortic flow patterns in patients with Marfan syndrome assessed by flow-sensitive four-dimensional MRI. J Magn Reson Imaging. 2012;35:594-600.

82. Barker AJ, Markl M, Bürk J, Lorenz R, Bock J, Bauer S, et al. Bicuspid Aortic Valve Is Associated With Altered Wall Shear Stress in the Ascending AortaClinical Perspective. Circ Cardiovasc Imaging. 2012;5:457-66.

83. Reiter G, Reiter U, Kovacs G, Kainz B, Schmidt K, Maier R, et al. Magnetic resonance-derived 3-dimensional blood flow patterns in the main pulmonary artery as a marker of pulmonary hypertension and a measure of elevated mean pulmonary arterial pressure. Circ Cardiovasc Imaging. 2008;1:23-30.

84. Geiger J, Markl M, Jung B, Grohmann J, Stiller B, Langer M, et al. 4D-MR flow analysis in patients after repair for tetralogy of Fallot. Eur Radiol. 2011;21:1651-7.

85. Francois CJ, Srinivasan S, Schiebler ML, Reeder SB, Niespodzany E, Landgraf BR, et al. 4D cardiovascular magnetic resonance velocity mapping of alterations of right heart flow patterns and main pulmonary artery hemodynamics in tetralogy of Fallot. J Cardiovasc Magn Reson. 2012;14:16.

86. Bachler P, Pinochet N, Sotelo J, Crelier G, Irarrazaval P, Tejos C, et al. Assessment of normal flow patterns in the pulmonary circulation by using 4D magnetic resonance velocity mapping. Magn Reson Imaging. 2013;31:178-88.

87. Harloff A, Albrecht F, Spreer J, Stalder A, Bock J, Frydrychowicz A, et al. 3D blood flow characteristics in the carotid artery bifurcation assessed by flow-sensitive 4D MRI at 3T. Magn Reson Med. 2009;61:65-74.

88. Markl M, Wegent F, Zech T, Bauer S, Strecker C, Schumacher M, et al. In Vivo Wall Shear Stress Distribution in the Carotid Artery: Effect of Bifurcation Geometry, Internal Carotid Artery Stenosis, and Recanalization Therapy. Circ Cardiovasc Imaging. 2010;3:647-55.

89. Meckel S, Leitner L, Bonati LH, Santini F, Schubert T, Stalder AF, et al. Intracranial artery velocity measurement using 4D PC MRI at 3 T: comparison with transcranial ultrasound techniques and 2D PC MRI. Neuroradiology. 2013;55:389-98.

90. Tao Y, Rilling G, Davies M, Marshall I. Carotid blood flow measurement accelerated by compressed sensing: Validation in healthy volunteers. Magn Reson Imaging. 2013;31:1485-91.

91. Bammer R, Hope TA, Aksoy M, Alley MT. Time-resolved 3D quantitative flow MRI of the major intracranial vessels: initial experience and comparative evaluation at 1.5T and 3.0T in combination with parallel imaging. Magn Reson Med. 2007;57:127-40.

92. Wetzel S, Meckel S, Frydrychowicz A, Bonati L, Radue EW, Scheffler K, et al. In vivo assessment and visualization of intracranial arterial hemodynamics with flow-sensitized 4D MR imaging at 3T. AJNR Am J Neuroradiol. 2007;28:433-8.

93. Boussel L, Rayz V, Martin A, Acevedo-Bolton G, Lawton MT, Higashida R, et al. Phase-contrast magnetic resonance imaging measurements in intracranial aneurysms in vivo of flow patterns, velocity fields, and wall shear stress: comparison with computational fluid dynamics. Magn Reson Med. 2009;61:409-17.

94. Isoda H, Ohkura Y, Kosugi T, Hirano M, Takeda H, Hiramatsu H, et al. In vivo hemodynamic analysis of intracranial aneurysms obtained by magnetic resonance fluid dynamics (MRFD) based on time-resolved three-dimensional phase-contrast MRI. Neuroradiology. 2010;52:921-8.

95. van Ooij P, Guedon A, Poelma C, Schneiders J, Rutten MC, Marquering HA, et al. Complex flow patterns in a real-size intracranial aneurysm phantom:

phase contrast MRI compared with particle image velocimetry and

computational fluid dynamics. NMR Biomed. 2011;25(1):14-26.

Schnell S, Ansari SA, Vakil P, Wasielewski M, Carr ML, Hurley MC, et al. Three-

dimensional hemodynamics in intracranial aneurysms: influence of size and

morphology. J Magn Reson Imaging. 2014;39:120-31.

Schrauben E, Johnson K, Huston J, Del Rio A, Reeder S, Field A, et al.

Reproducibility of cerebrospinal venous blood flow and vessel anatomy

with the use of phase contrast-vastly undersampled isotropic projection

reconstruction and contrast-enhanced MRA. Am J Neuroradiol.


Wählin A, Ambarki K, Birgander R, Wieben O, Johnson KM, Malm J, et al. Measuring pulsatile flow in cerebral arteries using 4D phase-contrast MR imaging. Am J Neuroradiol. 2013;34:1740-5.

Frydrychowicz A, Landgraf BR, Niespodzany E, Verma RW, Roldan-Alzate A, Johnson KM, et al. Four-dimensional velocity mapping of the hepatic and splanchnic vasculature with radial sampling at 3 tesla: A feasibility study in portal hypertension. J Magn Reson Imaging. 2011;34(3):577-84. Stankovic Z, Csatari Z, Deibert P, Euringer W, Blanke P, Kreisel W, et al. Normal and altered three-dimensional portal venous hemodynamics in patients with liver cirrhosis. Radiology. 2012;262:862-73. Stankovic Z, Csatari Z, Deibert P, Euringer W, Jung B, Kreisel W, et al. A feasibility study to evaluate splanchnic arterial and venous hemodynamics by flow-sensitive 4D MRI compared with Doppler ultrasound in patients with cirrhosis and controls. Eur J Gastroenterol Hepatol. 2013;25:669-75. Frydrychowicz A, Winterer JT, Zaitsev M, Jung B, Hennig J, Langer M, et al. Visualization of iliac and proximal femoral artery hemodynamics using time-resolved 3D phase contrast MRI at 3T. J Magn Reson Imaging. 2007;25:1085-92. Lum DP, Johnson KM, Paul RK, Turk AS, Consigny DW, Grinde JR, et al. Transstenotic pressure gradients: measurement in swine-retrospectively ECG-gated 3D phase-contrast MR angiography versus endovascular pressure-sensing guidewires. Radiology. 2007;245:751-60. BleyTA, Johnson KM, Francois CJ, Reeder SB, Schiebler ML BRL, Consigny D, et al. Noninvasive assessment of transstenotic pressure gradients in porcine renal artery stenoses by using vastly undersampled phase-contrast MR angiography. Radiology. 2011;261:266-73.

Frydrychowicz A, Landgraf B, Niespodzany E, Verma R, Roldan-Alzate A, Johnson K, et al. Four-dimensional velocity mapping of the hepatic and splanchnic vasculature with radial sampling at 3 tesla: A feasibility study in portal hypertension. J Magn Reson Imaging. 2011;34:577-84. Stankovic Z, Blanke P, Markl M. Usefulness of 4D MRI flow imaging to control TIPS function. Am J Gastroenterol. 2012;107:327-8. François CJ, Lum DP, Johnson KM, Landgraf BR, Bley TA, Reeder SB, et al. Renal arteries: isotropic, high-spatial-resolution, unenhanced MR angiography with three-dimensional radial phase contrast. Radiology. 2011;258:254. Hope MD, Dyverfeldt P, Acevedo-Bolton G, Wrenn J, Foster E, Tseng E, et al. Post-stenotic dilation: evaluation of ascending aortic dilation with 4D flow MR imaging. Int J Cardiol. 2012;156:e40-2.

Bissell MM, Hess AT, Biasiolli L, Glaze SJ, Loudon M, Pitcher A, et al. Aortic dilation in bicuspid aortic valve disease: flow pattern is a major contributor and differs with valve fusion type. Circ Cardiovasc Imaging. 2013;6:499-507. Mahadevia R, Barker AJ, Schnell S, Entezari P, Kansal P, Fedak PW, et al. Bicuspid Aortic Cusp Fusion Morphology Alters Aortic 3D Outflow Patterns, Wall Shear Stress and Expression of Aortopathy. Circulation. 2013;129(6):673-82. doi:10.1161/CIRCULATI0NAHA.113.003026. Uribe S, Bachler P, Valverde I, Crelier GR, Beerbaum P, Tejos C, et al. Hemodynamic assessment in patients with one-and-a-half ventricle repair revealed by four-dimensional flow magnetic resonance imaging. Pediatr Cardiol. 2013;34:447-51.

Geiger J, Hirtler D, Burk J, Stiller B, Arnold R, Jung B, et al. Postoperative

pulmonary and aortic 3D haemodynamics in patients after repair of

transposition of the great arteries. Eur Radiol. 2014;24:200-8.

Ebbers T. Flow imaging: cardiac applications of 3D cine phase-contrast MRI.

Current Cardiovascular Imaging Reports. 2011;4:127-33.

Markl M, Kilner PJ, Ebbers T. Comprehensive 4D velocity mapping of the

heart and great vessels by cardiovascular magnetic resonance. J Cardiovasc

Magn Reson. 2011;13:7.

Frydrychowicz A, Francois CJ, Turski PA. Four-dimensional phase contrast magnetic resonance angiography: Potential clinical applications. Eur J Radiol. 2011;80(1):24-35.

Markl M, Frydrychowicz A, Kozerke S, Hope M, Wieben O. 4D flow MRI. J Magn Reson Imaging. 2012;36:1015-36.

117. Hope MD, Sedlic T, Dyverfeldt P. Cardiothoracic magnetic resonance flow imaging. J Thorac Imaging. 2013;28:217-30.

118. Stankovic Z, Allen BD, Garcia J, Jarvis KB, Markl M. 4D flow imaging with MRI. Cardiovasc Diagn Ther. 2014;4:173-92.

119. Harloff A, Simon J, Brendecke S, Assefa D, Helbing T, Frydrychowicz A, et al. Complex plaques in the proximal descending aorta: an underestimated embolic source of stroke. Stroke. 2010;41:1145-50.

120. den Reijer PM, Sallee III D, van der Velden P, Zaaijer ER, Parks WJ, Ramamurthy S, et al. Hemodynamic predictors of aortic dilatation in bicuspid aortic valve by velocity-encoded cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2010;12:4.

121. Hope MD, Wrenn J, Sigovan M, Foster E, Tseng EE, Saloner D. Imaging biomarkers of aortic disease: increased growth rates with eccentric systolic flow. J Am Coll Cardiol. 2012;60:356-7.

122. Eriksson J, Dyverfeldt P, Engvall J, Bolger AF, Ebbers T, Carlhäll CJ. Quantification of presystolic blood flow organization and energetics in the human left ventricle. American Journal of Physiology-Heart and Circulatory Physiology. 2011;300:H2135-41.

123. Fredriksson AG, Zajac J, Eriksson J, Dyverfeldt P, Bolger AF, Ebbers T, et al. 4-D blood flow in the human right ventricle. American Journal of Physiology-Heart and Circulatory Physiology. 2011;301:H2344-50.

124. Hope MD, Sigovan M, Wrenn SJ, Saloner D, Dyverfeldt P. MRI hemodynamic markers of progressive bicuspid aortic valve-related aortic disease. J Magn Reson Imaging. 2013;40(1):140-5.

125. Stalder AF, Russe MF, Frydrychowicz A, Bock J, Hennig J, Markl M. Quantitative 2D and 3D phase contrast MRI: optimized analysis of blood flow and vessel wall parameters. Magn Reson Med. 2008;60:1218-31.

126. Bieging ET, Frydrychowicz A, Wentland A, Landgraf BR, Johnson KM, Wieben O, et al. In vivo three-dimensional MR wall shear stress estimation in ascending aortic dilatation. J Magn Reson Imaging. 2011;33:589-97.

127. Potters WV, Ooij P, Marquering H, vanBavel E, Nederveen AJ. Volumetric arterial wall shear stress calculation based on cine phase contrast MRI. J Magn Reson Imaging. 2015;41(2):505-16.

128. Yang GZ, Kilner PJ, Wood NB, Underwood SR, Firmin DN. Computation of flow pressure fields from magnetic resonance velocity mapping. Magn Reson Med. 1996;36:520-6.

129. Tyszka JM, Laidlaw DH, Asa JW, Silverman JM. Three-dimensional, time-resolved (4D) relative pressure mapping using magnetic resonance imaging. J Magn Reson Imaging. 2000;12:321-9.

130. Ebbers T, Wigstrom L, Bolger AF, Engvall J, Karlsson M. Estimation of relative cardiovascular pressures using time-resolved three-dimensional phase contrast MRI. Magn Reson Med. 2001;45:872-9.

131. Ebbers T, Wigstrom L, Bolger AF, Wranne B, Karlsson M. Noninvasive measurement of time-varying three-dimensional relative pressure fields within the human heart. J Biomech Eng. 2002;124:288-93.

132. Markl M, Wallis W, Brendecke S, Simon J, Frydrychowicz A, Harloff A. Estimation of global aortic pulse wave velocity by flow-sensitive 4D MRI. Magn Reson Med. 2010;63:1575-82.

133. Wentland AL, Wieben O, Francois CJ, Boncyk C, Munoz Del Rio A, Johnson KM, et al. Aortic pulse wave velocity measurements with undersampled 4D flow-sensitive MRI: comparison with 2D and algorithm determination.

J Magn Reson Imaging. 2013;37:853-9.

134. Dyverfeldt P, Sigfridsson A, Kvitting JP, Ebbers T. Quantification of intravoxel velocity standard deviation and turbulence intensity by generalizing phase-contrast MRI. Magn Reson Med. 2006;56:850-8.

135. Dyverfeldt P, Gardhagen R, Sigfridsson A, Karlsson M, Ebbers T. On MRI turbulence quantification. Magn Reson Imaging. 2009;27:913-22.

136. Dyverfeldt P, Kvitting JP, Sigfridsson A, Engvall J, Bolger AF, Ebbers T. Assessment of fluctuating velocities in disturbed cardiovascular blood flow: in vivo feasibility of generalized phase-contrast MRI. J Magn Reson Imaging. 2008;28:655-63.

137. Binter C, Knobloch V, Manka R, Sigfridsson A, Kozerke S. Bayesian multipoint velocity encoding for concurrent flow and turbulence mapping. Magn Reson Med. 2013;69:1337-45.

138. Sigovan M, Hope MD, Dyverfeldt P, Saloner D. Comparison of four-dimensional flow parameters for quantification of flow eccentricity in the ascending aorta.

J Magn Reson Imaging. 2011;34:1226-30.

139. Barker AJ, Ooij P, Bandi K, Garcia J, Albaghdadi M, McCarthy P, et al. Viscous energy loss in the presence of abnormal aortic flow. Magn Reson Med. 2013;72(3):620-8.

140. Heiberg E, Ebbers T, Wigstrom L, Karlsson M. Three-dimensional flow characterization using vector pattern matching. Visualization and Computer Graphics, IEEE Transactions on. 2003;9:313-9.

141. Töger J, Kanski M, Carlsson M, Kovacs SJ, Söderlind G, Arheden H, et al. Vortex ring formation in the left ventricle of the heart: analysis by 4D flow MRI and Lagrangian coherent structures. Ann Biomed Eng. 2012;40:2652-62.

142. Dumoulin CL, Souza SP, Walker MF, Wagle W. Three-dimensional phase contrast angiography. Magn Reson Med. 1989;9:139-49.

143. Bock J, Frydrychowicz A, Stalder AF, Bley TA, Burkhardt H, Hennig J, et al. 4D phase contrast MRI at 3 T: effect of standard and blood-pool contrast agents on SNR, PC-MRA, and blood flow visualization. Magn Reson Med. 2010;63:330-8.

144. Hess AT, Bissell MM, Ntusi NA, Lewis AJ, Tunnicliffe EM, Greiser A, et al. Aortic 4D flow: Quantification of signal-to-noise ratio as a function of field strength and contrast enhancement for 1.5 T, 3T, and 7T. Magn Reson Med. 2014;73(5):1864-71.

145. Pelc NJ, Bernstein MA, Shimakawa A, Glover GH. Encoding strategies for three-direction phase-contrast MR imaging of flow. J Magn Reson Imaging. 1991;1:405-13.

146. Dumoulin CL, Souza SP, Darrow RD, Pelc NJ, Adams WJ, Ash SA. Simultaneous acquisition of phase-contrast angiograms and stationary-tissue images with Hadamard encoding of flow-induced phase shifts.

J Magn Reson Imaging. 1991;1:399-404.

147. Johnson KM, Markl M. Improved SNR in phase contrast velocimetry with five-point balanced flow encoding. Magn Reson Med. 2010;63:349-55.

148. Walker PG, Cranney GB, Scheidegger MB, Waseleski G, Pohost GM, Yoganathan AP. Semiautomated method for noise reduction and background phase error correction in MR phase velocity data. J Magn Reson Imaging. 1993;3:521 -30.

149. Bernstein MA, Zhou XJ, Polzin JA, King KF, Ganin A, Pelc NJ, et al. Concomitant gradient terms in phase contrast MR: analysis and correction. Magn Reson Med. 1998;39:300-8.

150. Peeters JM, Bos C, Bakker CJ. Analysis and correction of gradient nonlinearity and B0 inhomogeneity related scaling errors in two-dimensional phase contrast flow measurements. Magn Reson Med. 2005;53:126-33.

151. Markl M, Bammer R, Alley MT, Elkins CJ, Draney MT, Barnett A, et al. Generalized reconstruction of phase contrast MRI: analysis and correction of the effect of gradient field distortions. Magn Reson Med. 2003;50:791 -801.

152. Bernstein MA, Grgic M, Brosnan TJ, Pelc NJ. Reconstructions of phase contrast, phased array multicoil data. Magn Reson Med. 1994;32:330-4.

153. Xiang QS. Temporal phase unwrapping for CINE velocity imaging. J Magn Reson Imaging. 1995;5:529-34.

154. Lorenz R, Bock J, KorvinkJ, Markl M. Importance of different correction methods for optimized 3D visualization of 3-directional MR velocity data. Proceedings 18th Scientific Meeting International Society for Magnetic Resonance in Medicine. 2010.

155. Gatehouse PD, Rolf MP, Graves MJ, Hofman MB, Totman J, Werner B, et al. Flow measurement by cardiovascular magnetic resonance: a multi-centre multi-vendor study of background phase offset errors that can compromise the accuracy of derived regurgitant or shunt flow measurements. J Cardiovasc Magn Reson. 2010;12:1-8.

156. Busch J, Vannesjo SJ, Barmet C, Pruessmann KP, Kozerke S. Analysis of temperature dependence of background phase errors in phase-contrast cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2014;16:97.

157. Dumoulin CL. Phase contrast MR angiography techniques. Magn Reson Imaging Clin N Am. 1995;3:399-411.

158. Markl M, Harloff A, Bley TA, Zaitsev M, Jung B, Weigang E, et al. Time-resolved 3D MR velocity mapping at 3T: improved navigator-gated assessment of vascular anatomy and blood flow. J Magn Reson Imaging. 2007;25:824-31.

159. Kheradvar A, Ebbers T, Tonti G, Fraser AG, Narula J. Emerging Trends in CV Flow Visualization. JACC Cardiovasc Imaging. 2012;5(3):305-16.

160. Rodriguez Munoz D, Markl M, Moya Mur JL, Barker A, Fernandez-Golfin C, Lancellotti P, et al. Intracardiac flow visualization: current status and future directions. European heart journal cardiovascular Imaging. 2013;14:1029-38.

161. Shadden SC, Arzani A. Lagrangian Postprocessing of Computational Hemodynamics. Annals ofbiomedical engineering. 2015;43(1):41-58.

162. Hsiao A, Tariq U, Alley MT, Lustig M, Vasanawala SS. Inlet and outlet valve flow and regurgitant volume may be directly and reliably quantified with accelerated, volumetric phase-contrast MRI. J Magn Reson Imaging. 2014;41(2):376-85.

163. Calkoen EE, Roest AA, van der Geest RJ, de Roos A, Westenberg JJ. Cardiovascular function and flow by 4-dimensional magnetic resonance imaging techniques: new applications. J Thorac Imaging. 2014;29:185-96.

164. Harloff A, Zech T, Wegent F, Strecker C, Weiller C, Markl M. Comparison of blood flow velocity quantification by 4D flow MR imaging with ultrasound at the carotid bifurcation. AJNR Am J Neuroradiol. 2013;34:1407-13.

165. Hofman MBM, Visser FC, Van Rossum AC, Vink GQM, Sprenger M, Westerhof N. In vivo validation of magnetic resonance blood volume flow measurements with limited spatial resolution in small vessels. Magn Reson Med. 1995;33:778-84.

166. van Ooij P, Potters WV, Collins J, Carr M, Carr J, Malaisrie SC, Fedak PW, McCarthy PM, Markl M, Barker AJ. Characterization of Abnormal Wall Shear Stress Using 4D Flow MRI in Human Bicuspid Aortopathy. Ann Biomed Eng. 2015;43(6):1385-97.

167. BockJ, Frydrychowicz A, Lorenz R, Hirtler D, Barker AJ, Johnson KM, et al. In vivo noninvasive 4D pressure difference mapping in the human aorta: Phantom comparison and application in healthy volunteers and patients. Magn Reson Med. 2011;66:1079-88.

168. Dyverfeldt P, Ebbers T, Länne T. Pulse Wave Velocity with 4D Flow MRI: Systematic Differences and Age-Related Regional Vascular Stiffness. Magn Reson Imaging. 2014;32(10):1266-71.

169. McGibney G, Smith M, Nichols S, Crawley A. Quantitative evaluation of several partial Fourier reconstruction algorithms used in MRI. Magn Reson Med. 1993;30:51 -9.

170. Chao TC, Chung HW, Hoge WS, Madore B. A 2D MTF approach to evaluate and guide dynamic imaging developments. Magn Reson Med. 2010;63:407-18.

171. Blaimer M, Ponce IP, Breuer FA, Jakob PM, Griswold MA, Kellman P. Temporal filtering effects in dynamic parallel MRI. Magn Reson Med. 2011;66:192-8.

172. Ramb R, Binter C, Schultz G, Assländer J, Breuer F, Zaitsev M, Kozerke S, Jung B. A g-factor metric for kt-GRAPPA-and PEAK-GRAPPA-based based parallel imaging. Magnetic Resonance in Medicine. 2015;74(1):125-35.

173. Petersson S, Dyverfeldt P, Ebbers T. Assessment of the accuracy of MRI wall shear stress estimation using numerical simulations. J Magn Reson Imaging. 2012;36:128-38.

174. Petersson S, Dyverfeldt P, Gärdhagen R, Karlsson M, Ebbers T. Simulation of phase contrast MRI of turbulent flow. Magn Reson Med. 2010;64:1039-46.

175. Jiang J, Johnson K, Valen-Sendstad K, Mardal K-A, Wieben O, Strother C. Flow characteristics in a canine aneurysm model: a comparison of 4D accelerated phase-contrast MR measurements and computational fluid dynamics simulations. Med Phys. 2011;38:6300-12.

176. Knobloch V, Binter C, Gülan U, Sigfridsson A, Holzner M, Lüthi B, et al. Mapping mean and fluctuating velocities by Bayesian multipoint MR velocity encoding-validation against 3D particle tracking velocimetry. Magn Reson Med. 2014;71:1405-15.

177. Dyverfeldt P, Sigfridsson A, Knutsson H, Ebbers T. A novel MRI framework for the quantification of any moment of arbitrary velocity distributions. Magn Reson Med. 2011;65:725-31.

178. Zwart NR, Pipe JG. Multidirectional high-moment encoding in phase contrast MRI. Magn Reson Med. 2013;69:1553-63.

179. Kellman P, McVeigh ER. Image reconstruction in SNR units: A general method for SNR measurement!. Magn Reson Med. 2005;54:1439-47.

180. Robson PM, Grant AK, Madhuranthakam AJ, Lattanzi R, Sodickson DK, McKenzie CA. Comprehensive quantification of signal-to-noise ratio and g-factor for image-based and k-space-based parallel imaging reconstructions. Magn Reson Med. 2008;60:895-907.

181. Giese D, Haeberlin M, Barmet C, Pruessmann KP, Schaeffter T, Kozerke S. Analysis and correction of background velocity offsets in phase-contrast flow measurements using magnetic field monitoring. Magn Reson Med. 2012;67:1294-302.

182. Gu T, Korosec FR, Block WF, Fain SB, Turk Q, Lum D, et al. PC VIPR: a high-speed 3D phase-contrast method for flow quantification and high-resolution angiography. AJNR Am J Neuroradiol. 2005;26:743-9.

183. Sigfridsson A, Petersson S, Carlhäll CJ, Ebbers T. Four-dimensional flow MRI using spiral acquisition. Magn Reson Med. 2012;68:1065-73.

184. Kecskemeti S, Johnson K, Wu Y, Mistretta C, Turski P, Wieben O. High resolution three-dimensional cine phase contrast MRI of small intracranial aneurysms using a stack of stars k-space trajectory. J Magn Reson Imaging. 2012;35:518-27.

185. Rolf MP, Hofman M, Kuijer J, van Rossum AC, Heethaar RM. 3D velocity quantification in the heart: Improvements by 3D PC-SSFP. J Magn Reson Imaging. 2009;30:947-55.

Nilsson A, Revstedt J, Heiberg E, Stählberg F, Bloch KM. Volumetric velocity measurements in restricted geometries using spiral sampling: a phantom study. Magnetic Resonance Materials in Physics, Biology and Medicine. 2015;28(2):103-18.

Baltes C, Kozerke S, Hansen MS, Pruessmann KP, Tsao J, Boesiger P.

Accelerating cine phase-contrast flow measurements using k-t BLAST and k-t

SENSE. Magn Reson Med. 2005;54:1430-8.

Knobloch V, Boesiger P, Kozerke S. Sparsity transform k-t principal

component analysis for accelerating cine three-dimensional flow

measurements. Magn Reson Med. 2013;70:53-63.

Liu J, Dyverfeldt P, Acevedo-Bolton G, Hope M, Saloner D. Highly

accelerated aortic 4D flow MR imaging with variable-density random

undersampling. Magn Reson Imaging. 2014;32:1012-20.

Jung B, Honal M, Ullmann P, Hennig J, Markl M. Highly k-t-space-accelerated

phase-contrast MRI. Magn Reson Med. 2008;60:1169-77.

Lustig M, Donoho D, Pauly JM. Sparse MRI: The application of compressed

sensing for rapid MR imaging. Magn Reson Med. 2007;58:1182-95.

Tsao J, Boesiger P, Pruessmann KP. k-t BLAST and k-t SENSE: Dynamic MRI

with high frame rate exploiting spatiotemporal correlations. Magn Reson

Med. 2003;50:1031-42.

van Ooij P, Potters WV, Collins J, Carr M, Carr J, Malaisrie SC, Fedak PW, McCarthy PM, Markl M, Barker AJ. Characterization of Abnormal Wall Shear Stress Using 4D Flow MRI in Human Bicuspid Aortopathy. Annals of biomedical engineering 2014, published online.

Submit your next manuscript to BioMed Central and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at