Here, we describe two measures of pulmonary function – barometric plethysmography, which allows the measurement of lung volume, and volumetric capnography, a tool to measure the anatomic dead space and airways uniformity. These techniques may be used independently or combined to assess airways function at different lung volumes.
Cite this ArticleCopy Citation | Download Citations
Seymour, M., Pritchard, E., Sajjad, H., Tomasson, E. P., Blodgett, C. M., Winnike, H., Paun, O. V., Eberlein, M., Bates, M. L. Combining Volumetric Capnography And Barometric Plethysmography To Measure The Lung Structure-function Relationship. J. Vis. Exp. (143), e58238, doi:10.3791/58238 (2019).
Translate text to:
Tools to measure lung and airways volume are critical for pulmonary researchers interested in evaluating the impact of disease or novel therapies on the lung. Barometric plethysmography is a classic technique to evaluate the lung volume with a long history of clinical use. Volumetric capnography utilizes the profile of exhaled carbon dioxide to determine the volume of the conducting airways, or dead space, and provides an index of airways homogeneity. These techniques may be used independently, or in combination to evaluate the dependence of airways volume and homogeneity on lung volume. This paper provides detailed technical instructions to replicate these techniques and our representative data demonstrates that the airways volume and homogeneity are highly correlated to lung volume. We also provide a macro for the analysis of capnographic data, which can be modified or adapted to fit different experimental designs. The advantage of these measures is that their advantages and limitations are supported by decades of experimental data, and they can be made repeatedly in the same subject without expensive imaging equipment or technically advanced analysis algorithms. These methods may be particularly useful for investigators interested in perturbations that change both the functional residual capacity of the lung and airways volume.
Gas washout techniques have been used for decades to provide important information about the structure and uniformity of the airway tree. The lung is classically described as having two compartments – a conducting zone that is comprised of the anatomic dead space and the respiratory zone where gas exchange occurs in the alveoli. The conducting airways are termed as “dead space” because they do not participate in the exchange of oxygen and carbon dioxide. In the single breath gas washout method, the concentration profile of an exhaled gas can be used to determine the volume of the anatomic dead space and to derive information about the uniformity of ventilation. Some methods rely on the breathing of inert gases to make these measures (N2, argon, He, SF6, etc.). The use of inert gas is well-established, supported by scientific consensus statements1, and there are available commercial equipment with user friendly interfaces. However, the exhaled profile of carbon dioxide (CO2) can be used to derive similar information. Evaluating the profile of CO2 as a function of the exhaled volume, or volumetric capnography, does not require the participant to breathe special gas mixtures and allows the investigator to gather additional information flexibly about metabolism and gas exchange with minimal adjustment to the technique.
During a controlled exhalation, the concentration of CO2 can be plotted against the total exhaled volume. At the beginning of an exhalation, the dead space is filled with atmospheric gas. This is reflected in Phase I of the exhaled CO2 profile where there is an undetectable amount of CO2 (Figure 1, top). Phase II marks the transition to the alveolar gas, where gas exchange occurs and CO2 is abundant. The volume at the midpoint of Phase II is the volume of the anatomic dead space (VD). Phase III contains alveolar gas. Because airways with different diameters empty at different rates, the slope (S) of Phase III provides information about airways uniformity. A steeper slope of Phase III suggests a less uniform airway tree proximal to the terminal bronchioles, or convection-dependent inhomogeneity2. In the case where a perturbation may change the rate of CO2 production, and to make comparisons between individuals, the slope can be divided by the area under the curve to normalize for differences in metabolism (NS or normalized slope). Volumetric capnography has been used previously to evaluate the changes in airways volume and uniformity following air pollutant exposure3,4,5,6.
Gas transport in the lung is governed by both convection and diffusion. Single breath washout measures are highly dependent on air flow and the measured value of VD occurs at the convection-diffusion boundary. Changing the flow rate of the exhalation or preceding inhalation changes the location of that boundary7. Capnography is also highly dependent on the volume of the lung immediately preceding the maneuver. Larger lung volumes distend the airways, resulting in larger values of VD8. One solution is to consistently make the measurement at the same lung volume – usually functional residual capacity (FRC). An alternative, described here, is to couple volumetric capnography with barometric plethysmography, in order to obtain the relationship between VD and lung volume. The participant then performs the maneuver at constant flow rates, while varying the lung volume. This still allows for classic capnographic measures to be made at FRC, but also for the relationship between the lung volume and dead space volume and between the lung volume and homogeneity to be derived. Indeed, the added value of coupling capnography with plethysmography comes from the ability to test hypotheses about the distensibility of the airways tree and the structure-function relationship of the lung. This may be a valuable tool for investigators aiming to quantify the influence of airways mechanics versus lung compliance and elastance on pulmonary function in healthy and diseased populations9,10,11. Furthermore, accounting for the absolute lung volume at which the volumetric capnographic measurements are being performed allows investigators to characterize the effects of conditions that can alter the inflation state of the lung, such as obesity, lung transplant, or interventions like chest wall strapping. Volumetric capnography may ultimately have clinical utility in the intensive care setting12,13.
This protocol has been previously approved by and follows the guidelines set by the University of Iowa Institutional Review Board. Data shown were collected as part of a project approved by the Institutional Review Board at the University of Iowa. Participants gave informed consent and the studies were performed in accordance with the Declaration of Helsinki.
- Check the equipment table to verify that all required equipment is available. Double check the configuration using the graphic depiction of the equipment in Figure 2.
NOTE: Barometric plethysmography is a well-described clinical tool and is performed using commercial equipment according to the consensus statements on standardizing lung volume measurements14,15. When necessary, lung flows and volumes are compared to predicted values from the NHANES data set and Goldman and Becklake16 that are included in the plethysmograph software.
- Perform calibration of the plethysmograph daily and prior to any experiments.
- Measure the temperature, barometric pressure and relative humidity using a standard barometer prior to the calibration and enter these values into the plethysmograph software as correction factors.
- Calibrate the flow sensor using a calibrated 3 L syringe at variable flow rates. Calibrate the box pressure using a precise 50 mL pump. Box pressure transducers should be checked monthly and re-calibrated as needed, per manufacturer’s recommendation.
- Immediately prior to the measurement, place the participant in the whole body plethysmograph and close the door. Make measurements after 30-60 s, which allows for thermal equilibration.
- Instruct the participant to place their mouth on the mouthpiece, put on nose clips, and place their hands on their cheeks. Preventing “puffing” of the cheeks during the maneuver minimizes changes in volume that result from changing the mouth volume.
- Instruct the participant to breathe normally, allowing at least four tidal breaths to be acquired and functional residual capacity (FRC) to be established.
- At the end of a normal exhalation (FRC), close the shutter. Coach the participant to pant lightly at 0.5-1 breaths/s for 3-4 s. Evaluate the relationship between the mouth pressure and plethysmograph pressure to ensure that it is a series of overlapping, straight lines without thermal drift.
- Open the shutter and allow the participant to take a normal breath. Coach the participant to exhale to residual volume (RV), followed by a maximal inspiratory maneuver to total lung capacity. Repeat at least three times until FRC values that agree within 5% are obtained
3. Volumetric Capnography
NOTE: Steps 3.1 – 3.4 are performed before the arrival of the research subject.
- Before proceeding, address the variables in Table 1 and modify if needed. It is important that these variables are adjusted during the study design phase and then held constant for the duration of the study.
- Before beginning a new experimental protocol, take care to accurately measure the time delay between the gas analyzer, which measures CO2 concentration, and the pneumotach, which measures flow. This allows for the CO2 and flow signals to be aligned.
- Measure the time delay experimentally with a stream of 5% CO2. Attach the gas line to a stopcock, followed by the mouthpiece.
- Open the stopcock, introducing the gas at a rate of 10 L/min. Determine the mean time delay between the response of the pneumotach and gas analyzer over 10 trials and enter into the macro.
- Maintain the time delay constant by maintaining the analyzer sampling rate. The time delay is highly dependent on the sampling rate of the gas analyzer and it is critical that this remain constant through the experiment and between participants.
- Define three “channels” for the collection of flow, exhaled CO2 (%), and volume. Flow and exhaled CO2 (%) are analog inputs and volume is the integral of flow.
- Confirm that flow and CO2 (%) are measured directly from the pneumotach and gas analyzer and that volume is calculated as the integral of flow. Figure 3 shows that these are being collected in channels 1,2, and 6.
- Calibrate the gas analyzer prior to each use. Include the O2 sensor if this is to be measured.
- Zero the analyzer with an inert gas. 100% calibration grade (<0.01% contaminant) N2 or He may be used, although helium is preferred because nitrogen may be contaminated with trace amounts of oxygen. Place the drying tube in a bag or connect to a mixing chamber. Flush the bag or chamber with inert gas at a rate of at least 10 L/min. Care should be taken not to pressurize the system as this can impact the calibration.
- Flood the bag or chamber with inert gas to displace O2 and the CO2. Once the displayed concentrations of CO2 and O2 stabilize, adjust the zero knobs until they both read zero.
- Repeat with 6% CO2 and room air (20.93% O2) as calibration gases. When the concentration of the desired gas stabilizes, adjust the span knob to match the concentration of the calibration gas.
- Recheck the inert gas and calibration gases and adjust the zero and span until both are accurate ±0.1%.
- Calibrate the heated pneumotach according to the manufacturer’s instructions.
- Briefly, allow the pneumotach to warm to 37 °C for at least 20 min prior to the study.
- Select the drop-down menu of the flow channel (Channel 1), select the Spirometer menu option, and click Zero to zero the pneumotach. Finish by selecting Okay.
- Directly connect a 3L syringe to the pneumotach using a flow head adapter. Highlight the calibration breath. Again, select the drop-down menu of the flow channel. Select Spirometer flow | Calibrate, type in 3L, and select Okay"\.
- Check the calibration by injecting 3L into the pneumotach at varying flow rates (0-4 L/s, 4-8 L/s, and 8-12 L/s). The difference from 3 L should be less than 5%.
- Collect the maneuver, ensuring that two sequential breaths are collected and that they are made at the same flow rate.
- Coach the subject to perform a single maneuver consisting of two pairs of breaths – a coaching breath and a breath for analysis. This is shown graphically in Figure 1 (bottom).
- During the maneuver, coach the participants to follow the flow guide on the computer monitor. The investigator may coach the subject by indicating “inhale now” or “exhale now”.
- Perform the maneuver so that there are two pairs of these breaths in a single maneuver. The first exhalation of the maneuver is 3 s and the second is 5 s. Consider adding a resistor in-line with the mouthpiece in order to make exhaled flow easier to control. A resistance with 5 cm H2O/L/s of resistance is generally well-tolerated.
NOTE: It is important that if a resistor is used, it is used throughout the study and for every participant because it increases mouth and airway pressure, which can change airway diameter. It is also important that participants not “puff out” their cheeks as this increases the dead space.
- Measurement protocol
- Instruct the participant to sit straight with both feet on the floor, put nose clips on their nose and place their mouth on the mouthpiece.
- Coach the participant to complete at least one minute of tidal breathing. This is for measures of metabolic function and allows the participant to familiarize themselves with the mouthpiece. After one minute, stop data collection.
- Next, coach the participants to vary their tidal volume, taking either normal, smaller- or larger than normal tidal breaths. This ensures that the capnograms are obtained at different lung volumes
- Coach the participant that they should transition to performing a capnogram maneuver as soon as they see the flow tracing appear on their screen.
- Resume data collection at a random point in the participant’s respiratory cycle. This allows for measurements to be made at different lung volumes.
- Finally, coach to perform a sigh at the end of each maneuver, completely relaxing the muscles of respiration. This allows for FRC to be determined.
- Stop data collection. Repeat Steps 3.6.3-3.6.5 until at least 6-8 maneuvers (12 -16 pairs of breaths for analysis) are completed.
4. Data Analysis
- Exporting Data. To run through the macro, each pair of breaths must be exported as a single text file that is then imported into the macro. Screen shots of this process are given in Supplemental Figure 1.
- Highlight each pair of breaths, taking to care to highlight a portion of the exhalation before the maneuver begins.
- Under the file menu, select Export, and name the subject's maneuver.
- Use the drop-down menu under Save As Type and save it as a data file. Then select Save.
- This will prompt an Export As Text box to appear. On the right deselect Block header Columns, Time, Date, Comments, and Event Markers.
- On the left, select Current Selection and Output NaN for Values. Select Downsample by and enter 10 into the box.
- Select the Flow Channel and the CO2 (%) Channel to be exported and click Okay. Consider making duplicates of these exported files as backups before beginning the analysis.
- Perform the macro analysis. The Step-by-step annotated screen shots of for analyzing exported maneuvers with the macro and comparing to lung volume are given in Supplemental Figure 2 and may be used as a guide.
- Open the macro, go to file, and select Open.
- Select the saved data file, saved with the .txt extension.
- A Text Import Wizard box will appear. In the upper left-hand corner, select Delimited and click Next. For step 2, select Tab under Delimiters and click Next. For step 3, select General under Column Data Format and click Finish.
- To run the macro, select View, Macro, View Macro, and Run in succession. Select Yes if there is a backup copy of the data.
- Allow the macro to run (approximately 90 s) and generate a workbook with four sheets. Of relevance to these measurements, Sheet 2 contains the numeric data and Chart 3 contains a plot of the capnogram.
- Return to the data and determine the volume for FRC. This is identified as the volume at the end of the sigh at which flow = 0 L/s.
- Determine the volume at which the second exhalation in each pair of breaths was begun. By subtracting this from the FRC volume, the starting volume above or below FRC can be determined for each breath.
Representative plethysmography results are given in Figure 4. This participant required four attempts in order to collect three FRC values with <5% variability from the mean.%Ref reflects the percent of the predicted value for each variable based on population regression equations that take into account sex, age, race, height and weight
Figure 1 (top) shows a representative single capnogram used in analysis and Figure 1 (bottom) shows the raw data of the entire sequence of the maneuver. In Figure 1 (bottom), the capnogram and flow tracing are not aligned to account for the time delay. Data generated from running a sequence of breaths through the macro are shown at the end of Supplemental Figure 2. This individual had a dead space of 0.266 L, a slope of 0.523% CO2/L and a normalized slope of 0.0826 L-1. Quality information about the maneuver are also given in columns F, G, I, J, and K. Column F gives the average exhaled flow rate, with the standard deviation in column G. The exhaled tidal volume is given in column J and the R-squared value for the slope is in column K.
Dead space and slope plotted as a function of lung volume are given in Figure 5. In the left panels, dead space and slope are plotted versus lung volume relative to FRC, where FRC=0 L. In the right panels, lung volume and slope are plotted versus absolute lung volume. In both cases, dead space and slope are significantly correlated to lung volume (p<0.05 for all four regression analyses). This suggests that dead space and airways homogeneity increase as lung volume increases, although little is known about this relationship in populations with lung disease or with bronchodilator therapy. The investigator may also choose to use these data to describe the numerical value of dead space and slope at specific lung volumes (FRC, residual volume, 50% of total lung capacity, etc.)3.
Figure 1. Sample capnogram (top), with exhaled CO2 (%) plotted as a function of the exhaled volume. I, II, and III indicate the three phases of the capnogram. The dotted line indicates the volume of the dead space and the solid line represents the slope of the alveolar plateau (Phase III). The slope can be divided by the area under the capnogram (shaded grey, labeled A) to yield the normalized slope. The four breath sequence is shown in the bottom panel, followed by a sigh breath to determine functional residual capacity. Each pair of breaths is analyzed as a single maneuver. Please click here to view a larger version of this figure.
Figure 2. Equipment setup for capnographic measurements. Shown in this figure are the pneumotach and gas analyzer required for capnographic measurements. The left monitor and tracing are used by the participant as a guide in generating the flow pattern while data are observed on the right monitor by the investigator. Please click here to view a larger version of this figure.
Figure 3. Channel settings for the acquisition of the volumetric capnogram. Flow is collected in Channel 1, CO2 concentration (%) is collected in Channel 2, and the tidal volume is calculated in Channel 3. Please click here to view a larger version of this figure.
Figure 4. Representative plethysmograph data from a healthy, male subject. Particularly relevant to the protocol reported here are the total lung capacity (TLC), residual volume (RV) and functional residual capacity (FRC). Please click here to view a larger version of this figure.
Figure 5. Dead space and alveolar slope plotted as a function of absolute lung volume (right panels) and as the volume relative to the functional residual capacity (volume-FRC, left). Note the dependence of the airways volume and lung heterogeneity on lung volume. Lung volume may be expressed as a function of FRC or absolute volume, depending on the experimental design. Please click here to view a larger version of this figure.
Figure 6. Factors impacting data accuracy. Data are given as the mean ± 95% confidence interval. Relationship between the CO2 sampling rate and the time delay between the gas analyzer and pneumotach (top). The time delay should accurately be determined before beginning the experiment. Measuring eight total maneuvers allows for the measurement of the dead space at a single lung volume with <5% variability (bottom). Please click here to view a larger version of this figure.
Here, a protocol for the measurement of VD and airways homogeneity (slope) is provided. These measurements can be made at FRC, or as a function of lung volume. Measuring FRC before the start of the experiment and after a perturbation allows VD and slope to be plotted as a function of lung volume and may provide useful information about the structure-function relationship of the lung that is not obtained from capnography at FRC alone.
Airways volume and high-resolution structure can be obtained from computed tomographic imaging17,18, but this requires exposure to radiation and expertise in image processing. With volumetric capnography, repeated measures can be made without increasing risk to the participant. It also does not require expensive equipment or advanced data processing capabilities. Volumetric capnography is an ideal method for experiments with multiple time points and multiple lung volumes and in-patient populations whose radiation exposure should be minimized.
With regard to the barometric plethysmography, care should be taken to perform the measurement according to consensus statements. When it is important to compare participant values to predicted population values, weight should be measured with a scale and height should be verified with a stadiometer. As noted in the protocol, the most critical component to measure before beginning volumetric capnography is the time delay between the pneumotach and the gas analyzer. The time delay is highly dependent on the analyzer sampling rate (Figure 5, top) and small changes in the sampling rate can have large influences on measured values. The analyzer flow rate should be checked at the beginning and throughout the experiment. Calibration of the analyzer and pneumotach are also critical and care should be taken to ensure their accuracy before beginning an experiment.
We have also determined the accuracy of the measurement at a single lung volume in 3 participants. Figure 5 (bottom) demonstrates that it is necessary to complete four maneuvers (8 total breaths) at a single lung volume to measure dead space so that the variation is <5%. Investigators should take care to make a sufficient number of measurements when having data at a particular lung volume is important. In a subset of 36 maneuvers analyzed in duplicate by two investigators, intra-investigator analysis variability was less than 0.5%.
These methods also require a technician or investigator that is skilled in coaching the participant to make the ventilatory maneuvers. A limitation in pulmonary function studies can be the participant’s ability to perform the maneuver. However, participants that are able to perform clinical pulmonary function are typically able to perform the capnographic maneuvers. If the study is designed such that capnography follows plethysmography and spirometry, participants that are unable to perform a coached spirometric or plethysmographic maneuver can be excluded. In 60 previous studies, one participant who performed clinical spirometry was excluded because they could not follow the capnographic breathing pattern. There are currently no consensus guidelines defining acceptable capnographic measurement criteria. However, intersubject variability is 8±1% of the target flow rate in our 10 most recent participants. Intrasubject (between maneuver) variability is 4±2%.
Issues relating to data accuracy and reproducibility are the result of errors in the time delay or the analyzer and pneumotach calibration. Before each experiment, take care to calibrate the analyzer with a set of known gases and generate a multi-point standard curve to confirm the analyzer’s accuracy.
Beyond the scope of the information provided here, the macro contains two additional calculations that may be of interest. When the maneuvers are made at FRC, the FRC column provides an estimate of FRC based on the Farmery method19. Calculation of the peripheral bronchial cross sectional area is based on the method described by Scherer, et al.20. Finally, if desired, the end tidal CO2 and average expired CO2 concentration can be used to calculate the physiological dead space for comparison to the anatomic dead space21,22.
The authors have nothing to disclose.
This work was funded by the Departments of Health and Human Physiology and Internal Medicine at the University of Iowa. This work was also supported by the Old Gold Fellowship (Bates) and Grant IRG-15-176-40 from the American Cancer Society, administered through The Holden Comprehensive Cancer Center at The University of Iowa (Bates)
|Computer with dual monitor||Dell Instruments|
|PowerLab 8/35*||AD Instruments||PL3508|
|LabChart Data Acquisition Software*||AD Instruments||Version 8|
|Gemini Respiratory Gas Analyzer* (upgraded option)||CWE, Inc||GEMINI 14-10000||*indicates that part is available in the Exercise Physiology package from AD Instruments|
|Heated Pneumotach with Heater Controller* (upgraded option)||Hans Rudolph, Inc||MLT3813H-V|
|3L Calibration Syringe||Vitalograph||36020|
|Nose Clip*||VacuMed||Snuffer 1008|
|Pulse Transducer*||AD Instruments||TN1012/ST|
|Flanged Mouthpiece*||AD Instruments||MLA1026|
|Nafion drying tube with three-way stopcock*||AD Instruments||MLA0343|
|Desiccant cartridge (optional for humid environments)*||AD Instruments||MLA6024|
|Resistor||Hans Rudolph, Inc||7100 R5|
|Flow head adapters*||AD Instruments||MLA1081|
|Modified Tubing Adapter (optional)||AD Instruments||SP0145|
|Two way non-rebreather valve (optional)*||AD Instruments||SP0146|
|High Purity Helium Gas||Praxair||He 4.8|
|6% CO2 and 16% O2 Calibration Gas||Praxair||Custom|
|Microsoft Excel||Microsoft||Office 365|
- Robinson, P. D., et al. Consensus statement for inert gas washout measurement using multiple- and single- breath tests. European Respiratory Journal. 41, (3), 507-522 (2013).
- Verbanck, S., Paiva, M. Gas mixing in the airways and airspaces. Comprehensive Physiology. 1, (2), 809-834 (2011).
- Bates, M. L., et al. Pulmonary function responses to ozone in smokers with a limited smoking history. Toxicology and Applied Pharmacology. 278, (1), 85-90 (2014).
- Bates, M. L., Brenza, T. M., Ben-Jebria, A., Bascom, R., Ultman, J. S. Longitudinal distribution of ozone absorption in the lung: comparison of cigarette smokers and nonsmokers. Toxicology and Applied Pharmacology. 236, (3), 270-275 (2009).
- Reeser, W. H., et al. Uptake of ozone in human lungs and its relationship to local physiological response. Inhalation Toxicology. 17, (13), 699-707 (2005).
- Taylor, A. B., Lee, G. M., Nellore, K., Ben-Jebria, A., Ultman, J. S. Changes in the carbon dioxide expirogram in response to ozone exposure. Toxicology and Applied Pharmacology. 213, (1), 1-9 (2006).
- Baker, L. G., Ultman, J. S., Rhoades, R. A. Simultaneous gas flow and diffusion in a symmetric airway system: a mathematical model. Respiration Physiology. 21, (1), 119-138 (1974).
- Fowler, W. S. Lung Function Studies. II. The Respiratory Dead Space. American Journal of Physiology-Legacy Content. 154, (3), 405-416 (1948).
- Eberlein, M., et al. Supranormal Expiratory Airflow after Bilateral Lung Transplantation Is Associated with Improved Survival. American Journal of Respiratory and Critical Care Medicine. 183, (1), 79-87 (2011).
- Eberlein, M., Schmidt, G. A., Brower, R. G. Chest wall strapping. An old physiology experiment with new relevance to small airways diseases. Annals of the American Thoracic Society. 11, (8), 1258-1266 (2014).
- Taher, H., et al. Chest wall strapping increases expiratory airflow and detectable airway segments in computer tomographic scans of normal and obstructed lungs. Journal of Applied Physiology. (2017).
- Verscheure, S., Massion, P. B., Verschuren, F., Damas, P., Magder, S. Volumetric capnography: lessons from the past and current clinical applications. Critical Care. 20, (1), 184 (2016).
- Suarez-Sipmann, F., Bohm, S. H., Tusman, G. Volumetric capnography: the time has come. Current Opinion in Critical Care. 20, (3), 333-339 (2014).
- Wanger, J., et al. Standardisation of the measurement of lung volumes. European Respiratory Journal. 26, (3), 511-522 (2005).
- Culver, B. H., et al. Recommendations for a Standardized Pulmonary Function Report. An Official American Thoracic Society Technical Statement. American Journal of Respiratory and Critical Care Medicine. 196, (11), 1463-1472 (2017).
- Goldman, H. I., Becklake, M. R. Respiratory function tests; normal values at median altitudes and the prediction of normal results. Am Rev Tuberc. 79, (4), 457-467 (1959).
- Shim, S. S., et al. Lumen area change (Delta Lumen) between inspiratory and expiratory multidetector computed tomography as a measure of severe outcomes in asthmatic patients. J The Journal of Allergy and Clinical. (2018).
- Smith, B. M., et al. Human airway branch variation and chronic obstructive pulmonary disease. Proceedings of the National Academy of Sciences of the United States of America. 115, (5), E974-E981 (2018).
- Farmery, A. D. Volumetric Capnography and Lung Growth in Children - a Simple-Model Validated. Anesthesiology. 83, (6), 1377-1379 (1995).
- Scherer, P. W., Neufeld, G. R., Aukburg, S. J., Hess, G. D. Measurement of Effective Peripheral Bronchial Cross-Section from Single-Breath Gas Washout. Journal of Biomechanical Engineering-Transactions of the Asme. 105, (3), 290-293 (1983).
- Sinha, P., Soni, N. Comparison of volumetric capnography and mixed expired gas methods to calculate physiological dead space in mechanically ventilated ICU patients. Intensive Care Medicine. 38, (10), 1712-1717 (2012).
- Bourgoin, P., et al. Assessment of Bohr and Enghoff Dead Space Equations in Mechanically Ventilated Children. Respiratory Care. 62, (4), 468-474 (2017).