Abstract. This study analyzes simulated regional-scale ozone burdens both near the surface and aloft, estimates process contributions to these burdens, and calculates the sensitivity of the simulated regional-scale ozone burden to several key model 15 inputs with a particular emphasis on boundary conditions derived from hemispheric or global scale models. The Community Multiscale Air Quality (CMAQ) model simulations supporting this analysis were performed over the continental U.S. for the Hypothetical bounding scenarios were performed to quantify the effects of emissions, boundary conditions, and ozone dry deposition on the simulated ozone burden. Analysis of these simulations confirms that the characterization of ozone outside the regional-scale modeling domain can have a profound impact on simulated regional-scale ozone. This was further investigated by using data from four hemispheric or global modeling systems (Chemistry -Integrated Forecasting Model (C-25 IFS), CMAQ extended for hemispheric applications (H-CMAQ), GEOS-Chem, and AM3) to derive alternate boundary conditions for the regional-scale CMAQ simulations. The regional-scale CMAQ simulations using these four different boundary conditions showed that the largest ozone abundance in the upper layers was simulated when using boundary conditions from GEOS-Chem, followed by the simulations using C-IFS, AM3, and H-CMAQ boundary conditions, consistent with the analysis of the ozone fields from the global models along the CMAQ boundaries. Using boundary conditions from 30 AM3 yielded higher springtime ozone columns burdens in the mid-and lower troposphere compared to boundary conditions from the other models. For surface ozone, the differences between the AM3-driven CMAQ simulations and the CMAQ simulations driven by other large-scale models are especially pronounced during spring and winter where they can reach more Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2017-676 Manuscript under review for journal Atmos. Chem. Phys. Discussion started: 8 August 2017 c Author(s) 2017. CC BY 4.0 License. 2 than 10 ppb for seasonal mean ozone mixing ratios and as much as 15 ppb for domain-averaged daily maximum 8-hr average ozone on individual days. In contrast, the differences between the C-IFS, GEOS-Chem, and H-CMAQ driven regional-scale CMAQ simulations are typically smaller. Comparing simulated surface ozone mixing ratios to observations and computing seasonal and regional model performance statistics revealed that boundary conditions can have a substantial impact on model performance. Further analysis showed that boundary conditions can affect model performance across the entire range of the 5 observed distribution, although the impacts tend to be lower during summer and for the very highest observed percentiles. The results are discussed in the context of future model development and analysis opportunities.