We determined the short-term variability of pulmonary function in infants recovering from bronchopulmonary dysplasia. Sixteen infants with birth weight of 1,231 +/- 929 grams (mean +/- SD) and gestational age of 29 +/- 4 weeks were studied twice at 17 +/- 8 weeks postnatally at intervals of 4 to 8 days during a period of clinical stability. The infants were still on supplemental oxygen but were off diuretics and bronchodilators. We used a modification of the rapid thoracic compression method to measure forced expiratory flow (Vmax FRC) and the time constant (tau) of expiratory flow at functional residual capacity. Thoracic gas volume (TGV), mean and total airway resistance (RawM and RawT), and mean and total specific airway conductance (SGawM and SGawT) were measured in a whole body pressure plethysmograph. An esophageal balloon was used to measure dynamic pulmonary compliance (Cdyn). Variabilities were defined as the standard deviation of percent changes between the first and second tests. They were 30% for VmaxFRC, 23% for tau, 12% for TGV, 20% for RawM, 35% for RawT, 25% for SGawM, 72% for SGawT, and 23% for Cdyn. All these tests are useful in assessing pulmonary function of infants with BPD; however, their variability must be taken into account when interpreting short-term changes.