The effect of changing ventilator settings on indices of ventilation inhomogeneity in small ventilated lungs
© Schmalisch et al; licensee BioMed Central Ltd. 2006
Received: 21 March 2006
Accepted: 18 August 2006
Published: 18 August 2006
In ventilated newborns the use of multiple breath washout (MBW) techniques for measuring both lung volume and ventilation inhomogeneity (VI) is hampered by the comparatively high dead space fraction. We studied how changes in ventilator settings affected VI indices in this particular population.
Using a computer simulation of a uniformly ventilated volume the interaction between VI indices (lung clearance index (LCI), moment ratios (M1/M0, M2/M0, AMDN1, AMDN2) of the washout curve) and tidal volume (VT), dead space (VD) and functional residual capacity (FRC) were calculated. The theoretical results were compared with measurements in 15 ventilated piglets (age <12 h, median weight 1135 g) by increasing the peak inspiratory pressure (PIP). FRC and VI indices were measured by MBW using 0.8% heptafluoropropane as tracer gas.
The computer simulation showed that the sensitivity of most VI indices to changes in VD/VT and VT/FRC increase, in particular for VD/VT > 0.5. In piglets, the raised PIP caused a significant increase of VT from 15.4 ± 9.5 to 21.9 ± 14.7 (p = 0.003) and of the FRC from 31.6 ± 14.7 mL to 35.0 ± 15.9 mL (p = 0.006), whereas LCI (9.15 ± 0.75 to 8.55 ± 0.74, p = 0.019) and the moment ratios M1/M0, M2/M0 (p < 0.02) decreased significantly. No significant changes were seen in AMDN1 and AMDN2. The within-subject variability of the VI indices (coefficient of variation in brackets) was distinctly higher (LCI (9.8%), M1/M0 (6.6%), M2/M0 (14.6%), AMDN1 (9.1%), AMDN2 (16.3%)) compared to FRC measurements (5.6%). Computer simulations showed that significant changes in VI indices were exclusively caused by changes in VT and FRC and not by an improvement of the homogeneity of alveolar ventilation.
In small ventilated lungs with a high dead space fraction, indices of VI may be misinterpreted if the changes in ventilator settings are not considered. Computer simulations can help to prevent this misinterpretation.
In ventilated newborns respiratory problems are often caused by impaired lung development and uneven alveolar ventilation. Surfactant deficiency or dysfunction may increase the ventilatory inhomogeneity by collapse or over distention of the alveoli. Thus, there is an increasing clinical interest in multiple breath washout techniques (MBW) to measure both the functional residual capacity (FRC) and ventilatory inhomogeneity (VI) indices [1–4]. Lung clearance index (LCI) and moment ratios of the wash out curve are the most frequently used indices in infancy . Commonly, the first and the second moment related to the zeroth moment (M1/M0, M2/M0) are calculated.
Most VI indices are easily calculated, however, a key disadvantage is their dependency on the breathing pattern . In 1975, Saidel et al.  suggested that this dependency can be reduced by performing a moment analysis of the washout curve plotted as a function of the cumulative exhaled volume related to the FRC. However, the dependency of the moment ratios on the ventilatory dead space (VD) remained. Therefore, Habib and Lutchen  replaced the cumulative exhaled volume by the cumulative alveolar volume to reduce the influence of VD. They referred to the first two moment ratios of the wash out curve as alveolar-based mean dilution numbers AMDN1 and AMDN2.
Dead space fractions (VD related to the tidal volume VT) in adults commonly lie between 0.05 to 0.2 . In ventilated newborns the dead space fraction VD/VT is often markedly higher  depending on the ventilator settings. Typical values lie between 0.4 and 0.6  and in preterm or surfactant-depleted lungs VD/VT can rise up to 0.7 . Such small lungs are ventilated with a relative low tidal volume to prevent volu-trauma. Mainstream flow sensors and gas analyzers considerably increase the apparatus dead space so that high dead space fractions are not uncommon. The effect of an increased VD/VT on the sensitivity of VI indices to parameter changes is not well known. Therefore, the aim of this study was to investigate how changing ventilator settings affect the different VI indices in this particular population by mathematical modelling and by measurements in newborn piglets using the MBW technique with heptafluoropropane (HFP) as tracer gas.
In patients ventilated with a constant VT the LCI is given by
where NLCI is the number of breaths required to lower the end tidal tracer gas concentration to 1/40th of the starting concentration . The ideal washout curve of an inert gas from a uniformly ventilated volume can be expressed as
where cn is the end-expiratory gas concentration of the nth breathing cycle and c0 is the initial gas concentration. Using computer simulations of the washout curve LCI can be calculated as a function of VT/FRC and VD/VT.
The moments M0, M1, M2 of the washout curve were calculated up to NLCI. For a constant VT the moment ratios are given by
where "i" is the number of the breathing cycle. For an infinite number of cycles M1/M0 has a fixed limit given by
Following Habib and Lutchen , the replacement of VT by the alveolar ventilation VT-VD yield the alveolar-based mean dilution numbers
For a well-mixed volume and an infinite number of breathing cycles, AMDN1 is equal one regardless of VD/VT or VT/FRC
and an AMDN1 > 1 implies inhomogeneity purely at alveolar level . For very low dead spaces (VD/VT≈0) the alveolar-based mean dilution numbers are equal the moment ratios (AMDN1 = M1/M0, AMDN2 = M2/M0). An important feature of all VI indices is that they rise with increasing inhomogeneity of the alveolar ventilation which can be shown easily by computer simulations of multi-compartment models.
Fifteen newborn piglets (age <12 h, median weight 1135 g) placed in supine position within a heated incubator were anesthetized (azaperon 8 mg·kg-1 and ketamin 10 mg·kg-1), intubated (shortened neonatal endotracheal tube (ETT) with 3.5 mm outer diameter, Vygon, Ecouen, France), paralyzed (pancuronium-bromide 0.2 mg·kg-1·hour-1) and mechanically ventilated with a neonatal ventilator (Babylog 8000, Draeger, Lübeck, Germany). During the study period, air flow (6 L/min), respiratory rate of 40/min and fraction of inspired oxygen (FiO2) of 1.0 were kept constant. Positive end-expiratory pressure (PEEP) was set to zero, peak inflation pressure (PIP) was set initially to 8.3 ± 3.1 cmH2O and elevated to 12.1 ± 5.0 cmH2O. Ventilatory parameters were taken from the Babylog 8000 and recorded continuously. Lung volume and VI indices were measured by heptafluoropropane (HFP) wash in and wash out as previously described . Briefly, a new infrared HFP sensor was sited between the flow sensor of the Babylog 8000 and the ETT. The total apparatus dead space of HFP sensor, flow sensor and ETT was 4.5 mL determined by water displacement. The constant flow of the ventilator was 8 L/min in all measurements. Using a mechanical valve to start wash in or wash out, a HFP flow from a gas cylinder (medical grade HFP, Solvay, Hannover, Germany) was fed into the inspiratory limb of the ventilatory circuit to achieve a constant HFP concentration of 0.8%. The flow signal of the Babylog 8000 and the concentration signal from the HFP sensor were used to calculate FRC and the VI indices from the wash-in or wash-out curve up to 1/40th of the starting concentration by an external computer. The Fowler dead space VD was determined from the first 5 cycles. The calculation was stopped automatically after N cycles when the total amount of alveolar turnovers exceed the tenfold of the calculated FRC (minimum number 40 cycles).
After instrumentation and onset of mechanical ventilation a stabilisation period of 15 minutes was allowed before the measurements were started with a HFP wash-in procedure (FRCwash-in) and a consecutive wash-out procedure (FRCwash-out). Such a cycle was accepted for evaluation if the deviation between FRCwash-in and FRCwash-out was lower than 20% and the VT was higher than 4.5 mL (VDapp). In order to investigate the effect of ventilator settings on VI indices PIP was increased by 4 cm H2O. After a stabilisation period of 15 minutes the MBW was repeated in the same manner.
A computer program written in Visual Basic (Microsoft Corpor., USA) as a macro of a EXCEL worksheet (Microsoft Office 2000, Microsoft Corpor., USA) was developed to compare the VI indices measured in the piglets with the VI indices of a uniformly ventilated volume using the same ventilator settings [see Additional file 1]. The program calculates the washout curve according equation 2 and the corresponding VI indices according equations 1, 3 and 5.
For each measurement in the piglets at least 5 wash-in and wash-out cycles were performed and averaged. Data are presented as mean ± SD and mean individual differences with 95% CI as appropriate. Differences in the animals were compared by the paired t-test. To assess the within-subject variability of repeated measurements the coefficient of variation (CV) was calculated for all parameters and compared by a rank test. A level of statistical significance of p < 0.05 was accepted.
Effect of increased peak inspiratory pressure on ventilation, FRC and ventilatory inhomogeneity indices in 15 ventilated piglets (Presented are groups means ± SD and the mean individual changes with 95%CI compared to the baseline)
Increased peak inspiratory pressure
Mean individual difference in % with 95%CI
V T (mL)
15.4 ± 9.5
21.9 ± 14.7
38.9 (20.4 to 57.45)
p = 0.003
31.6 ± 14.7
35.0 ± 15.9
11.1 (4.5 to 17.8)
p = 0.006
V T /FRC
0.48 ± 0.17
0.59 ± 0.20
25.0 (10.7 to 39.4)
p = 0.003
V D /V T
0.49 ± 0.15
0.43 ± 0.13
-10.2 (-21.0 to 0.5)
p = 0.006
9.17 ± 0.75
8.55 ± 0.74
-6.4 (-11.8 to -1.0)
p = 0.019
M 1 /M 0
2.03 ± 0.16
1.84 ± 0.16
-9.4 (-14.0 to -4.9)
p = 0.006
M 2 M 0
8.41 ± 1.60
7.22 ± 4.80
-12.1 (-23.1 to -1.4)
p = 0.017
1.034 ± 0.30
1.028 ± 0.19
-0.6 (-9.1 to 7.9)
p = 0.886
2.88 ± 1.56
3.01 ± 1.06
4.5 (-9.9 to 18.1)
p = 0.515
A significant decrease was also seen in the LCI (p = 0.019) and the moment ratios M1/M0 (p = 0.006) and M2/M0 (p = 0.017). No significant changes were seen in AMDN1 and AMDN2.
There was a strong correlation between LCI and the moment ratios (M1/M0, M2/M0) with r = 0.885 and r = 0.907, respectively, and independent of which PIP was used. Thus, it was not surprising that there was a similar effect of the increased PIP on LCI and the moment ratios (M1/M0, M2/M0). No statistically significant correlations were found between the alveolar-based mean dilution numbers (AMDN1, AMDN2) and LCI, M1/M0 and M2/M0.
The within-subject variability of the measured parameters showed considerable variations but it was not affected by the increase of the PIP. The median CV of all FRC measurements was 5.6%. The median CV of the LCI was significantly greater (9.8%, p = 0.0004). Compared with the LCI the CV of M1/M0 was significantly smaller (6.6%, p = 0.003), whereas the CV of M2/M0 was distinctly greater (14.6%, p = 0.004). The CVs of AMDN1 (9.1%) and AMDN2 (16.3%) were always greater than the CVs of the other moment ratios.
Modelling of ventilatory inhomogeneity indices of a uniformly ventilated volume using the measured VD/VT and VT/FRC of table 4
Difference in %
V D /V T
V T /FRC
M 1 /M 0
M 2 /M 0
The measurement of lung volume and ventilation inhomogeneity by MBW is a fascinating, non-invasive technique. It is relatively easily performed, even in ventilated patients. In a previous study , we have shown that by this technique the effect of surfactant-depletion by lung lavage on the FRC and the VI indices is reliably measured: before and after lavage VD/VT was not significant different, therefore, the significant increase of the VI indices has to be predominantly attributed to the effect of lung lavage. In the present study the measurements were performed in healthy lungs and VD/VT was distinctly changed by an increase of the PIP. As shown in Fig. 5 the changes of the VI indices are mainly caused by physical laws of gas mixing.
The interpretation of significant changes in VI indices may be misleading if their dependency on the ventilator settings is not considered. This is a particular problem in small lungs where the relatively high dead space fraction increases the sensitivity of VI indices to parameter changes. Any changes in VD (e.g. by applying of a new mainstream sensor) or changes in VT and FRC (e.g. by changing of ventilator settings or by surfactant substitution) will affect the VI indices and therefore hamper their comparability.
As shown by the computer simulation, most VI indices increase with increasing VD/VT. This may explain why in newborns much higher VI indices values were measured [4, 13, 14] than in spontaneously breathing children [15, 16]. These higher values in newborns are more likely an expression of functional dependencies than the result of impaired alveolar ventilation. The relatively good agreement between the VI indices measured in healthy piglets (Table 1) and the calculated VI indices of a uniformly ventilated volume (Table 2) was surprising. There was only a small difference in the VI indices between the animal measurements and the modelling which can be attributed to the more complex ventilation distribution in the lungs of the piglets.
In infancy the LCI is one of most frequently used VI index [16–19] and easy to comprehend. It describes the number of turnovers to lower the end tidal tracer gas concentration to 1/40th of the starting concentration. Theoretically, the LCI is a static value of the flat tail of the washout curve and may vary if the signal is noisy. This explains its relatively high within-subject variability. The limitation on NLCI breathing cycle seems to be useful to reduce the influence of the signal noise on the measured LCI. Its main disadvantage is its high dependency on VT/VD as shown in Fig. 1.
Moment ratios are more abstract mathematical measures considering the whole washout curve. Only for M1/M0 a theoretical value for a well ventilated volume exists (equation 4). M1/M0 reflect more the first part of the washout curve, whereas M2/M0 better describe the tail of the curve. Therefore the within-subject variability of M2/M0 is distinctly higher compared to M1/M0 and similar to the within-subject variability of the LCI. A moment analysis makes higher demands on the wash out curve compared to LCI measurements. It requires a rapid rise of the tracer gas after the switch-on so that the full tracers gas concentration is reached during the first inspiratory cycle. This is sometimes difficult to achieve, in particular if the tracer gas is fed into the inspiratory limb of the ventilator circuit, far from the ETT, which may delay such a swift rise. Such a delay is tolerable for FRC measurements but will affect the calculation of the moment ratios.
In contrast to LCI and the ratios M1/M0 and M2/M0, no significant effect of the increased PIP on AMDN1 and AMDN2 was seen, as predicted by the computer simulation. These parameters suggested by Habbib and Lutchen  seem indeed to be less sensitive to the changes in the breathing pattern than the other ones. This dos not necessarily mean that they have a higher predictive value: with the exception of the above authors , a higher diagnostic value of these corrected moments could not be demonstrated until now . The main problem with these parameters is that they need an exact determination of the Fowler dead space. This is often difficult to evaluate in small lungs because the three phases of the gas concentration-volume diagram of the exhaled air are often not well defined . This may explain why in animals AMDN1 and AMDN2 often showed a very high within-subject variability (>20%). This high variability may limit their diagnostic value.
A central problem of all moment ratios is their dependency on the number of evaluated breathing cycles . The computer simulation has already shown that the theoretical values for M1/M0 and AMDN1 were not reached due to the finite number of evaluated cycles (Fig. 2). This hampers the comparability of the data between different laboratories if the start and the end of the evaluated breathing cycles are not specified.
The within-subject variability of LCI, M1/M0 and M2/M0 in our study was similar to those measured by Shao et al.  in preterm infants. In both studies the variability of the VI indices was distinctly higher compared with the CV of the FRC. Thus, in small ventilated lungs the determination of VI indices needs a higher number of wash-in and washout cycles than for FRC measurements to obtain reproducible results.
With the availability of dead space-minimized mainstream gas analyzers there is an increasing interest to measure ventilation inhomogeneity by MBW techniques. However, the use of VI indices in small lungs needs particular attention. Especially in small ventilated lungs with a relatively high dead space fraction most indices are significantly affected by ventilator settings. Changes in tidal volume and lung volume, or changes in the apparatus dead space hamper their comparison. Model simulations of a uniformly ventilated volume can help to decide if the changes in the VI indices are caused by changing ventilator settings or whether they indicate any changes in the homogeneity of alveolar ventilation.
- cn :
End-expiratory tracer gas concentration of the nth breathing cycle
Coefficient of variation
- FIO2 :
Fraction of inspired oxygen
Functional residual capacity
- FRCwash-in :
Functional residual capacity measured by wash-in of the tracer gas
- FRCwash-out :
Functional residual capacity measured by wash-out of the tracer gas
Intensive care unit (for newborn infants)
Lung clearance index
- M1/M0 :
First-to-zeroth moment ratio
- M2/M0 :
Second-to-zeroth moment ratio
Multiple breath washout
Number of breathing cycles
- NLCI :
Number of breaths required to lower the tidal tracer gas concentration to 1/40th of the starting concentration
Positive end-expiratory pressure
Peak inflation pressure
- VD :
- VDapp :
Apparatus deadspace volume
- VT :
The authors thank Ariane Kusztrich for her assistance in the animal experiments and Jessica Blank for her support in data analysis.
- Vilstrup CT, Bjorklund LJ, Larsson A, Lachmann B, Werner O: Functional residual capacity and ventilation homogeneity in mechanically ventilated small neonates. J Appl Physiol. 1992, 73: 276-283.PubMedGoogle Scholar
- Sandberg KL, Lindstrom DP, Sjoqvist BA, Parker RA, Cotton RB: Surfactant replacement therapy improves ventilation inhomogeneity in infants with respiratory distress syndrome. Pediatr Pulmonol. 1997, 24: 337-343. 10.1002/(SICI)1099-0496(199711)24:5<337::AID-PPUL6>3.0.CO;2-F.View ArticlePubMedGoogle Scholar
- Edberg KE, Sandberg K, Silberberg A, Ekstrom-Jodal B, Hjalmarson O: Lung volume, gas mixing, and mechanics of breathing in mechanically ventilated very low birth weight infants with idiopathic respiratory distress syndrome. Pediatr Res. 1991, 30: 496-500.View ArticlePubMedGoogle Scholar
- Schibler A, Henning R: Positive end-expiratory pressure and ventilation inhomogeneity in mechanically ventilated children. Pediatr Crit Care Med. 2002, 3: 124-128. 10.1097/00130478-200204000-00006.View ArticlePubMedGoogle Scholar
- Pillow JJ, Frerichs I, Stocks J: Lung function tests in neonates and infants with chronic lung disease: global and regional ventilation inhomogeneity. Pediatr Pulmonol. 2006, 41: 105-121. 10.1002/ppul.20319.View ArticlePubMedGoogle Scholar
- Larsson A, Jonmarker C, Werner O: Ventilation inhomogeneity during controlled ventilation. Which index should be used?. J Appl Physiol. 1988, 65: 2030-2039.PubMedGoogle Scholar
- Saidel GM, Salmon RB, Chester EH: Moment analysis of multibreath lung washout. J Appl Physiol. 1975, 38: 328-334.PubMedGoogle Scholar
- Habib RH, Lutchen KR: Moment analysis of a multibreath nitrogen washout based on an alveolar gas dilution number. Am Rev Respir Dis. 1991, 144: 513-519.View ArticlePubMedGoogle Scholar
- Claure N, D'Ugard C, Bancalari E: Elimination of ventilator dead space during synchronized ventilation in premature infants. J Pediatr. 2003, 143: 315-320. 10.1067/S0022-3476(03)00299-3.View ArticlePubMedGoogle Scholar
- Wenzel U, Wauer RR, Schmalisch G: Comparison of different methods for dead space measurements in ventilated newborns using CO2-volume plot. Intensive Care Med. 1999, 25: 705-713. 10.1007/s001340050933.View ArticlePubMedGoogle Scholar
- Proquitte H, Krause S, Rudiger M, Wauer RR, Schmalisch G: Current limitations of volumetric capnography in surfactant-depleted small lungs. Pediatr Crit Care Med. 2004, 5: 75-80. 10.1097/01.PCC.0000102384.60676.E5.View ArticlePubMedGoogle Scholar
- Edelman NH, Mittman C, Norris AH, Shock NW: Effects of respiratory pattern on age differences in ventilation uniformity. J Appl Physiol. 1968, 24: 49-53.PubMedGoogle Scholar
- Proquitte H, Kusztrich A, Auwärter V, Pragst F, Wauer RR, Schmalisch G: Functional residual capacity measurement by heptafluoropropane in ventilated newborn lungs: in-vitro and in-vivo validation. Crit Care Med. 2006, 34: 1789-95. 10.1097/01.CCM.0000220065.93507.AB.View ArticlePubMedGoogle Scholar
- Shao H, Sandberg K, Sjoqvist BA, Hjalmarson O: Moment analysis of multibreath nitrogen washout in healthy preterm infants. Pediatr Pulmonol. 1998, 25: 52-58. 10.1002/(SICI)1099-0496(199801)25:1<52::AID-PPUL6>3.0.CO;2-Q.View ArticlePubMedGoogle Scholar
- Schibler A, Hall GL, Businger F, Reinmann B, Wildhaber JH, Cernelc M, Frey U: Measurement of lung volume and ventilation distribution with an ultrasonic flow meter in healthy infants. Eur Respir J. 2002, 20: 912-918. 10.1183/09031936.02.00226002.View ArticlePubMedGoogle Scholar
- Aurora P, Gustafsson P, Bush A, Lindblad A, Oliver C, Wallis CE, Stocks J: Multiple breath inert gas washout as a measure of ventilation distribution in children with cystic fibrosis. Thorax. 2004, 59: 1068-1073. 10.1136/thx.2004.022590.View ArticlePubMedPubMed CentralGoogle Scholar
- Kraemer R, Blum A, Schibler A, Ammann RA, Gallati S: Ventilation inhomogeneities in relation to standard lung function in patients with cystic fibrosis. Am J Respir Crit Care Med. 2005, 171: 371-378. 10.1164/rccm.200407-948OC.View ArticlePubMedGoogle Scholar
- Gustafsson PM, Aurora P, Lindblad A: Evaluation of ventilation maldistribution as an early indicator of lung disease in children with cystic fibrosis. Eur Respir J. 2003, 22: 972-979.View ArticlePubMedGoogle Scholar
- Hjalmarson O, Sandberg KL: Lung function at term reflects severity of bronchopulmonary dysplasia. J Pediatr. 2005, 146: 86-90. 10.1016/j.jpeds.2004.08.044.View ArticlePubMedGoogle Scholar
- Schibler A, Henning R: Measurement of functional residual capacity in rabbits and children using an ultrasonic flow meter. Pediatr Res. 2001, 49: 581-588.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2466/6/20/prepub