Heterogeneous reactions occurring at the surface of atmospheric aerosol particles regulate the production and lifetime of a wide array of atmospheric gases. Aerosol surface area plays a critical role in setting the rate of heterogeneous reactions in the atmosphere. Despite the central role of aerosol surface area, there are few assessments of the accuracy of aerosol surface area concentrations in regional and global models. In this study, we compare aerosol surface area concentrations in the EPA's Community Multiscale Air Quality (CMAQ) model with commensurate observations from the 2011 NASA flight-based DISCOVER-AQ (Deriving Information on Surface Conditions from COlumn and VERtically Resolved Observations Relevant to Air Quality) campaign. The study region includes the Baltimore and Washington, D.C. metropolitan area. Dry aerosol surface area was measured aboard the NASA P-3B aircraft using an ultra-high-sensitivity aerosol spectrometer (UHSAS). We show that modeled and measured dry aerosol surface area, Sa,mod and Sa,meas respectively, are modestly correlated (r2=0.52) and on average agree to within a factor of 2 (Sa,mod/Sa,meas=0.44) over the course of the 13 research flights. We show that Sa,mod/Sa,meas does not depend strongly on photochemical age or the concentration of secondary biogenic aerosol, suggesting that the condensation of low-volatility gas-phase compounds does not strongly affect model–measurement agreement. In comparison, there is strong agreement between measured and modeled aerosol number concentration (Nmod/Nmeas=0.87, r2=0.63). The persistent underestimate of Sa in the model, combined with strong agreement in modeled and measured aerosol number concentrations, suggests that model representation of the size distribution of primary emissions or secondary aerosol formed at the early stages of oxidation may contribute to the observed differences. For reactions occurring on small particles, the rate of heterogeneous reactions is a linear function of both Sa and the reactive uptake coefficient (γ). To assess the importance of uncertainty in modeled Sa for the representation of heterogeneous reactions in models, we compare both the mean and the variance in Sa,mod/Sa,meas to those in γ(N2O5)mod/γ(N2O5)meas. We find that the uncertainty in model representation of heterogeneous reactions is primarily driven by uncertainty in the parametrization of reactive uptake coefficients, although the discrepancy between Sa,mod and Sa,meas is not insignificant. Our analysis suggests that model improvements to aerosol surface area concentrations, in addition to more accurate parameterizations of heterogeneous kinetics, will advance the representation of heterogeneous chemistry in regional models.