# A Lumped-Parameter Subject-Specific Model of Blood Volume Response to Fluid Infusion

^{1}Department of Mechanical Engineering, University of Maryland, College Park, MD, USA^{2}Department of Emergency Medicine, Massachusetts General Hospital, Boston, MA, USA

This paper presents a lumped-parameter model that can reproduce blood volume response to fluid infusion. The model represents the fluid shift between the intravascular and interstitial compartments as the output of a hypothetical feedback controller that regulates the ratio between the volume changes in the intravascular and interstitial fluid at a target value (called “target volume ratio”). The model is characterized by only three parameters: the target volume ratio, feedback gain (specifying the speed of fluid shift), and initial blood volume. This model can obviate the need to incorporate complex mechanisms involved in the fluid shift in reproducing blood volume response to fluid infusion. The ability of the model to reproduce real-world blood volume response to fluid infusion was evaluated by fitting it to a series of data reported in the literature. The model reproduced the data accurately with average error and root-mean-squared error (RMSE) of 0.6 and 9.5% across crystalloid and colloid fluids when normalized by the underlying responses. Further, the parameters derived for the model showed physiologically plausible behaviors. It was concluded that this simple model may accurately reproduce a variety of blood volume responses to fluid infusion throughout different physiological states by fitting three parameters to a given dataset. This offers a tool that can quantify the fluid shift in a dataset given the measured fractional blood volumes.

## Introduction

The fluid shift between the intravascular and interstitial compartments is an essential mechanism of homeostasis, and a key determinant of the physiological response to circulatory pathology and medical therapy. The net fluid shift is determined by the summary action occurring across the body's massive network of microvasculature, and the determinants of flow for each microscopic segment are complex, including the permeability of the vessels and the local Starling forces (i.e., hydrostatic and oncotic pressure gradients). Each of these determinants can be, in turn, altered by a wide range of other factors, including the vasomotion of upstream and downstream vessels, lymphatic flow, and a myriad of endocrine and exocrine signals that affect the preceding determinants.

Given the undeniable complexity that underlies fluid shift between the intravascular and interstitial compartments, it can be challenging to mathematically model the fluid shift in an individual subject or patient. Existing models of fluid shift can reproduce the volume changes in the blood and interstitial fluid (Cervera and Moss, 1974; Pirkle and Gann, 1975; Hedlund et al., 1988; Mazzoni et al., 1988; Arturson et al., 1989; Carlson et al., 1996; Gyenge et al., 2003; Tatara et al., 2007) and even intracellular fluid (Hedlund et al., 1988; Mazzoni et al., 1988; Arturson et al., 1989; Carlson et al., 1996; Ursino and Innocenti, 1997; Gyenge et al., 2003; Fernandez de Canete and Del Saz Huang, 2010; Siam et al., 2013). However, there are so many disparate factors that, without exhaustive and impractical measurements, it is not possible to comprehensively characterize the fluid shift associated with an individual. A reasonable alternative is to use “typical” values for certain parameters (i.e., values that represent average values in a population) (Champion et al., 1975; Lewis, 1986; Wears and Winton, 1990; Mardel et al., 1995; Simpson et al., 1996; Hirshberg et al., 2006; Saito et al., 2013) but this is no longer an individualized model. Recently reported volume kinetic models (Ståhle et al., 1997; Svensén and Hahn, 1997; Drobin and Hahn, 1999, 2002; Hahn, 2010) may offer value via a tradeoff between simplicity and transparency. However, these models still have weakness in terms of physiological transparency. As an example, volume kinetic models dictate that the ratio between the changes in the intravascular and interstitial volumes depends only on their baseline values (see Appendix for details), which is an oversimplification of the known physiology that the ratio actually varies with the level of intravascular and interstitial volumes as well as physiological conditions (Guyton et al., 1975). To overcome the challenges in modeling the fluid shift in an individual, we sought to develop a lumped-parameter model that could reproduce the fluid shift at the macroscopic level without microscopic details between the intravascular and interstitial compartments. Specifically, we sought a model which is simple enough to be fitted to individual subjects with only a rudimentary set of measurements, and at the same time accurate enough to characterize that subject's response to fluid infusion. The ability of the model to reproduce real-world blood volume response to fluid infusion was evaluated by fitting it to a series of data.

## Materials and Methods

### Development of a Lumped-Parameter Model of Blood Volume Response to Fluid Infusion

The proposed lumped-parameter model involves only three parameters to describe an individual's fluid shift after fluid infusion: target volume ratio (i.e., ratio between the volume changes in the intravascular and interstitial fluids), α; feedback gain (specifying the speed of fluid shift), K; and initial blood volume, V_{B0}.

The first parameter of the proposed model is α, which varies depending on the overall physiological state of the subject. By way of background, the volume of fluid stored in the intravascular, interstitial and intracellular compartments is determined by the vessel permeability and the hydraulic and osmotic pressure gradients at the capillary walls and cell membranes. From microscopic standpoint, the kinetics of a number of ions and proteins as well as the pressure-volume relationships of the fluid compartments determine the hydraulic and osmotic pressure gradients, and therefore, the fluid volume stored in each compartment. However, from macroscopic standpoint, the consequence of the interaction among these complex mechanisms is that the ratio between the volume changes in the intravascular and interstitial fluid is summarized by a constant parameter value, denoted here as α (Guyton et al., 1975). Typically, interstitial fluid volume changes 2–3 times as much as intravascular volume changes (i.e., α = 2–3) when total blood volume is increased, up to a critical blood volume level. Beyond this level, however, majority of the fluid infused to the body is not stored in the intravascular compartment but is shifted to the interstitial compartment (i.e., α≫2–3) (Guyton et al., 1975; Bajwa and Kulshrestha, 2012). In sum, the fluid shift between intravascular and interstitial compartments, which is due to an array of complex physiological processes, may be viewed as the output of a hypothetical feedback controller that regulates the volume changes in the intravascular and interstitial fluid at a target ratio 1:α.

The second parameter of the model is K, which specifies the speed of fluid shift, i.e., the rate of fluid shift between the intravascular and interstitial compartments. K also varies depending on the physiological state of the subject. This parameter is predicated on a feedback control system analogy shown in Figure 1A. The left and right compartments represent the intravascular and interstitial compartments, respectively, while the valve represents a summary of all the physiological processes that produce the fluid shift between the two compartments. In this analogy, it is both fluid infusion (u) and loss (v; e.g., hemorrhage and urine) that act on the intravascular compartment to alter blood volume (V_{B}). The magnitude of the valve opening is a function of the discrepancy between the target versus actual changes in V_{B}: there is an increased rate of fluid shift,q, between the intravascular and interstitial compartments when the discrepancy between target versus actual V_{B} grows larger.

**Figure 1. Lumped-parameter model of blood volume response to fluid infusion**. **(A)** Feedback control system analogy via connected compartments. The left and right compartments represent the intravascular and interstitial compartments, while the valve represents a summary of all the physiologic processes that produce the fluid shift between the two compartments. The fluid infusion (u) and loss (v; e.g., hemorrhage and urine) act on the intravascular compartment to alter the blood volume (V_{B}), which in turn alters the interstitial fluid volume (V_{ISF}). The magnitude of valve opening is a function of the discrepancy between the target versus actual changes in blood volume. **(B)** Mathematical model. The fluid infusion (u) and loss (v; e.g., hemorrhage and urine) are the inputs to the model, while the change in blood volume (ΔV_{B}) is the output. The objective of the feedback controller is to retain the 1/(1+α) fraction of the inputted fluid volume in the intravascular compartment while shifting the remaining α/(1 + α) fraction to the interstitial compartment in the steady state. The fluid shift from the intravascular to interstitial compartment (q) acts as feedback control to steer ΔV_{B} to the target change in blood volume (r_{B}).

The third parameter of the model is V_{B0}. Its role is to normalize V_{B} to yield its fractional change from the initial V_{B} ((V_{B} − V_{B0})/V_{B0}) as the output of the model. Considering that fractional change in V_{B}, rather than its absolute change, is reported in many experimental protocols (Ståhle et al., 1997; Svensén and Hahn, 1997; Drobin and Hahn, 1999, 2002; Hedin and Hahn, 2005; Hahn, 2010; Hahn et al., 2013), the above normalization allows the model to readily analyze experimental data obtained from such protocols.

The representation in Figure 1A can be formalized into the mathematical model (or block diagram) shown in Figure 1B. Both u and v are the inputs to the model, while the change in V_{B} from V_{B0} (ΔV_{B} = V_{B} − V_{B0}) is the output. The objective of the feedback controller is to retain the 1/(1+α) fraction of the inputted fluid volume in the intravascular compartment while shifting the remaining α/(1+α) fraction to the interstitial compartment in the steady state. q acts as feedback control to steer ΔV_{B} to the target change in V_{B} (r_{B}). In this way, the steady-state volume changes associated with the intravascular and interstitial fluid achieve the 1:α ratio.

The governing equation associated with this mathematical model is derived as follows. The rate of change in V_{B} is given by the following ordinary differential equation:

q is formulated to eliminate the discrepancy between the target versus actual changes in V_{B}. As described above, r_{B} in response to u and v is the 1/(1 + α) fraction of the accumulated resultant fluid intake u − v:

To capture the macroscopic behavior of q with a possibly simplest mathematical expression, q is modeled to be proportional to e_{B} = r_{B} − ΔV_{B}:

This expression physically means that (1) q shifts from the intravascular to interstitial compartment if r_{B} < ΔV_{B} while from interstitial to intravascular compartment if r_{B} > ΔV_{B}, and that (2) q increases as the discrepancy between r_{B} and ΔV_{B} increases. Combining (Equations 1–3) yields the following ordinary differential equation as the lumped-parameter model that relates u − v to ΔV_{B}:

where $\frac{\text{d}}{\text{dt}}\left[\text{u}\left(\text{t}\right)-\text{v}\left(\text{t}\right)\right]=\stackrel{\cdot}{\text{u}}\left(\text{t}\right)-\stackrel{\cdot}{\text{v}}\left(\text{t}\right)$ was used. To rewrite (Equation 4) in terms of the fractional V_{B} response, ΔV_{B} must be normalized by V_{B0}. Dividing both sides of (Equation 4) by V_{B0} yields (where ${\stackrel{\u2323}{\text{V}}}_{\text{B}}(\text{t})=\Delta {\text{V}}_{\text{B}}\left(\text{t}\right)/{\text{V}}_{\text{B0}}$):

The parameters in the model (Equation 5) can be uniquely determined if the data are informative (Ljung, 1999). To illustrate, (Equation 5) can be transformed to the following linear regression:

The vector θ in a linear regression is identifiable (Ljung, 1999) (meaning its elements can be uniquely determined provided the data are informative). Hence, all the parameters can be uniquely determined from θ once it is determined via, e.g., the standard least-squares technique (Ljung, 1999) or even a numerical optimization technique (Storn and Price, 1997): K from the first element of θ; V_{B0} from the second element of θ; and α from the third element of θ using K and V_{B0} thus determined.

### Model Evaluation Method

The lumped-parameter model detailed above is intended to summarize V_{B} response to u and v in an individual subject or patient using just three parameters specific to the individual. However, it is possible that the underlying physiology is too complex to be adequately reproduced by such a simple model. Hence, we sought to evaluate the ability of this lumped-parameter model to reproduce real-world blood volume response to fluid infusion by fitting it to a multitude of different experimental datasets, to assess whether or not such a simple model is capable of accurately reproducing the experimental data.

#### Experimental Data

To analyze the model, a series of datasets reported in the literature were used that had ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response (versus time) to u and v. In all the datasets used in this study, ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ was measured using hemoglobin as the tracer substance. Assuming that the mass of the tracer substance remains constant, an increase in V_{B} results in a decrease in the concentration of the tracer substance in the blood. Hence, ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ can be computed from the change in the blood hemoglobin concentration as follows (Drobin and Hahn, 1999, 2002; Hahn, 2010):

where Hct(t) is the hematocrit at time t, and Hgb(t) is hemoglobin concentration at time t. The datasets used in this study included (1) ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response to bolus infusion of fluid under different V_{B} states (Dataset 1 Drobin and Hahn, 1999); (2) ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response to bolus infusion of crystalloid fluids (Dataset 2 Drobin and Hahn, 2002); and (3) ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response to bolus infusion of colloid fluids (Dataset 3 Hedin and Hahn, 2005).

From the aforementioned reports (Drobin and Hahn, 1999, 2002; Hedin and Hahn, 2005), the volume of fluid infused and ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response were extracted at 10 min intervals as the average of the maximum and minimum responses across all the subjects (although the model is intended for use in individual subjects, average responses were considered here because the response associated with each subject could not be extracted from the visually presented data in these reports). The urine excretion rate was estimated by dividing the total urine volume by the study duration (infusion time + post-infusion observation time) based on the simplifying assumption that it remained constant, because only the total urine volumes were available in these reports.

Overall, we studied a total of seven distinct protocols summarized in Table 1: in Dataset 1, there were three protocols (with 0, 450, and 900 ml pre-infusion hemorrhage); in Dataset 2, there were two protocols (with infusions of 0.9% saline and ringer's lactate); and in Dataset 3, there were two protocols (with infusions of 5% albumin and autologous plasma). The measured ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ responses associated with all the datasets are shown in Figure 2.

**Figure 2. Measured fractional blood volume responses to fluid used for evaluation of the lumped-parameter model**. Ringer's AC, ringer's acetate (Dataset 1; Drobin and Hahn, 1999); Saline,0.9% saline (Dataset 2; Drobin and Hahn, 2002); Ringer's LA, ringer's lactate (Dataset 2; Drobin and Hahn, 2002); Albumin, 5% albumin (Dataset 3; Hedin and Hahn, 2005); Autologous Plasma, autologous plasma (Dataset 3; Hedin and Hahn, 2005).

#### Model Fitting

Fitting the model to each dataset meant estimating the values of the three parameters α, K and **V**_{B0} associated with each dataset. First, Equation (5) was discretized into a difference equation using the Euler's method, where T_{S} is sampling interval (Nise, 2010):

Second, an optimization problem was formulated to derive the set of optimal parameters that minimize the error between ${\stackrel{\u2323}{\text{V}}}_{\mathbf{\text{B}}}$ and its model-reproduced counterpart ${\widehat{\stackrel{\u2323}{\text{V}}}}_{\mathbf{\text{B}}}$:

where Θ is a vector of (unknown) parameters, Θ* is the optimal Θ, ${\widehat{\stackrel{\u2323}{\text{V}}}}_{\text{B}}(\text{i})$ is model-reproduced ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$(i) at a discrete time instant i (i = 1⋯N) for a protocol. For a given set of Θ, ${\widehat{\stackrel{\u2323}{\text{V}}}}_{\text{B}}$ was computed by inputting u and v data to (Equation 8). The optimization problem (Equation 9) was solved using MATLAB's Optimization Toolbox (MathWorks, MA, USA).

#### Model Testing

First, the ability of the model to reproduce ${\widehat{\stackrel{\u2323}{\text{V}}}}_{\text{B}}$ data was assessed by analyzing the goodness of fit associated with each protocol, in terms of sample-by-sample error normalized by the average ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ (e) and root-mean-squared error normalized by the average ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ (RMSNE):

where $\overline{{\stackrel{\u2323}{\text{V}}}_{\text{B}}(\text{i})}$ is the average ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response over the study duration, and ${\widehat{\stackrel{\u2323}{\text{V}}}}_{\text{B}}(\text{i},\text{\hspace{0.17em}}{\Theta}^{*})$ is ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ reproduced by inputting u and v data to (Equation 8) characterized by Θ*.

Second, the intravascular and interstitial fluid volume responses to u and v were qualitatively compared across different protocols. $\Delta {\hat{\text{V}}}_{\text{B}}\left(\text{i},{\Theta}^{*}\right)$ associated with each protocol was estimated by $\Delta {\widehat{\text{V}}}_{\text{B}}(\text{i},\text{\hspace{0.17em}}{\Theta}^{*})={\text{V}}_{\text{B}0}^{*}\cdot {\widehat{\stackrel{\u2323}{\text{V}}}}_{\text{B}}(\text{i},\text{\hspace{0.17em}}{\Theta}^{*})$, while the change in interstitial fluid volume $\Delta {\hat{\text{V}}}_{\text{ISF}}\left(\text{i},{\Theta}^{*}\right)$ was estimated as the numerical integration of q(i, Θ*):

Third, the following hypotheses on the parameters were generated to gauge if they were meaningfully derived by the model fitting: (1) α* decreases as pre-infusion hemorrhage increases; (2) ${\text{V}}_{\text{B}0}^{*}$ decreases as pre-infusion hemorrhage increases; and (3) α* is larger in crystalloids than in colloids. These hypotheses were tested by examining the values of α* and ${\text{V}}_{\text{B}0}^{*}$.

### Parametric Sensitivity Analysis

The parametric sensitivity analysis was performed in the frequency domain. The sensitivity functions associated with the model (Equation 5) to V_{B0}, K, and α were derived as the normalized partial derivatives of G(jω), the frequency response function associated with (Equation 5), with respect to V_{B0}, K, and α:

Then, the frequency responses of these sensitivity functions were examined to elucidate how V_{B0}, K, and α are tuned to fit the model to the data on ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response to u and v.

## Results

Table 2 summarizes the results of data analysis, including the estimated parameters for each dataset and the associated error metrics. Overall, the model reproduced the data accurately with average normalized error and RMSNE of 0.6 and 9.5% across crystalloid and colloid fluids.

Figures 3–5 show ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ responses reproduced by the model, together with the measured data, associated with Datasets 1, 2, and 3, respectively. The upper and lower panels present the measured versus model-reproduced ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ responses as well as $\Delta {\hat{\text{V}}}_{\text{B}}$ and $\Delta {\hat{\text{V}}}_{\text{ISF}}$.

**Figure 3. Model evaluation for Dataset 1 ( Drobin and Hahn, 1999)**. Upper panel, Measured fractional blood volume responses to ringer's acetate (“Data”) under different pre-infusion hemorrhage states, versus response reproduced by the lumped-parameter model (“Model”). Lower panel, Model-reproduced changes in blood and interstitial fluid volumes. 1st column, 0 ml pre-infusion hemorrhage; 2nd column, 450 ml pre-infusion hemorrhage; 3rd column, 900 ml pre-infusion hemorrhage.

**Figure 4. Model evaluation for Dataset 2 ( Drobin and Hahn, 2002)**. Upper panel, Measured fractional blood volume responses to crystalloid fluids (“Data”), versus response reproduced by the lumped-parameter model (“Model”). Lower panel, Model-reproduced changes in blood and interstitial fluid volumes. 1st column, 0.9% saline; 2nd column, ringer's lactate.

**Figure 5. Model evaluation for Dataset 3 ( Hedin and Hahn, 2005)**. Upper panel, Measured fractional blood volume responses to colloid fluids (“Data”), versus response reproduced by the lumped-parameter model (“Model”). Lower panel, Model-reproduced changes in blood and interstitial fluid volumes. 1st column, 5% albumin; 2nd column, autologous plasma.

Figure 6 illustrates **(A)** the frequency responses of the model's sensitivity functions *(12)* to V_{B0}, K, and α and **(B)** the mechanism of how V_{B0}, K, and α are tuned to fit the model to the data on ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response to u and v.

**Figure 6. (A)** Frequency responses of the model's sensitivity functions (Equation 12) to V_{B0} and α. **(B)** Mechanism of how V_{B0} and α are tuned in fitting the model to fractional blood volume response data.

## Discussion

This paper presented a lumped-parameter model that can reproduce V_{B} response to u and v. This section considers the strengths and limitations of the model in terms of its parsimony, accuracy and physiological transparency.

### Lumped-Parameter Model

#### Parsimony

The lumped-parameter model presented in this paper is essentially a minimal model (Equation 5) characterized by only three parameters: α, K, and V_{B0}. That this simple model can be fully characterized just by these parameters and accurately reproduce subject-specific V_{B} responses is its key benefit compared with the existing physiology-based first principles models (Pirkle and Gann, 1975; Mazzoni et al., 1988; Arturson et al., 1989; Carlson et al., 1996; Gyenge et al., 2003) whose complexity makes it extremely challenging to reproduce subject-specific V_{B} responses to u and v without exhaustive and impractical measurements.

#### Accuracy

The results of data analysis (Table 2 and Figures 3–5) support the ability of the lumped-parameter model to accurately reproduce V_{B} responses to u and v. The goodness of fit between measured versus model-reproduced ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ responses is encouraging. These results, derived from five different fluids one of which was also associated with three different pre-infusion hemorrhage states (Dataset 1), suggest that the model may be suitable for summarizing V_{B} responses to a wide range of fluids over diverse V_{B} states.

#### Transparency

The proposed lumped-parameter model is a simple and highly structured model, constrained by the physiological principle that the ratio between the volume changes in the intravascular and interstitial fluid is regulated at a target value (Guyton et al., 1975). Thus, it is not obvious if the parameters estimated from the model fitting (Equation 9) would retain the intended physiological meanings. Since the actual values of the model parameters (V_{B0}, K, and α) were not directly measured, the physiological transparency of the estimated parameters could not be established in a strict sense (i.e., by comparing the actual versus estimated parameters). Yet, the physiological “plausibility” of the proposed model is supported by a few observations from the data analysis (Table 2), as detailed below.

First, the values of ${\text{V}}_{\text{B}0}^{*}$ and α* estimated by fitting the lumped-parameter model to measured ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ responses were comparable, at least in terms of the orders of magnitude, to the known nominal plasma volume (~3.2 l for an 70 kg male; Skillman et al., 1967; Gyenge et al., 1999; Hedin and Hahn, 2005) and the ratio between the volume changes in the intravascular and interstitial fluid due to a perturbation in the body fluid volume (~2.3; Guyton et al., 1975; Lewis, 1986; Wears and Winton, 1990; Simpson et al., 1996): the values of α* derived for the models associated with isotonic crystalloids (Datasets 1–2) ranged between 2.4 and 4.5, while the values of ${\text{V}}_{\text{B}0}^{*}$ derived for the models (Datasets 1–3) ranged between 3.1 l and 5.6 l.

Second, α* decreased with an increase in the pre-infusion hemorrhage volume (Dataset 1), which is consistent with the physiological principle that more fluid is retained in the intravascular compartment as its volume decreases (Guyton et al., 1975). In fact, comparing the lower panels in Figure 3, the model estimated that q decreases as pre-infusion hemorrhage increases via a decrease in α*. At the same time, ${\text{V}}_{\text{B}0}^{*}$ also showed physiologically plausible behavior, exhibiting a decreasing trend with an increase in the pre-infusion hemorrhage (Dataset 1), meeting the physiological expectation that pre-infusion hemorrhage would decrease V_{B0}.

Third, α* estimated for colloids (Dataset 3) was smaller than α* estimated for crystalloids (Datasets 1–2), which is consistent with the fact that colloid is better retained in blood than crystalloid due to its relatively large molecule size (Hahn, 2013) (note that this property of colloid should decrease α* since less fluid is shifted to the interstitial compartment). In fact, comparing the lower panels in Figures 3–5, the model estimated that the fluid volume shifted from the intravascular to interstitial compartment relative to the change in V_{B} is largely smaller for colloids than crystalloids (note that the colloid infusion dose is only 40% of the crystalloid infusion dose; Table 1).

In sum, the lumped-parameter model presented in this paper is (1) simple in that it can be characterized by just three parameters; (2) able to reproduce V_{B} responses to u and v; and (3) showed physiologically plausible behaviors in terms of the estimated model parameters. We speculate that the lumped-parameter model may have clinical uses. For example, each time it is fitted to the response of a subject to bolus infusion of fluid, it may (1) characterize V_{B} state of a subject (e.g., by way of ${\text{V}}_{\text{B}0}^{*}$); (2) estimate unmeasured variables of interest, such as the volume of inter-compartmental fluid shift and the change in interstitial fluid volume; and (3) predict future physiological responses of a subject, such as predicting V_{B} response to future u and v. These speculative applications will be evaluated in future work.

### Parametric Sensitivity Analysis: Mechanism Underlying Model Fitting

Parametric sensitivity analysis provided in-depth insight as to how the model parameters are tuned to fit the model to the data. The Bode magnitude plots of the sensitivity functions (Equation 12) to V_{B0}, K, and α suggest that the model's response is primarily affected by V_{B0} in the high frequency regime while comparably affected by V_{B0} and α in the low frequency regime, with the effect of K limited to mid-frequency regime (Figure 6A). Physically, this means that the model's transient response is largely affected by V_{B0} (and, to a modest extent, by K) while its steady-state response is affected comparably by both V_{B0} and α. Thus, in solving (Equation 9), the model tends to fit the transient portion of the response by primarily adjusting V_{B0}, while it tends to fit the steady-state portion of the response by primarily adjusting α once a value of V_{B0} that best fits the transient response is determined (Figure 6B).

The above mechanistic insight on model fitting may elucidate why V_{B0} was often estimated to be larger than what is regarded as nominal as detailed below (see, e.g., V_{B0} estimated for Dataset 2 in Table 2). In response to a fluid bolus, V_{B} first peaks and then decreases as the fluid is distributed to the extravascular compartment. Thus, the real α may first peak and then decrease to a steady-state value as well (Guyton et al., 1975). However, the model assumes that the value of α is fixed at its steady-state (i.e., the lowest) value. The model (Equation 5) dictates that this α, together with actual V_{B0}, yields an increase in the model-reproduced ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response. To compensate for this discrepancy, V_{B0} is tuned to be larger than its actual value to decrease the transient ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ response. This may be (at least) a mechanism underlying the large values of V_{B0} derived for the Datasets 1–2. Compared with crystalloid (Datasets 1–2), V_{B0} derived for colloid (Dataset 3) were physiologically more plausible, because the infusion dose of colloid was small (Table 1), resulting in less likelihood of volume overload, thereby making the assumption of constant α more realistic.

### Study Limitations

Despite promising initial results obtained for the proposed lumped-parameter model, there are a few limitations that must be investigated in more depth in future studies. First, the amount of data examined in this study was limited (seven datasets representing average subjects). Considering that the model is ultimately intended for use in individual subjects, the model must be assessed using datasets associated with individual (rather than average) subjects receiving fluid infusion under diverse physiological states. Noting that we only evaluated the model's ability to reproduce V_{B} response to fluid infusion (i.e., testing the model against the same data used for fitting), particular focus must be given to evaluate the model's ability to predict V_{B} response (i.e., testing the model against fresh data not used for fitting). Second, the physiological relevance of the model parameters was not fully established. Future work is required to compare estimated model parameters with their gold standard counterparts to assess the agreement between them. V_{B0} may be compared with its direct measurement, e.g., based on the spectrophotometric detection of indocyanine green (Henschen et al., 1993). α may be compared with the ratio between ΔV_{B} and ΔV_{ISF} in the steady state, which may be derived from the measurements of V_{B0} and ${\stackrel{\u2323}{\text{V}}}_{\text{B}}$ as well as u and v. It is noted that α varies with volume state and is also altered by many physiological factors (Guyton et al., 1975). Thus, α must be interpreted as a resultant (and perhaps) linearized ratio between ΔV_{B} and ΔV_{ISF} within the physiological regime presented by the data, especially under the presence of diverse physiological changes. Third, blood hemoglobin concentration measurement can be subject to inaccuracy under large volume changes (Dasselaar et al., 2012). Future work must investigate the influence of the measurement inaccuracy on the quality of the estimated model parameters. Finally, the amount and attributes (e.g., frequency contents) of data required to accurately determine subject-specific model are yet to be examined. In this study, we have shown via the identifiability analysis that the model can be uniquely determined if the dataset is informative, and that model parameters estimated from dataset containing both transient and steady state exhibit physiologically plausible behaviors. Yet, for the model to be clinically useful, model parameters must be reliably determined from minimal data. Future analysis on this aspect will thus be necessary.

## Author Contributions

RB, AR, and JH contributed to the development and validation of the lumped-parameter model as well as the manuscript preparation.

## Funding

This material is based on work supported by the Office of Naval Research (ONR) under Grant No. N000141410591 and N000141512018. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the ONR.

## Conflict of Interest Statement

The authors have a pending patent related to the lumped-parameter model developed in this study.

## References

Arturson, G., Groth, T., Hedlund, A., and Zaar, B. (1989). Computer simulation of fluid resuscitation in trauma: first pragmatic validation in thermal injury. *J. Burn Care Res.* 10, 292–299. doi: 10.1097/00004630-198907000-00002

Bajwa, S. S., and Kulshrestha, A. (2012). Diagnosis, prevention and management of postoperative pulmonary edema. *Ann. Med. Health Sci. Res.* 2, 180–185. doi: 10.4103/2141-9248.105668

Carlson, D. E., Kligman, M. D., and Gann, D. S. (1996). Impairment of blood volume restitution after large hemorrhage: a mathematical model. *Am. J. Physiol.* 270, R1163–R1177.

Cervera, A. L., and Moss, G. (1974). Crystalloid distribution following hemorrhage and hemodilution: mathematical model and prediction of optimum volumes for equilibration at normovolemia. *J. Trauma Acute Care Surg.* 14, 506–520. doi: 10.1097/00005373-197406000-00007

Champion, H. R., Sturdivan, L., Nolan, J., Stega, M., Adams Cowley, R., Sacco, W., et al. (1975). A mathematical model for volume replacement in bleeding patients. *J. Surg. Res.* 19, 297–302. doi: 10.1016/0022-4804(75)90053-0

Dasselaar, J. J., van der Sande, F. M., and Franssen, C. F. M. (2012). Critical evaluation of blood volume measurements during hemodialysis. *Blood Purif.* 33, 177–182. doi: 10.1159/000334142

Drobin, D., and Hahn, R. G. (1999). Volume kinetics of Ringer's solution in hypovolemic volunteers. *Anesthesiology* 90, 81–91. doi: 10.1097/00000542-199901000-00013

Drobin, D., and Hahn, R. G. (2002). Kinetics of isotonic and hypertonic plasma volume expanders. *Anesthesiology* 96, 1371–1380. doi: 10.1097/00000542-200206000-00016

Fernandez de Canete, J., and Del Saz Huang, P. (2010). First-principles modeling of fluid and solute exchange in the human during normal and hemodialysis conditions. *Comput. Biol. Med.* 40, 740–750. doi: 10.1016/j.compbiomed.2010.07.001

Guyton, A. C., Taylor, A. E., and Granger, H. J. (1975). *Dynamics and Control of Body Fluid*. Philadelphia, PA: W.B. Saunders Company.

Gyenge, C. C., Bowen, B. D., Reed, R. K., and Bert, J. L. (1999). Transport of fluid and solutes in the body I. Formulation of a mathematical model. *Am. J. Physiol.* 277, H1215–H1227.

Gyenge, C. C., Bowen, B. D., Reed, R. K., and Bert, J. L. (2003). Preliminary model of fluid and solute distribution and transport during hemorrhage. *Ann. Biomed. Eng.* 31, 823–839. doi: 10.1114/1.1581878

Hahn, R. G. (2010). Volume kinetics for infusion fluids. *Anesthesiology* 113, 470–481. doi: 10.1097/ALN.0b013e3181dcd88f

Hahn, R. G. (2013). Why are crystalloid and colloid fluid requirements similar during surgery and intensive care? *Eur. J. Anaesthesiol.* 30, 515–518. doi: 10.1097/EJA.0b013e328362a5a9

Hahn, R. G., Bergek, C., Gebäck, T., and Zdolsek, J. (2013). Interactions between the volume effects of hydroxyethyl starch 130/0.4 and Ringer's acetate. *Crit. Care* 17, R104. doi: 10.1186/cc12749

Hedin, A., and Hahn, R. G. (2005). Volume expansion and plasma protein clearance during intravenous infusion of 5% albumin and autologous plasma. *Clin. Sci. Lond.* 108, 217–224. doi: 10.1042/cs20040303

Hedlund, A., Zaar, B., Groth, T., and Arturson, G. (1988). Computer simulation of fluid resuscitation in trauma. I. Description of an extensive pathophysiological model and its first validation. *Comput. Methods Programs Biomed.* 27, 7–21. doi: 10.1016/0169-2607(88)90099-5

Henschen, S., Busse, M. W., Zisowsky, S., and Panning, B. (1993). Determination of plasma volume and total blood volume using indocyanine green: a short review. *J. Med.* 24, 10–27.

Hirshberg, A., Hoyt, D. B., and Mattox, K. L. (2006). Timing of fluid resuscitation shapes the hemodynamic response to uncontrolled hemorrhage: analysis using dynamic modeling. *J. Trauma* 60, 1221–1227. doi: 10.1097/01.ta.0000220392.36865.fa

Lewis, F. R. (1986). Prehospital intravenous fluid therapy: physiologic computer modelling. *J. Trauma Acute Care Surg.* 26, 804–811. doi: 10.1097/00005373-198609000-00005

Ljung, L. (1999). *System Identification: Theory for the User, 2 Edn*. Upper Saddle River, NJ: Prentice Hall.

Mardel, S. N., Simpson, S. H., Kelly, S., Wytch, R., Beattie, T. F., and Menezes, G. (1995). Validation of a computer model of haemorrhage and transcapillary refill. *Med. Eng. Phys.* 17, 215–218. doi: 10.1016/1350-4533(95)95712-J

Mazzoni, M. C., Borgström, P., Arfors, K. E., and Intaglietta, M. (1988). Dynamic fluid redistribution in hyperosmotic resuscitation of hypovolemic hemorrhage. *Am. J. Physiol.* 255, H629–H637.

Pirkle, J. C., and Gann, D. (1975). Restitution of blood volume after hemorrhage: mathematical description. *Am. J. Physiol.* 228, 821–827.

Saito, F., Shimazu, T., Miyamoto, J., Maemura, T., and Satake, M. (2013). Interstitial fluid shifts to plasma compartment during blood donation. *Transfusion* 53, 2744–2750. doi: 10.1111/trf.12120

Siam, J., Mandel, Y., and Barnea, O. (2013). Optimization of oxygen delivery in fluid resuscitation for hemorrhagic shock: a computer simulation study. *Cardiovasc. Eng. Technol.* 5, 82–95. doi: 10.1007/s13239-013-0169-z

Simpson, S. H., Menezes, G., Mardel, S. N., Kelly, S., White, R., and Beattie, T. (1996). A computer model of major haemorrhage and resuscitation. *Med. Eng. Phys.* 18, 339–343. doi: 10.1016/1350-4533(95)00044-5

Skillman, J. J., Awwad, H. K., and Moore, F. D. (1967). Plasma protein kinetics of the early transcapillary refill after hemorrhage in man. *Surg. Gynecol. Obstet.* 125, 983–996.

Ståhle, L., Nilsson, A., and Hahn, R. G. (1997). Modelling the volume of expandable body fluid spaces during i.v. fluid therapy. *Br. J. Anaesth.* 78, 138–143. doi: 10.1093/bja/78.2.138

Storn, R., and Price, K. (1997). Differential evolution – A simple and efficient heuristic for global optimization over continuous spaces. *J. Glob. Optim.* 11, 341–359. doi: 10.1023/A:1008202821328

Svensén, C., and Hahn, R. G. (1997). Volume kinetics of Ringer solution, dextran 70, and hypertonic saline in male volunteers. *Anesthesiology* 87, 204–212. doi: 10.1097/00000542-199708000-00006

Tatara, T., Tsunetoh, T., and Tashiro, C. (2007). Crystalloid infusion rate during fluid resuscitation from acute haemorrhage. *Br. J. Anaesth.* 99, 212–217. doi: 10.1093/bja/aem165

Ursino, M., and Innocenti, M. (1997). Modeling arterial hypotension during hemodialysis. *Artif. Organs* 21, 873–890. doi: 10.1111/j.1525-1594.1997.tb00247.x

Wears, R. L., and Winton, C. N. (1990). Load and go versus stay and play: analysis of prehospital i.v. fluid therapy by computer simulation. *Ann. Emerg. Med.* 19, 163–168.

## Appendix

### Volume Kinetics Model

For illustration purposes, a two-compartment volume kinetic model is considered. However, the finding applies to general volume kinetic models. The governing equation is given by (Hahn, 2010):

where v_{c} and v_{c0} are the blood volume and its baseline value, v_{t} and v_{t0} the interstitial fluid volume and its baseline value, R_{0} fluid infusion rate, CL_{0}, CL and CL_{d} the clearances associated with the baseline fluid loss, dilution-dependent elimination and inter-compartmental distribution, respectively. Denoting the changes in the intravascular and interstitial fluid volumes as state variables x_{1} and x_{2}, (*Eq. A1*) is reformulated into the following:

where K_{1} = CL/v_{c0}, K_{2} = CL_{d}/v_{c0}, K_{3} = CL_{d}/v_{t0} and u = R_{0} − CL_{0}. So, taking the Laplace transform and solving for x_{1} and x_{2} yields:

This results in the following:

Thus, the steady-state ratio between x_{1} and x_{2} (equivalent to α in our lumped-parameter model) reduces to lim_{s → 0}x_{2}(s)/x_{1}(s) = K_{2}/K_{3} = v_{t0}/v_{c0}, meaning that the ratio between the changes in the intravascular and interstitial fluid volumes is determined by the ratio of their initial values.

Keywords: fluid shifts, blood volume, fluid resuscitation, mathematical model, lumped-parameter model, subject-specific model

Citation: Bighamian R, Reisner AT and Hahn J-O (2016) A Lumped-Parameter Subject-Specific Model of Blood Volume Response to Fluid Infusion. *Front. Physiol*. 7:390. doi: 10.3389/fphys.2016.00390

Received: 16 May 2016; Accepted: 23 August 2016;

Published: 31 August 2016.

Edited by:

Zbigniew R. Struzik, University of Tokyo, JapanReviewed by:

Yossi Chait, University of Massachusetts, USATsuneo Tatara, Hyogo College of Medicine, Japan

Ofer Barnea, Tel Aviv University, Israel

John K-J. Li, Rutgers University, USA

Copyright © 2016 Bighamian, Reisner and Hahn. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jin-Oh Hahn, jhahn12@umd.edu