Articles | Volume 23, issue 20
https://doi.org/10.5194/acp-23-13061-2023
https://doi.org/10.5194/acp-23-13061-2023
Research article
 | 
17 Oct 2023
Research article |  | 17 Oct 2023

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

Yanfeng He and Kengo Sudo
Abstract

Lightning can cause natural hazards that result in human and animal injuries and fatalities, infrastructure destruction, and wildfire ignition. Lightning-produced NOx (LNOx), a major NOx (NOx=NO+NO2) 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 LNOx trends can provide important indicators for all lightning-related phenomena and for LNOx effects on atmospheric chemistry and global climate. Understanding how global warming and changes in AeroPEs influence historical lightning and LNOx trends can be helpful in providing a scientific basis for assessing future lightning and LNOx 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 LNOx 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 % yr−1 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 LNOx emissions in a positive direction (around 0.03 % yr−1 or 3 % K−1). However, past increases in AeroPEs exert an opposite effect on the lightning and LNOx trends (0.07 % to 0.04 % yr−1 for lightning and 0.08 % to 0.03 % yr−1 for LNOx) 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 LNOx 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.

1 Introduction

Lightning, an extremely energetic natural phenomenon, occurs at every moment somewhere on Earth: its average occurrence frequency is approximately 46 times per second (Cecil et al., 2014). Lightning generation is associated with electric charge separation, which is mainly realized by collisions between graupel and hail as well as other types of hydrometeors within convective clouds (Lopez, 2016). As a natural hazard, lightning can cause human and animal injuries and fatalities, infrastructure destruction, and wildfire ignition (Cerveny et al., 2017; Cooper and Holle, 2019; Jensen et al., 2022; Veraverbeke et al., 2022). Lightning-produced NOx (LNOx) accounts for around 10 % of the global tropospheric NOx (NOx=NO+NO2) source. It is regarded as the dominant NOx source in the middle to upper troposphere (Schumann and Huntrieser, 2007; Finney et al., 2016a). Moreover, LNOx plays a crucially important role in atmospheric chemistry and global climate by affecting the abundances of OH radicals, important greenhouse gases (GHGs) such as ozone and methane, and other trace gases (Labrador et al., 2005; Schumann and Huntrieser, 2007; Wild, 2007; Liaskos et al., 2015; Finney et al., 2016b; Murray, 2016; Tost, 2017; He et al., 2022b).

Reportedly, the lightning flash rate (LFR) is related to the stage of convective cloud development (Williams et al., 1989), convective available potential energy (CAPE) (Romps et al., 2014), cloud liquid–ice water content (Saunders et al., 1991; Finney et al., 2014), and even the convective precipitation volume (Buechler et al., 1990; McCaul et al., 2009; Romps et al., 2014). Long-term global warming is associated with changes in the overall temperature and relative humidity profiles in the atmosphere and global convective adjustment (Manabe and Wetherald, 1975; Del Genio et al., 2007), which can strongly affect the lightning-related factors described above. Consequently, long-term global warming can be a fundamentally important factor affecting long-term variations in global lightning activity. Findings from many earlier numerical simulation studies manifest that global lightning activities are sensitive to long-term global warming, with most studies showing 5 %–16 % (average around 10 %) increases in global lightning activities per 1 K global warming (Price and Rind, 1994; Zeng et al., 2008; Hui and Hong, 2013; Banerjee et al., 2014; Krause et al., 2014; Romps et al., 2014; Clark et al., 2017). However, other numerical simulation studies such as those using an ice-based lightning scheme or convective mass flux as a proxy to parameterize lightning have yielded opposite results, suggesting that global lightning activity will decrease under long-term global warming (Clark et al., 2017; Finney et al., 2018).

Aside from long-term global warming, changes in aerosol loading can also be responsible for long-term global lightning activity variations. Aerosols influence lightning activity through aerosol radiative and microphysical effects, but the degree to which the two distinct effects influence regional- or global-scale lightning activities remains unclear (Yuan et al., 2011; Yang et al., 2013; Tan et al., 2016; Altaratz et al., 2017; Wang et al., 2018; Liu et al., 2020). Further research is needed. It is urgently necessary to elucidate the effects of aerosol radiative and microphysical effects on lightning on a global scale. The aerosol radiative effects indicate that aerosols can heat the atmospheric layer and can cool the Earth's surface by absorbing and scattering solar radiation (Kaufman et al., 2002; Koren et al., 2004, 2008; Li et al., 2017). Thereby, convection and electrical activities are likely to be inhibited (Koren et al., 2004; Yang et al., 2013; Tan et al., 2016). The microphysical effects suggest that by acting as cloud condensation nuclei (CCN) or as ice nuclei, aerosols can reduce the mean size of cloud droplets, consequently suppressing the coalescence of cloud droplets into raindrops. As a result, more liquid water particles are uplifted to higher mixed-phase regions of the troposphere, where they invigorate lightning (Wang et al., 2018; Liu et al., 2020).

The Earth has experienced a considerable degree of global warming and changes in aerosol and aerosol precursor emissions (AeroPEs) since the 1960s (Hoesly et al., 2018; NOAA National Centers for Environmental Information (NCEI), 2022). However, how historical lightning has trended and how lightning has responded to historical global warming and changes in AeroPEs are not well examined. This topic is worth exploring because historical lightning densities are indicators of all lightning-related phenomena (Price and Rind, 1994). Exploring the historical global LNOx emission trend is also meaningful because it can indicate the effects of LNOx emissions on atmospheric chemistry and global climate. Furthermore, investigating the effects of historical global warming and increases in AeroPEs on historical lightning and LNOx trends can provide a basis for assessing future lightning and LNOx trends.

Large-scale volcanic eruptions such as the 1991 Pinatubo eruption inject tremendous amounts of sulfuric gas into the stratosphere, where it converts to H2SO4 aerosols. Consequently, the stratospheric aerosols have increased in abundance after the volcanic eruptions. The enhanced stratospheric aerosol layer can cool the Earth's surface heterogeneously and can decrease the total amount of water in the atmosphere (Soden et al., 2002; Boucher, 2015, p. 63). The near-global perturbations in the radiative energy balance and meteorological fields caused by such strong volcanic eruptions might influence global lightning activities. If so, there might be ramifications for all lightning-related phenomena. Nevertheless, they remain poorly understood.

In our earlier work, we developed a new process and ice-based lightning scheme called the ECMWF-McCAUL scheme (He et al., 2022b). This lightning scheme was developed by combining benefits of the lightning scheme used in the European Centre for Medium-Range Weather Forecasts (ECMWF) forecasting system (Lopez, 2016) and those presented in reports by McCaul et al. (McCaul et al., 2009). The ECMWF-McCAUL scheme simulated the best lightning density spatial distributions among four existing lightning schemes when compared against satellite lightning observations (Lightning Imaging Sensor (LIS) and Optical Transient Detector (OTD)) during 2007–2011. The sensitivity of global lightning activity to changes in surface temperature on a decadal timescale was estimated as 10.13 % K−1 using the ECMWF-McCAUL scheme (He et al., 2022b), which is close to most past estimates (average of around 10 % K−1).

Using the chemistry–climate model CHASER (MIROC) with two lightning schemes (the widely used cloud top height scheme and the ice-based ECMWF-McCAUL scheme), we quantitatively investigated historical lightning and LNOx trends and ascertained how global warming, increases in AeroPEs, and the Pinatubo eruption respectively influenced these trends. Using two lightning schemes, we demonstrated the sensitivities of different lightning schemes to historical global warming, increases in AeroPEs, and the Pinatubo eruption.

Research methods including the model description and experiment setup are described in Sect. 2. In Sect. 3.1, the simulated historical lightning distributions and trends are validated using LIS and OTD lightning observations. Section 3.2 presents the effects of global warming and increases in AeroPEs on historical lightning and LNOx trends. In Sect. 3.3, the Pinatubo volcanic eruption effects on historical lightning and LNOx trends are discussed. Section 3.4 elucidated model intercomparisons of LFR trends and variation between our study (CHASER) and other CMIP6 model outputs. Section 4 presents relevant discussions and conclusions based on these study findings.

2 Method

2.1 Chemistry–climate model

We used the CHASER (MIROC) global chemistry–climate model (Sudo et al., 2002; Sudo and Akimoto, 2007; Watanabe et al., 2011; Ha et al., 2021) for this study, which incorporated the consideration of detailed chemical and physical processes in the troposphere and stratosphere. The CHASER version adopted for this study simulates the distributions of 94 chemical species while reflecting the effects of 269 chemical reactions (58 photolytic, 190 kinetic, and 21 heterogeneous). As processes associated with tropospheric chemistry, non-methane hydrocarbon (NMHC) oxidation and the fundamental chemical cycle of OxNOxHOxCH4–CO are considered. CHASER simulates stratospheric chemistry involving the Chapman mechanisms and catalytic reactions associated with HOx, NOx, ClOx, and BrOx. Moreover, it simulates the formation of polar stratospheric clouds (PSCs) and heterogeneous reactions occurring on their surfaces. CHASER is online coupled to MIROC (Model for Interdisciplinary Research on Climate) AGCM (Atmospheric General Circulation Model) ver. 5.0 (Watanabe et al., 2011), which simulates cumulus convection (Arakawa–Schubert scheme) and grid-scale large-scale condensation to represent cloud and precipitation processes. The radiation flux is calculated using a two-stream k-distribution radiation scheme, which considers absorption, scattering, and emissions by aerosol and cloud particles as well as by gaseous species (Sekiguchi and Nakajima, 2008; Goto et al., 2015). The aerosol component in CHASER is coupled with the SPRINTARS aerosol model (Takemura et al., 2009), particularly for simulating primary organic carbon, sea salt, and dust, which is also based on MIROC. The aerosol radiation effects are considered in both large-scale condensation and cumulus convection schemes, although the aerosol microphysical effects are only reflected in the large-scale condensation scheme.

This study used a horizontal resolution of T42 (2.8×2.8), with a vertical resolution of 36 σ–p hybrid levels from the surface to approximately 50 km. Anthropogenic and biomass burning emissions were obtained from the CMIP6 forcing datasets (van Marle et al., 2017; Hoesly et al., 2018) for 1959–2014 (https://esgf-node.llnl.gov/search/input4mips/, last access: 19 September 2022). Interannual variations in biogenic emissions for isoprene, monoterpene, acetone, and methanol were considered using an offline simulation by the Vegetation Integrative Simulator for Trace Gases (VISIT) terrestrial ecosystem model (Ito and Inatomi, 2012). The residual biogenic emissions (ethane, propane, ethylene, propene) used are climatological values derived from the Model of Emissions of Gases and Aerosols from Nature (MEGAN) modeling system (Guenther et al., 2012).

The CHASER (MIROC) global chemistry–climate model originally parameterized lightning with the widely used cloud top height scheme (Price and Rind, 1992). Here, a newly developed ice-based lightning scheme called the ECMWF-McCAUL has been implemented into CHASER (MIROC) (He et al., 2022b). The ECMWF-McCAUL scheme computes LFRs as a function of CAPE and QRa (QRa represents the total volumetric amount of cloud ice, graupel, and snow in the charge separation region). Compared with the cloud top height, a salient advantage of the ECMWF-McCAUL scheme is that it has a direct physical link with the charging mechanism.

2.2 Lightning NOx emission parameterizations

We tested two lightning schemes for this study. The first lightning scheme is the widely used cloud top height (CTH) scheme (Price and Rind, 1992), which was originally used in CHASER (MIROC). This lightning scheme uses the following equations to calculate LFR.

(1)Fl=3.44×10-5H4.9(2)Fo=6.2×10-4H1.73

Herein, F represents the total flash frequency (fl.min-1), H stands for the cloud top height (km), and subscripts l and o respectively denote the land and ocean (Price and Rind, 1992). However, we realize the CTH scheme in CHASER using Eqs. 3 and 4 (Sudo et al., 2002). Each model layer's cumulus cloud fractions are used to weight the calculated lightning densities from that layer in the CTH scheme.

(3)Fl=i=1n=36adj_factor×Cu_CFi×(Hi-Hsurface)4.9(4)Fo=i=1n=36adj_factor×Cu_CFi×(Hi-Hsurface)1.73

In these equations, i represents the model layer index. In addition, adj_factor represents adjustment factors that differ for different model layers and model grids. Cu_CFi symbolizes the cumulus cloud fraction at model layer i. Hi and Hsurface respectively denote the altitude of model layer i and the altitude of the model's surface layer.

The second lightning scheme used for this study is a newly developed one named the ECMWF-McCAUL scheme (He et al., 2022b), which is based on the original ECMWF scheme and findings reported by McCaul et al. (2009). The ECMWF-McCAUL scheme calculates LFRs as a function of CAPE (m2 s−2) and QRa (QRa symbolizes the total volumetric amount of cloud ice, graupel, and snow in the charge separation region) as

(5)fl=αlQRaCAPE1.3,(6)fo=αoQRaCAPE1.3,

where fl and fo respectively symbolize the total flash density (fl.m-2s-1) over land and ocean. In addition, αl and αo are constants (fl.s1.6kg-1m-2.6) determined after calibration against LIS–OTD climatology respectively for land and ocean. For this study, αl and αo are set respectively as 2.67×10-16 and 1.68×10-17. In the charge separation region (from 0 to 25 C isotherm), QRa (kg m−2) is expressed as a proxy for the charging rate because of collisions between graupel and hydrometeors of other types (McCaul et al., 2009). Moreover, QRa represents the total volumetric amount of hydrometeors of three kinds (graupel, snow, and cloud ice) within the charge separation region, calculated as

(7) Q Ra = z 0 z - 25 ( q graup + q snow + q ice ) ρ d z ,

where qgraup, qsnow, and qice respectively represent the mass mixing ratios (kg kg−1) of graupel, snow, and cloud ice. In addition, qice was diagnosed using the Arakawa–Schubert cumulus parameterization. Then, qgraup and qsnow were computed at each vertical level of the model using the following equations.

(8) q graup = β P f ρ V graup
(9) q snow = ( 1 - β ) P f ρ V snow

In these equations, Pf represents the vertical profile of the frozen precipitation convective flux (kgm-2s-1), ρ denotes the air density (kg m−3), and Vgraup and Vsnow respectively express the typical fall speeds for graupel and snow set to 3.1 and 0.5 m s−1 for this study. For land, the dimensionless coefficient β is set to 0.7, whereas it is set to 0.45 for oceans to consider the observed lower graupel content over the oceans.

Based on the cold cloud depth, a fourth-order polynomial (Eq. 10) is used to calculate the proportion of total flashes that are cloud-to-ground lightning flashes (p). An earlier report of the literature describes the method (Price and Rind, 1993).

(10) p = 1 64.9 - 36.54 D + 7.493 D 2 - 0.648 D 3 + 0.021 D 4

The depth of the cloud above the 0 C isotherms is represented by D (km) in this equation.

According to recent studies, the intra-cloud (IC) lightning flashes are as efficient as cloud-to-ground (CG) lightning flashes at producing NOx. The lightning NOx production efficiency is estimated as 100–400 mol per flash (Ridley et al., 2005; Cooray et al., 2009; Ott et al., 2010; Allen et al., 2019). The LNOx production efficiencies for IC and CG are therefore set to the same value (250 mol per flash) in CHASER, which is the median of the commonly cited range of 100–400 mol per flash. Therefore, in this study, the distinctions between IC and CG do not affect the distribution or magnitude of LNOx emissions. It is noteworthy that marked uncertainties are involved in ascertaining the LNOx production efficiency (Allen et al., 2019; Bucsela et al., 2019). The choice of a different LNOx production efficiency might affect the simulation of LNOx emissions. Further research must be undertaken to implement and validate a more sophisticated parameterization of LNOx production efficiency in chemistry–climate models. The calculated total column LNOx for each grid was distributed into each model layer based on a prescribed “backward C-shaped” LNOx vertical profile (Ott et al., 2010).

2.3 Lightning observation data for model evaluation

We used LIS and OTD gridded climatology datasets for this study, consisting of climatologies of total LFRs observed using the Lightning Imaging Sensor (LIS) and Optical Transient Detector (OTD). The OTD is aboard the MicroLab-1 satellite, and LIS is aboard the Tropical Rainfall Measuring Mission (TRMM) satellite (Cecil et al., 2014). Both sensors detected lightning by monitoring pulses of illumination produced by lightning in the 777.4 nm atomic oxygen multiplet above background levels. In low Earth orbit, both sensors viewed Earth locations for approximately 3 min during the passing of the OTD or for 1.5 min during the passing of the LIS. Each day, OTD and LIS respectively orbited the globe 14 times and 16 times. OTD observed data between +75 and 75 latitude during May 1995–March 2000, whereas LIS collected data between +38 and 38 latitude during January 1998–April 2015. This study uses the LIS–OTD 2.5 Degree Low Resolution Time Series (LRTS), which provides daily LFRs on a 2.5 regular latitude–longitude grid for May 1995–April 2015.

2.4 CMIP6 model outputs for model comparison

For the comparison of different model outputs from our study (CHASER) and other Earth system models or chemistry–climate models, we used LFR and surface temperature data from the CMIP6 CMIP historical experiments from CESM2-WACCM (Danabasoglu, 2019), GISS-E2-1-G (Kelley et al., 2020), and UKESM1-0-LL (Tang et al., 2019). CESM2-WACCM uses the Community Earth System Model ver. 2 (Danabasoglu et al., 2020). The CESM2 is an open-source fully coupled Earth system model. The Whole Atmosphere Community Climate Model ver. 6 (WACCM6) is the atmospheric component coupled to the other components in CESM2. The GISS-E2-1-G is the NASA Goddard Institute for Space Studies (GISS) chemistry–climate model version E2.1 based on the GISS Ocean v1 (G01) model (Miller et al., 2014; Kelley et al., 2020). The UKESM1-0-LL is the UK's Earth system model, details of which were described by Sellar et al. (2019). We used 3 ensembles from CESM2-WACCM, 9 ensembles from GISS-E2-1-G, and 18 ensembles from UKESM1-0-LL. Table S1 in the Supplement presents all the ensemble members used for this study.

2.5 Experiment setup

We conducted six sets of experiments with each set of experiments conducted using both the ECMWF-McCAUL (abbreviated as F1) and the CTH (abbreviated as F2) schemes. Table 1 presents the major settings of all experiments with the relative explanations of those settings. STD-F1/F2 (STD-F1 and STD-F2) experiments are standard experiments with the simulation period of 1959–2014. They are intended to reproduce the historical trends of lightning and LNOx. Climate1959-F1/F2 experiments are experiments that keep the climate simulations fixed to 1959 to derive the effects of global warming on historical lightning trends. ClimateAero1959-F1/F2 experiments are intended to reflect the conditions with climate simulations as well as aerosol and aerosol precursor (BC, OC, NOx, SO2) emissions fixed to 1959. The Aero1959-F1/F2 experiments are the same as the STD-F1/F2 experiments, except for the fact that the AeroPEs are fixed to 1959. The fifth set of experiments (Volca-off-F1/F2) was intended to exclude the influences of the Pinatubo volcanic eruption in order to compare it to STD-F1/F2 and to evaluate the Pinatubo eruption effects on historical lightning and LNOx trends and variation.

Table 1All experiments conducted for this study.

a For the model simulations, the climate is simulated by the prescribed SSTs and sea ice fields and the prescribed varying concentrations of GHGs (CO2; N2O; methane; chlorofluorocarbons, CFCs; and hydrochlorofluorocarbons, HCFCs) used only in the radiation scheme. The SSTs and sea ice fields are obtained from the HadISST dataset (Rayner et al., 2003). The prescribed GHG concentrations are derived from CMIP6 forcing datasets (Meinshausen et al., 2017). b We use “F1” to stand for the ECMWF-McCAUL scheme; “F2” represents the CTH scheme. c Stratospheric aerosol radiative forcing is simulated using the prescribed stratospheric aerosol extinction, which is obtained from the NASA GISS (Sato et al., 1993) and CCMI (Arfeuille et al., 2013) stratospheric aerosol datasets. d The climate is fixed to 1959 for the whole simulation period using the 1959 SSTs and sea ice field and GHG concentrations during the simulation period. e Aerosol (BC, OC) and aerosol precursor (NOx, SO2) emissions (anthropogenic + biomass burning) are fixed to 1959 throughout the simulation period. f Several biogenic emissions are interannually varying, including isoprene, monoterpenes, acetone, and methanol, which were calculated using an offline simulation using the Vegetation Integrative Simulator for Trace Gases (VISIT) terrestrial ecosystem model (Ito and Inatomi, 2012). Some other reactive biogenic VOCs (ethane, propane, ethylene, propene) used are climatological data derived from the Model of Emissions of Gases and Aerosols from Nature (MEGAN) modeling system (Guenther et al., 2012). g Here the June 1991–May 1995 SSTs and sea ice data were replaced with HadISST SSTs and sea ice climatological data during 1985–1990.

Download Print Version | Download XLSX

We simulate volcanic aerosol forcing by considering the prescribed stratospheric aerosol extinction in the radiation scheme. We used the NASA Goddard Institute for Space Studies (GISS) (Sato et al., 1993) and the Chemistry–Climate Model Initiative (CCMI) (Arfeuille et al., 2013) stratospheric aerosol datasets as the stratospheric aerosol climate data. The NASA GISS dataset includes monthly zonal-mean stratospheric aerosol optical thickness in four spectral bands. The CCMI dataset for CHASER includes monthly zonal-mean stratospheric aerosol extinction coefficients in 20 spectral bands. To remove the volcanic perturbation while maintaining the stratospheric background aerosol in the Volca-off-F1/F2, we used the following equation to process the stratospheric aerosol climatology (SAC) during June 1991–May 1996.

(11) SAC no_pinatubo = SAC background , | SAC raw - SAC background | > 1.96 σ , SAC raw , | SAC raw - SAC background | 1.96 σ

In this equation, SACno_pinatubo denotes the stratospheric aerosol climatological data as input data for Volca-off-F1/F2 experiments, and SACbackground represents the stratospheric background aerosol climatological data. (For this study, SACbackground is the corresponding temporally averaged values of the NASA GISS and CCMI stratospheric aerosol datasets during June 1986–May 1991 and June 1996–May 2001, when the time was close to the eruption and the stratosphere was less affected by volcanic eruptions.) SACraw stands for the original values of the NASA GISS and CCMI stratospheric aerosol datasets during June 1991–May 1996. Moreover, σ symbolizes the standard deviations of stratospheric background aerosol climate data. (For this study, σ is the corresponding standard deviations of the NASA GISS and CCMI stratospheric aerosol datasets during June 1986–May 1991 and June 1996–May 2001.) As displayed in Eq. (11), when the absolute differences between SACraw and SACbackground are larger than 1.96σ, we replace the original values (June 1991–May 1996) of the SAC with the temporally averaged values of the NASA GISS and CCMI datasets during June 1986–May 1991 and June 1996–May 2001. When the absolute differences between SACraw and SACbackground are equal to or smaller than 1.96σ, we still use the original values (June 1991–May 1996) of the SAC for the Volca-off experiments. The value of 1.96σ corresponds to the 95 % confidence interval, which can remove the Pinatubo perturbation sufficiently while maintaining the background level of stratospheric aerosol during June 1991–May 1996. Furthermore, the influences of the Pinatubo eruption affected the HadISST SSTs and sea ice fields. To remove the Pinatubo eruption's influences on the SSTs and sea ice fields from the Volca-off experiments too, we replaced the June 1991–May 1995 SSTs and sea ice data with HadISST SSTs and sea ice climatological data during 1985–1990 when conducting the Volca-off experiments. The 1985–1990 period was chosen because it is approximately the period of June 1991–May 1995 and because the SSTs and sea ice fields were less affected by volcanic activity during 1985–1990.

All the experiments calculate the LNOx emissions rates interactively through LNOx emission parameterizations except STD-rVolcaoff experiments. The STD-rVolcaoff experiments are the same as the STD experiments except for reading the daily LNOx emission rates calculated from the Volca-off experiments. The STD-rVolcaoff experiments are conducted for comparison with the STD experiments to elucidate the effects of LNOx emissions changes caused by the Pinatubo eruption on atmospheric chemistry (typically methane lifetime).

3 Results and discussion

3.1 Validation of the simulated historical lightning distribution and trend

To increase the credibility of the conclusions obtained based only on the numerical simulations, the model calculations must be evaluated using observational data. We used the LIS and OTD observations to evaluate the spatial and temporal distribution and historical lightning trends simulated by CHASER (MIROC). Figure 1a–c show the annual mean spatial distributions of lightning observed by LIS and OTD and by model simulations using the ECMWF-McCAUL and CTH schemes. Both the ECMWF-McCAUL and the CTH schemes generally captured the hotspots of lightning (central Africa, the Maritime Continent, South America) well, with strong spatial correlations between observations and model simulations (R>0.75), even if the lightning distributions were not well captured over the ocean. Figure 1d exhibits a strong spatial correlation between observations and simulation results maintained throughout the simulation period (1959–2014).

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f01

Figure 1Annual mean lightning flash densities from (a) LIS and OTD satellite observations spanning 1996–2000, (b) the STD experiment (1960–2014) with the ECMWF-McCAUL scheme used, and (c) the STD experiment (1960–2014) with the CTH scheme used. R and RMSE as shown in the title of panel (b) are calculated between panels (b) and (a), while R and RMSE as shown in the title of panel (c) are calculated between panels (c) and (a). Panel (d) presents the spatial correlation coefficients between the modeled spatial lightning distribution of each year and LIS and OTD lightning climatologies during 1996–2000.

The LIS and OTD observations are also used to evaluate historical lightning trends simulated by CHASER (MIROC). We examined the ±37.5 latitude mean LFR anomaly (1996–2013) calculated from LIS and OTD observations and STD-F1/F2 numerical experiments (Fig. 2 and Table 2). We also note some missing values within the ±37.5 latitude in LIS and OTD observations. To constrain the comparisons between observations and simulations as like-for-like, when we encounter a missing value in the LIS and OTD observations during spatial averaging, we also treat the CHASER-simulated value at the same location as a missing value. As displayed in Fig. 2, we would not necessarily expect that interannual variations in the LFR anomaly can be captured because meteorological nudging was not applied, and the simulated LFRs were only controlled by the prescribed SSTs and sea ice data. Nevertheless, the overall trends in the LFR anomaly simulated using both schemes matched the LIS and OTD observations well, as portrayed in Fig. 2. We further investigated the trends shown in Fig. 2 by the Mann–Kendall rank statistic and Sen's slope estimator, and the statistical summary is displayed in Table 2 (Salmi et al., 2002; Hussain and Mahmud, 2019). Neither the LFR anomaly (within ±37.5 latitude) derived from LIS and OTD observations nor the simulations show a significant trend for 1996–2013 using the Mann–Kendall rank statistic test (significance inferred as 5 %). The global LFR anomaly during 1993–2013 obtained from simulations (STD-F1/F2) also shows no significant trend, which is consistent with the Schumann resonance (SR) intensity observations (1993–2013) at Rhode Island, USA (Williams, 2022). However, the SR observations in Rhode Island (USA) exclude a consideration of the influences of solar cycles, which makes it less appropriate for lightning trend evaluation.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f02

Figure 2LFR anomalies of 1996–2013 within ±37.5 latitude obtained from two numerical experiments (STD-F1/F2) and LIS and OTD satellite observations. The curves represent the monthly time series data of the ±37.5 latitude mean LFR anomalies with the 1-D Gaussian (denoising) filter applied. The lines are the fitting curves of the monthly time series data of the ±37.5 latitude mean LFR anomalies. The trends in the LFR anomalies in percent per year are also presented in the legends.

Download

Table 2A statistical summary of the trends shown in Fig. 2 by the Mann–Kendall rank statistic and Sen's slope estimator. The monthly time series data of the ±37.5 latitude mean LFR anomalies were estimated by the Mann–Kendall rank statistic and Sen's slope estimator. “Trend” shows whether these are significant trends with the significance set as 5 %, as well as the percentage trends in percent per year estimated by linear regression. The p value is calculated during the Mann–Kendall trend test. “Slope” shows Sen's slope of trend. Qmin and Qmax respectively denote the lower and upper limits of the 95 % confidence interval of Sen's slope.

Download Print Version | Download XLSX

We further investigated the seasonal variabilities in the simulated LFR and compared them against LIS and OTD observations. The results are depicted in Fig. 3. Both the CTH and the ECMWF-McCAUL schemes captured the peak during JJA, but the overestimation of LFR by F1/F2 during DJF is also noticeable. Figure S1 in the Supplement presents a comparison of the LFR global distribution in different seasons during 1996–2003 from LIS and OTD lightning observations and STD experiment outputs. Generally, CHASER captured the spatial distribution of LFR in all four seasons well when compared against LIS and OTD observations. The spatial correlation coefficients (R) between observations and simulations are the highest (R=0.80 for both lightning schemes) in DJF, indicating CHASER's considerable capability to reproduce the LFR spatial distribution in DJF. As displayed in the first row of Fig. S1, the overestimation of LFR by F1/F2 during DJF is primarily attributable to the overestimation of LFR within the Maritime Continent and South America, but this might also be attributable to the underestimation of LFR by LIS and OTD within these two regions. It is believed that the LIS and OTD lightning detection efficiency is highly sensitive to the characteristics of convective clouds (cloud albedo, cloud optical thickness, etc.) (Boccippio et al., 2002; Cecil et al., 2014). High cloud albedo and cloud optical thickness might engender the underestimation of LFR by LIS and OTD. It is also noteworthy that the seasonal variation and long-term trend of global lightning are strongly influenced by distinct factors. The seasonal variation in global lightning activities is the most strongly affected by the 23 obliquity of the Earth's orbit and the asymmetric distribution of the continent between the Northern Hemisphere and the Southern Hemisphere. However, the long-term global lightning trend we investigated for this study is mainly controlled by climate forcers such as aerosols and GHGs. To minimize the effects of LFR seasonal variation on our study's results, we deseasonalized the results shown in all figures and tables by calculating their anomaly based on raw data. The validation described above and the deseasonalization of our study's results justified that the LFR seasonal variation (and the uncertainties in the simulation of LFR seasonal variation) in our study has a limited effect on these study results.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f03

Figure 3The mean annual cycle in the area average LFR during 1996–2003. The area average was taken over the grid cells where valid LIS and OTD lightning observations exist. LFR is normalized by min–max normalization.

Download

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f04

Figure 4Monthly time series data of global mean surface temperature anomalies with the 1-D Gaussian (denoising) filter applied and their fitting curves calculated from the outputs of numerical experiments (a, b) and obtained from NCEI (c). Panels (d–f) are the same as panels (a–c), but the averaged surface temperature anomalies are only calculated within the global land regions. The trends of the fitting curves in kelvin per decade are also presented in the legends.

Download

3.2 Effects of global warming and increases in AeroPEs on historical lightning and LNOx trends

As introduced in Sect. 1, global warming and changes in AeroPEs are the two main factors which influence long-term (1960–2014) historical lightning trends. (Hereinafter, historical lightning trends represent lightning trends of 1960–2014.) Evidence shows that the Pacific Decadal Oscillation (PDO) can also affect lightning trends over decadal timescales (Macias Fauria and Johnson, 2006; Mallick et al., 2022), and further research is anticipated to verify it. To analyze the effects of global warming on historical lightning trends, we designed and conducted two sets of experiments: one set of experiments including global warming (STD-F1/F2) and another set of experiments excluding global warming (Climate1959-F1/F2). Figure 4a and b respectively depict the global surface temperature anomalies calculated using the ECMWF-McCAUL and CTH schemes. The STD and Aero1959 experiments show an increasing trend (around 0.11 K per decade) of global mean surface temperature anomalies, which closely approximates the trend (around 0.15 K per decade) obtained from NOAA's National Centers for Environmental Information (NCEI) (2022) (Fig. 4c and f). Global temperature change data from 1880 to the present are available from NCEI, which tracks variations in the Earth's temperature based on thousands of stations' observation data around the globe (Climate at a Glance, National Centers for Environmental Information (NCEI), 2022). When the prescribed SSTs and sea ice fields and GHG concentrations were fixed to 1959 throughout the simulation period, the simulated trends of global mean surface temperature anomalies turned out to be flat (Climate1959 and ClimateAero1959). To elucidate the effects of increases in AeroPEs on averaged surface temperature to the greatest extent possible, we also show the averaged surface temperature anomaly only over land regions (Fig. 4d–f). The simulated global mean land surface temperature anomalies are also well matched with the NCEI observational data. The aerosol cooling effect can be more evident when only examining surface temperature trends averaged over land (Fig. 4d and e).

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f05

Figure 5Panels (a) and (b) show monthly time series data of global mean LFR anomalies with the 1-D Gaussian (denoising) filter applied and their fitting curves of different experiments simulated respectively using the ECMWF-McCAUL and CTH schemes. Panels (c) and (d) are the same as panels (a) and (b), except for the fact that the averaged LFR anomalies are calculated only within global land regions. Trends of the fitting curves (% yr−1) are also shown in the legends.

Download

Figure 5a and b respectively portray the global mean LFR anomalies and their fitting curves obtained from the outputs of the ECMWF-McCAUL and CTH schemes. Moreover, Table 3 shows the statistical summary of the trends in Fig. 5 utilizing the Mann–Kendall rank statistic and Sen's slope estimator. The global lightning trend obtained from the STD-F1 experiment turned out to be statistically flat (0.0 % yr−1), whereas the outputs of the STD-F2 experiment exhibited a significant increasing global lightning trend (0.03 % yr−1) determined using the Mann–Kendall rank statistic (significance inferred as 5 %).

Table 3A statistical summary of the trends shown in Fig. 5 by the Mann–Kendall rank statistic and Sen's slope estimator. The monthly time series data of global or land mean LFR anomalies were estimated by the Mann–Kendall rank statistic and Sen's slope estimator. “Trend” shows whether these are significant trends with the significance set as 5 %, as well as the percentage trends in percent per year estimated by linear regression. The p value is calculated during the Mann–Kendall trend test. “Slope” shows Sen's slope of trend. Qmin and Qmax respectively denote the lower and upper limits of the 95 % confidence interval of Sen's slope.

Download Print Version | Download XLSX

A comparison of the lightning trends calculated from the STD and Climate1959 experiments showed that both lightning schemes demonstrated that historical global warming (1960–2014) enhances the global lightning trends toward positive trends (around 0.03 % yr−1 or 3 % K−1). Global warming effects on historical lightning trends were evaluated as significant using the Mann–Kendall rank statistic, with significance inferred as 5 %, when using the CTH scheme but not in the case of the ECMWF-McCAUL scheme (see rows STD-F1 – Climate1959-F1 (global) and STD-F2 – Climate1959-F2 (global) in Table 3). As shown in Table 3, the differences in global lightning trends simulated by the STD-F1/F2 and Aero1959-F1/F2 experiments indicate that the increases in AeroPEs during 1960–2014 significantly suppress the global lightning trends (0.07 % to 0.04 % yr−1). It is noteworthy that this suppression of lightning trends is only attributable to aerosol radiative effects. Further research must be conducted to elucidate the long-term effects of aerosols on lightning through aerosol microphysical effects. We also investigated lightning trends only over land regions (Fig. 5c and d and Table 3) to ascertain the effects of changes in AeroPEs to the greatest extent possible. When observing the lightning trends over land only, the degree of suppression of lightning trends attributable to increases in AeroPEs expands to 0.10 % to 0.05 % yr−1, which is attributable to most AeroPEs and their growth coming from land regions. It is noteworthy that we used the same SSTs and sea ice data in the Aero1959 as those used for STD experiments. The SSTs and sea ice data also reflected the effects of increases in AeroPEs. Therefore, we might underestimate the effects of increases in AeroPEs on lightning trends by comparing the results of STD and Aero1959 experiments.

For the ECMWF-McCAUL scheme, model outputs affirm that global warming can enhance the global mean CAPE anomaly slightly and suppress the global mean QRa anomaly (Fig. 6a and b). Earlier studies have also indicated that the total solid (cloud ice, snow, and graupel) mass mixing ratio within charge separation regions is lower under global warming. Moreover, possible explanations are given in those studies (Finney et al., 2018; Romps, 2019). Because global warming enhances global convection activities and because lightning formation is highly related to convection activity, global warming enhances the historical global lightning trend simulated using the ECMWF-McCAUL scheme, mainly as a result of the simulated CAPE trend, which is enhanced by global warming.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f06

Figure 6Panels (a) and (b) respectively show monthly time series data of the global mean CAPE and QRa anomalies with the 1-D Gaussian (denoising) filter applied and their fitting curves simulated using the ECMWF-McCAUL scheme. Panels (c) and (d) respectively show differences in the CAPE anomaly trend (Jkg-1yr-1) and QRa anomaly trend (kgm-2yr-1) of the STD-F1 and Aero1959-F1 experiments in the global map. Panels (e) and (f) respectively show monthly time series data of the ±35 latitude land region mean CAPE and QRa anomalies with the 1-D Gaussian (denoising) filter applied and their fitting curves simulated using the ECMWF-McCAUL scheme. Panel (g) portrays the vertical profiles of the trend of the global mean cumulus cloud fraction anomaly simulated by the CTH scheme. Panel (h) depicts the relative contributions of each layer's cumulus to total lightning density in 1960 and 2014, as calculated from the outputs of the STD-F2 experiment.

The past increases in AeroPEs exert negligible effects on the trends of global mean CAPE and QRa anomalies, as displayed in Fig. 6a and b. However, as also demonstrated in our study (see Fig. 1), most lightning flashes occur over tropical and subtropical land regions. It is displayed in Fig. 6c and d that the past increases in AeroPEs mostly suppress the CAPE and QRa absolute trends within regions with high lightning densities. We further investigated the trends of the ±35 latitude land region mean CAPE and QRa anomalies, and the results are portrayed in Fig. 6e and f. Figure 6e and f show that past increases in AeroPEs significantly suppress the QRa trend (0.08 % yr−1) and slightly suppress the CAPE trend (0.03 % yr−1) within ±35 latitude land regions. Weaker convection activities (smaller CAPE) and fewer hydrometeors (cloud ice, graupel, snow) in the charge separation regions (0 to 25 C isotherm) engender less lightning. In the case of the ECMWF-McCAUL scheme, CAPE and QRa trends were suppressed within ±35 latitude terrestrial regions. This constitutes the main reason for the suppression of the historical global lightning trends induced by increases in AeroPEs through aerosol radiative effects. It is noteworthy that, because the aerosol microphysical effects are only considered in the grid-scale large-scale condensation scheme, our study might underestimate the aerosol microphysical effects which can enhance the trends of QRa and LFR toward a positive direction.

To explain the results simulated by the CTH scheme, we investigated the vertical profiles of the trend of the global mean cumulus cloud fraction anomaly (Fig. 6g). Investigating cumulus cloud fraction is reasonable because each model layer's cumulus cloud fractions are used to weight the calculated lightning densities from that layer in the CTH scheme, as introduced in Eqs. (3) and (4). Figure 6h shows the relative contributions of each model layer's cumulus to the calculated global total lightning densities in 1960 and 2014 obtained using the CTH scheme. As is shown in Fig. 6h, the vertical profiles of the relative contribution in 1960 and 2014 are almost identical. Cumulus convection is positively correlated with lightning formation, which is the scientific basis of parameterizing lightning densities using the cumulus cloud top height: the CTH scheme. Historical global warming enhances the lightning trend simulated by the CTH scheme mainly because the simulated historical global warming increases the cumulus reaching 200 hPa, which contributes greatly to the simulated global total lightning density (Fig. 6g and h). The increases in the deep convective cloud are regarded as related to the increases in tropopause height attributable to global warming, as shown in Fig. S2 in the Supplement. The past increases in AeroPEs suppress the lightning trend simulated by the CTH scheme because increases in AeroPEs decrease the cumulus reaching 200 hPa as well as the cumulus within the lower to middle troposphere by aerosol radiative effects (Fig. 6g). In addition, in the Supplement, we present a figure (Fig. S3) resembling Fig. 6 but which includes only a consideration of land regions. The mechanisms of global warming and increases in AeroPEs affecting lightning trends over land regions are similar to those described above on a global scale. We do not discuss details of them here.

We also investigated lightning trends simulated in different experiments with the global map (Fig. 7). Both the ECMWF-McCAUL and the CTH schemes show that lightning increased significantly in most parts of the Arctic region and decreased in some parts of the Southern Ocean during 1960–2014 (Fig. 7a and e). The significant lightning trends presented in Fig. 7a became nearly nonexistent when the climate simulations were fixed to 1959 (Fig. 7b and f), indicating the considerable effects of global warming on the trend of global lightning activities. Furthermore, the effects of past global warming and increases in AeroPEs on the lightning trends on the global map are displayed in Fig. 8. Figure 8a and c show that past global warming enhances lightning activities within the Arctic region and Japan, which is consistent with findings of an earlier study from which thunder day data in Japan were reported (Fujibe, 2017). Figure 8a and c also show that historical global warming suppresses lightning activities around New Zealand and some parts of the Southern Ocean. Both lightning schemes demonstrated that the historical increases in AeroPEs suppress lightning activities in some parts of the Southern Ocean and South America. The ECMWF-McCAUL scheme also suggests that historical increases in AeroPEs suppress lightning activities by aerosol radiative effects in some parts of India and China, where AeroPEs increased dramatically during 1960–2014 because of rapid economic development and energy consumption. Many observation-based studies indicate that aerosols can invigorate lightning activities in some regions of China and India, typically under relatively clean conditions (e.g., AOD<1.0), which is attributable to the aerosol microphysical effects (Wang et al., 2011; Zhao et al., 2017, 2020; Lal et al., 2018; Liu et al., 2020; Shi et al., 2020). Therefore, a total positive effect of aerosol on historical lightning trends in China and India cannot be ruled out. We further provided the same figures as Figs. 7 and 8 but using different units (fl.km-2yr-2) in Figs. S4 and S5 in the Supplement. Figures S4 and S5 show that the absolute lightning trends (fl.km-2yr-2) and the effects of global warming and increases in AeroPEs on the absolute lightning trends are slight in high-latitude regions but prominent in tropical areas.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f07

Figure 7Trends of the LFR anomaly (% yr−1) during 1960–2014 on the 2-D map. The trend at every point was calculated from the function of the approximating curve for the 1960–2014 time series data (LFR anomaly) at each grid cell. The area in which the trend was found to be significant by the Mann–Kendall rank statistic test (significance inferred as 5 %) is marked with the hatched lines.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f08

Figure 8Differences in trends of the LFR anomaly during 1960–2014 on the global map. The area in which the trend of the differences in the LFR anomaly time series data was found to be significant by the Mann–Kendall rank statistic test (significance inferred as 5 %) is marked with the hatched lines.

Trends in historical annual global LNOx emissions for different scenarios are generally consistent with trends in historical global mean LFRs, as shown in Figs. 5a and b and 9. This finding is not surprising because, as the lightning NOx emission parameterizations introduced in Sect. 2.2 show, the simulated LFRs are linearly related to the simulated LNOx emissions in our study. A comparison of the LNOx trends calculated from the STD and Climate1959 experiments showed that both lightning schemes demonstrated that historical global warming (1960–2014) enhances the global LNOx trends toward positive trends (0.02 %–0.04 % yr−1). Global warming effects on historical LNOx trends were evaluated as significant using the Mann–Kendall rank statistic, with significance inferred as 5 %, when using the CTH scheme but not in the case of the ECMWF-McCAUL scheme (see rows STD-F1 – Climate1959-F1 and STD-F2 – Climate1959-F2 in Table 4). As shown in Table 4, the differences in global LNOx trends simulated by the STD and Aero1959 experiments indicate that the increases in AeroPEs during 1960–2014 significantly suppress the global LNOx trends (0.08 % to 0.03 % yr−1). The results presented in Fig. 9 and Table 4 imply that historical global warming and increases in AeroPEs can affect atmospheric chemistry and engender feedback by influencing LNOx emissions.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f09

Figure 9Time series data of the 1960–2014 annual global LNOx production anomalies (Tg N yr−1) and their fitting curves simulated using the ECMWF-McCAUL (a) and CTH schemes (b). Trends of the fitting curves in percent per year are presented in the legends.

Download

Table 4A statistical summary of the trends shown in Fig. 9 by the Mann–Kendall rank statistic and Sen's slope estimator. The time series data of the annual global LNOx production anomalies were estimated by the Mann–Kendall rank statistic and Sen's slope estimator. “Trend” shows whether these are significant trends with the significance set as 5 %, as well as the percentage trends in percent per year estimated by linear regression. The p value is calculated during the Mann–Kendall trend test. “Slope” shows Sen's slope of trend. Qmin and Qmax respectively denote the lower and upper limits of the 95 % confidence interval of Sen's slope.

Download Print Version | Download XLSX

3.3 Pinatubo volcanic eruption effects on historical lightning and LNOx trends

We estimate the Pinatubo eruption effects on historical lightning and LNOx trends and variation by comparing the simulation results of the STD and Volca-off experiments. The simulated global mean LFRs by the STD and Volca-off experiments are the same until April 1991. They then begin to show differences from May 1991 (the time series of the global mean LFRs is not shown). This result is reasonable because the Pinatubo volcanic perturbations are removed from the SAC during June 1991–May 1996 in the Volca-off experiments in Eq. (11) and because the SAC of May 1991 used in CHASER is interpolated between the SAC of April 1991 and June 1991.

Figure 10c and d portray the time series of the LFR anomalies and Relative_diff (values over the red lines) during 1991–1992. Relative_diff is the relative difference in the global mean LFR anomalies between the STD and Volca-off experiments calculated using the following equation.

(12) Relative_diff = 100 % × LFRA STD - LFRA Volca-off LFR Volca-off

In the equation, LFRASTD represents the global mean LFR anomalies simulated by the STD-F1/F2 experiments. LFRAVolca-off denotes the global mean LFR anomalies simulated by the Volca-off-F1/F2 experiments. LFRVolca-off symbolizes the global mean LFRs simulated by the Volca-off-F1/F2 experiments.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f10

Figure 10Time series of the LFR anomalies during 1990–1999 and 1991–1992. Panels (a) and (b) show the time series of the LFR anomalies and their smoothed curves by the 1-D Gaussian (denoising) filter for 1990–1999. Panels (c) and (d) present the time series of the LFR anomalies during 1991–1992. The values shown over the red lines in panels (c) and (d) are Relative_diff calculated using Eq. (12).

Download

The monthly time series data of Relative_diff for 1990–1999 for both lightning schemes are calculated. The probability density distributions (PDDs) of Relative_diff spanning 1990–1999 and June 1991–May 1992 are displayed in Fig. 11. The 1990–1999 Relative_diff values presented in Fig. 11 (colored blue) are all normally distributed as determined by the Kolmogorov–Smirnov test. The 95 % confidence interval of the 1990–1999 Relative_diff is calculated and shown in the titles of Fig. 11. As displayed in Fig. 10c and d, the underlined values (Relative_diff) exceeded the 95 % confidence interval, indicating significant differences in the calculated global mean LFR anomalies by the STD and Volca-off experiments. In other words, global lightning activities were suppressed significantly by the Pinatubo eruption during the first year after the eruption. The PDDs of the June 1991–May 1992 Relative_diff (colored green in Fig. 11) shifted to the left compared to the 1990–1999 PDDs, indicating that global lightning activities were suppressed in the first year after the eruption.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f11

Figure 11Probability density distributions (PDDs) of Relative_diff obtained from monthly time series data of Relative_diff during 1990–1999 and June 1991–May 1992 (a year after the Pinatubo eruption). The 1990–1999 Relative_diff for both lighting schemes is normally distributed, and N(μ,σ2) is displayed in the titles of this figure. The 95 % confidence interval of the 1990–1999 Relative_diff is also shown in the titles of this figure.

Download

Figure 12a and b show August 1991–April 1992 averaged LFR anomaly differences between the STD and Volca-off experiments on the global map. We found from Fig. 12a and b that lightning activities are suppressed significantly within the three hotspots of lightning activities (central Africa, the Maritime Continent, and South America) during August 1991–April 1992, when the global mean LFRs are found to be suppressed. To elucidate the potential reasons for the suppressed global lightning activities during the first year after the Pinatubo eruption, we first investigated the August 1991–April 1992 averaged differences in the CAPE and QRa anomalies between STD-F1 and Volca-off-F1 (Fig. 12c and d) because lightning densities are computed with CAPE and QRa by the ECMWF-McCAUL scheme. The results showed that the Pinatubo eruption can engender apparent reductions in CAPE and QRa within tropical and subtropical terrestrial regions (typically three hotspots of lightning activities) where lightning occurrence is frequent. These reductions constitute the main reason for the suppressed global lightning activities during the first year after the Pinatubo eruption simulated by the ECMWF-McCAUL scheme. We also examined the August 1991–April 1992 averaged differences of the 200–400 hPa averaged cumulus cloud fraction anomaly between STD-F2 and Volca-off-F2 on the global map (Fig. 12e). The cumulus cloud fractions of each model layer are used to weight the calculated lightning densities from that layer by the CTH scheme, as explained in Sect. 2.2. As depicted in Figs. 12e and S6 in the Supplement, the Pinatubo eruption led to marked reductions in the middle to upper tropospheric cumulus cloud fractions during August 1991–April 1992 over three hotspots of lightning activities (central Africa, the Maritime Continent, and South America). As displayed in Fig. 6h, the cumulus that reached the middle to upper troposphere is closely related to lightning formation. Consequently, the simulated global lightning activities by the CTH scheme were also suppressed considerably during the first year after the Pinatubo eruption.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f12

Figure 12The August 1991–April 1992 averaged LFR anomaly differences (a, b), CAPE anomaly differences (c), QRa anomaly differences (d), and differences of the 200–400 hPa averaged cumulus cloud fraction anomaly between the STD-F2 and Volca-off-F2 experiments (e) on the global map.

Aside from the global lightning activity suppression described earlier, the production of LNOx might also decrease after the Pinatubo eruption. To explore this conjecture, we compared the LNOx emissions in the STD and Volca-off experiments (Fig. 13). In the case of the ECMWF-McCAUL scheme, the reduction in LNOx emissions caused by the Pinatubo eruption started in 1991 (5.84 %) and continued until 1993, with the highest percentage reduction occurring in 1992 (8.47 %) (Fig. 13a). However, the CTH scheme showed a slightly different scenario of LNOx emissions reduction after the Pinatubo eruption. The LNOx emissions are almost evenly reduced during 1991–1994 in the case of the CTH scheme (Fig. 13b). In conclusion, our study indicates that the Pinatubo eruption can engender reductions in global LNOx emissions, which last 2–3 years. However, there exists some uncertainty in evaluating the magnitude of the reductions: from 1.99 % to 8.47 % for the annual percentage reduction found in our study.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f13

Figure 13The 1990–1999 annual global LNOx emissions calculated from the STD and Volca-off experiments' outputs simulated using the ECMWF-McCAUL (a) and CTH schemes (b). Values over the red lines represent the relative differences (%) between the red and blue lines, calculated with respect to the blue lines.

Download

The simulated reduced global LNOx emissions caused by the Pinatubo eruption might influence atmospheric chemistry significantly. Most importantly, the reduced global LNOx emissions might reduce OH radical production and extend the global mean tropospheric lifetime of methane against tropospheric OH radicals, abbreviated hereinafter as the methane lifetime. We investigated this point further by comparing the methane lifetime anomaly simulated by the STD and STD-rVolcaoff experiments. As introduced in Sect. 2.5, the settings of the STD-rVolcaoff experiments are the same as those use for the STD experiments, except that they use the daily LNOx emission rates calculated from the Volca-off experiments. We calculated the monthly CH4 lifetime anomalies during 1990–1999 and ΔτCH4 (the difference in the CH4 lifetime anomaly between the STD and STD-rVolcaoff experiments), which are shown in Fig. 14c and d. Figure 14a and b display the PDDs of ΔτCH4 monthly time series during 1990–1999. The ΔτCH4 values shown in Fig. 14a and b are all normally distributed, as determined using the Kolmogorov–Smirnov test. The 95 % confidence interval of ΔτCH4 is calculated and shown in the titles of Fig. 14a and b. The annual global LNOx production averaged during 1990–1999 is 3.56 Tg N yr−1 for STD-F1 and 4.79 Tg N yr−1 for STD-F2. At this level of annual global LNOx production, we found that within the first 2 years after the Pinatubo eruption, ΔτCH4 exceeded the 95 % confidence interval simulated by both lighting schemes (February 1992 and April 1992 in the case of the ECMWF-McCAUL scheme; December 1991 in the case of the CTH scheme). However, the widely cited range of annual global LNOx production is 2–8 Tg N yr−1 (Schumann and Huntrieser, 2007). Presuming that ΔτCH4 responds linearly to the LNOx emission level and that the annual global LNOx production is 8 Tg N yr−1, then the extension of the CH4 lifetime because of the reduced LNOx emissions can reach around 0.54 years for the ECMWF-McCAUL scheme. As a comparison, ultraviolet shielding effects caused by stratospheric aerosols after the Pinatubo eruption led to the maximum increase in the methane lifetime by about 0.6 years (Fig. 14c and d).

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f14

Figure 14Panels (a) and (b) show the probability density distributions (PDDs) of ΔτCH4 obtained from the monthly time series data of ΔτCH4 during 1990–1999. ΔτCH4 represents the difference in the CH4 lifetime anomaly between the STD and STD-rVolcaoff experiments. The 95 % confidence interval of ΔτCH4 is also presented in the titles of panels (a) and (b). Panels (c) and (d) show monthly time series of the CH4 lifetime anomalies simulated by the STD-F1/F2 and STD-rVolcaoff-F1/F2 experiments. The values over the red lines represent ΔτCH4.

Download

3.4 Model intercomparisons of LFR trends with CMIP6 model outputs

The historical lightning trends demonstrated in our study are undoubtedly worth comparing with the results of other chemistry–climate models or Earth system models. As introduced in Sect. 2.4, for comparison of the simulated LFR trends and variations in our study with those of other CMIP6 models' outputs, we used all available LFR data from the CMIP6 CMIP historical experiments from CESM2-WACCM (3 ensembles) (Danabasoglu, 2019), GISS-E2-1-G (9 ensembles) (Kelley et al., 2020), and UKESM1-0-LL (18 ensembles) (Tang et al., 2019). Table S1 presents a complete list of the ensemble members we used. It is noteworthy that the LFR data obtained from the three CMIP6 models described earlier are calculated using the CTH scheme. The results of model intercomparisons of LFR trends and variations are displayed in Fig. 15.

https://acp.copernicus.org/articles/23/13061/2023/acp-23-13061-2023-f15

Figure 15Comparisons of simulated global mean LFR anomalies found in our study (CHASER) and found using other CMIP6 models. All the figures are created based on the monthly time series data of global mean LFR anomalies with a 1-D Gaussian (Denoising) filter applied. For CMIP6 models, the ensemble mean is shown as the solid line, and the full ensemble range is shown as grey shading (c–e). Fitting curves and the trends of fitting curves (% yr−1) are also given in (a and b).

Download

As displayed in Fig. 15a and b and Table 6, both the ECMWF-McCAUL and the CTH schemes (STD-F1/F2) simulated almost flat global lightning trends (even the trend is estimated to be significant in the case of the CTH scheme (0.03 % yr−1)), but the ensemble means obtained from another three CMIP6 models exhibit much larger significant increasing global lightning trends (trends from 0.11 % to 0.25 % yr−1). Many reasons underlie the differences in the global lightning trends simulated by CHASER in our study and by the three CMIP6 models, including the use of different methods to determine SSTs and sea ice fields. Instead of using a coupled atmosphere–ocean general circulation model to calculate SSTs and sea ice fields dynamically in the three CMIP6 models, CHASER uses the prescribed HadISST data (Rayner et al., 2003), which are based on plenty of observational data. Changes in the global mean sea surface temperature anomaly during 1960–2014 (ΔSST) obtained from the STD-F1/F2 and CMIP6 model outputs are presented in Table 5. We also used the observation-based Extended Reconstructed SST (ERSST) dataset (Huang et al., 2017) constructed by NOAA to evaluate the ΔSST obtained from different models. The ΔSST calculated from ERSST during 1960–2014 is 0.549 C, which most closely approximates the ΔSST obtained from STD-F1/F2. Considered from the perspective of SSTs and sea ice fields alone, the results (global lightning trends) of our study are expected to be closer to the actual situation.

Table 5Changes in the global mean surface temperature anomaly (ΔTS), global mean sea surface temperature anomaly (ΔSST), and global mean lightning flash rate anomaly (ΔLFR), as well as the rate of change in the LFR anomaly corresponding to each degree Celsius increase in the global mean surface temperature anomaly (ΔLFR/ΔTS) obtained from the STD-F1/F2 and CMIP6 model outputs. The change in ΔSST obtained from the ERSST dataset is also shown in this table. Changes were obtained by calculating the difference between the rightmost and leftmost points of the approximating curve for the 1960–2014 time series data.

Download Print Version | Download XLSX

Actually, the three CMIP6 models simulated stronger global warming during 1960–2014 than CHASER in our study, as displayed in Fig. S7 in the Supplement and Table 5. The CTH scheme is reported to respond positively to simulated global warming (Price and Rind, 1994; Zeng et al., 2008; Hui and Hong, 2013; Banerjee et al., 2014; Krause et al., 2014; Clark et al., 2017). The simulated stronger global warming by the three CMIP6 models is regarded as being responsible for differences in the simulated global lightning trends between our study and the three CMIP6 models (Fig. 15a and b and Table 6). We further investigated the sensitivities of the global mean LFR anomaly change to the global mean surface temperature anomaly increase (%C-1) obtained from CHASER and the three CMIP6 models. The sensitivities in percentage per degrees Celsius are presented in Table 5. Overall, even when using the same CTH scheme, the sensitivities (ΔLFR/ΔTS) simulated by the three CMIP6 models are higher than those simulated by CHASER in our study. This different sensitivity might be partially attributable to the nonlinear relation between lightning response and climate change (Pinto, 2013; Krause et al., 2014). Compared to the CTH scheme, the ECMWF-McCAUL scheme simulated a statistically non-significant negative sensitivity (ΔLFR/ΔTS), which is attributable to the stronger suppression of positive global lightning trends caused by increases in AeroPEs simulated using the ECMWF-McCAUL scheme.

Table 6A statistical summary of the trends shown in Fig. 15a and b by Mann–Kendall rank statistic and Sen's slope estimator. The time series data of global mean LFR anomalies were estimated by Mann–Kendall rank statistic and Sen's slope estimator. “Trend” shows whether these are significant trends with the significance set as 5 %, as well as the percentage trends in % yr−1 estimated by linear regression. The p value is calculated during Mann–Kendall trend test. “Slope” shows Sen's slope of trend. Qmin and Qmax respectively denote the lower and upper limits of the 95 % confidence interval of Sen's slope.

Download Print Version | Download XLSX

Figure 15d and e affirm that the global lightning variation simulated by our study is basically within the full ensemble range of GISS-E2-1-G and UKESM1-0-LL. After the Pinatubo eruption, as described in Sect. 3.3 of this paper, the GISS-E2-1-G and UKESM1-0-LL models also manifest significant suppression of global lightning activities, but the CESM2-WACCM model shows no such phenomenon. The commonalities and differences in the global lightning trends and variations found in the model intercomparisons imply that great uncertainties existed in past (1960–2014) global lightning trend simulations. Such uncertainties deserve to be investigated further.

4 Discussion and conclusions

We used two lightning schemes (the CTH and ECMWF-McCAUL schemes) to study historical (1960–2014) lightning and LNOx trends and variations and their influencing factors (global warming, increases in AeroPEs, and the Pinatubo eruption) within the CHASER (MIROC) chemistry–climate model. The CTH scheme, which is the most widely used lightning scheme, nevertheless lacks a direct physical link with the charging mechanism. The ECMWF-McCAUL scheme is a newly developed process-based and ice-based lightning scheme with a direct physical link to the charging mechanism.

With only the aerosol radiative effects considered in the lightning–aerosol interaction, both lightning schemes simulated almost flat trends of global mean LFR during 1960–2014 (no trend is detected in the case of the ECMWF-McCAUL scheme, but a slightly significant increasing trend is detected in the case of the CTH scheme). Reportedly, because the aerosol microphysical effects can enhance lightning activities (Yuan et al., 2011; Wang et al., 2018; Liu et al., 2020), our study might underestimate the increasing trend of global mean LFR (our study only considered the aerosol radiative effects in aerosol–lightning interactions). Further research is anticipated, with consideration of the effects of aerosol microphysical effects on long-term lightning trends. Moreover, both lightning schemes manifest that past global warming enhances the historical trend of the global mean lightning density toward a positive direction (around 0.03 % yr−1 or 3 % K−1). However, past increases in AeroPEs exert the opposite effect on the lightning trend (0.07 % to 0.04 % yr−1). The effects of the increased AeroPEs on the lightning trend only over land regions expand to 0.10 % to 0.05 % yr−1, which implies that the effects are more significant over land regions. We obtained similar results for the historical global LNOx emissions trend, which indicates that historical global warming and increases in AeroPEs can affect atmospheric chemistry and engender feedback by influencing LNOx emissions. Although the CTH and ECMWF-McCAUL schemes use different parameters to simulate lightning, both lightning schemes indicate that the enhanced global convective activity under global warming is the main reason for the increase in lightning and LNOx emissions. In contrast, the increases in AeroPEs have decreased lightning and LNOx emissions by weakening the convective activity in the lightning hotspots. By analyzing the simulation results on the global map, we also found that the effects of historical global warming and increases in AeroPEs on lightning trends are heterogeneous across different regions. Our results indicate that historical global warming enhances lightning activities within the Arctic region and Japan but suppresses lightning activities around New Zealand and some parts of the Southern Ocean. Both lightning schemes demonstrated that the historical increases in AeroPEs suppress lightning activities in some parts of the Southern Ocean and South America. The ECMWF-McCAUL scheme also suggests that historical increases in AeroPEs suppress lightning activities in some parts of India and China when only the aerosol radiative effects are considered. This finding is plausible because both countries experienced dramatic increases in AeroPEs during 1960–2014 because of rapid economic growth.

Furthermore, this paper is the first describing significant suppression of global lightning activity during the first year after the Pinatubo eruption, which is indicated in both lightning schemes (global lightning activities decreased by up to 18.10 % simulated by the ECMWF-McCAUL scheme). This finding is mainly attributable to the Pinatubo eruption weakening of the convective activities within the hotspots of lightning, which in turn decreased QRa and middle-level to high-level cumulus cloud fractions in these regions. The simulation results also indicate that the Pinatubo eruption can engender reductions in global LNOx emissions, which last 2–3 years. However, some uncertainty exists in evaluating the magnitude of these reductions (from 1.99 % to 8.47 % for the annual percentage reduction in our study). The case study of the Pinatubo eruption in our research indicates that other large-scale volcanic eruptions can also engender significant reduction in global lightning activities and global-scale LNOx emissions.

Lastly, we compared the global lightning trends demonstrated in our study with the outputs of three CMIP6 models: CESM2-WACCM, GISS-E2-1-G, and UKESM1-0-LL. We used all available LFR data from the CMIP6 CMIP historical experiments from the three models described above. The three CMIP6 models suggest significant increasing trends in historical global lightning activities, which differs from the findings of our study with respect to the magnitude of lightning trends. Unlike the three CMIP6 models which use a coupled atmosphere–ocean general circulation model to calculate SSTs and sea ice fields dynamically, our study (CHASER) uses the prescribed HadISST SSTs and sea ice data, which more closely reflect the actual situation. Therefore, we believe that the results (the historical global lightning trends) obtained from our study (CHASER) more closely approximate the actual situation. However, model intercomparisons of global lightning trends still indicate that considerable uncertainties exist in historical (1960–2014) global lightning trend simulations and that such uncertainties deserve further investigation.

Code availability

The source code for CHASER to reproduce results obtained from this work is obtainable from the repository at https://doi.org/10.5281/zenodo.5835796 (He et al., 2022a).

Data availability

The LIS and OTD data used for this study are available from https://search.earthdata.nasa.gov/search?portal=ghrc (NASA, 2023). The CMIP6 model outputs (LFR and surface temperature) used for this study are available from https://aims2.llnl.gov/search (ESGF, 2023). The Extended Reconstructed SST data used for this study are available from https://www.ncei.noaa.gov/products/extended-reconstructed-sst (NCEI, 2023).

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/acp-23-13061-2023-supplement.

Author contributions

YH conducted all simulations, interpreted the results, and wrote the paper. KS developed the CHASER (MIROC) model code, conceived the presented idea, and supervised the findings of this work and the paper preparation.

Competing interests

The contact author has declared that neither of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

The authors would like to take this opportunity to thank the Interdisciplinary Frontier Next-Generation Researcher Program of the Tokai Higher Education and Research System. The simulations were completed using the supercomputer (NEC SX-Aurora TSUBASA) at NIES (Japan). We thank NASA scientists and staff for providing LIS and OTD lightning observation data. We acknowledge the World Climate Research Programme, which coordinated and promoted CMIP6 through its Working Group on Coupled Modelling. We extend our sincere gratitude to the climate modeling groups for producing and providing their model outputs, the Earth System Grid Federation (ESGF) for archiving the data and providing free downloads, and the multiple funding agencies that have supported CMIP6 and the Earth System Grid Federation. We also thank Do Thi Nhu Ngoc for her assistance in downloading the CMIP6 model outputs.

Financial support

This research has been supported by the Ministry of the Environment, government of Japan (grant nos. S–12 and S–20), the Japan Society for the Promotion of Science (grant nos. JP20H04320, JP19H05669, JP23H04971, and JP19H04235), and the Japan Science and Technology Agency (grant no. JPMJSP2125).

Review statement

This paper was edited by Kostas Tsigaridis and reviewed by three anonymous referees.

References

Allen, D. J., Pickering, K. E., Bucsela, E., Krotkov, N., and Holzworth, R.: Lightning NOx Production in the Tropics as Determined Using OMI NO2 Retrievals and WWLLN Stroke Data, J. Geophys. Res.-Atmos., 124, 13498–13518, https://doi.org/10.1029/2018JD029824, 2019. 

Altaratz, O., Kucienska, B., Kostinski, A., Raga, G. B., and Koren, I.: Global association of aerosol with flash density of intense lightning, Environ. Res. Lett., 12, 114037, https://doi.org/10.1088/1748-9326/aa922b, 2017. 

Arfeuille, F., Luo, B. P., Heckendorn, P., Weisenstein, D., Sheng, J. X., Rozanov, E., Schraner, M., Brönnimann, S., Thomason, L. W., and Peter, T.: Modeling the stratospheric warming following the Mt. Pinatubo eruption: uncertainties in aerosol extinctions, Atmos. Chem. Phys., 13, 11221–11234, https://doi.org/10.5194/acp-13-11221-2013, 2013. 

Banerjee, A., Archibald, A. T., Maycock, A. C., Telford, P., Abraham, N. L., Yang, X., Braesicke, P., and Pyle, J. A.: Lightning NOx, a key chemistry–climate interaction: impacts of future climate change and consequences for tropospheric oxidising capacity, Atmos. Chem. Phys., 14, 9871–9881, https://doi.org/10.5194/acp-14-9871-2014, 2014. 

Boccippio, D. J., Koshak, W. J., and Blakeslee, R. J.: Performance Assessment of the Optical Transient Detector and Lightning Imaging Sensor. Part I: Predicted Diurnal Variability, J. Atmos. Ocean. Tech., 19, 1318–1332, https://doi.org/10.1175/1520-0426(2002)019<1318:PAOTOT>2.0.CO;2, 2002. 

Boucher, O.: Atmospheric Aerosols, Springer Netherlands, Dordrecht, https://doi.org/10.1007/978-94-017-9649-1, 2015. 

Bucsela, E. J., Pickering, K. E., Allen, D. J., Holzworth, R. H., and Krotkov, N. A.: Midlatitude Lightning NOx Production Efficiency Inferred From OMI and WWLLN Data, J. Geophys. Res.-Atmos., 124, 13475–13497, https://doi.org/10.1029/2019JD030561, 2019. 

Cecil, D. J., Buechler, D. E., and Blakeslee, R. J.: Gridded lightning climatology from TRMM-LIS and OTD: Dataset description, Atmos. Res., 135–136, 404–414, https://doi.org/10.1016/j.atmosres.2012.06.028, 2014. 

Cerveny, R. S., Bessemoulin, P., Burt, C. C., Cooper, M. A., Cunjie, Z., Dewan, A., Finch, J., Holle, R. L., Kalkstein, L., Kruger, A., Lee, T., Martínez, R., Mohapatra, M., Pattanaik, D. R., Peterson, T. C., Sheridan, S., Trewin, B., Tait, A., and Wahab, M. M. A.: WMO Assessment of Weather and Climate Mortality Extremes: Lightning, Tropical Cyclones, Tornadoes, and Hail, Weather Clim. Soc., 9, 487–497, https://doi.org/10.1175/WCAS-D-16-0120.1, 2017. 

Clark, S. K., Ward, D. S., and Mahowald, N. M.: Parameterization-based uncertainty in future lightning flash density, Geophys. Res. Lett., 44, 2893–2901, https://doi.org/10.1002/2017GL073017, 2017. 

Cooper, M. A. and Holle, R. L.: Current Global Estimates of Lightning Fatalities and Injuries, in: Reducing Lightning Injuries Worldwide, edited by: Cooper, M. A. and Holle, R. L., Springer International Publishing, Cham, 65–73, https://doi.org/10.1007/978-3-319-77563-0_6, 2019. 

Cooray, V., Rahman, M., and Rakov, V.: On the NOx production by laboratory electrical discharges and lightning, J. Atmos. Sol.-Terr. Phy., 71, 1877–1889, https://doi.org/10.1016/j.jastp.2009.07.009, 2009. 

Danabasoglu, G.: NCAR CESM2-WACCM model output prepared for CMIP6 CMIP historical, https://doi.org/10.22033/ESGF/CMIP6.10071, 2019. 

Danabasoglu, G., Lamarque, J.-F., Bacmeister, J., Bailey, D. A., DuVivier, A. K., Edwards, J., Emmons, L. K., Fasullo, J., Garcia, R., Gettelman, A., Hannay, C., Holland, M. M., Large, W. G., Lauritzen, P. H., Lawrence, D. M., Lenaerts, J. T. M., Lindsay, K., Lipscomb, W. H., Mills, M. J., Neale, R., Oleson, K. W., Otto-Bliesner, B., Phillips, A. S., Sacks, W., Tilmes, S., van Kampenhout, L., Vertenstein, M., Bertini, A., Dennis, J., Deser, C., Fischer, C., Fox-Kemper, B., Kay, J. E., Kinnison, D., Kushner, P. J., Larson, V. E., Long, M. C., Mickelson, S., Moore, J. K., Nienhouse, E., Polvani, L., Rasch, P. J., and Strand, W. G.: The Community Earth System Model Version 2 (CESM2), J. Adv. Model. Earth Sy., 12, e2019MS001916, https://doi.org/10.1029/2019MS001916, 2020. 

Del Genio, A. D., Yao, M.-S., and Jonas, J.: Will moist convection be stronger in a warmer climate?, Geophys. Res. Lett., 34, L16703, https://doi.org/10.1029/2007GL030525, 2007. 

ESGF (Earth System Grid Federation): ESGF MetaGrid [data set], https://aims2.llnl.gov/search (last access: 1 February 2023), 2023. 

Finney, D. L., Doherty, R. M., Wild, O., Huntrieser, H., Pumphrey, H. C., and Blyth, A. M.: Using cloud ice flux to parametrise large-scale lightning, Atmos. Chem. Phys., 14, 12665–12682, https://doi.org/10.5194/acp-14-12665-2014, 2014. 

Finney, D. L., Doherty, R. M., Wild, O., and Abraham, N. L.: The impact of lightning on tropospheric ozone chemistry using a new global lightning parametrisation, Atmos. Chem. Phys., 16, 7507–7522, https://doi.org/10.5194/acp-16-7507-2016, 2016a. 

Finney, D. L., Doherty, R. M., Wild, O., Young, P. J., and Butler, A.: Response of lightning NOx emissions and ozone production to climate change: Insights from the Atmospheric Chemistry and Climate Model Intercomparison Project, Geophys. Res. Lett., 43, 5492–5500, https://doi.org/10.1002/2016GL068825, 2016b. 

Finney, D. L., Doherty, R. M., Wild, O., Stevenson, D. S., MacKenzie, I. A., and Blyth, A. M.: A projected decrease in lightning under climate change, Nat. Clim. Change, 8, 210–213, https://doi.org/10.1038/s41558-018-0072-6, 2018. 

Fujibe, F.: Long-term Change in Lightning Mortality and Its Relation to Annual Thunder Days in Japan, Journal of Natural Disaster Science, 38, 17–29, https://doi.org/10.2328/jnds.38.17, 2017. 

Buechler, D. E., Wright, P. D., and Goodman, S. J.: Lightning/rainfall relationships during COHMEX, 16th Conf. on Severe Local Storms, Kananaskis Park, Alberta, Canada, 22–26 October 1990, 714, http://www.jstor.org/stable/26227919 (last access: 13 October 2023), 1990. 

Goto, D., Nakajima, T., Dai, T., Takemura, T., Kajino, M., Matsui, H., Takami, A., Hatakeyama, S., Sugimoto, N., Shimizu, A., and Ohara, T.: An evaluation of simulated particulate sulfate over East Asia through global model intercomparison, J. Geophys. Res.-Atmos., 120, 6247–6270, https://doi.org/10.1002/2014JD021693, 2015. 

Guenther, A. B., Jiang, X., Heald, C. L., Sakulyanontvittaya, T., Duhl, T., Emmons, L. K., and Wang, X.: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions, Geosci. Model Dev., 5, 1471–1492, https://doi.org/10.5194/gmd-5-1471-2012, 2012. 

Ha, P. T. M., Matsuda, R., Kanaya, Y., Taketani, F., and Sudo, K.: Effects of heterogeneous reactions on tropospheric chemistry: a global simulation with the chemistry–climate model CHASER V4.0, Geosci. Model Dev., 14, 3813–3841, https://doi.org/10.5194/gmd-14-3813-2021, 2021. 

He, Y., Hoque, M. S. H., and Sudo, K.: Introducing new lightning schemes into the CHASER (MIROC) chemistry climate model, Zenodo [code], https://doi.org/10.5281/zenodo.5835796, 2022a. 

He, Y., Hoque, H. M. S., and Sudo, K.: Introducing new lightning schemes into the CHASER (MIROC) chemistry–climate model, Geosci. Model Dev., 15, 5627–5650, https://doi.org/10.5194/gmd-15-5627-2022, 2022b. 

Hoesly, R. M., Smith, S. J., Feng, L., Klimont, Z., Janssens-Maenhout, G., Pitkanen, T., Seibert, J. J., Vu, L., Andres, R. J., Bolt, R. M., Bond, T. C., Dawidowski, L., Kholod, N., Kurokawa, J.-I., Li, M., Liu, L., Lu, Z., Moura, M. C. P., O'Rourke, P. R., and Zhang, Q.: Historical (1750–2014) anthropogenic emissions of reactive gases and aerosols from the Community Emissions Data System (CEDS), Geosci. Model Dev., 11, 369–408, https://doi.org/10.5194/gmd-11-369-2018, 2018. 

Huang, B., Thorne, P. W., Banzon, V. F., Boyer, T., Chepurin, G., Lawrimore, J. H., Menne, M. J., Smith, T. M., Vose, R. S., and Zhang, H.-M.: Extended Reconstructed Sea Surface Temperature, Version 5 (ERSSTv5): Upgrades, Validations, and Intercomparisons, J. Climate, 30, 8179–8205, https://doi.org/10.1175/JCLI-D-16-0836.1, 2017. 

Hui, J. and Hong, L.: Projected Changes in NOx Emissions from Lightning as a Result of 2000–2050 Climate Change, Atmospheric and Oceanic Science Letters, 6, 284–289, https://doi.org/10.3878/j.issn.1674-2834.13.0042, 2013. 

Hussain, Md. M. and Mahmud, I.: pyMannKendall: a python package for non parametric Mann Kendall family of trend tests, Journal of Open Source Software, 4, 1556, https://doi.org/10.21105/joss.01556, 2019. 

Ito, A. and Inatomi, M.: Water-use efficiency of the terrestrial biosphere: A model analysis focusing on interactions between the global carbon and water cycles, J. Hydrometeorol., 13, 681–694, https://doi.org/10.1175/JHM-D-10-05034.1, 2012. 

Jensen, J. D., Thurman, J., and Vincent, A. L.: Lightning Injuries, in: StatPearls, StatPearls Publishing, Treasure Island (FL), https://www.ncbi.nlm.nih.gov/books/NBK431128/ (last access: Please provide your last access date.), 2022. 

Kaufman, Y. J., Tanré, D., Holben, B. N., Mattoo, S., Remer, L. A., Eck, T. F., Vaughan, J., and Chatenet, B.: Aerosol Radiative Impact on Spectral Solar Flux at the Surface, Derived from Principal-Plane Sky Measurements, J. Atmos. Sci., 59, 635–646, https://doi.org/10.1175/1520-0469(2002)059<0635:ARIOSS>2.0.CO;2, 2002. 

Kelley, M., Schmidt, G. A., Nazarenko, L. S., Bauer, S. E., Ruedy, R., Russell, G. L., Ackerman, A. S., Aleinov, I., Bauer, M., Bleck, R., Canuto, V., Cesana, G., Cheng, Y., Clune, T. L., Cook, B. I., Cruz, C. A., Del Genio, A. D., Elsaesser, G. S., Faluvegi, G., Kiang, N. Y., Kim, D., Lacis, A. A., Leboissetier, A., LeGrande, A. N., Lo, K. K., Marshall, J., Matthews, E. E., McDermid, S., Mezuman, K., Miller, R. L., Murray, L. T., Oinas, V., Orbe, C., García-Pando, C. P., Perlwitz, J. P., Puma, M. J., Rind, D., Romanou, A., Shindell, D. T., Sun, S., Tausnev, N., Tsigaridis, K., Tselioudis, G., Weng, E., Wu, J., and Yao, M.-S.: GISS-E2.1: Configurations and Climatology, J. Adv. Model. Earth Sy., 12, e2019MS002025, https://doi.org/10.1029/2019MS002025, 2020. 

Koren, I., Kaufman, Y. J., Remer, L. A., and Martins, J. V.: Measurement of the Effect of Amazon Smoke on Inhibition of Cloud Formation, Science, 303, 1342–1345, https://doi.org/10.1126/science.1089424, 2004. 

Koren, I., Martins, J. V., Remer, L. A., and Afargan, H.: Smoke Invigoration Versus Inhibition of Clouds over the Amazon, Science, 321, 946–949, https://doi.org/10.1126/science.1159185, 2008. 

Krause, A., Kloster, S., Wilkenskjeld, S., and Paeth, H.: The sensitivity of global wildfires to simulated past, present, and future lightning frequency, J. Geophys. Res.-Biogeo., 119, 312–322, https://doi.org/10.1002/2013JG002502, 2014. 

Labrador, L. J., von Kuhlmann, R., and Lawrence, M. G.: The effects of lightning-produced NOx and its vertical distribution on atmospheric chemistry: sensitivity simulations with MATCH-MPIC, Atmos. Chem. Phys., 5, 1815–1834, https://doi.org/10.5194/acp-5-1815-2005, 2005. 

Lal, D. M., Ghude, S. D., Mahakur, M., Waghmare, R. T., Tiwari, S., Srivastava, M. K., Meena, G. S., and Chate, D. M.: Relationship between aerosol and lightning over Indo-Gangetic Plain (IGP), India, Clim. Dynam., 50, 3865–3884, https://doi.org/10.1007/s00382-017-3851-2, 2018. 

Li, Z., Guo, J., Ding, A., Liao, H., Liu, J., Sun, Y., Wang, T., Xue, H., Zhang, H., and Zhu, B.: Aerosol and boundary-layer interactions and impact on air quality, Natl. Sci. Rev., 4, 810–833, https://doi.org/10.1093/nsr/nwx117, 2017. 

Liaskos, C. E., Allen, D. J., and Pickering, K. E.: Sensitivity of tropical tropospheric composition to lightning NOx production as determined by replay simulations with GEOS-5, J. Geophys. Res., 120, 8512–8534, https://doi.org/10.1002/2014JD022987, 2015. 

Liu, Y., Guha, A., Said, R., Williams, E., Lapierre, J., Stock, M., and Heckman, S.: Aerosol Effects on Lightning Characteristics: A Comparison of Polluted and Clean Regimes, Geophys. Res. Lett., 47, e2019GL086825, https://doi.org/10.1029/2019GL086825, 2020. 

Lopez, P.: A lightning parameterization for the ECMWF integrated forecasting system, Mon. Weather Rev., 144, 3057–3075, https://doi.org/10.1175/MWR-D-16-0026.1, 2016. 

Macias Fauria, M. and Johnson, E. A.: Large-scale climatic patterns control large lightning fire occurrence in Canada and Alaska forest regions, J. Geophys. Res.-Biogeo., 111, G04008, https://doi.org/10.1029/2006JG000181, 2006. 

Mallick, C., Hazra, A., Saha, S. K., Chaudhari, H. S., Pokhrel, S., Konwar, M., Dutta, U., Mohan, G. M., and Vani, K. G.: Seasonal Predictability of Lightning Over the Global Hotspot Regions, Geophys. Res. Lett., 49, e2021GL096489, https://doi.org/10.1029/2021GL096489, 2022. 

Manabe, S. and Wetherald, R. T.: The Effects of Doubling the CO2 Concentration on the climate of a General Circulation Model, J. Atmos. Sci., 32, 3–15, https://doi.org/10.1175/1520-0469(1975)032<0003:TEODTC>2.0.CO;2, 1975. 

McCaul, E. W., Goodman, S. J., LaCasse, K. M., and Cecil, D. J.: Forecasting lightning threat using cloud-resolving model simulations, Weather Forecast., 24, 709–729, https://doi.org/10.1175/2008WAF2222152.1, 2009. 

Meinshausen, M., Vogel, E., Nauels, A., Lorbacher, K., Meinshausen, N., Etheridge, D. M., Fraser, P. J., Montzka, S. A., Rayner, P. J., Trudinger, C. M., Krummel, P. B., Beyerle, U., Canadell, J. G., Daniel, J. S., Enting, I. G., Law, R. M., Lunder, C. R., O'Doherty, S., Prinn, R. G., Reimann, S., Rubino, M., Velders, G. J. M., Vollmer, M. K., Wang, R. H. J., and Weiss, R.: Historical greenhouse gas concentrations for climate modelling (CMIP6), Geosci. Model Dev., 10, 2057–2116, https://doi.org/10.5194/gmd-10-2057-2017, 2017. 

Miller, R. L., Schmidt, G. A., Nazarenko, L. S., Tausnev, N., Bauer, S. E., DelGenio, A. D., Kelley, M., Lo, K. K., Ruedy, R., Shindell, D. T., Aleinov, I., Bauer, M., Bleck, R., Canuto, V., Chen, Y., Cheng, Y., Clune, T. L., Faluvegi, G., Hansen, J. E., Healy, R. J., Kiang, N. Y., Koch, D., Lacis, A. A., LeGrande, A. N., Lerner, J., Menon, S., Oinas, V., Pérez García-Pando, C., Perlwitz, J. P., Puma, M. J., Rind, D., Romanou, A., Russell, G. L., Sato, M., Sun, S., Tsigaridis, K., Unger, N., Voulgarakis, A., Yao, M.-S., and Zhang, J.: CMIP5 historical simulations (1850–2012) with GISS ModelE2, J. Adv. Model. Earth Sy., 6, 441–478, https://doi.org/10.1002/2013MS000266, 2014. 

Murray, L. T.: Lightning NOx and Impacts on Air Quality, Current Pollution Reports, 2, 115–133, https://doi.org/10.1007/s40726-016-0031-7, 2016. 

NASA: Earthdata search :: GHRC Portal [data set], https://search.earthdata.nasa.gov/search?portal=ghrc (last access: 10 October 2023), 2023. 

NCEI (National Centers for Environmental Information): Extended Reconstructed SST, National Centers for Environmental Information (NCEI) [data set], https://www.ncei.noaa.gov/products/extended-reconstructed-sst (last access: 27 March 2023), 2023. 

NOAA National Centers for Environmental Information (NCEI), Climate at a Glance: Global Time Series, published June 2023, https://www.ncei.noaa.gov/access/monitoring/climate-at-a-glance/global/time-series/globe/land_ocean/3/8/1880-2020 (last access: 10 January 2023), 2022. 

Ott, L. E., Pickering, K. E., Stenchikov, G. L., Allen, D. J., DeCaria, A. J., Ridley, B., Lin, R. F., Lang, S., and Tao, W. K.: Production of lightning NOx and its vertical distribution calculated from three-dimensional cloud-scale chemical transport model simulations, J. Geophys. Res.-Atmos., 115, 4301, https://doi.org/10.1029/2009JD011880, 2010. 

Pinto, O.: Lightning and climate: A review, in: 2013 International Symposium on Lightning Protection (XII SIPDA), 2013 International Symposium on Lightning Protection (XII SIPDA), 402–404, https://doi.org/10.1109/SIPDA.2013.6729250, 7–11 October 2013, Belo Horizonte, Brazil, 2013. 

Price, C. and Rind, D.: A simple lightning parameterization for calculating global lightning distributions, J. Geophys. Res., 97, 9919–9933, https://doi.org/10.1029/92JD00719, 1992. 

Price, C. and Rind, D.: What determines the cloud-to-ground lightning fraction in thunderstorms?, Geophys. Res. Lett., 20, 463–466, https://doi.org/10.1029/93GL00226, 1993. 

Price, C. and Rind, D.: Possible implications of global climate change on global lightning distributions and frequencies, J. Geophys. Res., 99, 823–833, https://doi.org/10.1029/94jd00019, 1994. 

Rayner, N. A., Parker, D. E., Horton, E. B., Folland, C. K., Alexander, L. V., Rowell, D. P., Kent, E. C., and Kaplan, A.: Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century, J. Geophys. Res.-Atmos., 108, 4407, https://doi.org/10.1029/2002jd002670, 2003. 

Ridley, B. A., Pickering, K. E., and Dye, J. E.: Comments on the parameterization of lightning-produced NO in global chemistry-transport models, Atmos. Environ., 39, 6184–6187, https://doi.org/10.1016/j.atmosenv.2005.06.054, 2005. 

Romps, D. M.: Evaluating the Future of Lightning in Cloud-Resolving Models, Geophys. Res. Lett., 46, 14863–14871, https://doi.org/10.1029/2019GL085748, 2019. 

Romps, D. M., Seeley, J. T., Vollaro, D., and Molinari, J.: Projected increase in lightning strikes in the united states due to global warming, Science, 346, 851–854, https://doi.org/10.1126/science.1259100, 2014. 

Salmi, T., Määttä, A., Anttila, P., Ruoho-Airola, T., and Amnell, T.: Detecting Trends of Annual Values of Atmospheric Pollutants by the Mann–Kendall Test and Sen's Slope Estimates the Excel Template Application MAKESENS, Publications on Air Quality, 31, https://www.researchgate.net/publication/259356944_Detecting_Trends_of_Annual_Values_of_Atmospheric_Pollutants_by_the_Mann-Kendall_Test_and_Sen's_Solpe_Estimates_the_Excel_Template_Application_MAKESENS (last access: 13 October 2023), 2002. 

Sato, M., Hansen, J. E., McCormick, M. P., and Pollack, J. B.: Stratospheric aerosol optical depths, 1850–1990, J. Geophys. Res.-Atmos., 98, 22987–22994, https://doi.org/10.1029/93JD02553, 1993. 

Saunders, C. P. R., Keith, W. D., and Mitzeva, R. P.: The effect of liquid water on thunderstorm charging, J. Geophys. Res.-Atmos., 96, 11007–11017, https://doi.org/10.1029/91JD00970, 1991. 

Schumann, U. and Huntrieser, H.: The global lightning-induced nitrogen oxides source, Atmos. Chem. Phys., 7, 3823–3907, https://doi.org/10.5194/acp-7-3823-2007, 2007. 

Sekiguchi, M. and Nakajima, T.: A k-distribution-based radiation code and its computational optimization for an atmospheric general circulation model, J. Quant. Spectrosc. Ra., 109, 2779–2793, https://doi.org/10.1016/j.jqsrt.2008.07.013, 2008. 

Sellar, A. A., Jones, C. G., Mulcahy, J. P., Tang, Y., Yool, A., Wiltshire, A., O'Connor, F. M., Stringer, M., Hill, R., Palmieri, J., Woodward, S., de Mora, L., Kuhlbrodt, T., Rumbold, S. T., Kelley, D. I., Ellis, R., Johnson, C. E., Walton, J., Abraham, N. L., Andrews, M. B., Andrews, T., Archibald, A. T., Berthou, S., Burke, E., Blockley, E., Carslaw, K., Dalvi, M., Edwards, J., Folberth, G. A., Gedney, N., Griffiths, P. T., Harper, A. B., Hendry, M. A., Hewitt, A. J., Johnson, B., Jones, A., Jones, C. D., Keeble, J., Liddicoat, S., Morgenstern, O., Parker, R. J., Predoi, V., Robertson, E., Siahaan, A., Smith, R. S., Swaminathan, R., Woodhouse, M. T., Zeng, G., and Zerroukat, M.: UKESM1: Description and Evaluation of the U. K. Earth System Model, J. Adv. Model. Earth Sy., 11, 4513–4558, https://doi.org/10.1029/2019MS001739, 2019. 

Shi, Z., Wang, H., Tan, Y., Li, L., and Li, C.: Influence of aerosols on lightning activities in central eastern parts of China, Atmos. Sci. Lett., 21, e957, https://doi.org/10.1002/asl.957, 2020. 

Soden, B. J., Wetherald, R. T., Stenchikov, G. L., and Robock, A.: Global Cooling After the Eruption of Mount Pinatubo: A Test of Climate Feedback by Water Vapor, Science, 296, 727–730, https://doi.org/10.1126/science.296.5568.727, 2002. 

Sudo, K. and Akimoto, H.: Global source attribution of tropospheric ozone: Long-range transport from various source regions, J. Geophys. Res.-Atmos., 112, D12302, https://doi.org/10.1029/2006JD007992, 2007. 

Sudo, K., Takahashi, M., Kurokawa, J. I., and Akimoto, H.: CHASER: A global chemical model of the troposphere 1. Model description, J. Geophys. Res.-Atmos., 107, ACH 7-1–ACH 7-20, https://doi.org/10.1029/2001JD001113, 2002. 

Takemura, T., Egashira, M., Matsuzawa, K., Ichijo, H., O'ishi, R., and Abe-Ouchi, A.: A simulation of the global distribution and radiative forcing of soil dust aerosols at the Last Glacial Maximum, Atmos. Chem. Phys., 9, 3061–3073, https://doi.org/10.5194/acp-9-3061-2009, 2009. 

Tan, Y. B., Peng, L., Shi, Z., and Chen, H. R.: Lightning flash density in relation to aerosol over Nanjing (China), Atmos. Res., 174–175, 1–8, https://doi.org/10.1016/j.atmosres.2016.01.009, 2016. 

Tang, Y., Rumbold, S., Ellis, R., Kelley, D., Mulcahy, J., Sellar, A., Walton, J., and Jones, C.: MOHC UKESM1.0-LL model output prepared for CMIP6 CMIP historical, https://doi.org/10.22033/ESGF/CMIP6.6113, 2019. 

Tost, H.: Chemistry–climate interactions of aerosol nitrate from lightning, Atmos. Chem. Phys., 17, 1125–1142, https://doi.org/10.5194/acp-17-1125-2017, 2017. 

van Marle, M. J. E., Kloster, S., Magi, B. I., Marlon, J. R., Daniau, A.-L., Field, R. D., Arneth, A., Forrest, M., Hantson, S., Kehrwald, N. M., Knorr, W., Lasslop, G., Li, F., Mangeon, S., Yue, C., Kaiser, J. W., and van der Werf, G. R.: Historic global biomass burning emissions for CMIP6 (BB4CMIP) based on merging satellite observations with proxies and fire models (1750–2015), Geosci. Model Dev., 10, 3329–3357, https://doi.org/10.5194/gmd-10-3329-2017, 2017. 

Veraverbeke, S., Finney, D., van der Werf, G., van Wees, D., Xu, W., and Jones, M.: Global attribution of anthropogenic and lightning fires, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-1160, https://doi.org/10.5194/egusphere-egu22-1160, 2022. 

Wang, Q., Li, Z., Guo, J., Zhao, C., and Cribb, M.: The climate impact of aerosols on the lightning flash rate: is it detectable from long-term measurements?, Atmos. Chem. Phys., 18, 12797–12816, https://doi.org/10.5194/acp-18-12797-2018, 2018.  

Wang, Y., Wan, Q., Meng, W., Liao, F., Tan, H., and Zhang, R.: Long-term impacts of aerosols on precipitation and lightning over the Pearl River Delta megacity area in China, Atmos. Chem. Phys., 11, 12421–12436, https://doi.org/10.5194/acp-11-12421-2011, 2011. 

Watanabe, S., Hajima, T., Sudo, K., Nagashima, T., Takemura, T., Okajima, H., Nozawa, T., Kawase, H., Abe, M., Yokohata, T., Ise, T., Sato, H., Kato, E., Takata, K., Emori, S., and Kawamiya, M.: MIROC-ESM 2010: model description and basic results of CMIP5-20c3m experiments, Geosci. Model Dev., 4, 845–872, https://doi.org/10.5194/gmd-4-845-2011, 2011. 

Wild, O.: Modelling the global tropospheric ozone budget: exploring the variability in current models, Atmos. Chem. Phys., 7, 2643–2660, https://doi.org/10.5194/acp-7-2643-2007, 2007. 

Williams, E. R.: Earle Williams, https://web.mit.edu/earlerw/www/index.html, last access: 19 December 2022. 

Williams, E. R., Weber, M. E., and Orville, R. E.: The relationship between lightning type and convective state of thunderclouds, J. Geophys. Res.-Atmos., 94, 13213–13220, https://doi.org/10.1029/JD094iD11p13213, 1989. 

Yang, X., Yao, Z., Li, Z., and Fan, T.: Heavy air pollution suppresses summer thunderstorms in central China, J. Atmos. Sol.-Terr. Phy., 95–96, 28–40, https://doi.org/10.1016/j.jastp.2012.12.023, 2013. 

Yuan, T., Remer, L. A., Pickering, K. E., and Yu, H.: Observational evidence of aerosol enhancement of lightning activity and convective invigoration, Geophys. Res. Lett., 38, 4701, https://doi.org/10.1029/2010GL046052, 2011. 

Zeng, G., Pyle, J. A., and Young, P. J.: Impact of climate change on tropospheric ozone and its global budgets, Atmos. Chem. Phys., 8, 369–387, https://doi.org/10.5194/acp-8-369-2008, 2008. 

Zhao, P., Zhou, Y., Xiao, H., Liu, J., Gao, J., and Ge, F.: Total Lightning Flash Activity Response to Aerosol over China Area, Atmosphere, 8, 26, https://doi.org/10.3390/atmos8020026, 2017. 

Zhao, P., Li, Z., Xiao, H., Wu, F., Zheng, Y., Cribb, M. C., Jin, X., and Zhou, Y.: Distinct aerosol effects on cloud-to-ground lightning in the plateau and basin regions of Sichuan, Southwest China, Atmos. Chem. Phys., 20, 13379–13397, https://doi.org/10.5194/acp-20-13379-2020, 2020. 

Download
Short summary
Lightning has big social impacts. Lightning-produced NOx (LNOx) plays a vital role in atmospheric chemistry and climate. Investigating past lightning and LNOx trends can provide essential indicators of all lightning-related phenomena. Simulations show almost flat global lightning and LNOx trends during 1960–2014. Past global warming enhances the trends positively, but increases in aerosol have the opposite effect. Moreover, global lightning decreased markedly after the Pinatubo eruption.
Altmetrics
Final-revised paper
Preprint