To determine the impact of ventilation frequency, lung volume, and parenchymal stiffness on ventilation distribution, we developed an anatomically-based computational model of the canine lung. Each lobe of the model consists of an asymmetric branching airway network subtended by terminal, viscoelastic acinar units. The model allows for empiric dependencies of airway segment dimensions and parenchymal stiffness on transpulmonary pressure. We simulated the effects of lung volume and parenchymal recoil on global lung impedance and ventilation distribution from 0.1 to 100 Hz, with mean transpulmonary pressures from 5 to 25 cmH2O. With increasing lung volume, the distribution of acinar flows narrowed and became more synchronous for frequencies below resonance. At higher frequencies, large variations in acinar flow were observed. Maximum acinar flow occurred at first antiresonance frequency, where lung impedance achieved a local maximum. The distribution of acinar pressures became very heterogeneous and amplified relative to tracheal pressure at the resonant frequency. These data demonstrate the important interaction between frequency and lung tissue stiffness on the distribution of acinar flows and pressures. These simulations provide useful information for the optimization of frequency, lung volume, and mean airway pressure during conventional ventilation or high frequency oscillation (HFOV). Moreover our model indicates that an optimal HFOV bandwidth exists between the resonant and antiresonant frequencies, for which interregional gas mixing is maximized.