Significant increase in graupel and lightning occurrence in a warmer climate simulated by prognostic graupel ... - Nature.com

Significant increase in graupel and lightning occurrence in a warmer climate simulated by prognostic graupel parameterization | Scientific Reports

Abstract

There is little consensus among global climate models (CGMs) regarding the response of lightning flash rates to past and future climate change, largely due to graupel not being included in models. Here a two-moment prognostic graupel scheme was incorporated into the MIROC6 GCM and applied in three experiments involving pre-industrial aerosol, present-day, and future warming simulations. The new microphysics scheme performed well in reproducing global distributions of graupel, convective available potential energy, and lightning flash rate against satellite retrievals and reanalysis datasets. The global mean lightning rate increased by 7.1% from the pre-industrial period to the present day, which was attributed to increased graupel occurrence. The impact of future warming on lightning activity was more evident, with the rate increasing by 18.4\(\%\,\textrm{K}^{-1}\) through synergistic contributions of destabilization and increased graupel. In the Arctic, the lightning rate depends strongly on the seasonality of graupel, emphasizing the need to incorporate graupel into GCMs for more accurate climate prediction.

Introduction

Lightning plays a significant role in atmospheric chemistry and environmental processes by providing nitrogen oxides and initiating wildfires 1, 2. Despite the importance of climate prediction, the link between anthropogenic climate change and lightning variability remains controversial. Previous observation-based studies have shown that lightning activity has increased over the last decade 3, 4, particularly in the Arctic region 5, although some studies have reported insignificant trends 6. As lightning may trigger wildfires in high-latitude regions 7, 8, 9, its increasing frequency may exacerbate future warming through the release of \(\textrm{CO}_{2}\) and permafrost methane 10, 11.

Increased lightning in response to global warming is also indicated by modeling studies 12, although some simulations indicate decreasing trends with increasing warming 13. Lightning activity depends on various factors including atmospheric instability, aerosols, cloud depth, and the mass of ice hydrometeors 14, 15, making accurate simulation difficult 16, 17, 18, 19.

A limitation of global climate models (GCMs) is that they do not consider high-density ice hydrometeors (graupel and hail) in cloud microphysical parameterization. Many GCMs treat precipitation (rain, snow, and graupel/hail) diagnostically, with precipitating hydrometeors being removed from the atmosphere in a single model time-step 20. Although a few GCMs incorporate prognostic precipitation 21, 22, little consideration has been given to the link between graupel and lightning 23. An improved understanding of the effects of climate change on lightning activity, at a fundamental process level in GCMs, requires more sophisticated linking of cloud and precipitation parameterization to lightning.

To this end, prognostic precipitation in the Model for Interdisciplinary Research on Climate v.6 (MIROC6), was expanded here to allow for graupel hydrometeor triggering of lightning. The various lightning parameterizations applied have led to variability in model lightning evaluation 16, 17, so a function of convective available potential energy (CAPE) and column ice hydrometeors was employed, which has been shown to perform well against satellite observations 23. The present study considered how changes in graupel and meteorology affect lightning activity through a series of sensitivity experiments involving perturbations of temperature and aerosol conditions, examining regional and seasonal characteristics with a particular focus on the Arctic, where lightning effects may accelerate global warming.

Results

Global climatology for graupel and lightning

The graupel distribution introduced in the model and total ice hydrometeors should be verified before consideration of variations in lightning flash rate. Total precipitating ice hydrometeor (\(\textrm{Q}_{\textrm{frz}}\)) is widely distributed over the tropical eastern Indian Ocean, the Maritime Continent, and the western Pacific where deep convection occurs, and over mid-latitude oceanic regions such as the North Pacific, east coast of North America, and Southern Ocean where super-cooled liquid clouds are frequently observed (Fig. 1 a,b). The global annual mean graupel water path (GWP) is 3.8 \(\textrm{g}\,\textrm{m}^{-2}\), similar to that applied in other global embedded-cloud-resolving models such as Multiscale Modeling Framework models 24, 25.

Spatial distributions of (a) precipitating ice water paths (\(\textrm{Q}_{\textrm{frz}}\); \(\textrm{g}\,\textrm{m}^{-2}\)); (b) zonal mean graupel mass-mixing ratio (\(10^{-6}\,\textrm{kg}\,\textrm{kg}^{-1}\)); (c) graupel effective radius (\(\upmu \textrm{m}\)) at the surface; and (d) zonal mean graupel occurrence (%) simulated in the present-day experiment.

The graupel effective radius at the surface is mapped in Fig. 1 c. It is regime-dependent and varies widely. In general, graupel size is not simply correlated with the \(\textrm{Q}_{\textrm{frz}}\), as its growth is controlled in a complex manner by various microphysical processes (e.g., increasing through riming and collection among other hydrometeors, and decreasing through melting and/or production of smaller ice nuclei), depending on the cloud and environmental regimes 26, 27. Simulated graupel size and fall velocity (Supplementary Fig. 1) vary within a reasonable observed range 28. The modeled graupel occurrence frequency (Fig. 1 d), defined by the fractional coverage of precipitation, is reasonable when compared with the other GCM that includes prognostic graupel 29.

Spatial distributions of (top) the occurrence frequency of graupel (%) from (a) Global Precipitation Measurement Dual-frequency Precipitation Radar (GPM/DPR) retrieval for January 2021–December 2022, and (b) the MIROC6-CHIMERRA PD simulation; (middle) the lightning flash rate (\(\textrm{fl}.\,\textrm{km}^{-2}\,\textrm{yr}^{-1}\)) from (c) LIS/OTD satellite observations for July 1995–February 2014, and (d) MIROC6-CHIMERRA PD simulation; and (bottom) the convective available potential energy (CAPE; \(\textrm{J}\,\textrm{kg}^{-1}\)) from (e) ERA5 climatology for 1980–2021 30 and (f) MIROC6-CHIMERRA PD simulation.

Although there are no global observational products for graupel distribution, the Global Precipitation Measurement (GPM) Dual-frequency Precipitation Radar (DPR) provides the occurrence frequency of heavy ice precipitation (e.g., graupel and hail 31) at latitude of 65\(^\circ\) S to 65\(^\circ\) N 32, 33. Figure 2 a shows the GPM-observed annual occurrence frequency of graupel. Since the distribution is dependent on the definition of high-density ice hydrometeors in satellite retrieval algorithms, biases are introduced when evaluating the model performance against observations (Fig. 34. However, the model performs well in representing the overall graupel distribution.

The observed and simulated global lightning distributions are also shown in Fig. 2, in which the widely reported land–ocean contrast (Fig. 35. Thus, the spatial distribution of lightning is not determined by CAPE alone (Fig. 2 e,f).

The present study used the lightning scheme parameterized by both CAPE and precipitating ice hydrometeors coupled to prognostic graupel [Eq. ( 3)], and the simulated spatial pattern of flash rate is similar to that of LIS/OTD satellite observations (Fig. 2 c,d). The model captured higher lightning flash rates well, particularly over Central Africa, the Amazon, the eastern coast of North America, India, and Southeast Asia. The flash rate is somewhat underestimated over the tropical convective zones, partly because of the coarse vertical resolution and because GCMs cannot explicitly resolve convective clouds. These issues may be addressed in a future study.

The global mean lightning rate increased by 7.1% from the PI to PD simulations (Table 1), due mainly to increased ice hydrometeors over the tropics (Supplementary Fig. 2), implying a signal of ‘cloud invigoration hypothesis’ 36 linking changes in lightning, graupel, and aerosols 14, 37, 38, 39. The cloud liquid water path (CLWP), cloud ice water path (CIWP), and GWP increase in response to the change in aerosol from PI to PD, but the graupel is the most sensitive to this change (Table 1). This is likely because the increased cloud liquid and ice water contribute directly to the source of graupel. Given that CAPE is almost constant (Supplementary Table Methods” section) is − 0.62 \(\textrm{W}\,\textrm{m}^{-2}\) in the new MIROC6-CHIMERRA (Cloud-Hydrometeors Interactive Module with Explicit Rain and Radiation; see “ 34, which emphasizes the important role of the riming of underlying super-cooled liquid clouds by falling snowflakes and graupel at mid latitudes, effectively removing cloud water 40 even though anthropogenic aerosols are increased from PI levels.

Future changes in graupel and lightning

The response of lightning flash rate to a unit increase in global surface air temperature over higher northern latitudes (above 55\(^\circ\) N) is shown in Fig. 3 a. An increased flash rate is predicted for most regions, with a mean increase over the Arctic (66\(^\circ\) N–90\(^\circ\) N) of 12.2\(\%\,\textrm{K}^{-1}\). The spatial pattern is inherently linked to that of graupel (Fig.

Spatial distributions (above 55\(^\circ\) N) of future changes relative to PD (\(\%\,\textrm{K}^{-1}\)) in (a) lightning flash rate, (b) graupel occurrence, (c) CAPE, and (d) graupel effective radius. Dotted regions indicate a statistically significant difference at the 95% confidence level.

Although graupel and CAPE are the main controls on the lightning distribution in the model, the responses of meteorology and large-scale circulation are also important factors in lightning prediction because our simulations suggest that increased graupel is a result of incremental changes in CLWP and CIWP with increasing temperature (Table 1). Graupel size may be related to the cloud and environmental regimes (Fig. 3 d and Supplementary Fig. 41. Although beyond the scope of this study, an improved physical understanding of lightning occurrence and more elaborate parameterization may help to reduce model uncertainty and variability 12, 16.

The global mean sensitivities of graupel occurrence and CAPE to a unit temperature increase are comparable, + 7.9\(\%\,\textrm{K}^{-1}\) and + 9.7\(\%\,\textrm{K}^{-1}\), respectively (Table 1). Over polar regions, however, the contribution of change in graupel occurrence is more evident (Table 3 and 4). As a result of the combined effects of changes in graupel and CAPE, the global mean response of lightning flash rate is + 18.4\(\%\,\textrm{K}^{-1}\), higher than previous estimates of \(-5\) to + 12\(\%\,\textrm{K}^{-1}\) 13, 42.

This increase in lightning flash rate may be attributed mainly to different model performances in representing cloud vertical structures 17, 43, different treatments of convective clouds that are not resolved in GCMs 44, or different microphysics process representations 19. Although the model of the present study included the treatment of graupel in large-scale condensation, a fundamental physical consideration of future lightning changes with regard to convective clouds may have been overlooked 45. The present results indicate the need for a prognostic treatment of precipitation, including high-density ice crystallites, for both cumulus and large-scale condensation; such a GCM is not currently available.

Models and experimental configurations vary among studies in terms of, for example, future scenarios, \(\textrm{CO}_{2}\) and aerosol emissions, horizontal/vertical resolution, and microphysical complexity. Although the responses of atmospheric circulation, surface forcing, and consequent feedback also differ between atmosphere-only and coupled-ocean configurations, there is no detailed documentation. Further study of model perturbation of statistics should be undertaken within a single-model framework.

Change in seasonal variation under the warming

There is evident seasonal variation in lightning activity 46, 47, and the possible effects of future warming on such variation are of interest but poorly documented. In this context, particular attention should be paid to polar regions, as the distinct occurrence of day and night cycles throughout the year in such regions provides an excellent opportunity for study of seasonal changes 48.

Monthly mean annual cycles of (a,d) graupel occurrence (%); (b,e) CAPE (\(\textrm{J}\,\textrm{kg}^{-1}\)); and (c,f) lightning flash rate (\(\textrm{fl}.\,\textrm{km}^{-2}\,\textrm{yr}^{-1}\)) over the Arctic (a–c) and Antarctic (d–f). The PD simulation is in black and the SST + 4K simulation in red. The scale in (e) is one-tenth that in (b). The error bars represent ±1 standard deviation of interannual variability.

Monthly time series of simulated graupel, CAPE, and lightning rate over the Arctic and Antarctic are shown in Fig. 4. In the Arctic region, graupel occurrence is higher in summer, particularly during June–September, and is lower in winter (Fig. 47, 49.

Of note, although differences in CAPE between the PD and SST + 4K experiments were found only in their absolute values (Fig. 4 b), differences in graupel occurrence between these experiments were also found in their seasonality (Fig.

Such seasonal variation in graupel is also evident in the Antarctic (Fig. 4 d) during February–April, but lightning rarely occurs poleward of 65\(^\circ\) S 5, 50 due to the very stable environmental conditions in this region (Fig. 4 e,f). Nevertheless, the Antarctic may be exposed to increased lightning under global warming (+ 26.4\(\%\,\textrm{K}^{-1}\)).

The results indicate that a reliable representation of seasonality in lightning activity under a global warming scenario is difficult to achieve unless the model considers the link between graupel and lightning, thereby emphasizing the need for the inclusion of prognostic graupel in GCMs. Furthermore, lightning occurrence may also be controlled by aerosol loading and type 38, 49. An understanding of interactions between microphysical aspects and lightning morphology at a fundamental process level 39, 51 will be crucial in future studies.

Discussion

The response of lightning flash rate to aerosol perturbations and future global warming was considered using the MIROC6-CHIMERRA GCM. Most GCMs struggle to simulate lightning activity because graupel is ignored-a longstanding hindrance in lightning evaluation. Microphysical parameterization in CHIMERRA was therefore updated here by the inclusion of prognostic graupel.

The new microphysics scheme provided graupel distributions more consistent with those of cloud-resolving models and GPM/DPR satellite retrievals, and lightning flash rates consistent with LIS/OTD satellite observations (Figs. 1, 2). The global mean lightning flash rate is insensitive to PI through PD aerosol perturbations (+ 7.1%), but the impact of future warming on lightning activity is more evident. The global mean lightning rate is predicted to increase by 18.4\(\%\,\textrm{K}^{-1}\), more than in earlier studies 13, 16, 42.

The increased lightning rate in the present model is attributed to destabilization and more frequent graupel, with the latter predominating at higher latitudes (Fig. 3 and Supplementary Fig. 1). In particular, the lightning rate in the Arctic is strongly dependent on the seasonality of graupel (Fig. 4). The month of graupel maximum occurrence is delayed to late summer in response to future warming, resulting in a larger difference between PD and future climate during the autumn season (Fig. 21, 22, 52, cumulus precipitation is still treated as diagnostically, and graupel and hail are generally not represented in GCMs. Future work on improving lightning projections will also lead to an improved understanding of lightning-produced \(\textrm{NO}_{\textrm{x}}\) 53 and lightning-ignited wildfires 9.

Variations among models are likely attributable to model performance in simulating cloud vertical structures, different treatments of convective clouds, and different microphysics process representations 16, 17, 18, 54. Therefore, model and scheme intercomparisons should be undertaken to elucidate the sources of uncertainty at the fundamental process level, with fixed model configurations. The present study may provide a benchmark for further experiments with different model configurations (e.g., involving \(\textrm{CO}_{\textrm{2}}\) and aerosol emission scenarios, different horizontal/vertical resolutions, and atmosphere-only or ocean-coupled scenarios), microphysics, and lightning schemes. Overall, the degree of confidence in simulating future responses of lightning to global warming depends on prognostic graupel parameterization being included in microphysics–lightning interactions.

Methods

MIROC6 with CHIMERRA prognostic graupel parameterization

The latest version of the MIROC6 global aerosol–climate model was used. This differs from the same version in the Coupled Model Intercomparison Project Phase 6 (CMIP6) 55 in that it incorporates prognostic precipitation 22. The model includes a two-moment large-scale microphysics scheme, the Cloud-Hydrometeors Interactive Module with Explicit Rain and Radiation (CHIMERRA), which considers both the mass and number mixing ratios of rain and snow. Precipitating hydrometeors also impact the radiation budget, so precipitation was radiatively active in the present model 34. However, high-density ice hydrometeors (graupel and hail) are still not prognostic variables, which has been an obstacle to investigations of local orographic precipitation 29 and the link between graupel and lightning in GCMs.

Therefore, another version of the model was devised with prognostic graupel parameterization. The bulk graupel class was treated as a two-moment scheme, based on the following prognostic equations:

where \(q_{g}\) and \(N_{g}\) are the mass and number mixing ratios of graupel, respectively; \(\rho _{a}\) is air density (\(\textrm{kg}\,\textrm{m}^{-3}\)); \(\textbf{u}\) is the wind vector representing horizontal advection (\(\textrm{m}\,\textrm{s}^{-1}\)); \(v_{q_{g}}\) and \(v_{N_{g}}\) are mass- and number-weighted fall velocities (\(\textrm{m}\,\textrm{s}^{-1}\)), respectively, representing vertical sedimentation; and \(S_{q_{g}}\) and \(S_{N_{g}}\) are respectively the source and sink tendencies (\(\textrm{s}^{-1}\)) of graupel with respect to deposition and sublimation, collection among other hydrometeors (i.e., accretion and riming), and melting, based on Reisner et al. 56. Graupel was radiatively inactive.

The fall speed of larger graupel/hail may be \(>10\) \(\textrm{m}\,\textrm{s}^{-1}\), which is computationally expensive during simulations of vertical sedimentation; the Eulerian solver requires short time-steps consistent with Courant–Friedrichs–Lewy (CFL) criteria. Therefore, the vertical sedimentation scheme was updated to a time-implicit precipitation sedimentation scheme 52. The model was solved for microphysics with 60 s iterations 22, allowing the sedimentation process to be solved at the same time-step, reducing computational costs.

Lightning parameterization and satellite data

The new parameterization with prognostic graupel enabled an evaluation of the frequency of lightning activity on a global scale. Among the various types of lightning scheme, it was convenient to apply one that explicitly considers graupel variation in considering the link between lightning rate and graupel morphology.

A scheme was chosen that combines frozen ice hydrometeors in the air column, \(\textrm{Q}_{\textrm{frz}}\) (\(\textrm{kg}\,\textrm{m}^{-2}\)), and CAPE (\(\textrm{J}\,\textrm{kg}^{-1}\)) to parameterize the lightning flash rate, flight (\(\textrm{fl}.\,\textrm{m}^{-2}\,\textrm{s}^{-1}\)), following He et al. 23:

where \(\alpha\) is a constant, and set to \(2.67 \times 10^{-16}\) and \(1.68 \times 10^{-17}\) for land and ocean, respectively, as applied by He et al. 23 for best performance among the various schemes; and \(\textrm{Q}_{\textrm{frz}}\) is the column integrated precipitating ice hydrometeors (cloud ice, snow, and graupel) between the 0 \(^\circ\)C (\(z_{0}\)) and −25 \(^\circ\)C (\(z_{-25}\)) isotherm levels:

where \(q_{i}\), \(q_{s}\), and \(q_{g}\) are the mass mixing ratios (\(\textrm{kg}\,\textrm{kg}^{-1}\)) for cloud ice, snow, and graupel, respectively, which were prognosed in the model; and \(\textrm{Q}_{\textrm{frz}}\) is a proxy for the charging rate due to collisions of graupel with other hydrometeors within clouds 57. The parameterized lightning flash rate is highly variable, depending on CAPE and \(\textrm{Q}_{\textrm{frz}}\) (Supplementary Fig. 5). Although \(\textrm{Q}_{\textrm{frz}}\) includes cloud ice and snow in addition to graupel, the simulated lightning rate shows a good correlation with GWP (Supplementary Fig. 6).

For model validation, the Lightning Imaging Sensor and the Optical Transient Detector (LIS/OTD) are widely used 58, 59. This satellite dataset provides flash number density (\(\textrm{fl}.\,\textrm{km}^{-2}\)) with a \(2.5^\circ \times 2.5^\circ\) spatial resolution, and was used here in developing a climatology for the available period of July 1995–February 2014.

Experimental setup and evaluations

The following three sets of experiments were undertaken to investigate the sensitivity of lightning activity to aerosol perturbations and future warming: (1) a present-day (PD, year 2000) control experiment; (2) a pre-industrial (PI, year 1850) aerosol experiment; and (3) a future-warming experiment with a uniform 4 K increase in sea surface temperature (SST + 4K). Model resolution was \(1.4^\circ \times 1.4^\circ\) with 40 vertical levels, and an atmosphere-only configuration with year 2000 boundary conditions and prescribed climatological SST and sea ice. The model time-step was 12 min. A suite of simulations was undertaken for 11 years, and the subsequent 10 years used for analysis.

The impact of graupel treatment on the aerosol–cloud interaction (ACI) was assessed through \(\textrm{ERF}_{\textrm{aci}}\) 60 quantified in the PD and PI experiments as follows:

where CRE is cloud radiative effect under clean-sky conditions 61 for the removal of contamination of aerosols (i.e., aerosol–radiation interactions or direct effects).

Future warming effects on graupel and lightning were evaluated by differences between experiments (1) and (3). Cloud, precipitation, and meteorological representations are highly model-dependent, so here the percentage change is reported in response to a unit surface air temperature perturbation (\(\%\,\textrm{K}^{-1}\)), as reported in previous studies 13, 42.

Accession codes

LIS/OTD satellite data are available at https://cmr.earthdata.nasa.gov/search/concepts/C1995865015-GHRC_DAAC.ht…. ERA5 products can be obtained from https://doi.org/10.24381/cds.f17050d7. The GPM 2ADPR product is downloaded from JAXA Global Portal System ( https://gportal.jaxa.jp/gpr/). Simulation data is made publicly available through Zenodo ( https://doi.org/10.5281/zenodo.7988296).
Mao, J. et al. Global impact of lightning-produced oxidants. Geophys. Res. Lett. 48, 1–9.
https://doi.org/10.1029/2021GL095740(2021).
Hessilt, T. D. et al. Future increases in lightning ignition efficiency and wildfire occurrence expected from drier fuels in boreal forest ecosystems of western North America. Environ. Res. Lett.
https://doi.org/10.1088/1748-9326/ac6311(2022).
Saha, U. et al. On the association of lightning activity and projected change in climate over the Indian sub-continent. Atmos. Res. 183, 173–190.
https://doi.org/10.1016/j.atmosres.2016.09.001(2017).
Qie, X. et al. Significantly increased lightning activity over the Tibetan plateau and its relation to thunderstorm genesis. Geophys. Res. Lett. 49, e2022GL099894.
https://doi.org/10.1029/2022GL099894(2022).
Holzworth, R. H. et al. Lightning in the arctic. Geophys. Res. Lett. 48, 1–6.
https://doi.org/10.1029/2020GL091366(2021).
Cecil, D. J., Buechler, D. E. & 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).
Krause, A., Kloster, S., Wilkenskjeld, S. & Paeth, H. The sensitivity of global wildfires to simulated past, present, and future lightning frequency. J. Geophys. Res. Biogeosci. 119, 312–322.
https://doi.org/10.1002/2013JG002502(2014).
He, J., Loboda, T. V., Chen, D. & French, N. H. Cloud-to-ground lightning and near-surface fire weather control wildfire occurrence in Arctic Tundra. Geophys. Res. Lett. 49, e2021GL096814.
https://doi.org/10.1029/2021GL096814(2022).
Pérez-Invernón, F. J., Gordillo-Vázquez, F. J., Huntrieser, H. & Jöckel, P. Variation of lightning-ignited wildfire patterns under climate change. Nat. Commun.
https://doi.org/10.1038/s41467-023-36500-5(2023).
Chen, Y. et al. Future increases in Arctic lightning and fire risk for permafrost carbon. Nat. Clim. Chang. 11, 404–410.
https://doi.org/10.1038/s41558-021-01011-y(2021).
Finney, D. L. Lightning threatens permafrost. Nat. Clim. Chang. 11, 379–380.
https://doi.org/10.1038/s41558-021-01016-7(2021).
Romps, D. M., Seeley, J. T., Vollaro, D. & 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).
Finney, D. L. et al. African lightning and its relation to rainfall and climate change in a convection-permitting model. Geophys. Res. Lett.
https://doi.org/10.1029/2020GL088163(2020).
Thornton, J. A., Virts, K. S., Holzworth, R. H. & Mitchell, T. P. Lightning enhancement over major oceanic shipping lanes. Geophys. Res. Lett. 44, 9102–9111.
https://doi.org/10.1002/2017GL074982(2017).
Han, Y., Luo, H., Wu, Y., Zhang, Y. & Dong, W. Cloud ice fraction governs lightning rate at a global scale. Commun. Earth and Environ. 2, 1–9.
https://doi.org/10.1038/s43247-021-00233-4(2021).
Clark, S. K., Ward, D. S. & Mahowald, N. M. Parameterization-based uncertainty in future lightning flash density. Geophys. Res. Lett. 44, 2893–2901.
https://doi.org/10.1002/2017GL073017(2017).
Finney, D. L. et al. A projected decrease in lightning under climate change. Nat. Clim. Chang. 8, 210–213.
https://doi.org/10.1038/s41558-018-0072-6(2018).
Gordillo-Vázquez, F. J., Pérez-Invernón, F. J., Huntrieser, H. & Smith, A. K. Comparison of six lightning parameterizations in CAM5 and the impact on global atmospheric chemistry. Earth Space Sci. 6, 2317–2346.
https://doi.org/10.1029/2019EA000873(2019).
Charn, A. B. & Parishani, H. Predictive proxies of present and future lightning in a superparameterized model. J. Geophys. Res. Atmos. 126, e2021JD035461.
https://doi.org/10.1029/2021JD035461(2021).
Michibata, T. Aerosol–Cloud Interactions in the Climate System (Springer, 2023).
Gettelman, A., Morrison, H., Santos, S., Bogenschutz, P. & Caldwell, P. M. Advanced two-moment bulk microphysics for global models. Part II: Global model solutions and aerosol–cloud interactions. J. Clim. 28, 1288–1307.
https://doi.org/10.1175/JCLI-D-14-00103.1(2015).
Michibata, T., Suzuki, K., Sekiguchi, M. & Takemura, T. Prognostic precipitation in the MIROC6-SPRINTARS GCM: Description and evaluation against satellite observations. J. Adv. Model. Earth Syst. 11, 839–860.
https://doi.org/10.1029/2018MS001596(2019).
He, Y., Mohammed, H., Hoque, S. & 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(2022).
Wang, M. et al. Aerosol indirect effects in a multi-scale aerosol–climate model PNNL-MMF. Atmos. Chem. Phys. 11, 5431–5455.
https://doi.org/10.5194/acp-11-5431-2011(2011).
Chern, J. D. et al. Performance of the Goddard multiscale modeling framework with Goddard ice microphysical schemes. J. Adv. Model. Earth Syst. 8, 66–95.
https://doi.org/10.1002/2015MS000469(2016).
Fitch, K. E. & Garrett, T. J. Graupel precipitating from thin Arctic clouds with liquid water paths less than 50 g m−2. Geophys. Res. Lett. 49, e2021GL094075.
https://doi.org/10.1029/2021gl094075(2021).
Fitch, K. E. & Garrett, T. J. Measurement and analysis for the parameterization of Arctic snow with frequent occurrence of riming. J. Geophys. Res. Atmos. 127, e2021JD035980.
https://doi.org/10.1029/2021jd035980(2022).
Heymsfield, A., Szakáll, M., Jost, A., Giammanco, I. & Wright, R. A comprehensive observational study of graupel and hail terminal velocity, mass flux, and kinetic energy. J. Atmos. Sci. 75, 3861–3885.
https://doi.org/10.1175/JAS-D-18-0035.1(2018).
Gettelman, A., Morrison, H., Thayer-Calder, K. & Zarzycki, C. M. The impact of rimed ice hydrometeors on global and regional climate. J. Adv. Model. Earth Syst. 11, 1543–1562.
https://doi.org/10.1029/2018MS001488(2019).
Hersbach, H. et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 146, 1999–2049.
https://doi.org/10.1002/qj.3803(2020).
Iguchi, T., Kawamoto, N. & Oki, R. Detection of intense ice precipitation with GPM/DPR. J. Atmos. Oceanic Tech. 35, 491–502.
https://doi.org/10.1175/JTECH-D-17-0120.1(2018).
Mroz, K. et al. Hail-detection algorithm for the GPM core observatory satellite sensors. J. Appl. Meteorol. Climatol. 56, 1939–1957.
https://doi.org/10.1175/JAMC-D-16-0368.1(2017).
Le, M. & Chandrasekar, V. Graupel and hail identification algorithm for the dual-frequency precipitation radar (DPR) on the GPM core satellite. J. Meteorol. Soc. Jpn. 99, 49–65.
https://doi.org/10.2151/jmsj.2021-003(2021).
Michibata, T., Suzuki, K. & Takemura, T. Snow-induced buffering in aerosol–cloud interactions. Atmos. Chem. Phys. 20, 13771–13780.
https://doi.org/10.5194/acp-20-13771-2020(2020).
Deierling, W., Petersen, W. A., Latham, J., Ellis, S. & Christian, H. J. The relationship between lightning activity and ice fluxes in thunderstorms. J. Geophys. Res. Atmos. 113, D15210.
https://doi.org/10.1029/2007JD009700(2008).
Rosenfeld, D. et al. Global observations of aerosol–cloud–precipitation–climate interactions. Rev. Geophys.
https://doi.org/10.1002/2013RG000441.Received(2014).
Yuan, T., Remer, L. A., Pickering, K. E. & Yu, H. Observational evidence of aerosol enhancement of lightning activity and convective invigoration. Geophys. Res. Lett. 38, L04701.
https://doi.org/10.1029/2010GL046052(2011).
Pan, Z. et al. Coarse sea spray inhibits lightning. Nat. Commun.
https://doi.org/10.1038/s41467-022-31714-5(2022).
Sun, M. et al. Understanding the effects of aerosols on electrification and lightning polarity in an idealized supercell thunderstorm via model emulation. J. Geophys. Res. Atmos. 129, e2023JD039251.
https://doi.org/10.1029/2023JD039251(2024).
Glassmeier, F. & Lohmann, U. Constraining precipitation susceptibility of warm, ice- and mixed-phase clouds with microphysical equations. J. Atmos. Sci. 73, 5003–5023.
https://doi.org/10.1175/JAS-D-16-0008.1(2016).
Sun, M. et al. Aerosol effects on electrification and lightning discharges in a multicell thunderstorm simulated by the WRF-ELEC model. Atmos. Chem. Phys. 21, 14141–14158.
https://doi.org/10.5194/acp-21-14141-2021(2021).
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).
Heuscher, L., Liu, C., Gatlin, P. & Petersen, W. A. Relationship between lightning, precipitation, and environmental characteristics at mid-/high latitudes from a GLM and GPM perspective. J. Geophys. Res. Atmos. 127, e2022JD036894.
https://doi.org/10.1029/2022jd036894(2022).
Field, P. R., Roberts, M. J. & Wilkinson, J. M. Simulated lightning in a convection permitting global model. J. Geophys. Res. Atmos. 123, 9370–9377.
https://doi.org/10.1029/2018JD029295(2018).
Trapp, R. J., Hoogewind, K. A. & Lasher-Trapp, S. Future changes in hail occurrence in the United States determined through convection-permitting dynamical downscaling. J. Clim. 32, 5493–5509.
https://doi.org/10.1175/JCLI-D-18-0740.1(2019).
Price, C., Penner, J. & Prather, M. NOx from lightning 1. Global distribution based on lightning physics. J. Geophys. Res. Atmos. 102, 5929–5941.
https://doi.org/10.1029/96jd03504(1997).
Kaplan, J. O. & Lau, K. H. K. The WGLC global gridded lightning climatology and time series. Earth Syst. Sci. Data 13, 3219–3237.
https://doi.org/10.5194/essd-13-3219-2021(2021).
Goosse, H. et al. Quantifying climate feedbacks in polar regions. Nat. Commun.
https://doi.org/10.1038/s41467-018-04173-0(2018).
Wang, Q., Li, Z., Guo, J., Zhao, C. & 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).
Kaplan, J. O. & Lau, K.H.-K. World wide lightning location network (WWLLN) global lightning climatology (WGLC) and time series, 2022 update. Earth Syst. Sci. Data 14, 5665–5670.
https://doi.org/10.5194/essd-14-5665-2022(2022).
Tao, W. K., Chen, J. P., Li, Z., Wang, C. & Zhang, C. Impact of aerosols on convective clouds and precipitation. Rev. Geophys.
https://doi.org/10.1029/2011RG000369(2012).
Guo, H. et al. Two-moment bulk cloud microphysics with prognostic precipitation in GFDL’s Atmosphere Model AM4.0: configuration and performance. J. Adv. Model. Earth Syst. 13, e2020MS002453.
https://doi.org/10.1029/2020MS002453(2021).
Tost, H., Jöckel, P. & Lelieveld, J. Lightning and convection parameterisations—Uncertainties in global modelling. Atmos. Chem. Phys. 7, 4553–4568.
https://doi.org/10.5194/acp-7-4553-2007(2007).
Taylor, P. C., Boeke, R. C., Li, Y. & Thompson, D. W. J. Arctic cloud annual cycle biases in climate models. Atmos. Chem. Phys. 19, 8759–8782.
https://doi.org/10.5194/acp-19-8759-2019(2019).
Tatebe, H. et al. Description and basic evaluation of simulated mean state, internal variability, and climate sensitivity in MIROC6. Geosci. Model Dev. 12, 2727–2765.
https://doi.org/10.5194/gmd-12-2727-2019(2019).
Reisner, J., Rasmussen, R. M. & Bruintjes, R. T. Explicit forecasting of supercooled liquid water in winter storms using the MM5 mesoscale model. Q. J. R. Meteorol. Soc. 124, 1071–1107.
https://doi.org/10.1256/smsqj.54803(1998).
McCaul, E. W., Goodman, S. J., LaCasse, K. M. & Cecil, D. J. Forecasting lightning threat using cloud-resolving model simulations. Weather Forecast. 24, 709–729.
https://doi.org/10.1175/2008WAF2222152.1(2009).
Boccippio, D. J., Boeck, W. L. & Blakeslee, R. J. Performance assessment of the optical transient detector and lightning imaging sensor. Part I: Predicted diurnal variability. J. Geophys. Res. Atmos. 19, 1318–1332.
https://doi.org/10.1029/2006JD007787(2002).
Christian, H. J. et al. Global frequency and distribution of lightning as observed from space by the optical transient detector. J. Geophys. Res. Atmos.
https://doi.org/10.1029/2002JD002347(2003).
Bellouin, N. et al. Bounding global aerosol radiative forcing of climate change. Rev. Geophys. 58, e2019RG000660.
https://doi.org/10.1029/2019rg000660(2020).
Ghan, S. J. Technical note: Estimating aerosol effects on cloud radiative forcing. Atmos. Chem. Phys. 13, 9971–9974.
https://doi.org/10.5194/acp-13-9971-2013(2013).

Acknowledgements

This research was supported by the JST FOREST Program (Grant JPMJFR206Y); the Japan Society for the Promotion of Science KAKENHI (Grants JP23K13171 and JP19H05669); the Environment Research and Technology Development Fund (Grant JPMEERF21S12004) of the Environmental Restoration and Conservation Agency, provided by the Ministry of Environment of Japan; and the Advanced Studies of Climate Change Projection (SENTAN) of the Ministry of Education, Culture, Sports, Science, and Technology (Grant JPMXD0722680395). Simulations by MIROC were executed with the SX-Aurora TSUBASA supercomputer system of the National Institute for Environmental Studies, Japan. Finally, the author thanks two anonymous reviewers for providing constructive suggestions and comments, which have significantly improved the paper.

Ethics declarations

Competing interests

The author declares no competing interests.

Additional information

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

About this article

Cite this article

Michibata, T. Significant increase in graupel and lightning occurrence in a warmer climate simulated by prognostic graupel parameterization. Sci Rep 14, 3862 (2024). https://doi.org/10.1038/s41598-024-54544-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-54544-5

Keywords

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.