Analysis of Secondary Organic Aerosol Simulation Bias in the Community Earth System Model (CESM2.1)

Organic aerosol (OA) has been considered as one of the most important uncertainties in climate modeling due to the complexity in presenting its chemical production and depletion mechanisms. To better understand the capability of climate models and probe into the associated uncertainties in simulating OA, we evaluate the Community Earth System Model version 2.1 (CESM2.1) configured with the Community Atmosphere Model version 6 (CAM6) with comprehensive tropospheric and stratospheric chemistry representation (CAM6-Chem), through a long-term simulation (1988–2019) with observations collected from multiple datasets in the United States. We find that CESM generally reproduces the inter-annual variation and seasonal cycle of OA mass concentration at surface layer with correlation of 0.40 as compared to ground observations, and systematically overestimates (69 %) in summer and underestimates (-19 %) in winter. Through a series of sensitivity simulations, we reveal that modeling bias is primarily related to the dominant fraction of monoterpene-formed secondary organic aerosol (SOA), and a strong positive correlation of 0.67 is found between monoterpene emission and modeling bias in eastern US during summer. In terms of vertical profile, the model prominently underestimates OA and monoterpene concentrations by 37–99 % and 82–99 % respectively in the upper air (>500 m) as validated against aircraft observations. Our study suggests that the current Volatility Basis Set (VBS) scheme applied in CESM might be parameterized with too high monoterpene SOA yields which subsequently result in strong SOA production near emission source area. We also find that the model has difficulty in reproducing the decreasing trend of surface OA in southeast US, probably because of employing pure gas VBS to represent isoprene SOA which is in reality mainly formed through multiphase chemistry, thus the influence of aerosol acidity and sulfate particle change on isoprene SOA formation has not been fully considered in the model. This study reveals the urgent need to improve the SOA modeling in climate models.