Historical (1960–2014) lightning and LNO _x trends and their controlling factors in a chemistry–climate model

Lightning can cause natural hazards that result in human and animal injuries and fatalities, infrastructure destruction, and wildfire ignition. Lightning-produced inline-formulaNO_x (Linline-formulaNO_x), a major inline-formulaNO_x (inline-formula $M6inlinescrollmathmlchem{\mathrm{normal\; NO}}_x=chem\mathrm{normal\; NO}+chem{\mathrm{normal\; NO}}_{\mathrm{normal\; 2}}$ 85pt13ptsvg-formulamathimgf0add4bbe2151ecfa7cd944e28fa7e9e acp-23-13061-2023-ie00001.svg85pt13ptacp-23-13061-2023-ie00001.png ) source, plays a vital role in atmospheric chemistry and global climate. The Earth has experienced marked global warming and changes in aerosol and aerosol precursor emissions (AeroPEs) since the 1960s. Investigating long-term historical (1960–2014) lightning and inline-formulaLNO_x trends can provide important indicators for all lightning-related phenomena and for inline-formulaLNO_x effects on atmospheric chemistry and global climate. Understanding how global warming and changes in AeroPEs influence historical lightning and inline-formulaLNO_x trends can be helpful in providing a scientific basis for assessing future lightning and inline-formulaLNO_x trends. Moreover, global lightning activities' responses to large volcanic eruptions such as the 1991 Pinatubo eruption are not well elucidated and are worth exploring. This study employed the widely used cloud top height lightning scheme (CTH scheme) and the newly developed ice-based ECMWF-McCAUL lightning scheme to investigate historical (1960–2014) lightning and inline-formulaLNO_x trends and variations as well as their influencing factors (global warming, increases in AeroPEs, and the Pinatubo eruption) in the framework of the CHASER (MIROC) chemistry–climate model. The results of the sensitivity experiments indicate that both lightning schemes simulated almost flat global mean lightning flash rate anomaly trends during 1960–2014 in CHASER (the Mann–Kendall trend test (significance inferred as 5 %) shows no trend for the ECMWF-McCAUL scheme, but a 0.03 inline-formula% yr⁻¹ significant increasing trend is detected for the CTH scheme). Moreover, both lightning schemes suggest that past global warming enhances historical trends for global mean lightning density and global inline-formulaLNO_x emissions in a positive direction (around 0.03 inline-formula% yr⁻¹ or 3 inline-formula% K⁻¹). However, past increases in AeroPEs exert an opposite effect on the lightning and inline-formulaLNO_x trends (inline-formula−0.07 % to inline-formula−0.04 inline-formula% yr⁻¹ for lightning and inline-formula−0.08 % to inline-formula−0.03 inline-formula% yr⁻¹ for Linline-formulaNO_x) when one considers only the aerosol radiative effects in the cumulus convection scheme. Additionally, effects of past global warming and increases in AeroPEs in lightning trends were found to be heterogeneous across different regions when analyzing lightning trends on the global map. Lastly, this paper is the first of study results suggesting that global lightning activities were markedly suppressed during the first year after the Pinatubo eruption as shown in both lightning schemes (global lightning activities decreased by as much as 18.10 % as simulated by the ECMWF-McCAUL scheme). Based on the simulated suppressed lightning activities after the Pinatubo eruption, the findings also indicate that global inline-formulaLNO_x emissions decreased after the 2- to 3-year Pinatubo eruption (1.99 %–8.47 % for the annual percentage reduction). Model intercomparisons of lightning flash rate trends and variations between our study (CHASER) and other Coupled Model Intercomparison Project Phase 6 (CMIP6) models indicate great uncertainties in historical (1960–2014) global lightning trend simulations. Such uncertainties must be investigated further.