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 signifi-cant 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.

graph failures all over Northern America and Europe and produced visible aurorae as far south as Cuba and Hawaii. Another example is a flare which occurred in March 1989, when power was knocked out across most of the Quebec province in Canada. Recently, much attention has focused on increasing our understanding of the Carrington event, in order to better quantify what extreme space events can do to our current technological society [2]. Events of this magnitude can also lead to satellite breakdown, with high associated financial costs.
In addition, recent predictions indicate that such extreme events will be more likely over the coming decades [3].
There was a number of attempts to model the influence of large solar flares that took place in the past on the Earth [1] [2] [4] [5] [6] [7] [8]. In order to investigate their effect, certain increased ionization rates were assigned at different altitudes. For flares which occurred during the satellite era, the rates were calculated based on the proton flux or X-ray flux measurements (e.g. [4] [7]). In order to simulate solar flares which occurred before the spacecraft era, the rates of ionisation were either scaled up from the known events to an estimated intensity [1], or were taken from the ice core nitrate levels [2]. After assigning the ionisation rates the global climate model was run under specified conditions. The different chemical reactions which occur in the atmosphere and the transport processes are modelled.
Also, a number of papers attempted to model the effect of different forms of space activity, such as gamma-ray bursts [9] [10], cosmic rays [11] [12] [13] [14] and supernovas [15] [16] [17]. However, as for the modelling of solar flares, these activities were not modelled directly. Instead, their impact on various chemical species in the atmosphere was formulated and incorporated into global climate models.
In this paper, an overview of the approach for calculating insolation from the space perspective presented in the authors' work [18] is provided. The derivation of the model uses the methods of vector field theory and surface integrals. Next, three different approaches for explicit modelling of a solar flare striking the Earth are introduced. These approaches include modelling the solar flare as a vector, as well as modelling it as a cylinder with two different power distributions. These methods can also be used to model gamma-ray bursts or cosmic rays hitting the Earth. Finally, a formula for the total amount of radiation received during a solar flare event is derived.

An Overview of the Approach for Calculating Insolation from the Space Perspective
For an overall description, please see [18]. A summary of the approach along with a more detailed description of several aspects is given here.
In Figure 1, the radiation vector is   , , x y z coordinates. The reason for this choice is that the formulation of the radiation vector remains the same as the Earth rotates around the Sun. The projection of the latitudinal belt in the X 1 OZ 1 , X 1 OY 1 and Y 1 OZ 1 coordinate planes were obtained (see Figure 2).
We should note that the Earth was assumed to be exactly spherical in our formulation. The projection obtained in the Y 1 OZ 1 coordinate plane forms a nearly trapezoidal region, while the projections in the X 1 OY 1 and the X 1 OZ 1 coordinate plane form ellipses. Since the projections obtained in the Y 1 OZ 1 plane are the easiest for calculations, it was chosen for projecting a latitudinal belt in the ( ) , , x y z coordinate frame.
The amount of light which is received by any particular area changes throughout the year. In the current approach it was modelled by an ellipse which changes with time since it is determined by the angle of the Earth's rotation around the Sun (see Figure 2). An ellipse change affects the change of size of the area receiving sunlight (the black hatched area in Figure 2) as well as size of the Journal of Applied Mathematics and Physics area of a latitudinal belt receiving sunlight (the illumination area). For different latitudes the illumination areas will be different. In addition, they may vary with time for some latitudinal belts. For details on this, please see [18].
The equations of the ellipse, circle, lower boundary of the belt and upper boundary of the belt were derived for the Y 1 OZ 1 coordinate plane (see Figure 3).
In Figure 4, we show the direction of the radiation vector from the Earth's point of view and the surface orientation. The positively-oriented surface (S + ) in The lower ( For the derivation details, please see [18]. Based on the approach outlined above, the amount of insolation for a complete year cycle for any latitudinal belt was computed in Maple software with the time step of approximately 10 days. We have also calculated the annual average value for each latitudinal belt and compared it with satellite data from NASA (see [18]). The results show very good agreement for equatorial and middle latitudes region with an average accuracy of 98%. Slightly less agreement (94.5%) can be observed in the polar regions which are known to be difficult to model.
The average accuracy of the approach is 97%.

The Mathematical Formulation of Solar Flares and Other Space Activity
In this section, building on the work introduced in Section 1, several different

Modelling a Solar Flare as a Vector
Firstly, we aim to model a solar flare simply as an additional radiation vector coming from the Sun and entering the Earth at an angle φ . This, together with solar radiation vector, is illustrated in Figure 5.
The flare radiation vector f a was assigned as: where f N is the power of flare.
A flare radiation flux per unit time can be calculated as: where S is the surface area of the Earth and n is the outward unit normal vector.
This formulation of the flare radiation flux is similar to the computation of the regular radiation flux. Different types of radiation may also be included into the modelling in a similar way, i.e. gamma-ray bursts or cosmic rays.

Modelling a Solar Flare as a Cylinder with Power Distribution Depending on Magnitude of Position Vector
Next, we propose to model radiation flux released as the result of a solar flare via a cylinder with a certain radius a R (see Figure 6). However, the power of the flare   is different in different areas of the cylinder. In particular, it is the highest in the centre and decreases towards the edges of the cylinder. This is illustrated by arrows of a different length in Figure 6.
In order to assign a flare power at any point of the cylinder, a polar coordinate system is introduced (see Figure 7).
Since flare power depends on radial distance and does not depend on polar angle, power at any point on a cross section of the cylinder can be defined as: In order to include the whole area of the cross section, the following integral was calculated: Therefore, flare power distribution within a cylinder can be assigned as: The flare radiation vector then can be written as: This modelling approach allows us to take into account the variation of flare power within a flare flux.

Modelling a Solar Flare as a Cylinder with Power Distribution Depending on Magnitude of Position Vector and Time
In addition to the uneven distribution of flare power within a flux, it also varies as flux passes the Earth. This can be incorporated into the modelling by assigning a flare power with a certain distribution. if , , where 1 t corresponds to the start of the flare and 2 t corresponds to the end.
The distribution of flare power within a cylinder and the flare radiation vector then take the form: and respectively.
The proposed way of modelling solar flares allows for the incorporation of radiation dispersal at the beginning and end of a flare flux due to entropy.

Calculation of Total Amount of Radiation during Solar Flare
The total amount of radiation received per m 2 in the Earth is the sum of the amount of regular solar radiation and the amount of radiation received from a flare: . total main flare The amount of regular solar radiation during a flare can be calculated as: t t S main a n S t I S × = ∫ ∫∫ (15) The amount of radiation received from a flare can be found as: a t R f t S flare a n r S r t I S × = ∫ ∫ ∫∫ (16) where r is the Jacobian determinant.

Conclusions and Future Work
In this paper, built on the authors' earlier work, a number of different approaches for explicit modelling of a solar flare striking the Earth are presented.
These methods allow for taking into account the variation of flare power within a flare flux, as well as for the incorporation of radiation dispersal at the beginning and end of a flare flux due to entropy. In addition, it can be possible to model gamma-ray bursts or cosmic rays hitting the Earth in a similar way. Finally, the cumulative radiation received in a solar flare event is determined.
We should note that it was assumed that the flare flux is wide and covers the Earth completely. The cases when the flux is narrow and covers only a certain Journal of Applied Mathematics and Physics area of the Earth have not been considered so far. However, the most complex of the proposed approaches allows for modelling these cases in the future. In order to do this, the intersection of the surface area of a cylinder and the surface area of the Earth has to be calculated.
The future use of the proposed modelling approaches within global climate models will enable to explicitly model space activity within the models and to study the possible impact of this activity on our planet.