Stroke volume (SV) is defined as that quantity of blood ejected from the cardiac ventricles for each heart beat. More specifically, it relates to the volume of blood ejected from opening to closure of the semilunar valves, and in the case of the left ventricle, the aortic valve. The time interval over which left ventricular ejection occurs is known as left ventricular ejection time (T_{LVE}, s). For purposes of simplicity, and considering this is not an in-depth analytical treatise on pulsatile blood flow, per se, the following simplistic model is proposed as an operational tool for understanding the opposing hypotheses concerning SV determination by means of the transthoracic electrical bioimpedance technique (TEB).

In geometric terms, consider stroke volume V (mL) a cylinder of length S (cm) and cross-sectional area π^{2} (CSA, cm^{2}). In this model it is assumed that SV is surrounded by a thin-walled viscoelastic aorta, equal in length S to the SV contained within at end-systole. Hemodynamically, as a function of ventricular ejection, let S (cm) also represent the time velocity integral, also known as stroke distance (cm).

If the dependent variable, V, and independent variables r and S are continuously differentiable functions of time t, the rate of change of aortic volume dV/dt or flow Q(t) (mL·s^{–1}) is given as,

If 2π_{ISA} (cm^{2}) of an ascending aortic segment and dr(t)/dt aortic wall velocity (cm·s^{–1}), then

NB: Q(t) in lieu of

For derivative 2, if π^{2} is the CSA of the aorta at a discreet point of measurement and dS(t)/dt the axial blood velocity (^{–1}) measured at that point, then axial Q can be given as,

Inasmuch as A_{CSA} changes continuously over the cardiac cycle as a function of aortic valve radius dr, the following is an animated representation of derivative 2 of equation 2 and equation 4 [1].

At an operational level, derivative 2 of equation 2 and equations 4 and 5 represent the basis for SV determination by means of both the Doppler/echo method and electromagnetic flowmetry. To find SV simply requires integration of the velocity profile over the ejection period at a discreet point in the aortic root and multiplying this by the aortic root CSA at the point of velocity measurement. Common sites for measurement include the left ventricular outflow tract and aortic valve leaflets. Thus,

where _{0} to t_{1}, T_{LVE}. Extrapolating derivative 1 of equation 2 and equation 3 for determination of SV is something very much more complex, for it represents the geometric foundation of the _{s}. Upon the force of ventricular ejection in early systole, and because of distal vascular hindrance, the proximal aorta expands, storing a part of the inflow (Q(t)_{in}) as potential energy (Q(t)_{stored}). Simultaneously, due to open outflow, blood is expelled axially (Q(t)_{out}) to the periphery, where Q(t)_{out} is commonly referred to as “runoff”. When aortic pressure exceeds left ventricular pressure from mid-systole to aortic valve closure, the recoiling aorta converts potential energy to kinetic energy, expelling the stored blood to the periphery. From aortic valve closure (end-systole) to end-diastole, non-pulsatile flow dissipates along a quasi-mono-exponential pressure decay. Pulsatile blood flow is therefore buffered and converted from periodic oscillating flow to near-continuous flow for each cardiac cycle [2]. It is this aforementioned physical phenomenon that is the basis for conceptualization of the plethysmographic-based ICG-derived SV equations, and the mechanics underlying the ^{2}·s^{–1}), then Q(t) for derivative 1 can be given as,

In terms of mechanics, since dA(t) and dA(t)/dt are non-linear functions of vessel material composition, pressure P(t) (mmHg) and rate of change of aortic pressure dP(t)/dt (mmHg·s^{–1}), respectively, the following pertains [3]:

where the right hand side (rhs) of equation 8 defines the magnitude and shape of the flow wave, and C is the aortic compliance.

If aortic dP(t) is a function of force F and area A, both changing with time over the ejection interval, dP(t) (mmHg) and dP(t)/dt (mmHg ⋅ s^{–1}) are given as,

respectively, and compliance C is given as,

where dA/dP is a non-linear function of Pulse pressure (ΔP) and mean distending pressure P_{m}. Thus,

where the mechanical expression for dV/dt of equation 11 is equivalent to the geometrical definition for dV/dt of derivative 1 in equation 2. As an analytical expansion of equation 8, the rhs of equation 11 more precisely defines the temporal landmarks, shape and magnitude of the radially oriented flow wave.

Figure 1 shows the pressure curve P(t), the rate of change of aortic pressure, dP/dt, aortic blood flow (Q(t)) and the rate of change of flow, which is acceleration of flow, dQ(t)/dt. It is clear from equation 11 that dF/dt_{max}, the major derivative, will occur early in systole, when dA/dt is trivial. Likewise, dA/dt_{max} will occur when dF/dt = 0. When dP/dt is modulated by C, dA/dt_{max} will occur at the peak magnitude of the flow wave dV/dt_{max} (i.e. Q_{max}) and dF/dt_{max} will occur on the steepest portion of the flow wave, immediately following aortic valve opening in early systole.

When flow is differentiated with respect to time dQ/dt_{max} will occur when dA/dt and dV/dt are minimal. Also apparent is that when flow acceleration = 0, dV/dt is at maximum. Therefore, dQ(t)/dt_{max} and dV/dt_{max} are monotonically restricted to their respective time-domains, with acceleration of blood flow, dQ/dt_{max} peaking, on the mean, approximately 50ms before peak flow, dV/dt_{max}.

If S is a segment of aorta with closed outflow, then flow into the segment would be the total input flow from left ventricular ejection. If, however, segment S is an open outflow conduit, equation 8 would represent a net flow. For a simple 2-element _{in}) is given as [4],

where derivatives 1 and 2 on the rhs of (12) represent Q(t)_{stored} and Q(t)_{out}, respectively.

To obtain total flow Q(t)_{total} through the segment over a single systolic interval (T_{LVE}, s), Q(t)_{out} must be measured. With regard to a two-element _{out} is obtained by dividing the pressure change over the ejection interval P(t) by the distal, or systemic vascular resistance, R_{s}. R_{s} is defined as the quotient of mean pressure to mean flow P_{m}/Q_{m}. Q(t)_{total} is thus given as [4],

Integration of Q(t)_{strored} yields stored volume inflow, SV_{stored}:

where ΔP is the pulse pressure, which is aortic peak systolic pressure P_{s} minus aortic diastolic pressure P_{d} (P_{s}–P_{d}) [5].

Since SV = CΔP, where C is a non-constant time variable function of pressure P, the following expression also represents total flow into a segment of aorta of length S and internal surface area wall velocity dA/dt [6].

where derivatives 1 and 2 represent Q(t)_{stored} and Q(t)_{out}, respectively. Differentiating dC_{p}(t)/dt results in,

In terms of total flow into segment S, equations 12 and 15 are equivalent. Integration of equation 13 over an ejection interval yields SV and is thus given as [2],

where equation 17 represents a two-element _{stored} and V_{out}, respectively.

Although a theoretically correct model for SV, equation 17 is not easily solved. Clearly, P(t), dP(t)/dt, and P_{m} can be computed by direct measurement and signal processing, but C and R_{s} must be obtained independently by alternative means [2]. This is clear from equation 10 for the determination of C and the relationship of R_{s} to Q(t)_{mean} in equation 17. From equations 15 and 16, it is also clear that C is not a constant as equation 17 implies. Pressure dependent C and the rate of change of C with respect to time, dC/dt, are not accounted for in equation 17. To more accurately predict SV from the _{c}. Z_{c} is approximated by aortic pulse pressure divided by peak aortic blood flow (ΔP/Q_{max}). Incorporating characteristic aortic impedance into the

Impedance cardiography (ICG) is a branch of bioimpedance primarily concerned with the determination of left ventricular stroke volume (SV) [7–10]. It was the first and still the only truly noninvasive, continuous, operator-independent, hands-off, beat-to-beat method used in clinical practice. Despite these desirable attributes, general acceptance of the method as a replacement for invasive reference standards in critically-ill humans has not been realized. As summarized by Raaijmakers et al. and Moshkovitz et al., both correlation and agreement against standard reference methods have been far too variable to reliably apply the method clinically for modulation of patient therapy [7,12]. As emphasized by Kauppinen et al. [11], the major problem in cardiovascular impedance measurements is the inability to accurately correlate the magnitude of the impedance waveforms to the hemodynamic variables they pretend to mimic. The lack of consistently high correlation and agreement between ICG measurements and reference standards would suggest that, the equation models describing SV in ICG are not sufficiently robust and the physical acquisition of the impedance signals not sensitive and specific enough for clinical application in the critically ill [11].

While measurements can be obtained by the whole body technique, the method is usually implemented by means of the transthoracic approach (i.e. transthoracic electrical bioimpedance cardiography, TEB, ICG) [7]. The latter technique, and subject of this review, involves applying a current field longitudinally across a segment L (cm) of thorax by means of a constant magnitude, high frequency (50–100 kHz), low amplitude AC (I(t)) (1–4 mA). By Ohm’s Law, the voltage difference U(t) (volt, U) measured within the current field is proportional to the transthoracic electrical impedance Z.

where Z is the frequency-dependent AC analog of the static DC resistance R. The magnitude of the impedance, or modulus, │Z│, comprises a resistive component, also known as the real part, R, and an imaginary part X. X is further comprised of an inductive part (X_{L}) and capacitive part (X_{C}), both collectively known as the reactance. With regard to the magnitude of X, X_{L} is proportional and X_{C} inversely proportional (1/f) to the frequency of the applied AC. │Z│ is determined by the geometric addition of the resistive and reactive components.

where the reactance X is,

Thus,

where the magnitudes of R and X on Cartesian coordinates are determined by the phase angle φ, and thus, R

For the remainder of the discussion, let

For cylindrical electrical conductors, and as a companion equation to equation 19 [13],

where ρ = the resistivity (Ω ⋅ cm), L = the length of the conductor, and A and V are the CSA (cm^{2}) and volume (mL) of the conductor, respectively. Thus, substituting ρL^{2}/V in equation 23 for Z in equation (19) and rearranging,

For the remainder of the discussion, the character “U” is used in

If the thorax is considered a cylindrical bulk electrical conductor of length L between the voltage sensing electrodes, with ρT the transthoracic specific resistance and V_{T} the volume of thorax, then, without ventilatory or cardiac activity, Z in equations 19 and 23 represents the adynamic or static transthoracic base impedance Z_{0}. When cardiac and ventilatory activity are superimposed on Z_{0}, a time-variable transthoracic impedance is registered, Z(t). By eliminating the oscillating cardiac-asynchronous ventilatory component ΔZ_{vent}, Z(t) comprises, in parallel, a static DC component Z_{0} (22Ω–45Ω) and a dynamic AC component ΔZ(t) (0.1Ω–0.2Ω) [14,15] (Fig.2). It should be noted that the magnitude of Z_{0} not only varies between individuals and the frequency of the applied AC, but also with the electrode configuration used for signal acquisition. For four common electrode configurations a computer model predicts a Z_{0} range of 26.3Ω–34.3Ω [11].

When electrical resistances or impedances are added in parallel, they are summed as their reciprocals.

As a composite impedance with tissues of widely varying specific resistances, Z_{0} comprises, in parallel, a static, relatively non-conductive multi-compartmental tissue impedance Z_{t} (400Ω·cm–1020Ω·cm) [11,16], a highly conductive blood impedance Z_{b} (or resistance R_{b} 100Ω·cm–180Ω·cm) varying with hematocrit [17], and a very highly conductive interstitial extra-vascular lung water (EVLW) impedance Z_{e} (60Ω·cm–70Ω·cm).

Z_{0} is thus an intensity weighted mean of the reciprocal sum of all thoracic tissue impedances, the sum of the reciprocals being less than the lowest tissue impedance of the thorax.

In this model, the blood resistance R_{b} is considered a cylindrical tube of constant length L, surrounded by the highly conductive extravascular lung water impedance Z_{e} and both surrounded by a non-conductive thoracic encompassing cylinder Z_{t}. It should be noted that, at the frequencies used in ICG, blood is almost purely resistive with a trivial reactive component. Therefore, the term blood resistance, R_{b}, is justifiably used in lieu of blood impedance Z_{b}.

_{b}(t): The Transthoracic Cardiogenic Impedance Pulse Variation

The systolic portion of the cardiogenically-induced impedance pulse variation ΔZ(t), hereafter known as ΔZ_{b}(t), comprises two components, arguably of equal magnitude [14,18–21]. The first is a velocity-induced change in the specific resistance of axially directed flowing blood (Δρ_{b}(t), Ω·cm·s^{–1}), where, in end-diastole, the state of highest blood resistivity, red blood cells are randomly oriented. By contrast, during peak aortic blood acceleration and the rapid ejection phase of early systole, red blood cells become deformed and assume a well defined state of parallel orientation along their long axis of symmetry [19,22,23] (fig. 3).

Parallel alignment opens clear current pathways through the highly conductive plasma, causing a decrease in blood resistivity and transthoracic resistivity [19,22] (i.e.↑Conductivity). The rapidity of and degree to which attainment of complete parallel orientation is achieved determines the magnitude of the resistivity change Δρ_{b}(t)_{max} (Ω·cm·s^{–1}) [23]. The maximum resistivity change occurs when the long axis of the red cell is oriented within 20^{0} to the direction of blood flow [23]. Figure 4 shows that, in early ejection, maximum acceleration of red cell reduced average velocity (i.e. mean spatial velocity) parallels the change in blood resistivity (i.e. conductivity) (r = 0.99). By comparison, upon negative acceleration, the resistivity change does not parallel the reduced average velocity, with subsequent delay in reaching baseline [19,23]. This disparity is, no doubt, a consequence of the red cells’ inability to achieve complete randomization at end-systole.

The second component generating a change in transthoracic specific resistance is the transversely or laterally-oriented volume displacement of non-conductive alveolar gas (ρ =10^{20} Ω·cm) by stroke volume-induced expansion of the ascending aorta, principally [24,25], with highly conductive blood (ρ_{b}=100–180 Ω·cm) ΔV_{b}(t) (ΔΩ·cm(t)). As discussed under pulsatile blood flow, in addition to distal vascular hindrance R_{s}, the volumetric expansion of the aorta is due to pressure-induced, compliance modulated changes in cross sectional area, π^{r2} (cm^{2}), and, over segment L, a change in internal surface area, 2π·(dA/dr)L, (cm^{2}) [1,26]. For example, assume that equation 23 represents the impedance Z of a segment of very thin walled aorta embedded within the thorax. Assume that the aortic segment is of length L with open outflow, the content of which is blood of static specific resistance ρ_{b} and volume V_{b}. Thus, an increase in aortic volume ΔV_{b}(t) would result in a corresponding decrease in vessel and transthoracic impedance ΔZ_{b}(t)_{volume}.

where the rhs of equation 27 is the ohmic equivalent of equation 14, wherein ΔZ_{b}(t) and ΔV_{b}(t) represent net changes. The combined velocity and volume effects cause a steep drop in transthoracic impedance.

For the velocity-induced change in blood resistivity, the decrease in vessel and transthoracic impedance is given as,

where Δρ_{b}(t) is of diminishing value upon ejection. Combining equations 27 and 28,

Combining equations 24, 26, and 29, the following impedance model describes the static and dynamic impedances in parallel, when an AC field is applied longitudinally across the thorax (fig 2).

The original Nyboer equation [27] was specifically proposed for determination of segmental blood volume changes in the upper and lower extremities. It is based upon the assumption that the arteries of the extremities are rigidly encased in muscle and connective tissue, and that both the arteries and surrounding tissue can be approximated as cylindrical conductors placed in parallel alignment. Nyboer found that, by placing spaced-apart circumferential current injecting electrodes on a limb segment and voltage sensing electrodes within the current field proximate the current injectors, impedance changes proportional to strain gauge determined blood volume changes could be measured. He also determined that, in order to compensate for simultaneous volume outflow from the limb segment during arterial inflow, venous outflow occlusion was necessary. Thus, the ΔZ_{b}(t) waveform represents the sum of two signals: one caused by inflow and the other caused by outflow (runoff) between the voltage sensing electrodes [28]. Using venous outflow obstruction, absolute volume changes during each pressure pulse were accurately approximated. Nyboer referred to this method, and quite accurately as, electrical impedance plethysmography. The foundation of Nyboer’s equation, and for that matter the earliest SV equations, originate from equation 27, where, when solving for ΔV_{b}(t)_{net},

However, since a large extremity artery (i.e. brachial or femoral) is encased in muscle and connective tissue, ΔZ_{b}(t) had to be further refined, reflecting a parallel connection of the static tissue and dynamic blood impedances. Since the magnitude of Z(t) of an extremity is larger than extremity Z_{0} by a trivial factor of extremity ΔZ_{b}(t), the approximation that extremity Z(t) and Z_{0} are virtually identical is a plausible assumption. Without assumption, it also follows that Z(t)–Z_{0} = ΔZ_{b}(t). Thus, if equation 25 is solved for ΔZ_{b}(t), using the reciprocal rule for impedances added in parallel:

Substituting the right hand side of equation 32 for ΔZ_{b}(t) into equation 31 results in the operational model for measurement of blood volume changes in an extremity. For the maximum volume change ΔV_{b(max)} with venous outflow occlusion [28],

where ΔZ_{max} is the peak magnitude of the waveform and ρ_{b}L^{2}/Z_{0}^{2} is a constant. For an arterial conduit with closed outflow, equation 33 is the bioelectric equivalent of the hemodynamic expression given for SV in equation 14 with outflow obstruction. For an extremity vessel with open outflow, equation 33 is analogous to the integrated value for equation 12. Clearly, as applied to SV calculations, equation 33 is inadequate, because complete aortic outflow obstruction for calculation of total SV is impossible.

_{b}(t)

As applied to thoracic measurements of SV, Nyboer professed to have solved the outflow problem by manually determining the maximum systolic down-slope of ΔZ_{b}(t) and extrapolating it backward to the beginning of ejection [28]. Thus, the maximum impedance change resulting from the backward extrapolation method was believed to be equivalent to the maximum impedance change attained as if no arterial runoff occurred during ejection.

Thus, SV, or the maximum volume change ΔV_{b} measured between the voltage sensing electrodes is given as,

where Z_{0} is the static transthoracic base impedance and ΔZ_{max} is the maximum backward extrapolated value of the thoracic cardiogenic impedance pulse variation, which includes runoff. While credible SV values could be obtained with equation 34, it was not widely accepted, this being due to difficulty in manually determining the true maximum downslope.

_{b}(t)

As a consequence of the deficiencies of the backward extrapolation procedure, Kubicek et al. proposed a maximum systolic forward extrapolation of ΔZ_{b}(t), which is the basis for all subsequent plethysmographic conceptually-based SV equations using the transthoracic approach [29,30]. Kubicek et al. made the assumption that, if the maximum systolic upslope (i.e. ΔZ′ = ΔZ·s^{–1}, Ω·s^{–1}) is held constant throughout ejection, then compensation for outflow, before and after attainment of Q_{max} is achieved. Theory underlying Kubicek’s method assumes that little arterial runoff occurs during the inertial phase of rapid systolic ejection. Unlike Nyboer’s direct measurement of ΔZ_{max}, Kubicek’s method requires multiplying ΔZ′ (forward extrapolation) by left ventricular ejection time T_{LVE} (i.e. ΔZ′ × T_{LVE} = ΔZ_{max}). In the original description of the technique, T_{LVE} was determined by means of phono-cardiography. As opposed to manual extrapolation of ΔZ′, electronic differentiation of ΔZ_{b}(t) was implemented. The first time-derivative clearly defines dZ/dt_{max}, with the presumed units of Ω·s^{–1}, as a distinct point C and left ventricular ejection time (s) as the interval between points B (aortic valve opening) and X (aortic valve closure) (fig. 5). Thus, ΔZ_{max} = dZ/dt_{max} × T_{LVE} = Ω. Kubicek’s equation is thus obtained by substituting dZ/dt_{max} × T_{LVE} for ΔZ_{max} in equation 34. The Kubicek equation is given below and purportedly equivalent to equation 17. Thus [29,30],

where ρ_{b} is the static specific resistance of blood, which was initially fixed at 150 Ω·cm [29]. As concerns the appropriate value for ρ, Quail et al [31] rearranged the Kubicek equation (assuming it correct), solved for ρ as the dependent variable, and measured SV by EMF. By means of normovolemic hemodilution exchange transfusion, they showed that, over a wide range of hematocrit from 26%–66%, the value of ρ remained virtually constant about a mean of 135Ω·cm at hematocrit 40%. L (cm) is the longitudinal distance between the voltage sensing electrodes on the base of the neck and those on the lower thorax at the level of the xiphoid process. Z_{0} is the quasi-static base impedance measured between the voltage sensing electrodes. In an in vitro expansible tube model, when equation 35 is compared to equation 33 without outflow occlusion, equation 33 systematically underestimates SV obtained from 35 [26]. By elimination of a Z_{0} term, the volume conductor V_{C} in equation 35 is given as,

where, when Z_{0} varies with alterations in extravascular lung water, the V_{C} is non-constant. Equation 35 is, in theory, modeled as the ohmic equivalent of _{c}) and systemic vascular resistance R_{s} into the overall theoretical gestalt.

The transthoracic impedance Z(t) is considered the parallel connection of an aggregate of static cylindrical tissue impedances, Z_{t}, considered as one, and, encased within, a dynamic blood resistance, R_{b}, otherwise known as the blood impedance, Z_{b}.

The blood resistance is considered a homogeneously conducting blood-filled cylinder of constant length L, or a parallel connection of an aggregate of cylinders considered as one.

The current distribution in the blood resistance is uniform.

All current flows through the blood resistance.

The volume conductor V_{C} is homogeneously perfused with blood of specific resistance ρ_{b}.

The magnitude of stroke volume is directly proportional to power functions of measured distance L between the voltage sensing electrodes [29–30], or to height-based thoracic length equivalents [32].

All pulsatile impedance changes ΔZ_{b}(t) are due to vessel volume changes ΔV_{b}(t), and, in the context of assumption 2, ΔZ_{b}(t) is due exclusively to changes in vessel radius dr(t) and CSA dA(t).

In the context of assumption 7, dZ/dt_{max} is the bioelectric equivalent of dV/dt_{max},

Outflow, or runoff during ventricular ejection, can be compensated for by extrapolating the peak rate of change of ΔZ_{b}(t) over the ejection interval (i.e. dZ/dt_{max} × T_{LVE}) to obtain ΔZ_{max} [29–30].

The specific resistance (resistivity) of blood ρ_{b} is constant during ejection [29]

The transthoracic specific resistance, ρ_{T}, is constant [31].

_{EPT}

The assumptions of the Sramek-Bernstein equation are virtually identical to those of the Kubicek method, except for the physical definition and magnitude of the V_{C} [32]. In Sramek’s interpretation of Kubicek’s V_{C}, ρ and a Z_{0} variable are eliminated by mathematical substitution. This simplification assumes, as per the work of Quail et al. [29] that, the ρ_{T} equivalent, Z_{0}A/L can replace ρ_{T}. By this substitution, V_{C} is rendered a personal constant for each individual. Sramek [32] named the modification of Kubicek’s V_{C} as the volume of electrically participating thoracic tissue, V_{EPT}. Conceptually, as opposed to Kubicek’s cylindrical V_{C}, Sramek’s V_{EPT} is geometrically a frustum, or truncated cone. Full mathematical derivation and justification for Sramek’s model is provided elsewhere [32]. The original Sramek equation is given as,

where L is the measured distance (cm) between the voltage sensing electrodes and 4.25 an experimentally-derived constant. The Bernstein modification of V_{EPT} assumes that SV is not only a function of thoracic length, but also body weight and blood volume. Correcting for the magnitude of the blood resistance (i.e. intrathoracic blood volume, ITBV, mL) and using the work of Feldschuh and Enson [33], a factor δ (delta) was appended to the Sramek equation. The Sramek-Bernstein equation is given as [32],

where L is a thoracic length equivalent, equal to 17% of body height (i.e. 0.17⋅H) (cm), and delta δ is the weight correction for blood volume. δ is a dimensionless parameter, which corrects for deviation from ideal body weight (kg) at any given height and further modified for the indexed blood volume (mL·kg^{–1}) at that weight deviation [32,33].

In order for equations 35, 37 and 38 to yield equivalent results to equation 17, it must be demonstrated that dZ/dt_{max} is explicitly dependent on aortic systolic pressure P(t), the rate of change of aortic pressure dP/dt, mean arterial pressure P_{mean}, systemic vascular resistance R_{s} and aortic compliance C. In an early, well designed experimental study, Yamakoshi et al. [1] demonstrated that, in an _{max} and blood pressure levels, actually demonstrated that the correlations are weak, at best. Clinically, Brown et al. [35] demonstrated that, despite advancing age (20–80 years) and progressive stiffening of the aorta, ICG CO is nonetheless highly correlated with thermodilution CO (TDCO). It is also interesting to note that, as an alleged analog of dV/dt_{max}, dZ/dt_{max} supposedly provides an ohmic mean velocity analog necessary for SV calculation (_{EPT}) of the Kubicek equation (equation 36), it seems quite improbable that the transthoracic specific resistance, ρ_{T}, remains constant in the face of changing values of Z_{0}. By inspection of equation 23, L and V remaining constant, it is apparent that this is impossibility, because, ρ_{T} must vary directionally with its dependent variable Z_{0}. Other unanswered questions concerning the validity of the plethysmographic hypotheses include the following:

The proper value, physiologic definition, and theoretical basis for ρ: i.e. does ρ vary with hematocrit (ρ_{b}), or is it a constant based on thoracic resistivity (ρ_{T})?

The physiological relevance of measured L or thoracic length equivalents in the Kubicek and Sramek-Bernstein equations regarding SV.

The validity of the outflow extrapolation procedure [36].

A coherent physiologic basis and correlate for the empirically-derived volume conductors and their relevance to SV.

Lack of consideration for the changing transthoracic base impedance (ΔZ_{0}) in critical illness, typified by increased thoracic liquids (pulmonary edema) causing aberrant electrical conduction [37].

Lack of regard for the effect of the blood velocity-induced change in the transthoracic specific resistance Δρ_{b}(t) [37].

Origin of dZ/dt_{max} in the hemodynamic time domain [37].

Despite the objections to the plethysmographic hypothesis, correlation with reference standards is considered good (r^{2} = 0.67, range = 0.52–0.81 and, r = 0.82, range = 0.70–0.90 [7,12]. However, high correlation and reasonable agreement with invasive reference standards does not verify the plethysmographic hypothesis as correct. It may simply mean that the product of dZ/dt_{max}, T_{LVE}, and a best-fit volume conductor V_{C}, yield results that mimic reference method SV [14,15].

_{max}

If equation 23 is differentiated by parts with respect to time, the following results:

Furthermore, if all variables in equation 23 are continuously differentiable functions of time and are expressed within the respective dZ/dt derivatives of equation 39 [37],

Since dL and dL/dt are of trivial magnitude, dZ/dt comprises derivatives 2 and 3, the units of which are Ω·s^{–2} and Ω·s^{–1}, respectively.

Using older plethysmographic techniques, which are conceptually analogous to the simple two-element _{b}(t). The peak slope is conceptually believed to be the peak slope of ΔZ_{vol}(t), and, according to theoretical assumptions, represents the ohmic equivalent of the peak rate of change of aortic volume (peak flow, dV/dt_{max}, mL·s^{–1}), at and before which little outflow is thought to occur [36]. When extrapolated over the ejection interval, this convention is believed to be proportional to total SV [29,30]. To better define the peak slope, ΔZ_{b}(t) is electronically differentiated to dZ/dt (i.e. d[ΔZ_{b}(t)]/dt), the peak magnitude of which is dZ/dt_{max} and conceptually analogous to the peak value of derivative 3 of equation 40 (i.e. d[ΔZ_{vol}(t)]/dt_{max}). That is,

A recently introduced ICG method conceptualizes the peak value of derivative 2 of equation 40 to represent dZ/dt_{max}, which is the peak rate of change of the red cell velocity-induced blood resistivity variation dρ_{b}(t)/dt_{max} (Ω⋅cm·s^{–2}) (i.e d[ΔZ_{vel}(t)]/dt_{max}, Ω·s^{–2}) [37]. Thus,

While differential time domain analysis provides the two possible origins of dZ/dt_{max}, the definitive origin requires temporal correspondence of dZ/dt_{max} with either peak ascending aortic blood flow (Q(t)_{max}, dV/dt_{max}) or peak aortic blood acceleration (d_{max}, dQ/dt_{max}).

Comparative time domain analysis confirms that peak flow velocity occurs at ^{∼} 100±20ms after opening of the aortic valve [38], whereas dZ/dt, (as extrapolated from the data of Matsuda et al. [39] and Lozano et al. [40]) and acceleration d_{max}, data provided by Matsuda et al. [39] show that, when the time interval from the Q wave of the ECG to LVdP/dt_{max} is subtracted from the time interval between the Q wave and dZ/dt_{max}, TTP of dZ/dt_{max} = 60±15ms (i.e. [Q→dZ/dt_{max}(ms)] – [Q→LVdP/dt_{max} (ms)] = 60±15 ms). This analysis is valid, because, in the absence of aortic valve disease, LVdP/dt_{max} almost invariably occurs within 2–5 ms just prior to aortic valve opening. Corroborative evidence from Lozano et al., obtained by subtracting the R(ECG)→B(dZ/dt) interval (mean 68 ms, range 47–84 ms) from the R→dZ/dt_{max} interval (mean 120ms, range 80–160ms), shows that the mean rise time (TTP) from point B to dZ/dt_{max} = 52±20 ms (i.e. [R→dZ/dt_{max} ms] – [R→B ms] = 52±20 ms). These calculated rise times are corroborated by actual measurements obtained by Debski et al. [41]. It can also be shown that TTP of both d_{max} and dZ/dt_{max} are identical. Regression equations by Adler et al. [42] and graphic evidence interpolated from Lyssegen et al. [43] show that, at normocardia (70–90 bpm) TTP of LVdP/dt_{max} is 50–60 ms, which is equivalent to the rise time of dZ/dt_{max}. These data suggest that the inotropic forces governing the potential energy generated during isovolumic contraction time (IVCT) are transferred, unabated, converting potential energy to equivalent kinetic energy (as expressed as d_{max} or equivalently dZ/dt_{max}) during the earliest phase of ejection. Welham et al. [44] have definitively shown that, during halothane-induced myocardial depression LVdP/dt_{max}, d_{max}, and dZ/dt_{max} decrease and recover proportionately with the depth of anesthesia. With profound myocardial depression, it is easily demonstrated that dZ/dt_{max} peaks synchronously with d_{max} and appreciably before Q (t)_{max}.

Figure 6 shows that dZ/dt_{max} intersects the flow curve during peak acceleration and appreciably before Q(t)_{max}. Simultaneously obtained waveforms comparing TTP of peak flow velocity and dZ/dt_{max}, showing the latter temporally preceding the former, can be found elsewhere [30,45,46,47].

Figure 7 shows that dZ/dt_{max} peak synchronously with aortic dP/dt_{max}, where, as extrapolated from equation 9 and its following discussion (_{max} corresponds in time with dF/dt_{max} (major derivative) when dA/dt_{max} (minor derivative) is trivial. This indicates that the peak rate of change of force of ventricular ejection is manifest as dZ/dt_{max} in impedance cardiography. As equations 9 and 11 predict and figure 1 indicates, Q(t)_{max} will occur synchronously with dA/dt_{max}, when aortic dP/dt, dF/dt and dQ/dt=0.

Figure 6 indicates that, at Q(t)_{max}, dZ/dt = 0. These observations infer, if not prove that dZ/dt_{max} occurs in earliest systole, contemporaneously with aortic dP/dt_{max} and peak blood flow acceleration (d_{max}, dQ/dt_{max}). Further proof that dZ/dt_{max} peaks with and is the electrical analog of d_{max} is verified experimentally by their corresponding relationship with the “I” wave of the ultra-low frequency acceleration ballistocardiogram (aBCG). Winter et al. [48] showed that the “I” wave of the aBCG corresponds precisely in time with peak aortic blood acceleration d_{max}, while both Kubicek [49] and Mohapatra and Hill [50] demonstrated that the “I” wave corresponds precisely in time with dZ/dt_{max} (fig 8). Seitz and McIlroy [51] demonstrated that peak blood acceleration and the “I” wave of the HJ interval occur 40–50 ms after opening of the aortic valve, which is consistent with the TTP of dZ/dt_{max}. Consistent with figure 7, showing that dZ/dt_{max} peaks with aortic dP/dt_{max}, Reeves et al. [52] demonstrated that the “I” wave of the aBCG peaks with the second derivative of the carotid contour volume displacement curve (i.e. carotid dP/dt_{max}, pressure acceleration).

Gaw et al. [23] showed that the relative blood conductivity change (Δσ/σ (%), i.e. [Δρ_{b}(t)/ρ_{b}]^{–1}(%)) parallels the peak acceleration of the red cell reduced average velocity (mean spatial velocity, <^{–1})) with r = 0.99 (fig 4). This implies that the normalized peak rate of change of the blood resistivity variation dρ_{b}(t)/dt_{max}/ρ_{b(stat)} is hemodynamically equivalent to the peak rate of change of the red cell reduced average blood acceleration, which is concordant with equation 42. Thus, dZ/dt_{max} appears during the inertial phase of earliest ventricular ejection when little volume change occurs and represents the electrical analog of the most explosive phase of the initial ventricular impulse. Thus, dZ/dt_{max} is bioelectrically equivalent to equation 42 and possesses the units of Ω⋅s^{–2}.

As discussed by Gaw [23] and Visser [19], the relative change of blood resistivity (or equivalently blood conductivity) due to aortic blood flow is related, hemodynamically, to an exponential power (n) of the reduced average blood velocity (i.e. mean spatial velocity) [(<^{n}, (s^{–1})^{n}]. It thus follows that dZ/dt_{max}/Z_{0} (s^{–2}) is the ohmic analog of peak aortic reduced average blood acceleration (PARABA), this being the mean acceleration divided by the aortic valve radius.

Through the relationship for mean flow velocity given by Visser [19], PARABA can be reduced to mean blood flow velocity by square root transformation,

Preliminary evidence suggests that the exponent m is in the range of 1.15–1.25.

To obtain ohmic mean velocity (s^{–1}) from the mean acceleration analog, dZ/dt_{max}/Z_{0} must also undergo square root transformation.

It naturally follows that,

Inasmuch as dZ/dt_{max} is the electrodynamic equivalent of mean aortic blood acceleration, and SV is obtained from a mean velocity calculation, it is suggested that equations 35, 37, and 38 produce a mean acceleration surrogate of SV, which is impossible according to equation 6. Thus,

Based on the above discussion, this leads to the inevitable conclusion that,

In view of this discussion, the outflow correction factor, dZ/dt_{max} x T_{LVE}, is rendered moot, because dZ/dt_{max} represents axial blood acceleration and not radially-oriented rate of change of volume according to _{max}), dZ/dt_{max} × T_{LVE} provides acceleration facsimiles of ohmic mean velocity.

As per equation 46, the relationship between ohmic mean velocity and ohmic mean acceleration is parabolic (ohmic mean acceleration = (ohmic mean velocity)^{2} or y = x^{2}. This relationship predicts that, over wide ranges of dZ/dt_{max}, ohmic mean velocity will be overestimated and underestimated at its upper and lower extremes, respectively.

Indeed, Yamakoshi et al. [1] showed that ICG-derived SV overestimates reference method SV in healthy canines by as much as 70%, and correspondingly underestimates reference SV by as much as 25% when myocardial failure is induced. Similarly, Ehlert and Schmidt [56] could not find a linear relationship between ICG and EMF-derived SV over a wide range of hemodynamic perturbations.

_{max}and ICG-Derived SV and CO

As inferred, the static specific resistance of blood, ρ_{b}, nor the static specific resistance of the thorax, ρ_{T}, are included in the ohmic acceleration-based SV method. Despite Visser’s [19] and Hoetink’s [22] _{b}(t) (i.e ΔZ_{b}(t)_{max}) is hemocrit-dependent, its maximum upslope, dZ/dt_{max}, is not (figure 3 in Visser et al.). In a canine model, Quail et al. [31] showed that normovolemic hemodilution over a hematocrit range of 66%–26%, produced an appropriate increase in SV, which was not different from electromagnetic flow-meter-derived (EMF) SV. Corroboratively, an

_{EPT}, V

_{C})

In contrast to the earlier methods, where the swept volume of the thorax (or portion thereof) is considered the appropriate V_{C}, the V_{EPT} corresponding to the new method is conceptualized physiologically and by magnitude as the intrathoracic blood volume (ITBV, V_{ITBV}). As opposed to linear-based volume conductors, using thoracic length or a height-based equivalent, ITBV is biophysically assumption-free, inherently unambiguous in physiologic meaning and intuitively understood as the physical embodiment of the blood resistance, R_{b}. Physiologically, ITBV has been shown to be highly correlated with left ventricular preload, as expressed by left ventricular end-diastolic volume (LVEDV), and thus with absolute values for SV and directional changes thereof [58–61]. Computationally, V_{ITBV} is found through linear allometric equivalents of body mass (kg). Supporting this relationship, studies show that body mass correlates linearly, and much more closely with total blood volume (TBV), SV and CO than patient height [62–66]. By magnitude, the ITBV represents approximately 25% of TBV, or about 17.5 mL⋅ kg^{–1}, which results in V_{ITBV} = 17.5 x W_{kg}, or equivalently 16W_{kg}^{1.02} [37,67] By comparison, existing volume conductors associated with the plethysmographic hypothesis are modeled as simple geometric abstractions, which are firmly rooted in basic electrical theory. By virtue of their “best-fit” mathematical construction, they bare little relevance to, and have virtually no biophysical relationship with other commonly accepted physiologic, anatomic, or hemodynamic parameters.

_{max}

Evidence suggesting that dZ/dt_{max} varies inversely with aortic valve CSA, or radius r, as a function of PARABA, is inferentially demonstrated through the work of Sageman [68]. He clearly showed that an inversely proportional and highly negatively correlated (r = –0.75) relationship exists between dZ/dt_{max}/Z_{0} and body mass in healthy humans. Aortic valve CSA has been shown to correlate highly with body mass (kg) and body surface area (BSA, m^{2}). But by contrast with dZ/dt_{max}/Z_{0}, Doppler and EMF peak velocities and systolic velocity integrals are totally independent of body mass [69]. Thus, Newtonian-based peak velocities of equal magnitude, measured between age-matched individuals of different body mass and aortic valve CSAs, will produce correspondingly disparate values of dZ/dt_{max}/Z_{0}. It follows that, while there is no direct proportionality between hemodynamically-based and impedance-derived systolic velocity integrals within or between individuals, a linear equivalence exists through their respective mean flow values. Specifically, as a convoluted abstraction of the equation of continuity,

Because of the absolute dependency of dZ/dt_{max}/Z_{0} upon PARABA, this means that, for any given value of mean acceleration, the magnitude of dZ/dt_{max} is related to and explicitly dependent on aortic root CSA by its dependency on R. Since aortic root CSA is a function of body mass, age, and gender, dZ/dt_{max} will vary accordingly. Thus, the magnitude of dZ/dt_{max} is multi-factorial and not wholly dependent on the respective levels of myocardial contractility and Z_{0}. As a first order approximation,

One of the major drawbacks of the impedance technique has been its inability to correctly predict SV in the presence of excess EVLW [70–72]. Critchley et al. [73] tested the hypothesis that the poor agreement between invasive reference standards and ICG is due to excess EVLW. Typified by Sepsis, they were able to show that ICG CO underestimated its thermodilution counterpart and that the degree of underestimation was related to the degree of EVLW excess. They also noted that, in general, excess EVLW was associated with values of Z_{0}<20–22Ω. In a subsequent experimental model in canines, where pulmonary edema was induced by oleic acid infusion, they found a systematic progressive bias between flow probe and ICG CO [74]. As pulmonary edema progressively worsened, ICG CO progressively underestimated the flow probe estimate. In concert with the progressive divergence of the two methods, Z_{0} progressively decreased (fig. 10).

Expanding equation 30, the following is a useful analytical tool for studying the effect of excess EVLW on the various compartments of the transthoracic impedance with the AC field.

In this model it is assumed that AC (I(t)) flows exclusively through the blood resistance Z_{b} (100Ω·cm–180Ω·cm) (element 2), despite the fact that Z_{e} (element 3) is of lower specific resistance (60Ω·cm–70Ω·cm). This assumption is probably not accurate, because the blood resistance R_{b} is considered a cylindrical conductor surrounded by a more highly conductive EVLW impedance Z_{e}. By the reciprocal rule for parallel impedances, the resultant impedance (i.e. Z_{b}||Z_{e}) would be lower than either impedance alone. If Z_{e} is of variable magnitude and Z_{b} is held constant, then Z_{0} should vary with Z_{e}. This is precisely what is observed clinically. When Z_{e} reaches a critical volume (V_{e(CRIT)}), Z_{e} becomes the lowest impedance as relates to the components of Z_{0}. This causes electrical shunting, and in the extreme case (Z_{0} =10Ω–12Ω), a complete short circuiting of current away from the blood resistance. This results in preferential flow through Z_{e} at the expense of Z_{b}, with the result being a decrease in the magnitude of Z_{0} and the amplitude of ΔZ_{b}(t) [1]. At this critical level of Z_{0} (i.e. Z_{0(CRIT)}, Z_{C}), which in humans is 20Ω±3Ω, excess EVLW causes spuriously reduced values of dZ/dt_{max}. Thus, dZ/dt_{max} no longer parallels its hemodynamic equivalent, PARABA, resulting in systematic underestimation of ICG SV/CO. In order to compensate for pathologic conduction through excess EVLW, an index of transthoracic aberrant conduction has been derived (ζ, zeta), which, as its magnitude decreases, creates a larger V_{EPT}.

In patients without excess EVLW, V_{EPT} is equivalent to V_{ITBV}. For all values of Z_{0}<20 Ω, 0<ζ<1 and for Z_{0}≥20 Ω, ζ=1. For all values of Z_{0}<20 Ω, the following pertains:

where ζ is given as [37],

where Z_{C} = critical level of base impedance Z_{0(CRIT)} = 20Ω, corresponding to V_{e(critical)}. Z_{0} = measured transthoracic impedance ≤ Z_{C}, and K = a trivial constant →0.

where V_{ITBV} = 16W_{(kg)}^{1.02}. Equation 55 [37] has been prospectively tested in critical and non-critically-ill patients and been shown to provide SV and CO values comparable to standard reference methods [37,75–78]. As reviewed and compiled by Moshkovitz et al. [7] there are other whole body and transthoracic equations, but they have not been effectively prospectively tested.

From the results of this review, it is strongly suggested that the conceptually-based plethysmographic hypothesis for ICG-derived CO is flawed. The incontrovertible evidence provided herein shows quite conclusively that dZ/dt_{max} is a mean acceleration analog and not that of the peak rate of change of aortic volume. Therefore, to obtain ohmic mean velocity from dZ/dt_{max}/Z_{0}, square root transformation is obligatory. Because dZ/dt_{max} represents axial blood acceleration, the “outflow compensation” for runoff is entirely irrelevant as a theoretical problem to be proved or disproved, as per the work of Faes et al. (36). The fact that dZ/dt_{max} is a mean acceleration analog permits the square wave integration, dZ/dt_{max} x T_{LVE}, as a valid convention similar to equations 6 and 49 for the Doppler/EMF method of SV determination. As for the volume conductors implemented by both the Nyboer/Kubicek and Sramek/Bernstein models, there seems little physiologic justification for either approach. The volume of electrically participating thoracic tissue involved in dynamic conduction is clearly the blood resistance. By magnitude, the blood resistance translates directly into the intrathoracic blood volume. While basing the volume conductor on a fixed volume appears to be concordant with the Sramek/Bernstein model, this can only be true when dynamic conduction through the blood resistance is in parallel with peak aortic reduced average blood acceleration. This can only occur, when, by volume, the blood resistance is the compartment with the lowest impedance. Thus, when excess extravascular lung water exceeds some critical volume, current is diverted to the compartment of lowest impedance, which is now that of the extravascular lung water. With progressive diversion of AC away from the blood resistance, the magnitude of ΔZ_{b}(t) and dZ/dt_{max} diminish and therefore are not reflective of or proportional to the hemodynamic state. Clearly, a more robust mathematical solution for excess extravascular lung water is desirable; that is, if is at all possible.

Algorithms for reconstruction of impedance spectra from non-uniformly sampled step responses Monitoring the skin biophysical parameters among coronavirus patients for three days in a row: a preliminary study A high accuracy voltage approximation model based on object-oriented sensitivity matrix estimation (OO-SME model) in electrical impedance tomography Influence of acute water ingestion and prolonged standing on raw bioimpedance and subsequent body fluid and composition estimates Worth the wait? Time course of supine shifts in body water compartments on variables of bioelectrical impedance analysis Validation of non-empirical fat-free mass estimation model for a wrist-worn device Bioimpedance-based authentication of defrosted versus fresh pork at the end of refrigerated shelf life Bioelectrical impedance changes of the trunk are opposite the limbs following acute hydration change Probing deep lung regions using a new 6-electrode tetrapolar impedance method Intravascular, interstitial and intracellular volume changes during short term deep tissue massage of the calf: A case study The Russian bioimpedance database: an update Development of a small footprint device for measuring electrodermal activity in the palm of the hand Identification of contractions from Electrohysterography for prediction of prolonged labor Cardiovascular diseases diagnosis by impedance cardiography Opinion: The future of electrical impedance tomography Characteristics of Rheoencephalography and some associated factors on menopausal women Comprehensive biotechnical system for screening risk-based diagnosis of COVID-19 and post-COVID syndrome Diabetic foot assessment using skin impedance in a custom made sensor-sock Body composition changes in the immediate peri-operative period following total joint arthroplasty Correlation of the DNA concentration of human samples to electrical bioimpedance measurements: a pilot study A portable band-shaped bioimpedance system to monitor the body fat and fasting glucose level