Skip to main content

Volume kinetic evaluation of fluid turnover after oral intake of tap water, lemonade and saline in volunteers

Abstract

Background

Volume kinetic fluid turnover of three beverages was investigated for the purpose of estimating their rates of absorption and elimination as well as their maximum effect on the blood volume. The results were then used to simulate the effects of ingesting different combinations of these fluids.

Method

Ten healthy volunteers ingested 0.5 L of tap water, lemonade (90 g/L carbohydrates) and isotonic saline (9 g/L) on different occasions. Venous blood samples for measurement of the blood haemoglobin (Hb), haematocrit and glucose concentrations were collected on 10 occasions over 2 h. A kinetic model based on haemoglobin dilution and urinary excretion was used to estimate the rate of absorption, the blood volume expansion over time, and the rate of elimination. Obtained kinetic data was used to simulate combinations of the three beverages in order to reach a predetermined goal of a 1:1 hydration of the blood volume and peripheral tissues over 6 h.

Results

Tap water had the fastest absorption but primarily hydrated peripheral tissues. Maximum hydration was reached after 17 min. Lemonade effectively expanded the blood volume and was absorbed and excreted at a high rate. The maximum hydration from isotonic saline occurred 60 min after ingestion. Slow excretion could make it possible to use saline to prolong the effects of the other two beverages.

Conclusions

It is possible to use the kinetic model to evaluate fluid turnover and compartmental distribution. Composition and timing of fluid intake can be calculated mathematically to meet predetermined goals of hydration and distribution.

Trial registration

NCT01360333 Date of registration: 05/23/2011.

Peer Review reports

Background

Orally ingested fluids are absorbed and excreted differently, depending on the ingested volume, tonicity and composition of the fluid. When a certain effect is desired, the composition of the fluid has to be adjusted accordingly [1]. These effects may include treatment of dehydration, supportive treatment and hydration prior to sport activity. An important aspect of the choice of fluid is time dependency, that is, the fluid’s effect over time, which is governed by the absorption and compartmental distribution of the orally ingested fluid. Fluid kinetics of beverages has earlier been studied using single point samples, ultrasound and gastric tubes to measure gastric emptying, body weight or urine-specific gravity. However, these techniques provide only snapshots or simplified views of what happens with the fluid in the body from ingestion to excretion.

Volume kinetics is related to pharmacokinetics and has been developed to capture the effects of intravenous fluid therapy, but there is reason to believe that this methodology can be used to study oral ingestion of fluid as well. Instead of using indirect measures, such as the emptying of the ventricle, volume kinetics estimates the effects of ingested fluid on the blood volume and body fluids. This is of interest to sports medicine because the blood volume is a major determinant of cardiac output [28].

If the kinetics of a beverage is known, the structural parameters could possibly be used to tailor a desired effect by combining different fluids. In sports, for example such tailoring could be useful [1, 7].“Timing” of the intake is, in that case, essential. The beverage that offers optimal performance effects would contain glucose and/or salt in various combinations in order to support metabolic needs and replace fluid and salt losses by sweating. Further, the addition of sugar or salt reduces the volume of distribution of the fluid, which increases blood volume expansion.

The first aim of this study was to determine whether it is feasible to measure different effects of fluids with the volume kinetic method. The second aim was to create mathematical models that simulate what fluids and compositions are needed to achieve desired effects. The input for the analysis was serial measurements of haemodilution, which is the inverse of intravascular water volume [9, 10].

Methods

The study was approved by the regional ethics committee in Linköping, Sweden (Dnr 2010/241-31) and registered at ClincalTrials.com with identifier NCT01360333. Data was collected during 2011. Ten healthy volunteers, seven males and three females, aged 33 ± 10 years (mean ± standard deviation [SD]; range 21–48 years) and with a body weight of 77 ± 13 kg (range 58–95 kg), were recruited for the three experiments. After being informed about the study both orally and in writing, each volunteer gave his/her consent for participation.

The participants arrived at the Department of Intensive Care at Linköping University Hospital between 7:30 and 8:30 a.m. The volunteers had fasted since midnight, but to ensure that they were euhydrated prior to the study, they were told to eat one sandwich and drink one glass (2 dL) of clear liquid shortly prior to 6:00 a.m. All participants voided prior to beginning of the study and were then randomised to ingest one of three beverages. To achieve a haemodynamic steady state, volunteers rested for 30 min before the experiments began. During this time, a cannula was placed in the cubital vein of one arm to sample blood. In every case, the first blood sample was collected more than 2 h after any previous ingestion of liquid. The volunteers rested on a bed below an OPN Thermal Ceiling radiant warmer (Aragon Medical, River Vale, NJ) placed about 1 m above them, and the heat was adjusted to achieve optimal comfort.

Ingested fluids

Each volunteer underwent the following three experiments in random order, separated by at least seven days:

  1. A.

    Tap water, 0.5 L (orally).

  2. B.

    Carbohydrate-containing soft drink (90 g carbohydrates/L), 0.5 L (orally).

  3. C.

    Saline, NaCl (9 g/L), 0.5 L (orally).

The tap water contained an average of 0.009 g/L of sodium, 0.002 g/L of potassium and 0.019 g/L of calcium.

Measurements

Venous blood (1–3 mL per tube) was withdrawn from the venous cannula using a vacuum tube just before and at 5, 10, 15, 20, 30, 45, 60, 90 and 120 min after the ingestion of the beverage. The baseline sample was drawn in duplicate, and the mean was used in all calculations. Before each sampling, a small volume of blood was drawn, and the volume was replaced by 2 mL of 0.9 % saline to prevent clotting.

All blood samples were analysed in the hospital’s central laboratory. The blood haemoglobin (Hb) concentration and the haematocrit (Hct) were analysed in EDTA tubes on a Cell-Dyn Sapphire (Abbott Diagnostics, Abbott Park, IL). The coefficient of variation (CV) for these analyses, based on the duplicate baseline samples, was 0.8 %. Plasma was used for measurement of the sodium, potassium (litium-heparine plasma gel tube) and glucose concentrations (Na-fluorid/oxalat-tube) on an ADVIA® 1800 Chemistry System (Siemens, Eschborn, Germany) with a CV of 1.0, 1.5 and 1.5 %, respectively. At the end of the study, all participants voided and the volume was measured.

Calculations

The following volume kinetic calculations are derived from pharmacokinetics.

The absorption, volume of distribution and elimination of fluid were analysed by a one-volume kinetic model.

Fluid was absorbed at a rate determined by a rate constant ka to increase the volume of central body fluid space V to v. The rate of elimination was given as the product of the dilution V and the clearance CL. The differential equation is as follows:

$$ \mathrm{d}v/dt={k}_aD{e}^{-{k}_at}{\textstyle\ \hbox{-}\ }CL\left(v{\textstyle\ \hbox{-}\ }V\right)/V $$

where D is the volume of ingested drink. The Hb-derived fractional plasma dilution was used to indicate the volume expansion V resulting from the infusion [11]:

$$ \left(v{\textstyle\ \hbox{-}\ }V\right)/V=\left(\left[\mathrm{H}\mathrm{b}/\mathrm{h}\mathrm{b}\right]\ {\textstyle \hbox{-}\ }1\right)/\left(1{\textstyle\ \hbox{-}\ }\mathrm{H}\mathrm{c}\mathrm{t}\right) $$

Symbols in capital letters denote baseline values. The minor artificial dilution due to the blood sampling was corrected mathematically [11]. Elimination of fluid by evaporation was disregarded.

The structural parameters in the model (k a , V and CL) were estimated using population modelling in Phoenix software for non-linear mixed effects, NLME version 1.3 (Pharsight, St Louis, MO). The possibility to stabilise the analyses by using the renal clearance CLR as CL was used:

$$ CL=\frac{{\displaystyle \sum}\mathrm{urine}\ \mathrm{volume}\ }{AUC\kern0.5em for\ \left(v-V\right)/V} $$

where AUC = area under the curve. The half-life was obtained as ln 2 * V / Cl.

Simulations were performed with programs developed by us for the Matlab R2012b software (Math Works Inc., Natick, MA).

The plasma volume expansion was obtained as the product of the simulated plasma dilution, that is, ((v–V)/V) and the estimated plasma volume at baseline. The latter was obtained as (1–haematocrit) x 7 % of the body mass.

The amount of absorbed fluid that remained in the body at any time t was calculated as the product of the plasma dilution and V.

Statistics

Data are given as the mean, or best estimate, and SD. Differences were evaluated by the Wilcoxon matched-pair test.

Results

Haemoglobin dilution

All 10 volunteers ingested all three fluids. Tap water caused maximal haemodilution 20 min after ingestion (P < 0.05). Both saline and the sugar drink were followed by more pronounced haemodilution (P < 0.01 at 30 min for both), and saline continued to cause haemodilution up to the study’s conclusion at 120 min (P < 0.01) (Table 1 and Fig. 1).

Table 1 Series of blood haemoglobin (B-Hb) concentrations, baseline haematocrit (Hct) and baseline concentrations of serum sodium (S-Na), plasma glucose (P-glucose) and serum potassium (S-K). All values are presented as mean ± SD
Fig. 1
figure 1

Blood haemoglobin (B-Hb) concentration over the course of the study. Crude data before recalculation to plasma dilution. Each point is the mean of 10 measurements

Fluid balance

Results of the kinetic analysis are shown in Table 2. The differences between the ingested water volumes with regard to the rates of uptake and elimination are illustrated in Figs. 2 and 3. Tap water was absorbed at the highest rate but caused the least pronounced expansion of the blood volume. Most of the absorbed sugar drink volume remained in the blood, and minimal amounts were distributed to the tissues. The blood volume expansion yielded by saline was equally as pronounced as that yielded by the sugar drink, but it was more prolonged.

Table 2 Kinetic constants derived from kinetic analysis of the water volume content of three beverages taken by mouth on different occasions by 10 volunteers
Fig. 2
figure 2

Dilution of venous plasma after ingestion of three different beverages (Tap water (a), Sugar drink (b) and Normal saline (c)). Each point is one estimate of plasma dilution, and the solid line is the modelled average curve

Fig. 3
figure 3

Modelled blood volume expansion (a) and volume of infused fluid remaining in the body (b) after ingestion of three different beverages

Excretion of the ingested volume was fastest for tap water, intermediate for the sugar drink and slowest for saline. Their average half-lives were 13, 16 and 43 min, respectively.

Urinary excretion during the study period differed considerably between the three fluids, the most substantial being after ingestion of the sugar drink (567 ± 207 mL) and tap water (505 ± 194 mL). Both were significantly higher (P < 0.05) than excretion after the ingestion of saline, after which the urine volume amounted to only 266 ± 144 mL.

Glucose and electrolytes

The baseline concentrations of the serum sodium and potassium and the plasma glucose are shown at the bottom of Table 1.

The serum sodium concentration decreased after ingestion of both tap water (0 min vs 50 min; P < 0.01) and the sugar drink (P < 0.01; Fig. 4a).

Fig. 4
figure 4

Serum sodium (a), serum potassium (b) and plasma glucose (c) concentrations after ingestion of three different beverages. Each point is the mean of 10 measurements

Serum potassium increased slightly after tap water ingestion, but decreased after ingestion of the sugar drink (0 min vs 30 min; P < 0.05). Saline ingestion increased serum potassium from 3.68 ± 0.23 at baseline to 3.96 ± 0.31 mmol/L after 1 h (P < 0.05; Fig. 4b).

Blood glucose rose quickly from 5.3 ± 0.5 to a maximum of 8.3 ± 1.3 mmol/L at 30 min after ingestion of lemonade (P < 0.01; Fig. 4c). Blood glucose did not change in response to ingestion of the other two beverages.

Simulations

Computer simulations were performed with the aim of finding a pattern of fluid intake and a mixture of the three beverages that maintained a stable and marked blood volume expansion of up to 420 min (and also to hydrated, but not excessively, the peripheral tissues). A proportion of 1:1 between the blood and peripheral tissues was considered adequate.

To quickly hydrate both body fluid volumes, a combination of tap water and the sugar drink was found to be most effective. Later fluid intake should contain less volume but relatively more saline to prolong blood volume expansion (Fig. 5).

Fig. 5
figure 5

Blood volume expansion (a) and volume of infused fluid remaining in the body (b) when three beverages are combined to create a stable expansion of the blood and peripheral tissues in the approximate proportion 1:1. A low total (x1) and multiples thereof are shown. Computer simulation based on the kinetic constants from Table 2. Three intakes or doses (0, 60 and 240 min) of different fluid combinations and volumes are simulated up to 420 min. The first dose requires a different composition than the other two

Discussion

The kinetic analysis yielded the rates of uptake and elimination as well as the volume of distribution for the three liquids. The kinetic profiles of the tested fluids proved to be quite different, which is of interest for the timing and dosing of beverages [12]. Tap water was quickly absorbed but primarily hydrated the peripheral tissues, whereas the glucose drink more effectively expanded the blood volume. The downside of both of these fluids is that they were excreted quickly. Saline can be used to prolong the volume effect, as saline distributes to both the blood and to peripheral tissues, but is excreted more slowly.

Previous research contains some contradiction between the rate of gastric emptying and effect on the blood volume [13, 14]. Ingesting carbohydrate-rich fluid before a surgical procedure, in contrast to fasting, improves postoperative bowel function [15]. In the present study, the subjects were told to drink 2 dL of clear fluid and eat one sandwich 2 h prior to the beginning of the study. This may explain why the gastric emptying and uptake of carbohydrate fluid occurred quickly. Because the bowel was stimulated with fluid and food, it was better prepared to handle the carbohydrate fluid ingested during the study. Further, the subjects were at rest, thus optimising bowel function.

The computer simulation in Fig. 5 suggests how beverages may be combined to reach predefined goals. We assumed that increasing blood volume, which typically places an athlete higher on the Frank-Starling curve, would be such a goal. Pre-exercise hyperhydration [7, 16] reduces heart rate during exercise, indicating larger preload as well as stroke volume. Hydration of peripheral tissues for oxygenation and temperature regulation is also needed [7, 17, 18], but the focus would still fall primarily on expanding blood volume. In such simulations one may further consider losses of salt, water and glucose to make the situation more realistic with regard to various exercises, but such adaptations were not performed in this study.

The simulations highlight that a substantial amount of fluid must be ingested to expand the blood volume. To increase the blood volume by 400 mL within 30 min, one could ingest 900 mL of mostly lemonade and tap water with the addition of a small amount of saline, a combination that results in a glucose concentration of 4 %. The amount of further intake required to maintain a steady state is much lower. However, because our data are valid only for volunteers at rest, Fig. 5 should not be applied directly to athletes performing strenuous exercise.

The concepts used in the calculations were derived from pharmacokinetics, but their application to fluid volumes has been refined in anaesthesia research [19]. Haemodilution is calculated by a series of measurements of Hb that, together with the measured urine volume, is entered into a set of differential equations that include the unknown parameters in a kinetic model. In the present study, the three unknown parameters represented the rate of fluid uptake, the volume of distribution at baseline and the rate of elimination. The calculations of water volume distribution account only for the volume absorbed by the mechanism ka. Therefore, the fluid volume that remains in the gut at any time t is 500 ml minus the absorption function, which is the first part of the differential equation used to estimate the kinetic parameters. However, our calculations assumed complete absorption of the ingested water at the end of the 2-h period of the study.

The beverages caused changes in the blood concentration of some small ions. As expected, serum sodium decreased after ingestion of the two sodium-free fluids, and blood glucose increased markedly after intake of lemonade. The decrease in serum potassium following ingestion of lemonade is probably a co-phenomenon with the uptake of glucose by the body cells. Even though saline contains no potassium, it raised the serum potassium concentration. The reason for this effect is that isotonic saline contains 50 % more chloride ions than the extracellular fluid, which causes acidosis, which in turn promotes potassium translocation from the cells. The slow excretion of saline as compared to tap water and lemonade can be explained by a reduction of the glomerular filtration rate by 10 %–15 % [20] and the fact that the kidneys must concentrate the urine to excrete the sodium load [21]. The prolonged hydration capacity and content of salt makes saline a useful beverage in warm environments [1, 18].

The limitations of the present study include the low level of inflicted haemodilution, which increases the impact of sampling technique and measurement errors on the precision of the kinetic parameters. Therefore, a commonly used technique of pooling all Hb values for each fluid (the naïve approach) was used for the kinetic evaluations. The presented kinetic parameters do not capture between-subject variability but provide a representative picture of the overall effects of the beverages. To provide suggestions about the use of the three beverages during physical exercise, data should be collected in that setting.

Conclusions

It is possible to use the volume kinetic model to evaluate the fluid turnover and compartmental distribution of oral liquids on a group basis. The composition and timing of fluid intake can be calculated mathematically to meet predetermined goals of hydration and water volume distribution. If a rapid effect of short duration is desired, water or a carbohydrate-rich beverage should be chosen. If a more long-lasting effect is desired, beverages containing sodium should be used [4, 22, 23]. Further studies will focus on optimising the beverages with respect to exercise tasks.

References

  1. Baker LB, Jeukendrup AE. Optimal composition of fluid-replacement beverages (review article). Compr Physiol. 2014;4:575–620.

    Article  PubMed  Google Scholar 

  2. Convertino VA. Blood volume response to physical activity and inactivity. Am J Med Sci. 2007;334:72–9.

    Article  PubMed  Google Scholar 

  3. Coyle EF, Hemmert MK, Coggan AR. Effects of detraining on cardiovascular responses to exercise: Role of blood volume. Am J Physiol. 1986;60:96–9.

    Google Scholar 

  4. Greenleaf JE, Looft-Wilson R, Wisherd JL, McKenzie MA, Jensen CD, Whittam JH. Pre-exercise hypervolemia and cycle ergometer endurance in men. Biol Sport. 1997;14:103–14.

    CAS  PubMed  Google Scholar 

  5. Hagberg JM, Goldberg AP, Lakatta L, O’Connor FC, Becker LC, Lakatta EG, Fleg JL. Expanded blood volumes contribute to the increased performance of endurance-trained older men. J Appl Physiol. 1998;85:484–9.

    CAS  PubMed  Google Scholar 

  6. Horwitz LD, Atkins JM, Leshin SJ. Role of the Frank-Starling mechanism in exercise. Circ Res. 1972;31:868–75.

    Article  CAS  PubMed  Google Scholar 

  7. Gigou PY, Dion T, Asselin A, Berrigan F, Goulet ED. Pre-exercise hyperhydration-induced bodyweight gain does not alter prolonged treadmill running time-trial performance in warm ambient conditions. Nutrients. 2012;4:949–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Starling EH. The Linacre lecture on the law of the heart. London: Longmans, Green and Co.; 1918.

    Google Scholar 

  9. Jimenez C, Melin B, Koulmann N, Allevard AM, Launay JC, Savourey G. Plasma volume changes during and after acute variations of body hydration level in humans. Eur J Appl Physiol Occup Physiol. 1999;80:1–8.

    Article  CAS  PubMed  Google Scholar 

  10. Noakes TD, Wilson G, Gray DA, Lambert MI, Dennis SC. Peak rates of diuresis in healthy humans during oral fluid overload. S Afr Med J. 2001;91:852–7.

    CAS  PubMed  Google Scholar 

  11. Ewaldsson CA, Hahn RG. Kinetics and extravascular retention of acetated Ringer’s solution during isoflurane and propofol anesthesia for thyroid surgery. Anesthesiology. 2005;103:460–9.

    Article  CAS  PubMed  Google Scholar 

  12. Goulet EDB. Dehydration and endurance performance in competitive athletes. Nutr Rev. 2012;70 Suppl 2:S132–6.

    Article  PubMed  Google Scholar 

  13. Evans GH, Shirreffs SM, Maughan RJ. Acute effects of ingesting glucose solution on blood and plasma volume. Br J Nutr. 2009;101:1503–8.

    Article  CAS  PubMed  Google Scholar 

  14. Evans GH, Shirreffs SM, Maughan RJ. The effects of repeated ingestion of high and low glucose-electrolyte solutions on gastric emptying and blood 2H2O concentration after an overnight fast. Br J Nutr. 2011;106:1732–9.

    Article  CAS  PubMed  Google Scholar 

  15. Luttikhold J, Oosting A, van den Braak CC, van Norren K, Rijna H, van Leeuwen PA, Bouritius H. Preservation of the gut by preoperative carbohydrate loading improves postoperative food intake. Clin Nutr. 2013;32:556–61.

    Article  CAS  PubMed  Google Scholar 

  16. Goulet ED, Rousseau SF, Lamboley CR, Plante GE, Dionne IJ. Pre-exercise hyperhydration delays dehydration and improves endurance capacity during 2 h of cycling in a temperate climate. J Physiol Anthropol. 2008;27:263–71.

    Article  PubMed  Google Scholar 

  17. Sims ST. Rehrer, NJ, Bell ML, Cotter JD: Preexercise sodium loading aids fluid balance and endurance for women exercising in the heat. J Appl Physiol. 2007;103:534–41.

    Article  CAS  PubMed  Google Scholar 

  18. Sims ST, van Vliet L, Cotter JD, Rehrer NJ. Sodium loading aids fluid balance and reduces physiological strain of trained men exercising in the heat. Med Sci Sports Exerc. 2007;39:123–30.

    Article  CAS  PubMed  Google Scholar 

  19. Hahn RG. Volume kinetics of infusion fluids (review). Anesthesiology. 2010;113:470–81.

    Article  PubMed  Google Scholar 

  20. Chowdbury AH, Cox EF, Francis ST, Lobo DN. A randomized, controlled, double-blind crossover study on the effects of 2-L infusions of 0.9% saline and Plasma-Lyte 148 on renal blood flow velocity and renal cortical tissue perfusion in healthy volunteers. Ann Surg. 2012;256:18–24.

    Article  Google Scholar 

  21. Hahn RG, Drobin D. Rapid water and slow sodium excretion of Ringer’s solution dehydrates cells. Anesth Analg. 2003;97:1590–4.

    Article  CAS  PubMed  Google Scholar 

  22. Greenleaf JE, Looft-Wilson R, Wisherd JL, Jackson CG, Fung PP, Ertl AC, Whittam JH. Hypervolemia in men from fluid ingestion at rest and during exercise. Aviat Space Environ Med. 1998;69:374–86.

    CAS  PubMed  Google Scholar 

  23. Frey MA, Riddle J, Charles JB, Bungo MW. Blood and urine responses to ingesting fluids of various salt and glucose concentrations. J Clin Pharmacol. 1991;31:880–7.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Assistance was given by nurse anaesthetist Susanne Lind.

Funding

This work was supported by a grant from Östergötland County Council.

Availability of data and materials

All relevant data are contained within the manuscript and its additional files.

Authors’ contributions

JZ contributed to conception, design and interpretation of the study. AM took part in design of the study and the acquisition of data. RH contributed to design of the study as well as interpretation of the data. All three authors were involved in the drafting and revising of the manuscript and all three authors have given their approval to the final version of the article.

Authors’ information

Joachim Zdolsek MD, Assoc prof, works as a consultant at the Department of Anaesthesiology and Intensive care at Vrinnevi Hospital in Norrköping, Sweden.

Annika Metander MD works as a consultant at the Department of Thoracic Anaesthesiology at Linköping University Hospital, Sweden.

Robert Hahn MD, Prof, works as Research Director at Södertälje Hospital, Sweden.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The study was approved by the regional ethics committee in Linköping, Sweden (Dnr 2010/241-31) and registered at ClincalTrials.com with identifier NCT01360333. All participants gave their written informed consent to participate in the study.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Joachim Zdolsek.

Additional file

Additional file 1:

Complete laboratory results. (XLS 48 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zdolsek, J., Metander, A. & Hahn, R. Volume kinetic evaluation of fluid turnover after oral intake of tap water, lemonade and saline in volunteers. BMC Sports Sci Med Rehabil 8, 22 (2016). https://doi.org/10.1186/s13102-016-0045-x

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13102-016-0045-x

Keywords