We investigated the millennial variability (1000 A.D.-2000 A.D.) of global biogenic volatile organic compound (BVOC) emissions by using two independent numerical models: The Model of Emissions of Gases and Aerosols from Nature (MEGAN), for isoprene, monoterpene, and sesquiterpene, and Lund-Potsdam-Jena-General Ecosystem Simulator (LPJ-GUESS), for isoprene and monoterpenes. We found the millennial trends of global isoprene emissions to be mostly affected by land cover and atmospheric carbon dioxide changes, whereas monoterpene and sesquiterpene emission trends were dominated by temperature change. Isoprene emissions declined substantially in regions with large and rapid land cover change. In addition, isoprene emission sensitivity to drought proved to have significant short-term global effects. By the end of the past millennium MEGAN isoprene emissions were 634 TgC yr-1 (13% and 19% less than during 1750-1850 and 1000-1200, respectively), and LPJ-GUESS emissions were 323 TgC yr-1(15% and 20% less than during 1750-1850 and 1000-1200, respectively). Monoterpene emissions were 89 TgC yr-1(10% and 6% higher than during 1750-1850 and 1000-1200, respectively) in MEGAN, and 24 TgC yr-1 (2% higher and 5% less than during 1750-1850 and 1000-1200, respectively) in LPJ-GUESS. MEGAN sesquiterpene emissions were 36 TgC yr-1(10% and 4% higher than during 1750-1850 and 1000-1200, respectively). Although both models capture similar emission trends, the magnitude of the emissions are different. This highlights the importance of building better constraints on VOC emissions from terrestrial vegetation.