The study, conducted by Thomas Murchie, tackles the question with analytical mathematics rather than new animal or human experiments. Its starting point is the three-element Windkessel model, in which the arterial system is represented by a peripheral resistance, a characteristic impedance of the proximal aorta, and a compliance that stores blood during systole and releases it during diastole. In that framework, Ea depends on total systemic resistance, on arterial compliance through the time constant tau, and critically on the durations of systole and diastole. Clinicians, however, rarely measure Ea that way. Instead they use a convenient surrogate, the ratio of end-systolic pressure to stroke volume, which tracks the Windkessel value closely in living subjects but hides the separate contributions of resistance and compliance beneath a single ratio.

Over the years, several laboratories have proposed a compromise: a linear approximation of the form Ea equals R plus k divided by C, where R is a resistance term, C is compliance, and k is a coefficient that must be fitted to data. Empirical studies produced values of k around 0.314 in work by Segers and colleagues in 2002 and 0.42 in work by Chemla and colleagues in 2003. More provocatively, other analyses concluded that compliance contributes almost nothing at all, effectively setting k near zero and reducing Ea to roughly total resistance divided by cycle length. The physiological determinants of k were never described analytically, leaving the field with a scatter of numbers and no principled way to choose among them.

The new paper resolves the ambiguity with a higher-order approximation. By replacing the exponential term in the Windkessel equation with a polynomial expansion and carrying the algebra one step beyond the truncation that yields the k-equals-zero result, Murchie derives a closed-form expression: Ea is approximately total resistance divided by the cardiac cycle length, plus a compliance term whose coefficient k equals one half multiplied by the square of the ratio of diastolic duration to total cycle duration. In other words, k is not a universal constant at all. It is a timing factor, growing as diastole occupies a larger share of the heartbeat, which happens characteristically when the heart rate is low. At high heart rates, where diastole is squeezed short, the compliance term shrinks toward irrelevance, which is precisely why some studies concluded that compliance barely matters.

To test the approximation, the analysis first compared it against the parent Windkessel equation across a deliberately broad parameter space: 693 evenly distributed combinations of total resistance, compliance, and heart rate, spanning values seen in healthy adults and in patients treated for septic and cardiogenic shock. Heart rates of 60, 80, and 100 beats per minute were used, with ejection times adjusted according to the relationship observed during exercise by Mertens and colleagues. The agreement was essentially exact for practical purposes, with a coefficient of determination above 0.9999 and a mean absolute percentage error of just 0.51 percent. The worst-case error, arising only when heart rate, resistance, and compliance were simultaneously low, was 3.96 percent, a deviation traced to the limits of the truncated Taylor series at large values of the ratio of diastolic duration to the Windkessel time constant.

The second validation was more demanding, because it pitted the formula against a completely independent model of the circulation. Using a publicly available database of 4,374 simulated adults generated by Charlton and colleagues with a one-dimensional vascular model reproducing age-appropriate cardiovascular distributions from 25 to 75 years, the study compared the new approximation against the clinical surrogate of end-systolic pressure divided by stroke volume, sampled carefully at the end of the small end-systolic flow reversal in the aortic root. The new timing-dependent formula achieved a coefficient of determination of 0.9980 with a mean absolute percentage error of 0.87 percent. Segers’ fixed-coefficient regression, by contrast, yielded 0.9903 with errors up to nearly 15 percent, and the k-equals-zero simplification fell to 0.9610 with maximum errors above 25 percent. The ordering was unambiguous: retaining compliance with a timing-scaled coefficient beats ignoring it outright.

Beyond the curve-fitting, the paper offers a physical intuition for what the coefficient k actually represents. Since Ea approximates end-systolic pressure over stroke volume, and resistance over cycle length approximates mean arterial pressure over stroke volume, the compliance term corresponds to the excess of end-systolic pressure over mean pressure, per unit of stroke volume. That excess depends on two timing-dependent factors. The first is the fraction of the stroke volume still stored in the arterial compartment at the end of systole rather than already run off through the peripheral resistance, which vanishes as diastole shrinks. The second is the fraction of the resulting arterial pressure change that sits above mean arterial pressure, which depends on how the arterial pressure decay is weighted in time. Both factors converge to zero as diastole occupies an ever-smaller slice of the cycle, which is exactly the behavior captured by the squared ratio of diastolic to total cycle duration.

The implications reach beyond theory into the design of empirical studies. Because the ratio of diastolic duration to cycle length varies substantially with heart rate, the apparent relationship between compliance and Ea measured across a population with natural heart-rate variation may not reflect the true contribution of compliance within any single individual. If diastolic fraction happens to covary systematically with compliance across a cohort, the compliance-dependent term could appear nearly constant even though Ea physically depends on compliance in every subject. That mechanism, the paper argues, offers a coherent explanation for the divergent empirical estimates of k that have accumulated since the early 2000s: analyses that reduce Ea to resistance over cycle length and analyses that fit a constant k are simply working at different orders of approximation of the same underlying exponential expression.

The author is careful to delineate the limits of the work. The analysis addresses the interpretation of Windkessel parameters in the systemic circulation only, does not tackle the uncertainty inherent in estimating those parameters from clinical data, excludes the pulmonary circulation, and does not examine how the derived load expressions relate to pressure-flow behavior in living vessels. No new human or animal data were collected, so ethics approval was not required. Even so, the central message is likely to resonate widely among physiologists and clinicians who use Ea to gauge ventricular afterload: the stiffness of the arteries and the rhythm of the heart are not independent contributors to cardiac load but are woven together by the clock of the cardiac cycle, and any formula that ignores that coupling risks misreading how hard the heart truly works.

Subject of Research: The dependence of effective arterial elastance on arterial compliance and cardiac cycle timing in a Windkessel model of the systemic circulation

Article Title: Cardiac timing modulates the contribution of compliance to effective arterial elastance

Article References: Murchie, T. (2026). Cardiac timing modulates the contribution of compliance to effective arterial elastance. Physiological Reports, 14(19), Article e71112. https://doi.org/10.14814/phy2.71112

Image Credits: AI Generated

DOI: 10.14814/phy2.71112

Keywords: effective arterial elastance, arterial compliance, Windkessel model, cardiac timing, diastole, heart rate, ventricular afterload, end-systolic pressure, stroke volume, hemodynamics, cardiovascular physiology, mathematical modeling