Simulation of radon-222 with the GEOS-Chem global model: emissions, seasonality, and convective transport

Radon-222 (Rn-222) is a short-lived radioactive gas naturally emitted from land surfaces and has long been used to assess convective transport in atmospheric models. In this study, we simulate Rn-222 using the GEOS-Chem chemical transport model to improve our understanding of Rn-222 emissions and surface concentration seasonality and characterize convective transport associated with two Goddard Earth Observing System (GEOS) meteorological products, the Modern-Era Retrospective analysis for Research and Applications (MERRA) and GEOS Forward Processing (GEOS-FP). We evaluate four global Rn-222 emission scenarios by comparing model results with observations at 51 surface sites. The default emission scenario in GEOS-Chem yields a moderate agreement with surface observations globally (68.9 % of data within a factor of 2) and a large underestimate of winter surface Rn-222 concentrations at Northern Hemisphere midlatitudes and high latitudes due to an oversimplified formulation of Rn-222 emission fluxes (1 atom cm(-2) s(-1) over land with a reduction by a factor of 3 under freezing conditions). We compose a new global Rn-222 emission scenario based on Zhang et al. (2011) and demonstrate its potential to improve simulated surface Rn-222 concentrations and seasonality. The regional components of this scenario include spatially and temporally varying emission fluxes derived from previous measurements of soil radium content and soil exhalation mod- els, which are key factors in determining Rn-222 emission flux rates. However, large model underestimates of surface Rn-222 concentrations still exist in Asia, suggesting unusually high regional Rn-222 emissions. We therefore propose a conservative upscaling factor of 1.2 for Rn-222 emission fluxes in China, which was also constrained by observed deposition 210pb fluxes of (a progeny of Rn-222). With this modification, the model shows better agreement with observations in Europe and North America (> 80 % of data within a factor of 2) and reasonable agreement in Asia (close to 70 %). Further constraints on Rn-222 emissions would require additional concentration and emission flux observations in the central United States, Canada, Africa, and Asia. We also compare and assess convective transport in model simulations driven by MERRA and GEOS-FP using observed Rn-222 vertical profiles in northern midlatitude summer and from three short-term airborne campaigns. While simulations with both GEOS products are able to capture the observed vertical gradient of Rn-222 concentrations in the lower troposphere (0-4 km), neither correctly represents the level of convective detrainment, resulting in biases in the middle and upper troposphere. Compared with GEOS-FP, MERRA leads to stronger convective transport of 222 Rn, which is partially compensated for by its weaker large-scale vertical advection, resulting in similar global vertical distributions of Rn-222 concentrations between the two simulations. This has important implications for using chemical transport models to interpret the transport of other trace species when these GEOS products are used as driving meteorology.