# Understanding Volume Kinetics: The Role of Pharmacokinetic Modeling and Analysis in Fluid Therapy

^{1}Department of Clinical Studies, Ontario Veterinary College, University of Guelph, Guelph, ON, Canada^{2}Research Unit, Södertälje Hospital, Södertälje, Sweden^{3}Karolinska Institutet, Danderyds Hospital (KIDS), Stockholm, Sweden^{4}College of Veterinary Medicine, Lincoln Memorial University, Harrogate, TN, United States

Fluid therapy is a rapidly evolving yet imprecise clinical practice based upon broad assumptions, species-to-species extrapolations, obsolete experimental evidence, and individual preferences. Although widely recognized as a mainstay therapy in human and veterinary medicine, fluid therapy is not always benign and can cause significant harm through fluid overload, which increases patient morbidity and mortality. As with other pharmaceutical substances, fluids exert physiological effects when introduced into the body and therefore should be considered as “drugs.” In human medicine, an innovative adaptation of pharmacokinetic analysis for intravenous fluids known as volume kinetics using serial hemoglobin dilution and urine output has been developed, refined, and investigated extensively for over two decades. Intravenous fluids can now be studied like pharmaceutical drugs, leading to improved understanding of their distribution, elimination, volume effect, efficacy, and half-life (duration of effect) under various physiologic conditions, making evidence-based approaches to fluid therapy possible. This review article introduces the basic concepts of volume kinetics, its current use in human and animal research, as well as its potential and limitations as a research tool for fluid therapy research in veterinary medicine. With limited evidence to support our current fluid administration practices in veterinary medicine, a greater understanding of volume kinetics and body water physiology in veterinary species would ideally provide some evidence-based support for safer and more effective intravenous fluid prescriptions in veterinary patients.

## Introduction

Fluid therapy paradigms are constantly changing (1–3) due to new discoveries (4, 5) and ongoing debate on the ideal fluid choice, dose, rate, and efficacy in different patient populations (6–8). The limited and poor quality scientific literature on IV fluid therapy in veterinary medicine (9), particularly in the feline species, has given rise to empirical fluid therapy recommendations (10, 11) that are based on broad assumptions, outdated physiological principles, clinician's anecdotal experiences, or extrapolation from human clinical trials and canine experimental models. Extrapolation of data from human and canine studies are not ideal as there may be critical species variation and different fluid kinetics associated with various IV solutions. As a routine treatment of hospitalized patients, fluid therapy has the potential to cause fluid overload which increases morbidity and mortality (12), especially in cats (13), yet we have little evidence to support our current fluid administration practices in this species (9).

The recent discovery by the “Fluid Expansion as Supportive Therapy” (FEAST) study in African children whereby standard fluid boluses significantly increased 48 h mortality (14) took the medical profession by surprise and subsequently raised many questions regarding our understanding of fluid pharmacokinetics. Since then, clinicians are more receptive to the idea that fluids should be considered “drugs” and prescribed according to the four D's (drug, dosing, duration, de-escalation) and four phases (resuscitative, optimization, stabilization, evacuation) of fluid therapy (2, 3, 15). The FEAST study is an important reminder that IV fluid therapy can exert varying physiologic effects dependent upon the context in which they are administered (16) and can be detrimental if administered inappropriately (6, 7).

For decades, research in fluid therapy was constrained by the lack of effective methods to assess important effects and outcomes. Fluid dynamics and volume expansion effects have been investigated using conventional indicator dilution techniques (17, 18), blood volume monitoring (18–20), hemoglobin dilution (21), and bioelectrical impedance analysis (22–24). In 1997, several authors (25–28) pioneered a novel and innovative pharmacokinetic (PK) model adapted for body fluid spaces that could provide descriptive data on the physiologic behavior of infused IV fluids in terms comparable to those employed in conventional PK practice (29).

For the past 20 years, Robert G. Hahn and his research group have been studying and refining the concept of volume kinetics (VK), providing a research platform to understand how administration of various IV fluids affect body fluid compartments under different physiologic conditions. Intravenous fluids can now be studied like pharmaceutical drugs, leading to improved understanding of their time-volume effects or elimination half-life on plasma and interstitial fluid compartments (30), making evidence-based approaches to fluid therapy possible. Although this concept originally served as a research model for humans (30–33), a handful of experimental studies have been performed in animals including rabbits (34), pigs (35), and sheep (36–45).

The study of VK in companion animal species is unprecedented, and there is no scientific literature currently available. In response to the call for more scientific research to build an evidence-based foundation for veterinary fluid therapy, the purpose of this review is to introduce the basic concepts of VK, its current use in human and animal research, as well as its potential and limitations as a research tool for fluid therapy research in veterinary medicine.

## Basic Pharmacokinetic Concepts and Models

Pharmacokinetics and pharmacodynamics are complex disciplines that utilize advanced mathematical principles that practitioners sometimes find challenging to comprehend. However, these two branches of pharmacology are pivotal in understanding ways to enhance efficacy and decrease harm for any exogenous substance prescribed (46), thus critical to clinical drug delivery. Pharmacokinetic analysis and simulation is the established approach to guide pharmaceutical dosing and frequency since it outlines a substance's concentration time course within the body (46). Pharmacodynamics, on the other hand, illustrate the dose-related effect of a substance in the body over time (46).

Absorption, distribution, metabolism, and excretion form the classic principles which make up the “ADME” acronym that describes PK (47). Absorption is the process where a substance travels from the site of administration to the site of measurement (46), while bioavailability refers to the degree and rate that an active substance is absorbed and present at the site of action (48). During and after the absorption phase, distribution occurs when the substance travels to and from the site of measurement (46). Volume of distribution (*V*_{d}) is a proportionality constant between measured plasma concentration and the corresponding amount of substance within the body according to the state of substance disposition, i.e., at initial time 0 (*V*_{d(c)}), under steady-state conditions (*V*_{d(ss)}), pseudo-equilibrium conditions (*V*_{d(area)}), or at any given time (*V*_{d(t)}) (49). Volume of distribution is primarily designed to determine the appropriate loading dose for rapid achievement of targeted therapeutic plasma concentration, but often interpreted indirectly to encompass the extent of substance distribution following recourse to physiological models involving plasma protein and tissue drug binding (49). Finally, elimination is the sum of metabolism and excretion (46). The term excretion, defined as the irreversible loss of chemically unchanged compound from the body, is routinely confused and used interchangeably with the term elimination, which is the irreversible loss of substance from the site of measurement (50). Most substances are eliminated by a first-order process and occurs when the amount of substance eliminated at any time is directly proportional to the amount of substance in the body, such that the fraction of elimination over time remains constant (51). Zero-order elimination, on the other hand, occurs when the amount of substance eliminated for each time interval is constant regardless of the amount of substance in the body (51).

The management of a substance by the body is an intricate process as the various components of “ADME” occur simultaneously, thereby constantly altering substance concentration in tissues and fluids (51). Application of mathematical principles simplify and help to describe these complex body processes and facilitate the anticipation of a substance's behavior within the body (51). The most fundamental model utilized in PK is the compartmental model, often designated by the number of compartments needed to depict the substance's behavior in the body, e.g., one-compartment, two-compartment, and multi-compartment models (51). The compartmental model is known as a mammillary or catenary model and is dependent on whether peripheral compartments are linked in parallel or in series to the central compartment, respectively (52). In a mammillary compartmental model, the compartments do not represent a specific organ or fluid space, but rather a group of body tissues or fluid spaces that possess similar substance distribution patterns. Blood (plasma) and highly perfused tissues such as the heart, lungs, liver, and kidneys are often collectively regarded as the central compartment (i.e., vessel rich group), while moderately (i.e., muscle) and poorly perfused tissues (i.e., adipose tissue, cerebrospinal fluid) are considered the peripheral compartment (53).

Compartmental models are said to be deterministic because the model that best describes the PK of a substance is determined by the substance's concentrations in the blood (51). If a substance's plasma concentration-time profile best fits a one-compartment model, this denotes that the substance is dispersed instantaneously and rapidly throughout its volume of distribution (53). The plasma concentration-time profile displays a mono-exponential decline as only a single process, i.e., elimination, contributes to the decrease in plasma concentration, while the distribution phase occurs too rapidly to be charted (53, 54). If a substance's PK pattern is better fitted by a bi-exponential model, this suggests that the substance first distributes in the central compartment and then more slowly to a peripheral compartment (53). The substance distributes back and forth between these two compartments and leaves the body from the central compartment (53). This leads to a plasma-concentration time profile with a bi-exponential decline characterized by two distinct phases: (1) the distribution phase which is the initial steep decline in plasma concentration, and (2) the elimination phase which is the subsequent slower decline in plasma concentration sustained by redistribution of substance from the peripheral compartment (54).

The merit of any model depends on how well it mimics the substance's concentration in the body. Typically, the simplest model that sufficiently predicts changes in substance concentration over time is selected (51). A concentration-time curve can be plotted following serial substance plasma concentration measurement and logarithmically transformed to demonstrate a mono-phasic or bi-phasic time profile (53). Subsequently, any unknown concentrations can be predicted and simulated (51). In PK studies, appropriate doses and dosing intervals are determined using computer simulation after repeated experimental measurements of plasma drug concentration are conducted and compared with theoretical computer-generated concentration-time data using non-linear least squares regression (55).

The conventional way to present PK output is to analyze each subject individually and report the mean parameter values for a group as the “typical parameter values” for that cohort. An alternative approach is to analyze all subjects in one single analysis while keeping track of the between-individual variability. The widespread adoption of such population modeling in human pharmacology, also known as non-linear mixed effects (NLME) modeling, provides a framework to characterize sources of variability in drug disposition and response using statistical models that account for between- and within- individual and experimental variation (56–58). Identification and quantification of primary determinants of between-individual variability allow for improved and tailored use of therapeutic drugs, including individualization of dosage regimes and optimization of dosing schedules (58). Population models are referred to as NLME models due to the inclusion of fixed (i.e., constant within the population) and random (i.e., unexplained and varied within the population) effects parameters that possess a non-linear relationship (58). Between-individual variability can be explained by population characteristics or covariates (e.g., species, breed, age, sex, body weight, disease status, etc.) that are included additively or proportionally to the vector of population parameters or fixed effects (58). Random effect parameters represent variability in PK parameter estimates secondary to between- and within- individual variability (58). NLME modeling can be used for analysis of sparse data collected in scenarios where frequent blood sampling is challenging, dense data collected from a small study population, separate plasma and tissue PK studies, and meta-analysis of published data across disparate study designs (57, 58). These attractive advantages of NLME are recognized to be of value with potential widespread applications in veterinary pharmacology (57, 58).

## Basic Concepts of Volume Kinetics

Volume kinetics (VK) is an adaptation of PK theory allowing insights into the effects of IV fluid on theoretical body fluid compartments (25, 30–33, 55). This approach describes changes in plasma volume (PV) during and after IV fluid infusion over time (27). Volume kinetics can also be used to explore the distribution-elimination pattern of IV fluids under various physiological conditions and simulate any infusion fluid choice or rate to achieve desirable PV expansion (32, 33). Visual inspection of the plasma dilution-time curves provides imperative information on the magnitude and time course of PV expansion following actual IV fluid infusion (26). Meanwhile, computer simulation of plasma dilution-time profiles for any hypothetical fluid infusion using the generated kinetic parameter estimates allows insight into how IV fluid therapy should be planned. Figure 1 illustrates an example of information attainable from visual inspection of plasma dilution-time curves (59).

**Figure 1**. Principle of volume kinetics. Measured plasma dilution-time curves of individual subjects (thin lines) and VK modeled plasma dilution-time curve (thick line) during and after a 30 min IV infusion of 10 mL/kg of 6% HES 130/0.4/9:1 (Voluven^{®}) in 10 healthy male volunteers with mean body weight of 79 kg (59). Serial Hb measurements were obtained throughout the study; hemodilution was corrected to baseline HCT to express plasma dilution. Dividing the total infused volume (790 mL) by the plasma dilution extrapolated to 0.27 at time 0 yielded a *V* of 2.9 L, which approximates the expected size of PV in healthy humans, implying that Voluven^{®} only distributes in the plasma space in health. Plotting the plasma dilution-time curve on a logarithmic paper suggested that Voluven^{®} intravascular half-life is 120 min in healthy humans.

When “ADME” concepts are put into the context of IV fluid administration, the absorption is instantaneous with a 100% bioavailability, metabolism is irrelevant, therefore leaving only distribution and elimination to be considered. The central concepts of both PK and VK analyses are therefore dependent upon the volume of distribution and clearance of an administered substance (31). As in PK, a theoretical mathematical model that corresponds reasonably well to the physiological dynamics of IV fluid is constructed (32, 33). The basic mathematical models developed in VK are known as the volume of fluid space (VOFS) kinetic models. These are designated by the number of fluid spaces necessary that mimic the kinetics of IV fluid within the body akin to traditional PK analysis, e.g., one-volume fluid space (1-VOFS), two-volume fluid space (2-VOFS) kinetic models, etc. (25). Both PK and VK models are admittedly challenging to conceptualize, given that they are mathematical derivatives that do not possess precise tangible anatomical or physiological connections (31).

Unlike pharmaceuticals, water is the main component of both infusion fluids and plasma (30, 31), therefore plasma concentration of IV fluids has to be expressed differently. Blood hemoglobin (Hb) has become an attractive endogenous dilution indicator for blood volume (BV) estimation (21) and correlates inversely with blood water concentration (26). Therefore, serial measurement of acute changes in Hb concentration, which is confined to the intravascular space, can be used as a surrogate marker for IV fluid concentration in VK (32, 55). Serum albumin has been explored as a potential endogenous dilution indicator for VK analysis, however possible loss of albumin into the interstitium negated its use (26). The main assumption of VK modeling and analysis is that Hb remains constant and distributes homogenously with the infused fluid in the circulation, such that serial Hb changes are reflective of IV fluid kinetics.

Volume of distribution (*V*_{d}) in conventional PK analysis is not easily applied to fluid therapy. Unlike conventional PK, where *V*_{d} is assumed constant as the volume of administered drug is often negligible, the volume of distribution of IV fluid (denoted by the symbol *v*) continually changes as a direct consequence of IV fluid administration. Therefore, the baseline output parameter, volume of expandable body fluid space (*V*), representing the volume of distribution of initial fluid infusion is estimated instead in VK (34). Depending on the kinetic model used, *V* can be further categorized as the volume of expandable central body fluid space (*V*_{c}) and the volume of expandable peripheral body fluid space (*V*_{p}). *V* can only be estimated from plasma dilution measurement and represents the functional fluid space actually expanded by IV fluid (26). Therefore, *V* may not always equate to the size of body fluid spaces estimated by external indicator-dilution techniques that disperse into the extracellular fluid spaces of minimally expandable structures such as bone, cartilage, dense connective tissue, and organs confined by tight fibrous capsules (34).

Aside from the above interpretation differences, the mathematical solutions and data analysis are analogous to PK compartmental modeling. Volume kinetic models are based on the assumption that plasma dilution follows a mono-exponential (1-VOFS kinetic model) or bi-exponential (2-VOFS kinetic model) curve. The dilution-time profiles, sometimes augmented by urinary excretion (60), are entered into commercial computer software and fitted to the differential equations describing the fluid shifts in these kinetic models (25). Non-linear least-squares regression is used to generate the best kinetic parameter estimates from the differential equations relevant to each kinetic model (25). Subsequently, analytical or matrix solutions to the differential equations are solved to compute *V* as well as the kinetic constants governing fluid distribution (25). Similar to PK modeling, the kinetic model with the simplest solution should be chosen (31) and statistical tests can be applied to select the best fitted model for presentation (61, 62). Agreement between model-predicted renal clearance (*Cl*_{R}) and actual urine output measurement may also assist in final model selection (32).

Following that, computer simulation of dilution-time profiles for hypothetical fluid infusion can be accomplished by inputting the mean parameter estimates obtained from non-linear least-squares regression into the differential equation solution describing the kinetic model (25). Various commercial mathematical programs such as MATLAB^{®} (MathWorks^{®} Inc, Natick, MA, USA) (63, 64) or specialized PK modeling and simulation software such as Monolix (Lixoft, Antony, France), Phoenix WinNonlin^{®} or Phoenix^{®} NLME^{TM} (Certara USA, Inc., Princeton, NJ, USA) (65) can be utilized for this purpose.

## History and Evolution of Volume Kinetic Modeling

The original VK models were first developed in 1997 based on the concept of clearance and plasma dilution (25–27). The most commonly encountered problems with earlier VK analyses included overestimation and poor precision when estimating the size of the peripheral fluid space (*V*_{p}), which became an issue with slow fluid elimination (25, 26, 66, 67). Over the years, efforts to refine and improve the robustness of these models were made, including the incorporation of urine output as an input variable to VK analysis (60).

Parameter estimates from earlier VK studies are challenging to interpret and compare because they are derived using different kinetic methodologies. This eventually led to a major modification to the VK approach in 2006 which allows for a single comprehensive model solution for isotonic crystalloid solutions (67). In this modified “micro-constant” analysis, the fractional plasma dilution is substituted with absolute volume expansion, flow rate (unit: mL/min) parameters are substituted with rate parameters (unit: /min), and the existence of *V*_{p} is acknowledged but its size is not routinely estimated (67). Despite these changes, this modified model still describes the same functional system as the antecedent models (27, 60, 66, 68) and is able to address bi-directional fluid flux, i.e., outflow when plasma dilution is high, inflow when distributed fluid exerts greater force than *V*_{c} (67). Since then, progressively sophisticated extensions have been explored to address specific clinical research questions including three-volume kinetic modeling which accounts for hypertonic osmotic fluid shifts (37, 68), volume turnover kinetics to predict fluid dynamics following hemorrhage and fluid resuscitation (41), and population kinetic modeling that enables the inclusion of various physiological covariates (69–74).

Over the past 10 years, VK has gradually moved away from the clearance model to the micro-constant model where clearance variables (*Cl*, *Cl*_{d}, *Cl*_{0}) are substituted by rate constants (*k*_{10}, *k*_{12}, *k*_{21}, *k*_{0}) and fractional plasma dilution ($\frac{({v}_{c}-{V}_{c})}{{V}_{c}}$) is replaced by absolute volume expansion (*v*_{c} − *V*_{c}) (32, 33) in the kinetic equations. The clearance and micro-constant models are essentially one and the same, except that the former deals with fractional plasma dilution and flow rate (unit: mL/min), while the latter deals with absolute volume expansion and rate (unit: /min). The same Hb-derived plasma dilution is still used to indicate the volume expansion of *V*_{c} resulting from IV infusions (59). The micro-constant model is preferred over the clearance model as fluid volumes within the different body fluid spaces are now directly apparent and unequal fluid distribution across *V*_{c} and *V*_{p} compartments can be further investigated (30, 33).

With the evolution of VK models over time, many parameters have been represented differently, and so it is important to review the literature with an eye to this evolution to not be confused by the changes in symbols and parameters. This is further complicated by the use of proprietary symbols and expressions between PK computer software and PK studies (75). This has created a great degree of confusion and ambiguity, making PK and VK difficult to understand. Consistency in the use of symbols, units, and nomenclature would hopefully clarify the precise meaning of the term or concept as defined, curtail erroneous interpretation, and allow meaningful comparison between studies.

The following symbols and abbreviations are deduced to be equivalent dependent on whether a clearance model or a micro-constant model is presented: *V*_{1} = *V*_{c}, *V*_{2} = *V*_{t} = *V*_{p}, *v*_{1} = *v*_{c}, *v*_{2} = *v*_{t} = *v*_{p}, *k*_{i} = *R*_{0}, *k*_{b} = *Cl*_{0} = *k*_{0}, *k*_{r} = *Cl* ≅ *k*_{10}, and *k*_{t} = *Cl*_{d} ≅ *k*_{12} and *k*_{21} (25, 32, 59, 67). To improve uniformity and clarity, the parameters, symbols, and abbreviations analogous to those used in PK (75, 76) should be used. Comparison of conventional PK and analogous VK parameters are summarized in Table 1.

## Two-Volume Fluid Space Kinetic Model

Most IV fluids can be modeled using a 2-VOFS kinetic model (33, 55), also referred to as the bi-exponential model. The 2-VOFS kinetic model has been reported to best describe the effects of fluid therapy in human patients receiving crystalloid infusions during surgery, general anesthesia, dehydration, and hypovolemia (32), but does not always seem to describe the situation in healthy human volunteers (26, 27, 68) likely due to the rapid elimination of crystalloids in healthy awake states.

The 2-VOFS model, as illustrated by Figure 2, suggests that IV fluid administered at a certain rate (*R*_{0}) would expand the volume of a central body fluid space (*V*_{c}) to a larger volume (*v*_{c}) (32). As transcapillary hydrostatic and colloid osmotic pressure changes, the fractional plasma dilution ($\frac{{v}_{c}-{V}_{c}}{{V}_{c}}$) (26, 32) or absolute volume expansion (*v*_{c} − *V*_{c}) (67, 77, 78) indicate distribution of fluid to a peripheral body fluid space (*V*_{p}) which in turn expands to a greater volume (*v*_{p}). These central and peripheral body fluid spaces attempt to retain their respective initial volumes to maintain homeostasis (26, 32). The net rate of fluid flux between the central and peripheral body fluid spaces is governed by the volume expansion difference between *V*_{c} and *V*_{p} multiplied by the distribution clearance (*Cl*_{d}) in the clearance model (32, 79), or the intercompartmental rate constants (*k*_{12}, *k*_{21}) in the micro-constant model (67, 77, 78). In the micro-constant model, intercompartmental rate constants can be further categorized into central to peripheral intercompartmental rate constant (*k*_{12}) and peripheral to central intercompartmental rate constant (*k*_{21}) akin to PK analysis (67, 77).

**Figure 2**. Two-volume fluid space (2-VOFS) kinetic model. Adapted from previous VK work (25, 26, 32, 67, 77–79). The symbols in black represent the current micro-constant model (67, 77, 78), whereas the symbols within red parentheses represent the original model (25, 26) and the symbols within blue parentheses represent the clearance model (32, 79).

Fluid elimination occurs from the central body fluid space via two mechanisms: (a) baseline water loss which is accounted for by zero-order elimination (*Cl*_{0} or *k*_{b}) often pre-set between 0.3 and 0.5 mL/min in humans (66, 68, 80), and (b) dilution-dependent, primarily renal excretion governed by first-order elimination (*Cl* or *k*_{10}) (32). The rate of elimination is given by the product of fractional plasma dilution ($\frac{{v}_{c}-{V}_{c}}{{V}_{c}})$ (26, 32) or absolute volume expansion (*v*_{c} − *V*_{c}) (67, 77, 78), and the first-order elimination clearance (*Cl*) or elimination rate constant (*k*_{10}) (32).

The unknown parameters in the 2-VOFS kinetic model that are estimated by non-linear least-squares regression are *V*_{c}, *Cl*_{d} (or *k*_{12} and *k*_{21}), and *Cl* (or *k*_{10}) (31). If urine output is measured, *k*_{0} can even be estimated and includes all other body fluid not allocated by the VK model along with the basal fluid loss (32, 60). The first-order elimination (*Cl* or *k*_{10}) is then set to a fixed value determined by the total urinary excretion obtained. In this scenario, *Cl* is analogous to the renal clearance which is calculated by dividing the urine output by the integral of the plasma dilution-time curve, also known as the area under concentration-time curve (*AUC*), following the assumption that half of the *Cl*_{0} or *k*_{0} due to insensible fluid losses appears as urine (60). The following differential equations describe the clearance model for the 2-VOFS kinetic model (30):

The differential equations describing the micro-constant model for the 2-VOFS kinetic model are as follows (30):

By comparing Equation (1) for the clearance model with Equation (3) for the micro-constant model, net *Cl*_{d} is obtained as the product of *V*_{c} × *k*_{12} and *V*_{p} × *k*_{21} (30). Similarly, *Cl* is the product of *V*_{c} × *k*_{10} (30).

## One-Volume Fluid Space Kinetic Model

Intravenous fluid plasma dilution-time profiles may not always display the bi-exponential form of a 2-VOFS kinetic model and sometimes statistically fit the 1-VOFS kinetic model better (32). A common scenario where this occurs is when crystalloids are eliminated rapidly from the body (27, 68), leading to an increased *Cl* : *Cl*_{d} (or *k*_{10} : *k*_{12} and *k*_{21}) ratio. This reduces the time for full *V*_{p} expansion prior to elimination, therefore *V*_{c} and the partially expanded *V*_{p} fuse into a moderate size single body fluid space, allowing 1-VOFS kinetic model to sufficiently describe these dilution-time profiles (32). In addition, colloidal solutions with high molecular weight polymers that do not distribute extensively into the extravascular tissues should theoretically be described by the 1-VOFS kinetic model (31).

The 1-VOFS kinetic model, also known as the mono-exponential model, is illustrated by Figure 3. As the name suggests, IV fluid administered at a given rate (*R*_{0}) is thought to expand the volume of a single body fluid space (*V*) to a larger volume (*v*) (32). Elimination of fluid still occurs via basal fluid losses (*Cl*_{0} or *k*_{b}) and a dilution-dependent, primary renal excretion governed by a first-order elimination (*Cl* or *k*_{10}) (32). The unknown parameters in the 1-VOFS kinetic model that are estimated by non-linear least-squares regression include *V* and *Cl* (or *k*_{10}) (31). The fractional plasma dilution ($\frac{(v\text{}-\text{}V)}{V}$) in the clearance model or absolute volume expansion (*v* − *V*) in the micro-constant model, is obtained from plasma dilution computed from serial Hb concentration. Similar to the 2-VOFS kinetic model, *Cl*_{0} and *Cl* can be estimated through VK modeling or calculated if urine output is available. The following differential equation describes the clearance model for the 1-VOFS kinetic model (32, 79):

In comparison, the differential equation that describes the micro-constant model for the 1-VOFS kinetic model is as below (67, 77, 78):

**Figure 3**. One-volume fluid space (1-VOFS) kinetic model. Adapted from previous VK work (25, 26, 32, 67, 77–79). The symbols in black represent the current micro-constant model (67, 77, 78), whereas the symbols within red parentheses represent the original model (25, 26) and the symbols within blue parentheses represent the clearance model (32, 79).

## Mass Balance Principles

The volumes of individual body fluid compartments are conventionally estimated using external indicator-dilution techniques. Some examples include but are not limited to, deuterium oxide (81) or radioactive tritium oxide (82) for total body water determination; bromide (83, 84) for extracellular fluid volume determination; and radio-iodinated serum albumin (18), Evans Blue dye (18, 36), or indocyanine green dye (17, 18, 38) for plasma volume determination. Based on the law of mass conservation, indicator-dilution techniques fulfill the mass balance principle at steady-state, whereby the total indicator mass following distribution in the fluid compartment will remain the same as the total mass injected into the compartment (80, 85). Hence, the volume of a body fluid compartment can be estimated provided that the indicator's concentration following homogenous mixing can be accurately determined (84).

Although exogenous indicator-dilution techniques are considered the conventional “gold standard” method for volume determination (17, 86), they can result in hypersensitivity reactions and potential errors as a result of indicator extravasation, plasma turbidity, Hb contamination, or controversial back-extrapolation methods (18, 87–90). Radioactive isotopes are cumbersome to use and there is a significant time delay for exogenous indicators to equilibrate to a steady-state distribution (80, 84). This precludes the ability to distinguish between the distribution and elimination phases (34, 67) and does not permit analysis of acute dynamic fluid shifts (36).

In 1987, an endogenous indicator-dilution technique using blood Hb concentration was introduced as an accurate method for estimation of BV changes following validation with the radio-iodinated serum albumin technique (21). Endogenous Hb is confined solely within the plasma space, not associated with any time-delay for steady-state equilibrium, much easier to use, and widely measurable in routine clinical laboratories (21, 67). Therefore, serial measurement of Hb concentrations following IV fluid administration can be used to quantify plasma dilution in the blood and represent acute dynamic fluid shifts within the body (21). The mass balance technique using serial Hb dilution requires the assumption that total body Hb remains constant and homogenously distributed within the circulation at all times, except during sampling losses or hemorrhage to which correction for Hb losses can be made (21, 38).

## Hemoglobin-Derived Plasma Dilution

Plasma and Hb are the main components of the central body fluid space, and the change in *V* or *V*_{c} is directly proportional to plasma dilution (26). Therefore, serial measurement of Hb concentrations following IV fluid administration can be used to quantify plasma dilution in the blood. Hemoglobin dilution must be converted to the corresponding plasma dilution considering that extracellular fluid, and not red blood cells (RBC), are expanded by the infused fluid. Moreover, PV is the actual entity that distributes across body fluid spaces and from which water is ultimately eliminated from the body. Plasma dilution, which equals the dilution of *V* or *V*_{c}, is used to quantitate fluid volume load. As such, the concentration of IV fluid, or volume expansion, is expressed in terms of Hb-derived plasma dilution as a function of time (32), given by the reference Equation (7). This reference equation can be applied directly into the non-linear least square regression curve-fitting procedure if blood sampling and hemorrhage are negligible.

Simultaneously measured RBCs are diluted similarly as Hb during IV fluid loading but quantified using a different technique by the automated hematology analyzer. Therefore, if available, RBC count and RBC mass could also be included in the curve-fitting procedure to improve precision and accuracy (32), although the author's (RGH) experience with this adaptation yielded negligible effect. The same calculations using the reference Equation (7) are performed for RBC count and RBC mass over time, and the mean value for Hb-derived plasma dilution and RBC-derived plasma dilution can be used as the overall plasma dilution (91). Calculated hematocrit value obtained from automated hematology analyzers is preferred over packed cell volume for plasma dilution determination as the latter frequently overestimates true RBC fraction due to incomplete microhematocrit centrifugation compaction (80).

## Corrections for Hemoglobin Loss

Serial blood sampling can result in sufficient loss of Hb creating a “false” dilution that is unrelated to IV fluid therapy. Therefore, blood loss and sampled volume should be corrected prior to the curve-fitting procedure. The calculations to correct for Hb loss and alteration in erythrocyte size consider the PV expansion instead of plasma dilution. To account for sampling loss, change in BV is calculated by estimating the total Hb mass in the circulation from which blood loss is subtracted (21). An assumption of initial BV has to be preset to time 0 (31). The total baseline Hb mass (*MHb*_{(0)}) is first obtained, then the total Hb mass corrected for sampling at any time (*t*), followed by the expanded BV corrected for sampling at any time (*t*) (32):

Changes in BV at any time (*t*) can be provided by Equation (11), and the BV expression can be easily transformed to PV using Equation (12):

Corrected plasma dilution,

The relationship between baseline Hb (*Hb*_{(0)}) and diluted Hb at any time (*t*) is written as $\frac{H{b}_{(0)}}{H{b}_{(t)}}$ in the reference Equation (7), whereas the inverse relationship is used in Equations (13, 14) (32). The lower Hb concentration must be placed as the denominator of the ratio in the reference Equation (7) to arrive at a correct proportion between changes in Hb and water volume.

The degree of error introduced by the wide range of initial BV has been shown to be forgiving (66), however the error related to frequent blood sampling and the use of successive Hb ratio for volume estimation should not be ignored as the volume-estimation errors can increase exponentially (92). Given the importance of an accurate baseline value, averaging duplicate and triplicate Hb measurements on the initial blood sample is strongly recommended, and the percent coefficient of variation for the entire sampling and analysis should be determined (92). Duplicate Hb measurements for subsequent blood samples throughout the studies have been attempted by the author (RGH) with marginal improvement to VK analysis.

## Corrections for Alternation in Mean Corpuscular Volume and Osmotic Fluid Shift

When investigating IV fluids that alter plasma osmolality, for example hypertonic saline (HS) or glucose infusions, a correction for changes in mean corpuscular volume should also be considered. Mean corpuscular volume can be altered in humans following infusion of HS because of osmotically induced decreases in intracellular water, this has also been reported in dogs and cats (93). Changes in mean corpuscular volume can be considered by the addition of a term for the relationship between the baseline (*MCV*_{(0)}) and at any time (*t*) into Equation (13) above:

If HS or glucose solutions are infused, an osmotic fluid shift theoretically occurs across the cell membrane and water moves from the intracellular fluid space (~40% of BW) to the extracellular fluid space (~20% of BW) (80). Using species-specific baseline serum osmolality (*Osm*_{(0)}), which is ~295 mOsm/kg for humans (32), ~317 mOsm/kg for cats (94), and ~302 mOsm/kg for dogs (95), the translocated volume (*f*_{t}) can be obtained from Equation (16) (37, 68):

The osmotic force diminishes progressively with subsequent amount of infused fluid; therefore, *f*_{t} should be entered as a linear function in the analysis process such that *f*_{t} at each point in time is governed by the total amount of infused fluid (32).

## Volume Kinetic Models in Physiological Contexts

Recall that PK compartment models and VK kinetic models are mathematical concepts that describe body spaces which a substance or fluid appears to occupy, but may not correspond to an actual anatomical space or physiological volume (31, 54). In PK analysis, *V*_{d} is theoretical and does not equate to a real volume (53, 54). A substance could have a large *V*_{d} that substantially exceeds the total body volume (54), implying that it is highly distributed in tissues (53). The smallest *V*_{d} is the PV (51) which suggests that the substance is poorly distributed and confined to the plasma (53). Similarly, this concept applies to *V* and *V*_{c} which substitute *V*_{d} in VK analysis.

In healthy human VK studies, the size of *V*_{c} for crystalloids using a 2-VOFS kinetic model has been reported as 3–4 L (4–6% of BW) (26, 31, 60, 68), which is similar to the expected (80) and measured human PV (69). However, when a 1-VOFS kinetic model is used, the size of *V* is usually twice the size of PV (59, 60, 68) as it combines both the plasma volume and approximately half of the interstitial fluid volume. As is expected for colloids, the size of *V* of a 1-VOFS kinetic model was found to be analogous to the PV (59).

The size of *V*_{p} in healthy adult humans was discovered to be 6–8 L (8–11% of BW) (26, 68, 96), which is smaller than the expected size of the interstitial fluid space (11 L or 15% of BW) (80). The *V*_{p} encompasses about 66.6% of the interstitial fluid space since bones, cartilage, dense connective tissues, and organs with a tight fibrous capsule are either not expandable or minimally expandable by fluid (26, 34, 96). Unlike the conventional exogenous indicator-dilution techniques, VK only highlights expandable body fluid spaces (26). Some interstitial tissues have increased compliance for volume expansion while others require very elevated fluid pressure before expansion can occur (97, 98). The *V*_{p} was discovered to be larger when massive fluid infusions beyond what is routinely safe in clinical practice were administered to sheep (39). The precision of *V*_{p} is often lower than *V*_{c} (26, 32, 60) but is less crucial since the size of *V*_{p} does not contribute significantly to interpretable knowledge (67). So far, estimates of *V*_{p} that are significantly smaller than *V*_{c} have not been reported (67).

In healthy humans, there should be no restriction on fluid movement into the peripheral body fluid space or return to the central body fluid space. Distribution clearance (*Cl*_{d}) or the intercompartmental rate constants (*k*_{12}, *k*_{21}) can be appointed the same value for bi-directional flow based on the assumption that fluid is equally prone to flow in either direction (32). Alternatively, uneven distribution can be analyzed using the micro-constant model by splitting the intercompartmental rate constants into *k*_{12} and *k*_{21}, which could serve to quantify the accumulation of peripheral edema (30, 36, 38). Fluid exchange is considered normal with *k*_{12} that is twice as high as *k*_{21} (33, 71, 77).

According to the model, the parameter *Cl* (or *k*_{10}) that is estimated by non-linear least square regression should approximate the renal clearance (*Cl*_{R}), since the ultimate destiny of IV fluid is elimination via the kidneys (31). It is, however, important to remember that parameter estimates generated from VK modeling are not absolute measurements of physiologic body fluid variables, but rather represent estimates of abstract concepts which provide insights into the body's handling of IV fluid (31). The first-order elimination rate constant has been found to correlate well with urine output in human volunteers receiving isotonic crystalloids (26, 27, 60, 66) and normovolemic conscious sheep (36). However, a discrepancy in measured and estimated values resulting in underestimation of *V*_{p} was reported in sheep during general anesthesia receiving positive pressure ventilation (36). The most common way to calculate the elimination half-life (*t*_{1/2}) of infused fluid is by utilizing *Cl*_{R} which is obtained by dividing urine output with the AUC for the fractional plasma dilution or absolute volume expansion in the clearance or micro-constant model, respectively (30).

## Clinical and Research Applications

The use of VK overcomes the limitations of external indicator-dilution techniques for physiologic spaces estimation or measurements of hemodynamic end-points (37). An external indicator distributes within conventional physiological body fluid spaces but may not reflect the effects of IV fluid under dynamic, non-steady state circumstances. Hemodynamic end-points undeniably provide a different set of useful information, but fundamental information about fluid shifts, functional body fluid space volumes, or mechanisms behind differences in fluid dynamics are not assessable. Although VK fluid space models and parameter estimates are less intuitive and describe fluid handling without precisely measuring physiologic spaces, interpreting them in light of mass balance analysis may improve understanding and describe changes in traditional physiologic compartments that are more familiar to both scientist and clinicians (38). The clinical implications of VK analysis in healthy and sick animals are summarized in Table 2.

### Published Work in Human Literature

There are now more than 60 research publications in human medicine describing the VK of isotonic crystalloids, balanced crystalloids, hypertonic solutions, and colloids under diverse surgical and physiological conditions (30–33). Although first introduced two decades ago, VK studies continue to appear in the scientific literature (70–74, 99–102) and the research method is starting to gain recognition (15, 16), especially in light of the increasing debate on the safety and changing practices of fluid therapy.

The first *in vivo* VK study described the effects of commonly used equipotent IV fluid boluses on the expandable fluid spaces in healthy male human volunteers (mean 80 kg) (26). Balanced crystalloid (25 mL/kg acetated Ringer's solution over 30 min) was found to generate the largest plasma dilution, while the dilutions corresponding to colloidal solution (5 mL/kg 6% Dextran 70 in 0.9% NaCl over 30 min) and HS (3 mL/kg 7.5% NaCl over 30 min) were comparable but with a longer duration of effect (26). Urinary excretion was not measured, and *k*_{10} (reported as *k*_{r} in the original model) was computed by VK modeling (26). This work was further extended to investigate the effects of different balanced crystalloid infusion rates (25 mL/kg acetated Ringer's solution over 15, 30, 45, and 80 min) and volumes (12.5 mL/kg acetated Ringer's solution over 30 min) in healthy female human volunteers (mean 60 kg) (27). Acetated Ringer's solution was discovered to have a more effective and extended plasma expansion effect when infused over a longer time period, and reached a maximum effect of 36% volume expansion (550 mL in an adult human) regardless of infusion rates (27).

Over the years, VK analysis has provided some mechanistic explanations as to why and when patients are sensitive to large volumes of fluid (103). Slow fluid distribution to the peripheral compartment has resulted in a 50–75% larger plasma dilution during a crystalloid infusion than would be expected if the distribution had been immediate (32). In healthy conscious volunteers, isotonic balanced crystalloids are reported to take up to 25–30 min to distribute following rapid IV infusion (25 mL/kg over 30 min) (32). Due to this lag time for crystalloids to equilibrate between plasma and interstitium (66), PV expansion during the actual infusion is much larger than the commonly suggested expansion of 20–25% of the infused volume (32). Several studies have shown that the immediate PV expansion can be substantial. In a group of normovolemic human volunteers, 50% of a 2 L acetated Ringer's solution given over 20 min was retained within the intravascular compartment at the end of the infusion (66). In a study where acetated Ringer's solution was infused continuously throughout transurethral resection of the prostate in male adults under general anesthesia (104), fluid retention averaged 60% of the infused volume. The fraction of infused fluid that persists in the intravascular compartment is higher for slower infusions (27). Therefore, crystalloids may be better PV expanders than currently acknowledged, provided that the infusion is continuous and not administered as a bolus.

Volume kinetic analysis has also been used to describe the distribution and elimination of isotonic and hypertonic fluids within the body (68). When IV administration of 25 mL/kg of 0.9% NaCl was used as a reference fluid and compared to an equal volume of lactated Ringer's solution, acetated Ringer's solution, 5 mL/kg of 7.5% NaCl, and 3 mL/kg of 7.5% NaCl in 6% dextran solution (HSD), all infused over 30 min, plasma dilution efficiency according to the AUC of dilution-time profiles was found to be 0.88, 0.91, 3.97, and 7.22, respectively in 10 healthy male volunteers (mean 81 kg) (68). Based on VK analysis and simulation, the strength of these respective fluids to dilute the plasma by 20% within 30 min was 0.94, 0.97, 4.44, and 6.15 times greater than of 0.9% NaCl, respectively (68). Urinary excretion was found to be 1.8 times and 2.7 times larger than the infused volume of HS and HSD, which the authors attributed to natriuresis induced by sodium load (68). In this study, comparison between fluids was complicated by the need for several models, including a three-volume fluid space kinetic model for hypertonic fluid, to account for possible osmotic fluid shift (*f*_{t}).

In a more recent study (59), the volume effects of 10 mL/kg of 6% HES 130/0.4/9:1 (Voluven^{®}), 20 mL/kg of acetated Ringer's solution, and a combination of HES and acetated Ringer's solution administered 75 min apart was explored using VK in 10 healthy male volunteers (mean 79 kg). The kinetic models were successfully fitted to all experiments using the 1-VOFS kinetic model for HES and the 2-VOFS kinetic model for acetated Ringer's solution. The *V*_{c} for HES in both series of experiments was 3.14 L (~5% of BW), close to the PV estimated by anthropometry, with a *t*_{1/2} of 2 h (59) similar to the *t*_{1/2} reported in the product monograph (105). Hydroxyethyl starch was found to induce diuresis with 85% of the infused volume excreted as urine when administered alone (59). The *V*_{c} for acetated Ringer's solution averaged 4.88 L (~6.2% of BW) with a *t*_{1/2} of 88 min (59), similar to the *t*_{1/2} of 82 min reported with 2% dehydration induced by furosemide (79). The authors thus attributed the prolonged *t*_{1/2} to mild dehydration following an overnight fast. The *t*_{1/2}of acetated Ringer's solution was otherwise reported to average 21 min in euhydrated conscious volunteers (79). When combined with HES, the distribution and elimination of acetated Ringer's solution occurred more slowly than in the single-infusion experiments (59).

Volume kinetic modeling has been used to evaluate IV fluid dynamics using different types of IV fluid solutions (26, 59, 68, 74) and rates (71) in various age groups (106, 107), as well as fluid shifts during inhalation anesthesia (99, 108), epidural anesthesia (109, 110), subarachnoid block (106, 108, 111–113), dehydration (78, 79, 114), hemorrhage (66, 73, 114), acute systemic inflammatory states (72), adrenergic influence (77), perioperative period (115), surgical trauma (96, 116), glucose supplementation (117–119), and pre-eclampsia (120). Through the various VK studies accumulated over the years, the *t*_{1/2} of various IV solutions can now be summarized and were discovered to be significantly variable governed by multiple factors such as the type of fluid and the patient's physiological conditions, further supporting the importance of clinical context (30). Of note, this body of work has been conducted exclusively by the same group of researchers, perhaps due to the intricacy of these analyses, and VK has not found widespread utilization by other researchers. Several international research teams have published on the development of locally modified mathematical fluid kinetic models (121–124) and external validation attempts of these models (125–127). Microcirculatory exchange models that predict fluid, protein, and small ion distribution in the vascular, interstitial, lymphatic, and intracellular compartments using mass balance equations for fluid and individual solutes along with auxiliary transport equations, instead of plasma dilution, have also been explored (22, 128). These model-predicted fluid volume changes were comparable to published experimental data (129) and clinical data obtained using segmental bioelectrical impedance analysis (22).

### Previous Work in Animal Research Models

To date, VK analysis has only been utilized in rabbits (34), pigs (35), and sheep (36–45) which served as experimental research models for human medicine. Sheep and pigs are common biomedical research models for the study of major human physiological systems including the cardiovascular, respiratory, and renal systems due to their well-defined anatomy, physiology, and large body size which permits frequent blood sampling as well as instrumentation with monitoring and sampling devices (130, 131). Translation of experimental work in animal models to humans is feasible, although some species differences exist, such as the reservoir function of the spleen. Thus far, VK modeling and analysis in sheep, rabbits, and pigs appears robust with clinically relevant results that are physiologically sound, providing credibility for wider application to other mammalian species that share similar cardiovascular-renal anatomy and physiologic characteristics.

The distribution and elimination of IV fluids have been investigated under the influence of inhalant anesthesia (36, 38, 42), various fluid infusion rates and duration (37, 39), hemorrhage (41), hypoproteinemia (44), sepsis (34, 40), and use of vasoactive agents (43, 45) in sheep, rabbit, and pig models. Aside from the aforementioned studies performed by the same group of researchers, there is only one veterinary study that investigated real-time IV fluid dynamics using a different method in healthy anesthetized dogs (20). Studies exploring the use of VK in companion animal species could not be identified following a search of the English literature from the MEDLINE^{®} database using search terms such as “volume kinetic,” “fluid kinetic,” “fluid dynamics,” “fluid behavior,” “fluid distribution,” “fluid elimination,” “fluid half-life,” “canine,” “feline,” “dog,” “cat,” and “veterinary” via the PubMed^{®}, Europe PubMed Central^{®} (PMC), and Ovid^{®} search engines. Therefore, VK studies in animal research models and the sole canine study form the only basis for additional discussion of fluid dynamics in veterinary medicine.

In 2002, Brauer et al. (36) investigated fluid dynamics of isotonic crystalloid infusion (25 mL/kg 0.9% NaCl over 20 min) in 6 normovolemic splenectomised sheep (mean 42 kg) while under isoflurane anesthesia and while conscious. Using a cross-over experimental design, the 2-VOFS kinetic model revealed that isotonic crystalloid solution was rapidly eliminated from *V*_{c} via urinary excretion (median 863 mL, range 604–1122 mL) in conscious sheep; however when anesthetized and mechanically ventilated, urinary excretion was markedly reduced (median 9 mL, range 4–150 mL) resulting in peripheral fluid accumulation (36). Unlike indicator-dilution techniques, VK analysis was able to demonstrate how fluid elimination from *V*_{c} transpired more rapidly than accounted for by urinary excretion, thus rather than being excreted, the fluid distributed and accumulated within the peripheral space (38). Plasma volume expansion was otherwise similar for both groups and reached 40–50% fluid efficacy at the immediate end of the infusion despite the marked differences in fluid distribution and elimination (36). Estimated VK parameters correlated well with parameters measured using the indicator-dilution technique with Evans blue dye in the conscious sheep and *V*_{c} was 1.6 L (~4% of BW) (36). However, clearance (reported as *k*_{r}) predicted by VK analysis (58.5 mL/min) was significantly different from the manual calculation using median urinary excretion (0.6 mL/min) when sheep were anesthetized and mechanically ventilated (36). Estimates of *k*_{r} by VK analysis erroneously underestimated peripheral fluid accumulation by corresponding to the sum of urinary output and extravascular fluid accumulation (36). This led to the discovery that *k*_{r} simply reflects net fluid outflow from *V*_{c} and does not approximate urinary excretion of infused fluid as previously presumed (26, 27, 60, 109). Unfortunately, the combination of mechanical ventilation and isoflurane anesthesia in the study group precluded separation of individual intervention effects.

A subsequent experimental cross-over study by Connolly et al. (38) in 7 normovolemic splenectomised sheep (mean 28 kg) was designed to discern the effects of isoflurane anesthesia from mechanical ventilation. The study effectively demonstrated that isoflurane anesthesia alone was responsible for significantly decreased urinary excretion and promoting ECF accumulation during isotonic crystalloid volume loading (38). This study was conducted very thoroughly using 4 protocols: conscious spontaneous ventilation, conscious mechanical ventilation through tracheostomy tube, anesthetized spontaneous ventilation, and anesthetized mechanical ventilation. In this study, baseline PV (~5% BW) increased rapidly during the infusion of 25 mL/kg 0.9% NaCl over 20 min, achieving a 40% fluid efficacy at the immediate end of the infusion (38). This degree of immediate fluid expansion corroborated previous study findings in normovolemic sheep (36) and human volunteers (66). However, despite the initial PV expansion, rapid decline of volume expansion immediately ensued for a duration of 30 min (14% fluid efficacy) followed by a slower phase of decline until the end of the 3 h experiment (4% fluid efficacy) (38). Interestingly, the major findings of Connolly's study were contradictory to a common expectation that positive pressure ventilation alters circulating blood volume and peripheral fluid accumulation by impeding venous return and changing cardiac output (132). The absence of positive end-expiration pressure along with the use of normal tidal volumes (10–15 mL/kg) in euvolemic sheep with healthy compliant lungs may have caused inadequate dynamic changes in intrathoracic pressure and lung volume to induce appreciable cardiovascular effects and secondary fluid shifts. The physiologic mechanism by which isoflurane inhibits diuretic response to volume load and increases extravascular fluid retention has not been determined but reduced glomerular filtration rate and renal blood flow, as well as the involvement of antidiuretic hormone and atrial natriuretic peptide have been speculated to be possible contributing factors (38). A recent matched case-control study of 23 dogs, anesthetized for an elective orthopedic procedure receiving 10 mL/kg/hr of intravenous lactated Ringer's solution for 4 h, reported similar findings of decreased urine output production (<0.5 mL/kg/hr) and substantial peripheral fluid retention as evidenced by significant BW gain, positive fluid balance, as well as increased TBW and ECF volume measured using bioimpedance spectroscopy (133).

Brauer et al. (37) also studied the impact of infusion duration of 0.9% NaCl (6 mL/kg over 5 min, 24 mL/kg over 20 min) and 7.5% NaCl in 6% dextran solution (HSD) (4 mL/kg over 2 min and 20 min) in 6 conscious splenectomised sheep (mean 36 kg). The maximum arterial plasma dilution at the end of the 5 and 20 min 0.9% NaCl infusion were 10 and 22%, respectively. Meanwhile, maximum arterial plasma dilution after 2 and 20 min of HSD infusions were 24 and 21%, respectively. Therefore, the authors concluded that VK variables obtained during a short infusion can be used to predict the outcome of longer infusions, even if the longer infusion also delivers a larger volume (37). This work involving clinically relevant IV fluid doses suggests that VK modeling conforms to linearity such that kinetic parameters obtained can be used to simulate the outcome of other experiments. Subsequently, Svensen et al. (39) published their findings on the VK effects of various isotonic crystalloid infusion volumes and rates (25, 50, and 100 mL/kg of 0.9% NaCl over 20 min) in 6 conscious, splenectomised sheep. Elimination of isotonic crystalloid solution from *V*_{c} was found to be proportional to the magnitude of plasma dilution regardless of infused volumes and rates, and elimination occurs via expansion of *V*_{p} when renal excretion fails to increase in proportion to the volume of infused fluid (39). Therefore, large and rapid fluid boluses that exceed short-term renal excretory capacity contribute to peripheral fluid accumulation. This study also revealed that markedly supraphysiologic infusion doses exceed the limits of linearity of a VK model.

In 2010, Brauer et al. (44) found that severe acute hypoproteinemia induced by plasmapheresis does not reduce the PV expansion of isotonic crystalloid infusion (27 mL/kg 0.9% NaCl over 20 min) in non-hemorrhaged, non-splenectomised sheep. The depletion of mean total protein concentration from 5.4 to 2.5 g/dL following plasmapheresis resulted in a parallel reduction in mean plasma COP from 20 to 9.6 mmHg (44). Plasma volume expansion reached ~20% at the end of the infusion and stayed at 10–15% during the experiment (44). No difference in the PV expansion and cumulative urinary output was found between hypoproteinemic and normal sheep. This is contrary to the conventional physiological reasoning which predicts that a reduction in plasma oncotic pressure will increase fluid filtration, resulting in accumulation of interstitial fluid, reduced plasma protein concentration in capillary filtrate, and increased lymph flow. Splenic contraction and release of RBCs into the circulation affecting the accuracy of plasma dilution calculations were speculated to contribute to the findings however no evidence of hemoconcentration was documented (44), making this explanation less likely. It is intriguing to know if a plasma COP of 50% is sufficient to maintain normal fluid filtration.

The effects of systemic illness have also been explored using VK modeling. Svensen et al. (34) investigated the volume effect of a balanced crystalloid solution (25 mL/kg acetated Ringer's solution over 30 min) in 10 conscious rabbits (mean 4.4 kg) and found that early endotoxemia alters the body's handling of crystalloid solution (34). The expandable volume, *V*, which represents distribution to the ECF space decreased from 473 ± 37 mL (10% of BW) to 327 ± 54 mL (7.5% of BW), and the rate of elimination (reported as *k*_{r} in the original model) increased from 2.9 ± 0.5 to 5.9 ± 2.8 mL/min during early normotensive endotoxemia (34). Visual comparison of individual plasma dilution-time curves demonstrated markedly variable volume effect of acetated Ringer's solution following endotoxin administration. This study provided a glimpse into the net effects of various pathophysiological responses to endotoxemia on the volume effect of IV crystalloid solution, however individual responses that could account for the observed changes could not be isolated. The *V* of 10% BW obtained in healthy rabbits was in agreement with the size of *V* obtained in healthy humans (26, 27) but smaller than the anatomical ECF volume (~20% of BW) to which crystalloids were expected to distribute across.

Almost a decade later, Svensen et al. (40) found contradictory results in a population of septic sheep. The authors found that the distribution and elimination of isotonic crystalloid infusion (25 mL/kg 0.9% NaCl over 20 min) in 6 splenectomised sheep (mean 42 kg) were unchanged by early or late sepsis induced by *Pseudomonas aeruginosa* bacteremia (40). Plasma volume expansion was 312 ± 50 mL (~29.7% fluid efficacy), 386 ± 34 mL (~36.7% fluid efficacy), and 400 ± 51 mL (~38% fluid efficacy) in the control, early sepsis, and late sepsis group (40). Similar peak PV expansion, plasma dilution-time curves, and volume kinetic parameters were obtained for both the control and septic group. These results were contrary to clinical impressions that crystalloid fluid requirements are greatly increased in septic patients due to increased vascular permeability and rapid loss from the intravascular space. The elimination pattern of isotonic crystalloid in septic sheep was similar to control animals and septic sheep also were able to maintain similar levels of PV expansion in response to fluid infusion. The difference between the septic rabbit model (34) and septic sheep model (40) was postulated to be due to the different hypodynamic and hyperdynamic septic shock states that may have been encountered by the models.

Following these studies, fluid dynamics were explored subsequent to administration of vasoactive agents (dopamine, isoprenaline, phenylephrine) in a population of healthy awake sheep (43). Vasoactive drugs were discovered to markedly change the distribution and elimination of crystalloid fluid, thereby altering the PV expansion, urinary excretion, and the risk of peripheral edema (43). An α_{1}-adrenergic receptor agonist (phenylephrine) was found to promote renal excretion of infused fluid at the expense of fluid distribution to the periphery thus limiting the volume expanding effects of the fluid infusion, while β_{1}-adrenergic receptor agonist (isoprenaline) had the opposite effect. The interaction between fluid administration and vasoactive drug (norephinephrine, phenylephrine, dopamine, and esmolol) delivery was recently further explored in an experimental sepsis model in sheep (45). Results were similar in most respects to the healthy sheep model (43) in that α_{1}-adrenergic receptor stimulation with vasoactive drugs accelerated the distribution and elimination of infused fluid, while β_{1}-adrenergic receptor stimulation retarded the distribution and elimination of infused fluid. Having said that, α_{1}-adrenergic receptor stimulation did not increase urinary excretion (elimination) to the same extent as it did in healthy awake sheep (43) due to the compounding inhibitory effects of sepsis, general anesthesia, and hypovolemia on diuresis, therefore worsening overall peripheral fluid accumulation (45). In addition, the tendency for peripheral fluid accumulation in the septic sheep model was also pronounced due to the virtual absence of fluid redistribution from peripheral tissues back to the central compartment, in particular when phenylephrine was given (45). The results of this study suggest that crystalloid fluids given in the early phase of sepsis have a marked tendency to accumulate in extravascular peripheral tissues, and drugs that exert a strong stimulating effect on β_{1}-adrenergic receptors help to limit this aberrant fluid distribution thereby improving hemodynamics through more effective PV expansion.

Aside from the aforementioned animal research models performed by the same research group, there is only one veterinary study that investigated fluid dynamics of various IV fluid solutions in a prospective crossover experimental study involving 4 healthy anesthetized dogs (range 23–25 kg) (20). Silverstein and colleagues utilized a non-invasive continuous in-line hematocrit monitor (Crit-Line^{TM} IIR hematocrit monitor, In-Line Diagnostics, Kaysville, UT) that optically measured hematocrit along an extracorporeal circuit established between a central venous jugular catheter and a cephalic venous catheter to describe the real-time BV changes over 4 h, following rapid IV fluid administration (20). This monitor provides hematocrit value and calculates percent change in BV every 20 s. Each dog served as its own control and received the following treatments on separate occasions every week: 80 mL/kg of 0.9% NaCl at 150 mL/min, 20 mL/kg of Dextran 70 at 150 mL/min, 30 mL/kg of 6% hetastarch at 150 mL/min, 4 mL/kg of 7.5% NaCl at 1 mL/kg/min, and no IV fluid as control (20). In this study, immediate and rapid increases in BV were described during the infusion of each IV fluid bolus. The volume expansion effect of each IV fluid bolus was directly influenced by the volume of fluid administered, with 0.9% NaCl resulting in the greatest increase in BV (76.4 ± 10.0% change) followed by synthetic colloids (35.9 ± 7.3% change for Dextran 70; 27.2 ± 6.4% change for HES), and finally 7.5% NaCl (17.1 ± 3.2% change) immediately at the end of the infusions (20). Although the volume expanding effect of HS was significantly less than that of any other fluid, its efficiency ratio was the greatest (2.7 ± 0.5) of all the fluids while the efficiency ratio for 0.9% NaCl was the smallest (0.8 ± 0.1). The efficiency ratio for Dextran 70 and HES were 0.9 ± 0.4 and 1.1 ± 0.3 respectively. Upon discontinuing both 0.9% NaCl and HS infusions, the rise in BV ceased immediately with a steep decline in BV for 10 min followed by a more gradual decline thereafter. The fall in BV in the 0.9% NaCl group fell below that of HES by 30 min post-infusion. In contrast, the rise in BV continued for 10 min at the end of colloid infusions and a plateau was observed for the remainder of the experiment. By 240 min, colloidal solutions sustained the greatest volume expansion (25.6 ± 16.1 and 26.6 ± 8.6% for Dextran 70 and HES, respectively) compared to the crystalloid solutions (18.0 ± 9.7 and 2.9 ± 6.1% for 0.9% NaCl and HS, respectively) (20). Similar to VK studies of healthy conscious human volunteers (32), isotonic crystalloid resulted in a substantial BV increase during the infusion period (76.4 ± 10.0% change) and distributed into the ECF compartment by 30 min post-infusion with only ~25% of the delivered volume remaining within the intravascular space (20).

Following publication of this study by Silverstein and colleagues, IV fluid expansion and retention of differing resuscitative fluids have been widely disseminated in the veterinary literature and textbooks and generalized to the feline species. The use of a non-invasive continuous in-line hematocrit monitor in Silverstein's study was an innovative and unique method which obviated the need for repeated invasive blood sampling, the use of sophisticated mathematical modeling, or isotope administration to estimate BV. However, this method requires placement of a central venous catheter, heavy sedation or general anesthesia, and the use of an extracorporeal circuit which is not practical in the clinical setting, except for patients receiving renal replacement therapies.

## Limitations of Volume Kinetics

As mentioned previously, VK has not found widespread adoption in human medicine due to the complexity of these analyses requiring some mathematical and statistical background, combined with an understanding of biology, pharmacology, and physiology. Volume kinetic modeling and analysis involves specialized PK modeling software that utilizes distinctive PK nomenclature and programming languages which may not be intuitive to the general scientific audience. Data collection is simple and straightforward; however, data analysis and interpretation involve a steep learning curve for those who are unfamiliar with these concepts. Training courses and first-hand utilization of the PK modeling software would be able to accelerate this learning curve.

The main restriction to extensive clinical use of VK is the requirement for repetitive invasive Hb sampling (134). Successive blood sampling over a period of time, usually 3–4 h, is necessary to capture and understand the dynamic fluid flux within the body (32). This presents a greater challenge in veterinary patients who are much smaller with unforgiving blood volumes. A non-invasive modality that could continuously trend a patient's Hb level with accuracy and precision is essential for the safe translation of VK investigation to the clinical population. A continuous non-invasive Hb monitoring device (Radical-7^{®} Pulse CO-Oximeter^{®}, SET V7.4.0.9 and SET V.6.0.1, Handheld R.7.7.1.0, D-station R5.1.2.7, Masimo Corp., Irvine, CA, USA) successfully generated plasma dilution curves and useful kinetic data for group-level human VK analyses, but was unsuccessful in replacing invasive blood sampling due to wide between-subject variation and exaggeration of plasma dilution (134–136). This technology, in terms of hardware and proprietary software algorithms, has been continuously improved over the years; thus, it is worthwhile to reexplore its potential in VK analysis. Unfortunately, this non-invasive Hb monitoring technology is not ready for veterinary use following unsuccessful validation in anesthetized dogs (137) and the lack of published studies in cats.

Mass balance technique using laboratory blood Hb changes has been widely used to estimate intravascular volume status (21, 138–140) but its accuracy has recently been criticized (92, 141) despite past validation against the radio-iodinated serum albumin external indicator-dilution technique (21). Given the possible ripple effect of measurement errors on volume estimation (92), collection of duplicate or triplicate Hb samples, adequate mixing of sampled blood (142), and high-precision laboratory analysis are important to reduce between-sample variability that would affect serial plasma dilution determination and subsequent VK analysis (32).

Recommendations from published VK studies include utilizing a safe yet adequate IV infusion volume, e.g., 20–25 mL/kg crystalloid fluid bolus over 30 min, to minimize “noisy” data and estimation errors associated with very small infusion volumes (32, 92). Many fluid dynamic studies have infused fluids over 30 min to several hours (17, 143) and these prolonged infusion times preclude the differentiation of actual volume effects from distribution-excretion effects (20). Lastly, physiologic variations as a result of a change in body position, splenic contraction and sequestration, dehydration, hypotension, and the use of diuretic or adrenergic drugs, can affect extrapolation of RBC mass to vascular volume changes (20, 32) and should be minimized during the study period. Although of little importance in humans (144), catecholamine-mediated splenic contraction is known to cause a transient increase in hematocrit and BV in dogs, cats, and sheep (31, 145–148) making splenectomy an important consideration in the experimental research setting.

## Conclusions

Volume kinetics is an innovative research method that is gaining recognition for its wealth of accumulated evidence in this new era where clinicians are searching for context-sensitive fluid therapy paradigms. Despite its limitations, the pharmacokinetics of IV fluid therapy is still an appealing concept that has the potential to serve as a new research tool in veterinary medicine to provide insights on the distribution and elimination of commonly prescribed IV fluids. This research method is feasible and deserves a thorough investigation in the companion animal species. The detailed time course of IV fluids within the body, volume expansion effect, efficacy, half-life (duration of effect), and body water physiology in different patient populations under various clinical conditions may allow for more evidence-based IV fluid therapy prescriptions for our veterinary patients. Until a less invasive method of serial Hb monitoring has been validated for veterinary patients, VK will likely remain as a fundamental research tool for modeling and simulation of IV fluid therapy.

## Author Contributions

XY wrote the first draft of the manuscript. SB, AB, WM, and RH provided critical revision of the manuscript. All authors read and approved the final version of the manuscript to be published.

## Conflict of Interest

The authors declare that the manuscript was written in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## References

1. Woodcock TE, Woodcock TM. Revised Starling equation and the glycocalyx model of transvascular fluid exchange: an improved paradigm for prescribing intravenous fluid therapy. *Br J Anaesth*. (2012) 108:384–94. doi: 10.1093/bja/aer515

2. Hoste EA, Maitland K, Brudney CS, Mehta R, Vincent JL, Yates D, et al. Four phases of intravenous fluid therapy: a conceptual model. *Br J Anaesth*. (2014) 113:740–7. doi: 10.1093/bja/aeu300

3. Malbrain MLNG, Van Regenmortel N, Owczuk R. It is time to consider the four D's of fluid management. *Anaesthesiol Intensive Ther*. (2015) 47:1–5. doi: 10.5603/AIT.a2015.0070

4. Chappell D, Jacob M. Role of the glycocalyx in fluid management: small things matter. *Best Pract Res Clin Anaesthesiol*. (2014) 28:227–34. doi: 10.1016/j.bpa.2014.06.003

5. Monnet X, Marik PE, Teboul J-L. Prediction of fluid responsiveness: an update. *Ann Intensive Care*. (2016) 6:111. doi: 10.1186/s13613-016-0216-7

6. Raghunathan K, Shaw AD, Bagshaw SM. Fluids are drugs: type, dose and toxicity. *Curr Opin Crit Care*. (2013) 19:290–8. doi: 10.1097/MCC.0b013e3283632d77

7. McDermid RC, Raghunathan K, Romanovsky A, Shaw AD, Bagshaw SM. Controversies in fluid therapy: type, dose and toxicity. *World J Crit Care Med*. 2014 Feb (4) 3:24–33. doi: 10.5492/wjccm.v3.i1.24

8. Cazzolli D, Prittie J. The crystalloid-colloid debate: consequences of resuscitation fluid selection in veterinary critical care. *J Vet Emerg Crit Care*. (2015) 25:6–19. doi: 10.1111/vec.12281

9. Muir WW, Ueyama Y, Noel-Morgan J, Kilborne A, Page J. A systematic review of the quality of IV fluid therapy in veterinary medicine. *Front Vet Sci*. (2017) 4:127. doi: 10.3389/fvets.2017.00127

10. Davis H, Jensen T, Johnson A, Knowles P, Meyer R, Rucinsky R, et al. 2013 AAHA/AAFP fluid therapy guidelines for dogs and cats. *J Am Anim Hosp Assoc*. (2013) 49:149–59. doi: 10.5326/JAAHA-MS-5868

11. Glover PA, Rudloff E, Kirby R. Hydroxyethyl starch: a review of pharmacokinetics, pharmacodynamics, current products, and potential clinical risks, benefits, and use. *J Vet Emerg Crit Care*. (2014) 24:642–61. doi: 10.1111/vec.12208

12. Cavanagh AA, Sullivan LA, Hansen BD. Retrospective evaluation of fluid overload and relationship to outcome in critically ill dogs. *J Vet Emerg Crit Care*. (2016) 26:578–86. doi: 10.1111/vec.12477

13. Ostroski CJ, Drobatz KJ, Reineke EL. Retrospective evaluation of and risk factor analysis for presumed fluid overload in cats with urethral obstruction: 11 cases (2002-2012). *J Vet Emerg Crit Care*. (2017) 27:561–8. doi: 10.1111/vec.12631

14. Maitland K, Kiguli S, Opoka RO, Engoru C, Olupot-Olupot P, Akech SO, et al. Mortality after fluid bolus in African children with severe infection. *N Engl J Med*. (2011) 364:2483–95. doi: 10.1056/NEJMoa1101549

15. Malbrain MLNG, Van Regenmortel N, Saugel B, De Tavernier B, Van Gaal P-J, Joannes-Boyau O, et al. Principles of fluid management and stewardship in septic shock: it is time to consider the four D's and the four phases of fluid therapy. *Ann Intensive Care*. (2018) 8:66. doi: 10.1186/s13613-018-0402-x

16. Tatara T. Context-sensitive fluid therapy in critical illness. *J Intensive Care*. (2016) 4:20. doi: 10.1186/s40560-016-0150-7

17. Lamke LO, Liljedahl SO. Plasma volume changes after infusion of various plasma expanders. *Resuscitation*. (1976) 5:93–102. doi: 10.1016/0300-9572(76)90029-0

18. Jones JG, Wardrop CA. Measurement of blood volume in surgical and intensive care practice. *Br J Anaesth*. (2000) 84:226–35. doi: 10.1093/oxfordjournals.bja.a013407

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

20. Silverstein DC, Aldrich J, Haskins SC, Drobatz KJ, Cowgill LD. Assessment of changes in blood volume in response to resuscitative fluid administration in dogs. *J Vet Emerg Crit Care*. (2005) 15:185–92. doi: 10.1111/j.1476-4431.2005.00138.x

21. Hahn RG. A haemoglobin dilution method (HDM) for estimation of blood volume variations during transurethral prostatic surgery. *Acta Anaesthesiol Scand*. (1987) 31:572–8. doi: 10.1111/j.1399-6576.1987.tb02623.x

22. Tatara T, Tashiro C. Quantitative analysis of fluid balance during abdominal surgery. *Anesth Analg*. (2007) 104:347–54. doi: 10.1213/01.ane.0000253031.70916.27

23. da Silva AT, Hauschild DB, de Almeida Oliveira LD, de Fragas Hinnig P, Franco Moreno YM, Wazlawik E. Association of hyperhydration evaluated by bioelectrical impedance analysis and mortality in patients with different medical conditions: systematic review and meta-analyses. *Clin Nutr ESPEN*. (2018) 28:12–20. doi: 10.1016/j.clnesp.2018.08.022

24. Asklöf M, Kjølhede P, Wodlin NB, Nilsson L. Bioelectrical impedance analysis; a new method to evaluate lymphoedema, fluid status, and tissue damage after gynaecological surgery - a systematic review. *Eur J Obstet Gynecol Reprod Biol*. (2018) 228:111–9. doi: 10.1016/j.ejogrb.2018.06.024

25. Stahle L, Nilsson A, Hahn RG. Modelling the volume of expandable body fluid spaces during I.V. fluid therapy. *Br J Anaesth*. (1997) 78:138–43. doi: 10.1093/bja/78.2.138

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

27. Hahn RG, Drobin D, Stahle L. Volume kinetics of Ringer's solution in female volunteers. *Br J Anaesth*. (1997) 78:144–8. doi: 10.1093/bja/78.2.144

28. Hahn RG, Svensén C. Plasma dilution and the rate of infusion of Ringer's solution. *Br J Anaesth*. (1997) 79:64–7. doi: 10.1093/bja/79.1.64

29. Stanski DR. The pharmacokinetics of intravenous fluids. *Anesthesiology*. (1997) 87:200–1. doi: 10.1097/00000542-199708000-00003

30. Hahn RG, Lyons G. The half-life of infusion fluids: an educational review. *Eur J Anaesthesiol*. (2016) 33:475–82. doi: 10.1097/EJA.0000000000000436

31. Svensen CH, Rodhe PM, Prough DS. Pharmacokinetic aspects of fluid therapy. *Best Pract Res Clin Anaesthesiol*. (2009) 23:213–24. doi: 10.1016/j.bpa.2008.11.003

32. Hahn RG. Volume kinetics for infusion fluids. *Anesthesiology*. (2010) 113:470–81. doi: 10.1097/ALN.0b013e3181dcd88f

33. Hahn RG. Understanding volume kinetics. *Acta Anaesthesiol Scand*. (2020) 64:570–8. doi: 10.1111/aas.13533

34. Svensen C, Hjelmqvist H, Hahn RG. Volume kinetics of Ringer solution during endotoxinemia in conscious rabbits. *J Endotoxin Res*. (1997) 4:425–30. doi: 10.1177/096805199700400606

35. Hahn RG, Gebäck T. Fluid volume kinetics of dilutional hyponatremia; a shock syndrome revisited. *Clinics*. (2014) 69:120–7. doi: 10.6061/clinics/2014(02)08

36. Brauer KI, Svensen C, Hahn RG, Traber LD, Prough DS. Volume kinetic analysis of the distribution of 0.9% saline in conscious versus isoflurane-anesthetized sheep. *Anesthesiology*. (2002) 96:442–9. doi: 10.1097/00000542-200202000-00033

37. Brauer LP, Svensen CH, Hahn RG, Kilicturgay S, Kramer GC, Prough DS. Influence of rate and volume of infusion on the kinetics of 0.9% saline and 7.5% saline/6.0% dextran 70 in sheep. *Anesth Analg*. (2002) 95:1547–56. doi: 10.1097/00000539-200212000-00014

38. Connolly CM, Kramer GC, Hahn RG, Chaisson NF, Svensen CH, Kirschner RA, et al. Isoflurane but not mechanical ventilation promotes extravascular fluid accumulation during crystalloid volume loading. *Anesthesiology*. (2003) 98:670–81. doi: 10.1097/00000542-200303000-00015

39. Svensen CH, Brauer KP, Hahn RG, Uchida T, Traber LD, Traber DL, et al. Elimination rate constant describing clearance of infused fluid from plasma is independent of large infusion volumes of 0.9% saline in sheep. *Anesthesiology*. (2004) 101:666–74. doi: 10.1097/00000542-200409000-00015

40. Svensen CH, Clifton B, Brauer KI, Olsson J, Uchida T, Traber LD, et al. Sepsis produced by *Pseudomonas bacteremia* does not alter plasma volume expansion after 0.9% saline infusion in sheep. *Anesth Analg*. (2005) 101:835–42. doi: 10.1213/01.ane.0000166763.24782.26

41. Norberg A, Brauer KI, Prough DS, Gabrielsson J, Hahn RG, Uchida T, et al. Volume turnover kinetics of fluid shifts after hemorrhage, fluid infusion, and the combination of hemorrhage and fluid infusion in sheep. *Anesthesiology*. (2005) 102:985–94. doi: 10.1097/00000542-200505000-00018

42. Hahn RG, Brauer L, Rodhe P, Svensen CH, Prough DS. Isoflurane inhibits compensatory intravascular volume expansion after hemorrhage in sheep. *Anesth Analg*. (2006) 103:350–8. doi: 10.1213/01.ane.0000222469.94134.ba

43. Ewaldsson C-A, Vane LA, Kramer GC, Hahn RG. Adrenergic drugs alter both the fluid kinetics and the hemodynamic responses to volume expansion in sheep. *J Surg Res*. (2006) 131:7–14. doi: 10.1016/j.jss.2005.09.012

44. Brauer KI, Brauer LP, Prough DS, Rodhe P, Hahn RG, Traber DL, et al. Hypoproteinemia does not alter plasma volume expansion in response to a 0.9% saline bolus in awake sheep. *Crit Care Med*. (2010) 38:2011–5. doi: 10.1097/CCM.0b013e3181ef45e2

45. Li Y, Xiaozhu Z, Guomei R, Qiannan D, Hahn RG. Effects of vasoactive drugs on crystalloid fluid kinetics in septic sheep. *PLoS ONE*. (2017) 12:e0172361. doi: 10.1371/journal.pone.0172361

46. Loucks J, Yost S, Kaplan B. An introduction to basic pharmacokinetics. *Transplantation*. (2015) 99:903–7. doi: 10.1097/TP.0000000000000754

47. Nelson E. Kinetics of drug absorption, distribution, metabolism, and excretion. *J Pharm Sci*. (1961) 50:181–92. doi: 10.1002/jps.2600500302

48. Committee for Proprietary Medicinal Products (CPMP). *Note for Guidance on the Investigation of Bioavailability and Bioequivalence [Internet]*. London: The European Agency for the Evaluation of Medicinal Products (2001). p. 1–19. Report No.: CPMP/EWP/QWP/1401/98. Available online at : https://www.ema.europa.eu/en/documents/scientific-guideline/note-guidance-investigation-bioavailability-bioequivalence_en.pdf (accessed September 26, 2020).

49. Toutain PL, Bousquet-Mélou A. Volumes of distribution. *J Vet Pharmacol Ther*. (2004) 27:441–53. doi: 10.1111/j.1365-2885.2004.00602.x

50. Rowland M, Tozer TN. Fundamental concepts and terminology. In: *Clinical Pharmacokinetics and Pharmacodynamics: Concepts and Applications*. 4th ed. Baltimore, MD: Wolters Kluwer Health; Lippincott William & Wilkins (2011). p. 17–45.

51. DiPiro JT, Spruill WJ, Wade WE, Blouin RA, Pruemer JM. Introduction to pharmacokinetics and pharmacodynamics. In: *Concepts in Clinical Pharmacokinetics*. 5th ed. Bethesda, MD: American Society of Health-System Pharmacists (2010). p. 1–18.

52. Atkinson AJ, Ruo TI, Frederiksen MC. Physiological basis of multicompartmental models of drug distribution. *Trends Pharmacol Sci*. (1991) 12:96–101. doi: 10.1016/0165-6147(91)90515-T

53. Dhillon S, Gill K. Basic pharmacokinetics. In: Dhillon S, Kostrzewski A, editors. *Clinical Pharmacokinetics*. 1st ed. London: Pharmaceutical Press (2006). p. 1–44.

54. Peck TE, Williams M. Mathematics and pharmacokinetics. In: *Pharmacology for Anaesthesia and Intensive Care*. 1st ed. London: Cambridge University Press (2000). p. 47–66.

55. Hahn RG. The use of volume kinetics to optimize fluid therapy. *J Trauma*. (2003) 54:155–8. doi: 10.1097/01.TA.0000064515.51896.DE

56. Mould DR, Upton RN. Basic concepts in population modeling, simulation, and model-based drug development. *CPT Pharmacometrics Syst Pharmacol*. (2012) 1:e6. doi: 10.1038/psp.2012.4

57. Riviere JE, Gabrielsson J, Fink M, Mochel J. Mathematical modeling and simulation in animal health. Part I: moving beyond pharmacokinetics. *J Vet Pharmacol Ther*. (2016) 39:213–23. doi: 10.1111/jvp.12278

58. Bon C, Toutain PL, Concordet D, Gehring R, Martin-Jimenez T, Smith J, et al. Mathematical modeling and simulation in animal health. Part III: using nonlinear mixed-effects to characterize and quantify variability in drug pharmacokinetics. *J Vet Pharmacol Ther*. (2018) 41:171–83. doi: 10.1111/jvp.12473

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

60. Hahn RG, Drobin D. Urinary excretion as an input variable in volume kinetic analysis of Ringer's solution. *Br J Anaesth*. (1998) 80:183–8. doi: 10.1093/bja/80.2.183

61. Boxenbaum HG, Riegelman S, Elashoff RM. Statistical estimations in pharmacokinetics. *J Pharmacokinet Biopharm*. (1974) 2:123–48. doi: 10.1007/BF01061504

62. Motulsky HJ, Ransnas LA. Fitting curves to data using nonlinear regression: a practical and nonmathematical review. *FASEB J*. (1987) 1:365–74. doi: 10.1096/fasebj.1.5.3315805

63. Wada DR, Stanski DR, Ebling WF. A PC-based graphical simulator for physiological pharmacokinetic models. *Comput Methods Programs Biomed*. (1995) 46:245–55. doi: 10.1016/0169-2607(95)01624-3

64. Easterling MR, Evans MV, Kenyon EM. Comparative analysis of software for physiologically based pharmacokinetic modeling: simulation, optimization, and sensitivity analysis. *Toxicol Methods*. (2000) 10:203–29. doi: 10.1080/10517230050121615

65. Kiang TKL, Sherwin CMT, Spigarelli MG, Ensom MHH. Fundamentals of population pharmacokinetic modelling: modelling and software. *Clin Pharmacokinet*. (2012) 51:515–25. doi: 10.1007/BF03261928

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

67. Drobin D. A single-model solution for volume kinetic analysis of isotonic fluid infusions. *Acta Anaesthesiol Scand*. (2006) 50:1074–80. doi: 10.1111/j.1399-6576.2006.01117.x

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

69. Norberg A, Hahn RG, Li H, Olsson J, Prough DS, Børsheim E, et al. Population volume kinetics predicts retention of 0.9% saline infused in awake and isoflurane-anesthetized volunteers. *Anesthesiology*. (2007) 107:24–32. doi: 10.1097/01.anes.0000268387.34758.6d

70. Hahn RG. The elimination half-life of crystalloid fluid is shorter in female than in male volunteers: a retrospective population kinetic analysis. *Biol Sex Differ*. (2016) 7:54. doi: 10.1186/s13293-016-0105-7

71. Hahn RG, Drobin D, Zdolsek J. Distribution of crystalloid fluid changes with the rate of infusion: a population-based study. *Acta Anaesthesiol Scand*. (2016) 60:569–78. doi: 10.1111/aas.12686

72. Li Y, Yi S, Zhu Y, Hahn RG. Volume kinetics of Ringer's lactate solution in acute inflammatory disease. *Br J Anaesth*. (2018) 121:574–80. doi: 10.1016/j.bja.2018.04.023

73. Nyberg J, Li H, Wessmark P, Winther V, Prough DS, Kinsky MP, et al. Population kinetics of 0.9% saline distribution in hemorrhaged awake and isoflurane-anesthetized volunteers. *Anesthesiology*. (2019) 131:501–11. doi: 10.1097/ALN.0000000000002816

74. Hahn RG, Zdolsek M, Hasselgren E, Zdolsek J, Björne H. Fluid volume kinetics of 20% albumin. *Br J Clin Pharmacol*. (2019) 85:1303–11. doi: 10.1111/bcp.13897

75. Baggot JD. Pharmacokinetic terms: symbols and units. *J Vet Pharmacol Ther*. (2001) 24:81–2. doi: 10.1046/j.1365-2885.2001.00340.x

77. Li YH, Zhu HB, Zheng X, Chen HJ, Shao L, Hahn RG. Low doses of esmolol and phenylephrine act as diuretics during intravenous anesthesia. *Crit Care*. (2012) 16:R18. doi: 10.1186/cc11175

78. Hahn RG, Bahlmann H, Nilsson L. Dehydration and fluid volume kinetics before major open abdominal surgery. *Acta Anaesthesiol Scand*. (2014) 58:1258–66. doi: 10.1111/aas.12416

79. Zdolsek J, Li Y, Hahn RG. Detection of dehydration by using volume kinetics. *Anesth Analg*. (2012) 115:814–22. doi: 10.1213/ANE.0b013e318261f6ba

80. Hall JE. The body fluid compartments: extracellular and intracellular fluids; edema. In: *Guyton and Hall Textbook of Medical Physiology*. 13th ed. Philadelphia, PA: Elsevier (2016). p. 305–20.

81. Zanghi BM, Cupp CJ, Pan Y, Tissot-Favre DG, Milgram NW, Nagy TR, et al. Noninvasive measurements of body composition and body water via quantitative magnetic resonance, deuterium water, and dual-energy x-ray absorptiometry in cats. *Am J Vet Res*. (2013) 74:721–32. doi: 10.2460/ajvr.74.5.721

82. Foy JM, Schnieden H. Estimation of total body water (virtual tritium space) in the rat, cat, rabbit, guinea-pig and man, and of the biological half-life of tritium in man. *J Physiol*. (1960) 154:169–76. doi: 10.1113/jphysiol.1960.sp006571

83. Zdolsek JH, Lisander B, Hahn RG. Measuring the size of the extracellular fluid space using bromide, iohexol, and sodium dilution. *Anesth Analg*. (2005) 101:1770–7. doi: 10.1213/01.ANE.0000184043.91673.7E

84. Finch NC, Heiene R, Elliott J, Syme HM, Peters AM. Determination of extracellular fluid volume in healthy and azotemic cats. *J Vet Intern Med*. (2015) 29:35–42. doi: 10.1111/jvim.12506

85. Geankoplis CJ, Hersel AA, Lepek DH. Introduction to engineering principles and units. In: *Transport Processes and Separation Process Principles*. 5th ed. Hoboken, NJ: Prentice Hall (2018). p. 13 (Conservation of mass and material balances).

86. Rehm M, Haller M, Orth V, Kreimeier U, Jacob M, Dressel H, et al. Changes in blood volume and hematocrit during acute preoperative volume loading with 5% albumin or 6% hetastarch solutions in patients before radical hysterectomy. *Anesthesiology*. (2001) 95:849–56. doi: 10.1097/00000542-200110000-00011

87. Gore CJ, Hopkins WG, Burge CM. Errors of measurement for blood volume parameters: a meta-analysis. *J Appl Physiol*. (2005) 99:1745–58. doi: 10.1152/japplphysiol.00505.2005

88. Jacob M, Conzen P, Finsterer U, Krafft A, Becker BF, Rehm M. Technical and physiological background of plasma volume measurement with indocyanine green: a clarification of misunderstandings. *J Appl Physiol*. (2007) 102:1235–42. doi: 10.1152/japplphysiol.00740.2006

89. Polidori D, Rowley C. Optimal back-extrapolation method for estimating plasma volume in humans using the indocyanine green dilution method. *Theor Biol Med Model*. (2014) 11:33. doi: 10.1186/1742-4682-11-33

90. Hahn RG. Must hypervolaemia be avoided? A critique of the evidence. *Anaesthesiol Intensive Ther*. (2015) 47:449–56. doi: 10.5603/AIT.a2015.0062

91. Ewaldsson C-A, Hahn RG. Kinetics and extravascular retention of acetated ringer's solution during isoflurane or propofol anesthesia for thyroid surgery. *Anesthesiology*. (2005) 103:460–9. doi: 10.1097/00000542-200509000-00006

92. Wolf MB. Hemoglobin-dilution method: effect of measurement errors on vascular volume estimation. *Comput Math Methods Med*. (2017) 2017:3420590. doi: 10.1155/2017/3420590

93. Boisvert AM, Tvedten HW, Scott MA. Artifactual effects of hypernatremia and hyponatremia on red cell analytes measured by the Bayer H^{*}1 analyzer. *Vet Clin Pathol*. (1999) 28:91–6. doi: 10.1111/j.1939-165X.1999.tb01055.x

94. Dugger DT, Mellema MS, Hopper K, Epstein SE. Comparative accuracy of several published formulae for the estimation of serum osmolality in cats. *J Small Anim Pract*. (2013) 54:184–9. doi: 10.1111/jsap.12046

95. Dugger DT, Epstein SE, Hopper K, Mellema MS. A comparison of the clinical utility of several published formulae for estimated osmolality of canine serum. *J Vet Emerg Crit Care*. (2014) 24:188–93. doi: 10.1111/vec.12153

96. Svensén C, Ponzer S, Hahn RG. Volume kinetics of Ringer solution after surgery for hip fracture. *Can J Anaesth*. (1999) 46:133–41. doi: 10.1007/BF03012547

97. Aukland K, Reed RK. Interstitial-lymphatic mechanisms in the control of extracellular fluid volume. *Physiol Rev*. (1993) 73:1–78. doi: 10.1152/physrev.1993.73.1.1

98. Bhave G, Neilson EG. Body fluid dynamics: back to the future. *J Am Soc Nephrol*. (2011) 22:2166–81. doi: 10.1681/ASN.2011080865

99. Hahn RG. Arterial pressure and the rate of elimination of crystalloid fluid. *Anesth Analg*. (2017) 124:1824–33. doi: 10.1213/ANE.0000000000002075

100. Hahn RG. Influences of red blood cell and platelet counts on the distribution and elimination of crystalloid fluid. *Medicina*. (2017) 53:233–41. doi: 10.1016/j.medici.2017.07.005

101. Atlas G, Li JK-J, Amin S, Hahn RG. Development and retrospective clinical assessment of a patient-specific closed-form integro-differential equation model of plasma dilution. *Biomed Eng Comput Biol*. (2017) 8:1–20. doi: 10.1177/1179597217730305

102. Zdolsek M, Hahn RG, Zdolsek JH. Recruitment of extravascular fluid by hyperoncotic albumin. *Acta Anaesthesiol Scand*. (2018) 62:1255–60. doi: 10.1111/aas.13150

103. Hahn RG. Fluid therapy might be more difficult than you think. *Anesth Analg*. (2007) 105:304–5. doi: 10.1213/01.ane.0000270218.31147.67

104. Hahn RG. Volume effect of Ringer's solution in the blood during general anaesthesia. *Eur J Anaesthesiol*. (1998) 15:427–32. doi: 10.1046/j.1365-2346.1998.00313.x

105. *Voluven ^{®} 6% Hydroxyethyl Starch 130/0,.4 in 0.9% Sodium Chloride Injection*. [Internet]. Toronto, ON: Fresenius Kabi Canada Ltd. (2020). Available online at: https://www.fresenius-kabi.com/en-ca/documents/Product-Monograph-Voluven-022520.pdf (accessed May 10, 2020).

106. Li Y, Zhu S, Hahn RG. The kinetics of Ringer's solution in young and elderly patients during induction of general anesthesia with propofol and epidural anesthesia with ropivacaine. *Acta Anaesthesiol Scand*. (2007) 51:880–7. doi: 10.1111/j.1399-6576.2007.01351.x

107. Li Y, Hahn RG, Hu Y, Xiang Y, Zhu S. Plasma and renal clearances of lactated Ringer's solution in pediatric and adult patients just before anesthesia is induced. *Paediatr Anaesth*. (2009) 19:682–7. doi: 10.1111/j.1460-9592.2009.03047.x

108. Ewaldsson C-A, Hahn RG. Bolus injection of Ringer's solution and dextran 1 kDa during induction of spinal anesthesia. *Acta Anaesthesiol Scand*. (2005) 49:152–9. doi: 10.1111/j.1399-6576.2004.00561.x

109. Drobin D, Hahn RG. Time course of increased haemodilution in hypotension induced by extradural anaesthesia. *Br J Anaesth*. (1996) 77:223–6. doi: 10.1093/bja/77.2.223

110. Hahn RG. Increased haemodilution in hypotension induced by epidural anaesthesia. *Acta Anaesthesiol Scand*. (1993) 37:357–60. doi: 10.1111/j.1399-6576.1993.tb03728.x

111. Hahn RG, Resby M. Volume kinetics of Ringer's solution and dextran 3% during induction of spinal anaesthesia for caesarean section. *Can J Anaesth*. (1998) 45:443–51. doi: 10.1007/BF03012580

112. Ewaldsson C-A, Hahn RG. Volume kinetics of Ringer's solution during induction of spinal and general anaesthesia. *Br J Anaesth*. (2001) 87:406–14. doi: 10.1093/bja/87.3.406

113. Hahn RG, Lindahl CC, Drobin D. Volume kinetics of acetated Ringer's solution during experimental spinal anaesthesia. *Acta Anaesthesiol Scand*. (2011) 55:987–94. doi: 10.1111/j.1399-6576.2011.02493.x

114. Hahn RG, Drobin D, Li Y, Zdolsek J. Kinetics of Ringer's solution in extracellular dehydration and hemorrhage. *Shock*. (2020) 53:566–73. doi: 10.1097/SHK.0000000000001422

115. Holte K, Hahn RG, Ravn L, Bertelsen KG, Hansen S, Kehlet H. Influence of “liberal” versus “restrictive” intraoperative fluid administration on elimination of a postoperative fluid load. *Anesthesiology*. (2007) 106:75–9. doi: 10.1097/00000542-200701000-00014

116. Olsson J, Svensen CH, Hahn RG. The volume kinetics of acetated Ringer's solution during laparoscopic cholecystectomy. *Anesth Analg*. (2004) 99:1854–60. doi: 10.1213/01.ANE.0000134809.07605.3C

117. Sjöstrand F, Hahn RG. Validation of volume kinetic analysis of glucose 2.5% solution given by intravenous infusion. *Br J Anaesth*. (2003) 90:600–7. doi: 10.1093/bja/aeg102

118. Sjöstrand F, Hahn RG. Volume kinetics of glucose 2.5% solution during laparoscopic cholecystectomy. *Br J Anaesth*. (2004) 92:485–92. doi: 10.1093/bja/aeh095

119. Hahn RG, Edsberg L, Sjöstrand F. Volume kinetic analysis of fluid shifts accompanying intravenous infusions of glucose solution. *Cell Biochem Biophys*. (2003) 39:211–22. doi: 10.1385/CBB:39:3:211

120. Drobin D, Hahn RG. Distribution and elimination of crystalloid fluid in pre-eclampsia. *Clin Sci*. (2004) 106:307–13. doi: 10.1042/CS20030349

121. Lee EH, Lee JY, Yeo YG, Choi KT. Modelling the expandable body fluid space after I.V. fluid infusion. *Korean J Anesthesiol*. (2006) 51:36. doi: 10.4097/kjae.2006.51.1.36

122. Choi KT, Yeo YK. Kinetics of volume expansion during infusion of Ringer's solution based on single volume model. *Korean J Chem Eng*. (2006) 23:610–6. doi: 10.1007/BF02706803

123. Choi KT, Lee EH, Yeo YK. Kinetics of volume expansion during infusion of Ringer's solution based on two volume model. *Korean J Chem Eng*. (2006) 23:867–73. doi: 10.1007/s11814-006-0001-7

124. Lee J-H, Choo Y-J, Lee Y-H, Rhim J-H, Lee S-H, Choi BM, et al. Population-based volume kinetics of Ringer's lactate solution in patients undergoing open gastrectomy. *Acta Pharmacol Sin*. (2019) 40:710–6. doi: 10.1038/s41401-018-0179-3

125. Kang SJ, Choi KT, Lee JH, Lee DH. Assessment of fluid space after infusion of Hartmann's solution and hydroxethyl starch solutions during general anesthesia. *Korean J Anesthesiol*. (2008) 54:300. doi: 10.4097/kjae.2008.54.3.300

126. Lee EH, Kim S-K, Yeo YG, Choi KT. Effects of anesthesia on fluid volume kinetics after infusion of colloid solution during blood donation. *Korean J Anesthesiol*. (2010) 58:514–20. doi: 10.4097/kjae.2010.58.6.514

127. Choi BM, Karm MH, Jung KW, Yeo YG, Choi KT. The predictive performance of infusion strategy nomogram based on a fluid kinetic model. *Korean J Anesthesiol*. (2015) 68:128–35. doi: 10.4097/kjae.2015.68.2.128

128. Gyenge CC, Bowen BD, Reed RK, Bert JL. Transport of fluid and solutes in the body I. Formulation of a mathematical model. *Am J Physiol*. (1999) 277:1215–27. doi: 10.1152/ajpheart.1999.277.3.H1215

129. Gyenge CC, Bowen BD, Reed RK, Bert JL. Transport of fluid and solutes in the body II. Model validation and implications. *Am J Physiol*. (1999) 277:1228–40. doi: 10.1152/ajpheart.1999.277.3.H1228

130. Scheerlinck J-PY, Snibson KJ, Bowles VM, Sutton P. Biomedical applications of sheep models: from asthma to vaccines. *Trends Biotechnol*. (2008) 26:259–66. doi: 10.1016/j.tibtech.2008.02.002

131. Swindle MM, Makin A, Herron AJ, Clubb FJ, Frazier KS. Swine as models in biomedical research and toxicology testing. *Vet Pathol*. (2012) 49:344–56. doi: 10.1177/0300985811402846

132. Pinsky MR. Cardiopulmonary interactions: physiologic basis and clinical applications. *Ann Am Thorac Soc*. (2018) 15:45–8. doi: 10.1513/AnnalsATS.201704-339FR

133. Boscan P, Pypendop BH, Siao KT, Francey T, Dowers K, Cowgill L, et al. Fluid balance, glomerular filtration rate, and urine output in dogs anesthetized for an orthopedic surgical procedure. *Am J Vet Res*. (2010) 71:501–7. doi: 10.2460/ajvr.71.5.501

134. Hahn RG, Li Y, Zdolsek J. Non-invasive monitoring of blood haemoglobin for analysis of fluid volume kinetics. *Acta Anaesthesiol Scand*. (2010) 54:1233–40. doi: 10.1111/j.1399-6576.2010.02321.x

135. Sjöstrand F, Rodhe P, Berglund E, Lundström N, Svensen C. The use of a noninvasive hemoglobin monitor for volume kinetic analysis in an emergency room setting. *Anesth Analg*. (2013) 116:337–42. doi: 10.1213/ANE.0b013e318277dee3

136. Lian Q, Li H, Zeng R, Lang J, Shangguan W, Liu H, et al. The use of a noninvasive hemoglobin monitor for determining fluid distribution and elimination in pediatric patients undergoing minor surgery. *J Clin Monit Comput*. (2015) 29:41–6. doi: 10.1007/s10877-014-9566-6

137. Read MR, Rondeau J, Kwong GPS. Validation of noninvasive hemoglobin measurements using co-oximetry in anesthetized dogs. *Can Vet J*. (2016) 57:1161–5.

138. Dill DB, Costill DL. Calculation of percentage changes in volumes of blood, plasma, and red cells in dehydration. *J Appl Physiol*. (1974) 37:247–8. doi: 10.1152/jappl.1974.37.2.247

139. Brecher ME, Monk T, Goodnough LT. A standardized method for calculating blood loss. *Transfusion*. (1997) 37:1070–4. doi: 10.1046/j.1537-2995.1997.371098016448.x

140. Dasselaar JJ, Lub-de Hooge MN, Pruim J, Nijnuis H, Wiersum A, de Jong PE, et al. Relative blood volume changes underestimate total blood volume changes during hemodialysis. *Clin J Am Soc Nephrol*. (2007) 2:669–74. doi: 10.2215/CJN.00880207

141. Meunier A, Petersson A, Good L, Berlin G. Validation of a haemoglobin dilution method for estimation of blood loss. *Vox Sang*. (2008) 95:120–4. doi: 10.1111/j.1423-0410.2008.01071.x

142. Narayanan S. *Preanalytical Issues Related to Blood Sample Mixing [Internet]*. Brønshøj: acutecaretesting.org (2005). p. 1–6. Available online at: https://acutecaretesting.org/-/media/acutecaretesting/files/pdf/preanalytical-issues-related-to-blood-sample-mixing.pdf (accessed Apr 21, 2019).

143. Solanke TF, Khwaja MS, Madojemu EI. Plasma volume studies with four different plasma volume expanders. *J Surg Res*. (1971) 11:140–3. doi: 10.1016/0022-4804(71)90038-2

144. Stewart IB, McKenzie DC. The human spleen during physiological stress. *Sports Med*. (2002) 32:361–9. doi: 10.2165/00007256-200232060-00002

145. Carneiro JJ, Donald DE. Blood reservoir function of dog spleen, liver, and intestine. *Am J Physiol*. (1977) 232:67–72. doi: 10.1152/ajpheart.1977.232.1.H67

146. Breznock EM, Strack D. Blood volume of nonsplenectomized and splenectomized cats before and after acute hemorrhage. *Am J Vet Res*. (1982) 43:1811–4.

147. Breznock EM, Strack D. Effects of the spleen, epinephrine, and splenectomy on determination of blood volume in cats. *Am J Vet Res*. (1982) 43:2062–6.

148. Ojiri Y, Noguchi K, Chibana T, Sakanashi M. Effects of adrenergic stimulants on the splenic diameter, haemoglobin content and haematocrit in anaesthetized dogs: determination of the adrenoceptor subtype responsible for changes in the splenic diameter. *Acta Physiol Scand*. (1993) 149:31–9. doi: 10.1111/j.1748-1716.1993.tb09589.x

## Nomenclature

Keywords: fluid therapy, pharmacokinetics, volume kinetics, hemoglobin dilution, urine output, distribution, elimination, half-life

Citation: Yiew XT, Bateman SW, Hahn RG, Bersenas AME and Muir WW (2020) Understanding Volume Kinetics: The Role of Pharmacokinetic Modeling and Analysis in Fluid Therapy. *Front. Vet. Sci.* 7:587106. doi: 10.3389/fvets.2020.587106

Received: 24 July 2020; Accepted: 12 October 2020;

Published: 20 November 2020.

Edited by:

Pablo Martín-Vasallo, University of La Laguna, SpainReviewed by:

Jonathan Paul Mochel, Iowa State University, United StatesPierre-Louis Toutain, Ecole Nationale Vétérinaire de Toulouse, France

Copyright © 2020 Yiew, Bateman, Hahn, Bersenas and Muir. 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) and the copyright owner(s) 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: Xiu Ting Yiew, xyiew@uoguelph.ca