Abstract-Ozone is the most important plant-damaging air pollutant in the United States today, causing annual crop losses estimated at greater than two billion dollars. The atmospheric ozone concentration that surrounds a plant is not the concentration that actually impinges upon the plant cells, because the plant's cuticle acts as a barrier to direct diffusion of ozone into cells for much of the plant surface. The primary avenue for ozone entry is via the stomata, which are adjustable pores in the epidermis. Ozone production and vapor pressure deficit (VPD) exhibit nonlinear diurnal cycles, which vary from location to location and, at a given location, vary with the seasons. VPD, a measure of the joint effect of temperature and humidity on the water potential gradient between a plant surface and the air, is highly correlated with stomatal closing and can be calculated from atmospheric data.Each plant species has its own threshold set of VPD and atmospheric ozone concentration, below which stomates are fully open and above which stomates are closed. Hence, combining ozone and VPD data with knowledge of how these variables affect the stomata allows one to model the effective ozone dose reaching the plant cells, a dose which typically cannot be measured directly.In this paper, we derive a diffusion model, consistent with Fick's first law, that predicts the ozone concentration reaching the cells adjacent to the substomatal cavity at a given time, using the concurrent ambient ozone concentration and VPD together with species-specific thresholds of ozone concentration and VPD below which stomates are fully open and above which stomates are closed.-1--2-Combining this diffusion model with particular functional forms for the daily curves of atmospheric ozone concentration and VPD allows one to calculate the expected daily ozone dose that the plant's cells receive and the variance of that dose. In addition, this methodology can be modified using hierarchical models to provide realistic regional estimates of the effective daily ozone dose for a species and the variance of the dose, which reflect the regional variation in the diurnal cycles of both atmospheric ozone concentration and VPD. The ozone dose and its variance predicted by this model can be used to assess ozone impact on red spruce in the northeastern US.