Comprehensive Climate Models Research Articles

Revisiting the reanalysis-model discrepancy in Southern Hemisphere winter storm track trends

Southern Hemisphere (SH) storminess has increased in the satellite era and recent work suggests comprehensive climate models significantly underestimate the trend. Here, we revisit this reanalysis-model trend discrepancy to better understand the mechanisms underlie it. A comprehensive like-for-like analysis shows reanalysis trends exhibit large uncertainty, and coupled climate model simulations exhibit weaker trends than most but not all reanalyses. However, simulations with prescribed sea surface temperature (SST) exhibit significantly greater storminess trends, particularly in the South Pacific, implying SST trend discrepancies in coupled simulations impact storminess trends. Using pacemaker simulations that correct Southern Ocean and tropical east Pacific SST trend discrepancies, we show that storminess trends in coupled simulations are underestimated because they do not capture the enhanced storminess resulting from Southern Ocean cooling and La-Nina-like teleconnection trends. Our findings emphasize large reanalysis uncertainty in SH circulation trends and the impact of regional SST trend discrepancies on circulation trends.

Exploring Mechanisms for Model‐Dependency of the Stratospheric Response to Arctic Warming

AbstractThe Arctic is estimated to have warmed up to four times faster than the rest of the globe since the 1980s. There is significant interest in understanding the mechanisms by which such warming may impact weather and climate at lower latitudes. One such mechanism is the “stratospheric pathway”; Arctic warming is proposed to induce a wave‐driven weakening of the stratospheric polar vortex, which may subsequently impact large‐scale tropospheric circulation. However, recent comprehensive model studies have found systematic differences in both the magnitude and sign of the stratospheric response to Arctic warming. Using a series of idealized model simulations, we show that this response is sensitive to characteristics of the warming and mean polar vortex strength. In all simulations, imposed polar warming amplifies upward wave propagation from the troposphere, consistent with comprehensive models. However, as polar warming strength and depth increases, the region through which waves can propagate is narrowed, inducing wave breaking and deceleration of the flow in the lower stratosphere. Thus, the mid‐stratosphere is less affected, with reduced sudden stratospheric warming frequency for stronger and deeper warming compared to weaker and shallower warming. We also find that the sign of the stratospheric response depends on the mean strength of the vortex, and that the stratospheric response in turn plays a role in the magnitude of the tropospheric jet response. Our results help explain the spread across multimodel ensembles of comprehensive climate models.

Replicating the Hadley cell edge and subtropical jet latitude disconnect in idealized atmospheric models

Abstract. Recent work has shown that variability in the subtropical jet's (STJ) latitude, ϕSTJ, is not coupled to that of the Hadley cell (HC) edge, ϕHC, but the robustness of this disconnect has not been examined in detail. Here, we use meteorological reanalysis products, comprehensive climate models, and an idealized atmospheric model to determine the necessary processes for a disconnect between ϕHC and ϕSTJ in the Northern Hemisphere's December–January–February season. We find that a decoupling can occur in a dry general circulation model, indicating that large-scale dynamical processes are sufficient to reproduce the metrics' relationship. It is therefore not reliant on explicit variability in the zonal structure, convection, or radiation. Rather, the disconnect requires a sufficiently realistic climatological basic state. Further, we confirm that the robust disconnect between ϕSTJ and ϕHC across the model hierarchy reveals their differing sensitivities to midlatitude eddy momentum fluxes; ϕHC is consistently coupled to the latitude of maximum eddy momentum flux, but ϕSTJ is not.

State‐Dependence of the Equilibrium Climate Sensitivity in a Clear‐Sky GCM

AbstractThe climate sensitivity peaks around 310 K in a wide variety of climate models, ranging from idealized single column models to fully comprehensive climate models. Here, we increase CO2 using a clear‐sky three‐dimensional atmospheric model with a radiation scheme which maintains accuracy for high CO2 and temperature levels. In contrast, the Equilibrium Climate Sensitivity (ECS) of our model plateaus around 310 K. We show that this is due to the moistening of the subtropical regions caused by a slowdown in atmospheric circulation, which increases the ECS at very high CO2 values. When relative humidity is fixed, the ECS peak is consistent with single column model results. This work does not rule out that clouds or other complex processes impact the ECS in comprehensive climate models. Though the changes in CO2 here are extreme, this study underlines the importance of changes in atmospheric circulation and relative humidity in quantitative assessments of climate sensitivity.

A Semi‐Analytical Model for Water Vapor, Temperature, and Surface‐Albedo Feedbacks in Comprehensive Climate Models

AbstractRadiative feedbacks govern the Earth's climate sensitivity and elucidate the geographic patterns of climate change in response to a carbon‐dioxide forcing. We develop an analytical model for patterned radiative feedbacks that depends only on changes in local surface temperature. The analytical model combines well‐known moist adiabatic theory with the radiative‐advective equilibrium that describes the energy balance in high latitudes. Together with a classic analytical function for surface albedo, all of the non‐cloud feedbacks are represented. The kernel‐based analytical feedbacks reproduce the feedbacks diagnosed from global climate models at the global, zonal‐mean, and seasonal scales, including in the polar regions, though with less intermodel spread. The analytical model thus provides a framework for a quantitative understanding of radiative feedbacks from simple physics, independent of the detailed atmospheric and cryospheric responses simulated by comprehensive climate models.

Southern Ocean Surface Temperatures and Cloud Biases in Climate Models Connected to the Representation of Glacial Deep Ocean Circulation

Abstract Simulating and reproducing the past Atlantic meridional overturning circulation (AMOC) with comprehensive climate models are essential to understanding past climate changes as well as to testing the ability of the models in simulating different climates. At the Last Glacial Maximum (LGM), reconstructions show a shoaling of the AMOC compared to modern climate. However, almost all state-of-the-art climate models simulate a deeper LGM AMOC. Here, it is shown that this paleodata–model discrepancy is partly related to the climate model biases in modern sea surface temperatures (SST) over the Southern Ocean (70°–45°S). Analysis of model outputs from three phases of the Paleoclimate Model Intercomparison Project shows that models with warm Southern Ocean SST biases tend to simulate a deepening of the LGM AMOC, while the opposite is observed in models with cold SST biases. As a result, a positive correlation of 0.41 is found between SST biases and LGM AMOC depth anomalies. Using sensitivity experiments with a climate model, we show, as an example, that changes in parameters associated with the fraction of cloud thermodynamic phase in a climate model reduce the biases in the warm SST over the modern Southern Ocean. The less biased versions of the model then reproduce a colder Southern Ocean at the LGM, which increases formation of Antarctic Bottom Water and causes shoaling of the LGM AMOC, without affecting the LGM climate in other regions. The results highlight the importance of sea surface conditions and clouds over the Southern Ocean in simulating past and future global climates. Significance Statement To test the ability of comprehensive climate models, simulations of the Last Glacial Maximum (LGM) have been conducted. However, most models simulated a deeper Atlantic meridional overturning circulation (AMOC) in the LGM, which contradicts paleodata suggesting a shallower AMOC. Here, using multimodel analysis and sensitivity experiments with a climate model, we show that this paleodata–model discrepancy is partly related to model biases in the modern Southern Ocean. Improvements in Southern Ocean surface temperatures and clouds reproduce a colder climate over the Southern Ocean at the LGM, which causes an intense shoaling of the AMOC due to increased formation of Antarctic Bottom Water. These results demonstrate the important effect of model biases over the Southern Ocean on simulating past climates.

North‐West Europe Hottest Days Are Warming Twice as Fast as Mean Summer Days

AbstractEurope has seen a rapid increase in the frequency and intensity of hot extremes in recent decades. In this study it is shown, using ERA5 reanalysis data 1960–2021, that the hottest summer days in North‐West Europe are warming approximately twice as fast as mean summer days. Moreover, this pattern stands out as relatively unusual across the Northern Hemisphere. It is also shown that comprehensive climate models fail to capture this difference in trends. A hypothesis is suggested to explain the differential rate of warming between the mean and hottest days, namely that the hottest days are often linked to warm advection from Iberia and North Africa, areas that are warming faster than North‐West Europe. This hypothesis can account for about 25% of the difference between ERA5 and a climate model ensemble and hence further research is needed to understand the drivers of the differing trends in mean and extreme temperature.

Energetic Constraints on the Pattern of Changes to the Hydrological Cycle under Global Warming

Abstract The response of zonal-mean precipitation minus evaporation (P − E) to global warming is investigated using a moist energy balance model (MEBM) with a simple Hadley cell parameterization. The MEBM accurately emulates zonal-mean P − E change simulated by a suite of global climate models (GCMs) under greenhouse gas forcing. The MEBM also accounts for most of the intermodel differences in GCM P − E change and better emulates GCM P − E change when compared to the “wet-gets-wetter, dry-gets-drier” thermodynamic mechanism. The intermodel spread in P − E change is attributed to intermodel differences in radiative feedbacks, which account for 60%–70% of the intermodel variance, with smaller contributions from radiative forcing and ocean heat uptake. Isolating the intermodel spread of feedbacks to specific regions shows that tropical feedbacks are the primary source of intermodel spread in zonal-mean P − E change. The ability of the MEBM to emulate GCM P − E change is further investigated using idealized feedback patterns. A less negative and narrowly peaked feedback pattern near the equator results in more atmospheric heating, which strengthens the Hadley cell circulation in the deep tropics through an enhanced poleward heat flux. This pattern also increases gross moist stability, which weakens the subtropical Hadley cell circulation. These two processes in unison increase P − E in the deep tropics, decrease P − E in the subtropics, and narrow the intertropical convergence zone. Additionally, a feedback pattern that produces polar-amplified warming partially reduces the poleward moisture flux by weakening the meridional temperature gradient. It is shown that changes to the Hadley cell circulation and the poleward moisture flux are crucial for understanding the pattern of GCM P − E change under warming. Significance Statement Changes to the hydrological cycle over the twenty-first century are predicted to impact ecosystems and socioeconomic activities throughout the world. While it is broadly expected that dry regions will get drier and wet regions will get wetter, the magnitude and spatial structure of these changes remains uncertain. In this study, we use an idealized climate model, which assumes how energy is transported in the atmosphere, to understand the processes setting the pattern of precipitation and evaporation under global warming. We first use the idealized climate model to explain why comprehensive climate models predict different changes to precipitation and evaporation across a range of latitudes. We show this arises primarily from climate feedbacks, which act to amplify or dampen the amount of warming. Ocean heat uptake and radiative forcing play secondary roles but can account for a significant amount of the uncertainty in regions where ocean circulation influences the rate of warming. We further show that uncertainty in tropical feedbacks (mainly from clouds) affects changes to the hydrological cycle across a range of latitudes. We then show how the pattern of climate feedbacks affects how the patterns of precipitation and evaporation respond to climate change through a set of idealized experiments. These results show how the pattern of climate feedbacks impacts tropical hydrological changes by affecting the strength of the Hadley circulation and polar hydrological changes by affecting the transport of moisture to the high latitudes.

Variations on a Pathway to an Early Eocene Climate

AbstractThe climate of the early Eocene was characterized by much higher temperatures and a smaller equator‐to‐pole surface temperature gradient than today. Comprehensive climate models have been reasonably successful in simulating that climate in the annual average. However, good simulations of the seasonal variations, and in particular much warmer Arctic winters over land, have proven more difficult. Further, while increased greenhouse gases seems necessary to achieve an Eocene climate, it is unclear whether there is a unique combination of factors that leads to agreement with all available proxies. Here we use a very flexible General Circulation Model to examine the sensitivity of the modeled climate to differences in CO2 concentration, land surface properties, ocean heat transport, and cloud extent and thickness. Even in the absence of ice or changes in cloudiness, increasing the CO2 concentration leads to a polar‐amplified surface temperature change because of increased water vapor levels combined with the lack of convection at high latitudes, with the nonlinear dependence of longwave radiation on temperature amplifying the increase in winter over land. Additional low clouds over Arctic land generally decrease summer temperatures and further increase winter temperatures (except at very high CO2 levels). An increase in the land surface heat capacity, plausible given large changes in vegetation, also decreases the Arctic land seasonality. Thus, different combinations of factors—high CO2 levels, changes in low‐level clouds, and an increase in land surface heat capacity—can lead to a simulation within the proxy uncertainty range of the majority of proxy data.

Blocking and General Circulation in GFDL Comprehensive Climate Models

Abstract This study examines the climatology and dynamics of atmospheric blocking, and the general circulation features that influence blocks in GFDL’s atmosphere-only (AM4) and coupled atmosphere–ocean (CM4) comprehensive models. We compare AM4 and CM4 with reanalysis, focusing on winter in the Northern Hemisphere. Both models generate the correct blocking climatology and planetary-scale signatures of the stationary wave. However, at regional scales some biases exist. In the eastern Pacific and over western North America, both models generate excessive blocking frequency and too strong of a stationary wave. In the Atlantic, the models generate too little blocking and a weakened stationary wave. A block-centered compositing analysis of block-onset dynamics reveals that the models 1) produce realistic patterns of high-frequency (1–6-day) eddy forcing and 2) capture the notable differences in the 500-hPa geopotential height field between Pacific and Atlantic blocking. However, the models fail to reproduce stronger wave activity flux convergence in the Atlantic compared to the Pacific. Overall, biases in the blocking climatology in terms of location, frequency, duration, and area are quite similar between AM4 and CM4 despite the models having large differences in sea surface temperatures and climatological zonal circulation. This could suggest that other factors could be more dominant in generating blocking biases for these GCMs. Significance Statement Atmospheric blocks are persistent high pressure systems that can lead to hazardous weather. Historically, climate models have had trouble capturing blocks, but recent changes in the models might lead to improvements. As such, the work herein investigates the spatial distribution, prevalence, duration, size, and dynamics of wintertime blocking in recent NOAA climate models. Overall, these models capture the long-term-average spatial pattern of blocking, and properly reproduce key dynamical features. However, the models produce too much blocking in the western United States, and too little over the northern Atlantic Ocean and Europe. These blocking biases are consistent with atmospheric stationary waves biases, but not jet stream bias. This downplays the role of jet biases in the models being responsible for blocking biases.

On Timescales and Reversibility of the Ocean's Response to Enhanced Greenland Ice Sheet Melting in Comprehensive Climate Models

AbstractWarming of the North Atlantic region in climate history often was associated with massive melting of the Greenland Ice Sheet. To identify the meltwater's impacts and isolate these from internal variability and other global warming factors, we run single‐forcing simulations including small ensembles using three complex climate models differing only in their ocean components. In 200‐year‐long preindustrial climate simulations, we identify robust consequences of abruptly increasing Greenland runoff by 0.05 Sv: sea level rise of 44 ± 10 cm, subpolar North Atlantic surface cooling of 0.7°C, and a moderate AMOC decline of 1.1–2.0 Sv. The latter two emerge in under three decades—and reverse on the same timescale after the perturbation ends in year 100. The ocean translates the step‐change perturbation into a multidecadal‐to‐centennial signature in the deep overturning circulation. In all simulations, internal variability creates notable uncertainty in estimating trends, time of emergence, and duration of the response.

Robust atmospheric river response to global warming in idealized and comprehensive climate models

AbstractAtmospheric rivers (ARs), narrow intense moisture transport, account for much of the poleward moisture transport in midlatitudes. While studies have characterized AR features and the associated hydrological impacts in a warming climate in observations and comprehensive climate models, the fundamental dynamics for changes in AR statistics (e.g., frequency, length, width) are not well understood. Here we investigate AR response to global warming with a combination of idealized and comprehensive climate models. To that end, we developed an idealized atmospheric GCM with Earth-like global circulation and hydrological cycle, in which water vapor and clouds are modeled as passive tracers with simple cloud microphysics and precipitation processes. Despite the simplicity of model physics, it reasonably reproduces observed dynamical structures for individual ARs, statistical characteristics of ARs, and spatial distributions of AR climatology. Under climate warming, the idealized model produces robust AR changes similar to CESM large ensemble simulations under RCP8.5, including AR size expansion, intensified landfall moisture transport, and an increased AR frequency, corroborating previously reported AR changes under global warming by climate models. In addition, the latitude of AR frequency maximum shifts poleward with climate warming. Further analysis suggests the thermodynamic effect (i.e., an increase in water vapor) dominates the AR statistics and frequency changes while both the dynamic and thermodynamic effects contribute to the AR poleward shift. These results demonstrate that AR changes in a warming climate can be understood as passive water vapor and cloud tracers regulated by large-scale atmospheric circulation, whereas convection and latent heat feedback are of secondary importance.

Persistent Model Biases in the CMIP6 Representation of Stratospheric Polar Vortex Variability

AbstractSudden stratospheric warmings (SSWs) can have major impact on surface wintertime weather, especially at mid‐high latitudes. We do not yet have a complete understanding of why some of these events influence our weather more than others, but one factor may be the dynamical nature of the SSW; whether it involves a split or a displacement of the polar vortex, and one way to explore this is through comprehensive climate models. Here, we analyze the stratospheric dynamics of SSWs within models from the sixth Coupled Model Intercomparison Project (CMIP6). All CMIP6 models simulate SSWs to some degree, although we find a persistent bias in the relative underrepresentation of split vortex events. When comparing with CMIP5 models, large biases persist despite significant model improvements in resolution and in representing atmospheric processes. We show that the simulated displacement frequency is strongly related to climatological lower stratospheric eddy heat flux. The split frequency, on the other hand, is not related to lower stratospheric eddy heat flux, but is strongly related to both the vortex geometry (aspect ratio) and lower stratospheric zonal winds. This suggests that those models with a large positive bias in zonal winds may inhibit the propagation of zonal wavenumber 2 planetary waves from the troposphere, which are associated with split events. Our results suggest how future model development may address these longstanding biases.

Stringent mitigation substantially reduces risk of unprecedented near-term warming rates

Following the Paris Agreement, many countries are enacting targets to achieve net-zero GHG emissions. Stringent mitigation will have clear societal benefits in the second half of this century by limiting peak warming and stabilizing climate. However, the near-term benefits of mitigation are generally thought to be less clear because forced surface temperature trends can be masked by internal variability. Here we use observationally constrained projections from the latest comprehensive climate models and a simple climate model emulator to show that pursuing stringent mitigation consistent with holding long-term warming below 1.5 °C reduces the risk of unprecedented warming rates in the next 20 years by a factor of 13 compared with a no mitigation scenario, even after accounting for internal variability. Therefore, in addition to long-term benefits, stringent mitigation offers substantial near-term benefits by offering societies and ecosystems a greater chance to adapt to and avoid the worst climate change impacts.

Changes in Local and Global Climate Feedbacks in the Absence of Interactive Clouds: Southern Ocean–Climate Interactions in Two Intermediate-Complexity Models

AbstractThe global-mean climate feedback quantifies how much the climate system will warm in response to a forcing such as increased CO2 concentration. Under a constant forcing, this feedback becomes less negative (increasing) over time in comprehensive climate models, which has been attributed to increases in cloud and lapse-rate feedbacks. However, out of eight Earth system models of intermediate complexity (EMICs) not featuring interactive clouds, two also simulate such a feedback increase: Bern3D-LPX and LOVECLIM. Using these two models, we investigate the causes of the global-mean feedback increase in the absence of cloud feedbacks. In both models, the increase is predominantly driven by processes in the Southern Ocean region. In LOVECLIM, the global-mean increase is mainly due to a local longwave feedback increase in that region, which can be attributed to lapse-rate changes. It is enhanced by the slow atmospheric warming above the Southern Ocean, which is delayed due to regional ocean heat uptake. In Bern3D-LPX, this delayed regional warming is the main driver of the global-mean feedback increase. It acts on a near-constant local feedback pattern mainly determined by the sea ice–albedo feedback. The global-mean feedback increase is limited by the availability of sea ice: faster Southern Ocean sea ice melting due to either stronger forcing or higher equilibrium climate sensitivity (ECS) reduces the increase of the global mean feedback in Bern3D-LPX. In the highest-ECS simulation with 4 × CO2 forcing, the feedback even becomes more negative (decreasing) over time. This reduced ice–albedo feedback due to sea ice depletion is a plausible mechanism for a decreasing feedback also in high-forcing simulations of other models.

Linearity of Outgoing Longwave Radiation: From an Atmospheric Column to Global Climate Models

AbstractThe linearity of global‐mean outgoing longwave radiation (OLR) with surface temperature is a basic assumption in climate dynamics. This linearity manifests in global climate models, which robustly produce a global‐mean longwave clear‐sky (LWCS) feedback of 1.9 W/m2/K, consistent with idealized single‐column models (Koll & Cronin, 2018, https//:doi.org/10.1073/pnas.1809868115). However, there is considerable spatial variability in the LWCS feedback, including negative values over tropical oceans (known as the “super‐greenhouse effect”) which are compensated for by larger values in the subtropics/extratropics. Therefore, it is unclear how the idealized single‐column results are relevant for the global‐mean LWCS feedback in comprehensive climate models. Here we show with a simple analytical theory and model output that the compensation of this spatial variability to produce a robust global‐mean feedback can be explained by two facts: (1) When conditioned upon free‐tropospheric column relative humidity (RH), the LWCS feedback is independent of RH, and (2) the global histogram of free‐tropospheric column RH is largely invariant under warming.

A Concise and Effective Expression Relating Subsurface Temperature to the Thermocline in the Equatorial Pacific

AbstractThe temperature of the subsurface water entrained into the surface mixed layer plays a key role in controlling the sea surface temperature (SST) and its interannual variability in the equatorial Pacific. In this paper, we combine a hyperbolic tangent function bounded by the warm pool SST and centered at the thermocline depth with a variable sharpness parameter to describe the time‐space evolutions of the subsurface temperature. Under simple approximations of the sharpness parameter, this concise expression becomes remarkably efficient in capturing the observed and climate‐model simulated subsurface temperature variability in terms of anomalies of the thermocline depth and SST of the El Niño‐Southern Oscillation (ENSO) phenomenon. The formulations for the subsurface temperature and thermocline sharpness developed in this work should be useful tools for evaluating and understanding the role of the thermocline feedback in ENSO behaviors in both theoretical and comprehensive climate models.

Estimating Impacts and Trade‐offs in Solar Geoengineering Scenarios With a Moist Energy Balance Model

AbstractThere are large uncertainties in the potential impacts of solar radiation modification (SRM) and in how these impacts depend on the way SRM is deployed. One open question concerns trade‐offs between latitudinal profiles of insolation reduction and climate response. Here, a moist energy balance model is used to evaluate several SRM proposals, providing fundamental insight into how the insolation reduction profile affects the climate response. The optimal SRM profile is found to depend on the intensity of the intervention, as the most effective profile for moderate SRM focuses the reduction at high latitudes, whereas the most effective profile for strong SRM is tropically amplified. The effectiveness of SRM is also shown to depend on when it is applied, an important factor to consider when designing SRM proposals. Using an energy balance model allows us to provide physical explanations for these results while also suggesting future avenues of research with comprehensive climate models.

Dynamics and Predictability of Hemispheric-Scale Multidecadal Climate Variability in an Observationally Constrained Mechanistic Model

AbstractThis paper addresses the dynamics of internal hemispheric-scale multidecadal climate variability by postulating an energy-balance (EBM) model comprising two deep-ocean oscillators in the Atlantic and Pacific basins, coupled through their surface mixed layers via atmospheric teleconnections. This system is linear and driven by the atmospheric noise. Two sets of the EBM model parameters are developed by fitting the EBM-based mixed-layer temperature covariance structure to best mimic basin-average North Atlantic/Pacific sea surface temperature (SST) covariability in either observations or control simulations of comprehensive climate models within the CMIP5 project. The differences between the dynamics underlying the observed and CMIP5-simulated multidecadal climate variability and predictability are encapsulated in the algebraic structure of the two EBM model versions so obtained: EBMCMIP5 and EBMOBS. The multidecadal variability in EBMCMIP5 is overall weaker and amounts to a smaller fraction of the total SST variability than in EBMOBS, pointing to a lower potential decadal predictability of virtual CMIP5 climates relative to that of the actual climate. The EBMCMIP5 decadal hemispheric teleconnections (and, by inference, those in CMIP5 models) are largely controlled by the variability of the Pacific, in which the ocean, due to its large thermal and dynamical memory, acts as a passive integrator of atmospheric noise. By contrast, EBMOBS features a stronger two-way coupling between the Atlantic and Pacific multidecadal oscillators, thereby suggesting the existence of a hemispheric-scale and, perhaps, global multidecadal mode associated with internal ocean dynamics. The inferred differences between the observed and CMIP5 simulated climate variability stem from a stronger communication between the deep ocean and surface processes implicit in the observational data.

A New Pathway for Incorporating the Impact of Solar Flares and Other Space Activity into Climate Modelling

A recent development in climate research is the recognition that space activity such as solar flares, gamma ray bursts, supernovas etc. can have a significant impact on the Earth. The existing methods of modelling space activity do not model these phenomena explicitly. Instead, their effect on ionization rates is formulated and then built into existing climate models. In this paper, based on the authors’ earlier work, several analytical approaches for modelling a solar flare hitting the Earth are presented. These methods can also be used for modelling gamma-ray bursts or cosmic rays striking our planet. In addition, a method of calculating the total amount of radiation received during a solar flare is proposed. The future use of the proposed modelling approaches within comprehensive global climate models allows for explicit modelling of space activity and provides a new pathway for studying the possible effects of this activity on the Earth.