Abstract. Biogenic volatile organic compounds (BVOCs) emitted by terrestrial vegetation significantly influence the atmospheric chemistry and global climate. Previous studies calculated BVOC emissions outside the comprehensive terrestrial ecosystem processes and oversimplified the representation of canopy environments which exert substantial impacts on BVOC emissions. The sophisticated land surface model CLM (Community Land Model) reproduces essential ground and canopy characteristics and can calculate BVOC emissions as a step of integrated biogeochemical processes. In this study, the land surface scheme CLM version4 (CLM4) of WRF (Weather Research and Forecasting model) was used to estimate BVOC emissions in China. Based on highly-resolved meteorological outputs derived from WRF, CLM4 computed real-time physical and biological variables to drive MEGAN (Model of Emissions of Gases and Aerosols from Nature), a BVOC estimation algorithm embedded within biogeochemistry component of CLM4, to calculate plant emission flux. MODIS (Moderate Resolution Imaging Spectroradiometer) land-use data with high resolutions was introduced in WRF to replace the outdated land surface parameters. An emission inventory of isoprene and monoterpenes (including α-pinene, β-pinene, 3-carene, t-β-ocimene, limonene, sabinene and myrcene) with high spatiotemporal resolution (12 × 12 km, hourly) was established for the year 2018. The annual BVOC emission in China was 14.7 Tg C, in which isoprene contributes about 78.3 % (11.5 Tg C), followed by α-pinene (1.2 Tg C) and β-pinene (0.7 Tg C). Due to the strong emission capacity and large areas, broadleaf forests contribute to 76.8 % of total isoprene emission and 72.1 % of monoterpenes emission, respectively. BVOC emissions showed marked seasonal and diurnal patterns with the peak emission occurring in summer and midday. Spatially, high emissions of BVOC were mainly concentrated in southern and northeastern China, as well as the Qinling Mountains in central China, accounting for 91.4 % of national emission. Guangxi, Yunnan and Hunan provinces are significant emitters due to large area of vegetation with high emission rate and favored environmental conditions. The emission estimates are compared to past modeling results, field measurements and further evaluated against top-down isoprene emission estimates. Generally, the coupled mode produced a reasonable simulation in both emission amounts and the spatiotemporal distribution of BVOCs. The WRF–CLM–MEGAN coupling framework could be further integrated with atmospheric chemistry model to investigate BVOC chemistry and their effects on regional pollution and climate.