Abstract. Emissions of biogenic volatile organic compounds (BVOCs) have changed in the past millennium due to changes in land use, temperature, and CO 2 concentrations. Recent reconstructions of BVOC emissions have predicted that global isoprene emissions have decreased, while monoterpene and sesquiterpene emissions have increased; however, all three show regional variability due to competition between the various influencing factors.In this work, we use two modeled estimates of BVOC emissions from the years 1000 to 2000 to test the effect of anthropogenic changes to BVOC emissions on secondary organic aerosol (SOA) formation, global aerosol size distributions, and radiative effects using the GEOS-Chem-TOMAS (Goddard Earth Observing System; TwO-Moment Aerosol Sectional) global aerosol microphysics model. With anthropogenic emissions (e.g., SO 2 , NO x , primary aerosols) turned off and BVOC emissions changed from year 1000 to year 2000 values, decreases in the number concentration of particles of size D p > 80 nm (N80) of > 25 % in year 2000 relative to year 1000 were predicted in regions with extensive landuse changes since year 1000 which led to regional increases in the combined aerosol radiative effect (direct and indirect) of > 0.5 W m −2 in these regions. We test the sensitivity of our results to BVOC emissions inventory, SOA yields, and the presence of anthropogenic emissions; however, the qualitative response of the model to historic BVOC changes remains the same in all cases. Accounting for these uncertainties, we estimate millennial changes in BVOC emissions cause a global mean direct effect of between +0.022 and +0.163 W m −2 and the global mean cloud-albedo aerosol indirect effect of between −0.008 and −0.056 W m −2 . This change in aerosols, and the associated radiative forcing, could be a largely overlooked and important anthropogenic aerosol effect on regional climates.