KR102422718B1 - Method for Calculating Mass of Halo Coronal Mass Ejection based on Synthetic Coronal Mass Ejection and Apparatus thereof - Google Patents

Abstract

A method for calculating the mass of synthetic coronal mass ejection (CME)-based halo CME and an apparatus thereof are disclosed. In the method, first, a three-dimensional factor is calculated from observation data of a single satellite using a full ice cream cone model, and then an observation density distribution at a position angle of a leading edge with a fastest projection speed in the observation data is calculated. Next, a synthetic CME is generated based on the full ice cream cone structure to which the 3-dimensional factor is applied and a power-law distribution following the electron density of a CME structure, and a synthetic CME density distribution is calculated at the position angle of the leading edge in the generated synthetic CME. Then, the observed density distribution and the synthetic CME density distribution are compared to determine a power-function density distribution of the synthetic CME when an error value is a minimum value, and the mass of the synthetic CME is calculated using the power-function density distribution of the synthetic CME. The present invention can more accurately calculate the mass of the Halo CME based only on brightness data measured from a single satellite.

Description

Method for Calculating Mass of Halo Coronal Mass Ejection based on Synthetic Coronal Mass Ejection and Apparatus thereof

The present invention relates to a method and apparatus for calculating the mass of a synthetic CME-based halo CME.

Coronal Mass Ejection (CME) is a very explosive phenomenon in the Sun, well known as the cause of geomagnetic storms.

The white light coronagraph used for these CME observations measures the brightness due to Thomson scattering of light by electrons distributed along the line of sight.

As such, the mass of the CME can be calculated if the brightness due to Thomson scattering and the distribution of electrons along the line-of-sight direction are known. However, in the case of white light coronagraph observation, the distribution of electrons along the gaze direction cannot be known due to the projection effect, so the exact mass of the CME cannot be measured or an underestimated value is obtained from the actual mass. In particular, in the case of halo CMEs propagating toward the Earth, the observed mass is even more underestimated because it is more affected by projection effects.

On the other hand, if CME observation is possible at multiple points at the same time, the mass of CME can be measured more accurately using the brightness measured at different locations. However, in the case of two satellites observing the sun at different points, There is a problem in that it is difficult to accurately measure the mass because simultaneous observation data cannot be continuously obtained due to the difference in .

An object of the present invention is to provide a method and apparatus for calculating the mass of a halo CME based on synthetic CME, which can more accurately calculate the mass of the halo CME based on only brightness data measured from a single satellite.

A characteristic configuration of the present invention for achieving the above-described object of the present invention and realizing the characteristic effects of the present invention described later is as follows.

According to one aspect of the present invention, there is provided a synthetic CME-based halo CME mass calculation method, the method comprising:

A halo coronal mass ejection (CME) mass ejection device calculates the mass of halo CMEs based on synthetic CMEs. It uses a full ice-cream cone model in which the CME structure is a flat cone combined with a hemisphere. Calculating a three-dimensional factor from observation data of a satellite, calculating an observation density distribution at a position angle of a leading edge having the fastest projection speed in the observation data, a full ice cream cone to which the three-dimensional factor is applied generating a composite CME based on a power-law distribution followed by a structure and the electron density of the CME structure, calculating a composite CME density distribution at the position angle of a leading edge in the composite CME, the determining a power density distribution of the synthetic CME when an error value is a minimum by comparing the observed density distribution with the synthetic CME density distribution, and using the power density distribution of the synthetic CME to determine the mass of the synthetic CME It includes the step of calculating.

According to another aspect of the present invention, there is provided a synthetic CME-based halo CME mass calculating device, the device comprising:

a processor, a memory, a communicator and an input/output device, wherein the communication device communicates with an external device, the input/output device displays information to the outside or receives information or commands input from the outside, the memory is a set of codes A process for calculating a three-dimensional factor from observation data of a single satellite input through the communicator or the input/output device, using a full ice cream cone model, the code having the fastest projection speed in the observation data A process of calculating an observed density distribution at a position angle of a leading edge, a process of generating a synthetic CME based on a power function distribution following a full ice cream cone structure to which the three-dimensional factor is applied and an electron density of the CME structure, the synthetic CME a process of calculating a composite CME density distribution at a position angle of a leading edge in a process of determining a power function density distribution of the composite CME when an error value is a minimum value by comparing the observed density distribution and the composite CME density distribution; and controlling the processor to execute a process of calculating the mass of the synthetic CME by using the power function density distribution of the synthetic CME.

According to the present invention, the mass of the halo CME can be more accurately calculated based on only the brightness data measured from a single satellite by generating the synthetic CME.

Due to this, it is possible to accurately measure the mass and energy of the halo CME, so that better research results can be obtained on the relationship between the CME and other elements of the space environment.

In addition, it can be utilized for the development and improvement of physical models using the mass of the CME and the energy calculated using it as an input factor.

1 is a schematic flowchart of a method for calculating the mass of a synthetic CME-based halo CME according to an embodiment of the present invention.
2 is a diagram illustrating the geometric structure of a full ice cream cone model that is a CME structure according to an embodiment of the present invention.
3 is a diagram illustrating a schematic configuration of an apparatus for calculating a mass of a synthetic CME-based halo CME according to an embodiment of the present invention.
4 is a diagram illustrating a basic difference image of CME observed on September 10, 2014 and a corresponding composite CME image.
5 (a) and (b) are diagrams showing the density distribution in MPA of the observed and synthetic CMEs together with the projected heights on November 17, 2001 and September 10, 2014, respectively.
6 (a) and (b) are plots showing the observed mass and two synthetic CME masses for event 1 and event 2, respectively, as a function of radius height.

Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those of ordinary skill in the art to which the present invention pertains can easily implement them. However, the present invention may be embodied in several different forms and is not limited to the embodiments described herein. And in order to clearly explain the present invention in the drawings, parts irrelevant to the description are omitted, and similar reference numerals are attached to similar parts throughout the specification.

Throughout the specification, when a part "includes" a certain component, it means that other components may be further included, rather than excluding other components, unless otherwise stated. In addition, terms such as “…unit”, “…group”, and “module” described in the specification mean a unit that processes at least one function or operation, which may be implemented as hardware or software or a combination of hardware and software. have.

It should be understood that throughout the specification, the term "and/or" describes only an association relationship for describing a related object and indicates that three relationships may exist. For example, A and/or B may represent three cases in which only A exists, both A and B exist, and only B exists. Also, throughout the specification, the character "/" generally indicates an "or" relationship between related objects.

Hereinafter, a method for calculating the mass of a synthetic CME-based halo CME according to an embodiment of the present invention will be described with reference to the drawings.

(여기서,

는 상수임)를 따른다는 것이다.Before the description, two assumptions are preset in order to calculate the mass of the halo CME based on the synthetic CME according to an embodiment of the present invention. The first assumption is that the structure of the CME is a hemispherical ice cream cone (or full ice-cream cone model) shape, and the second assumption is that the electron density of the CME has a power-law distribution, for example

(here,

is a constant).

1 is a schematic flowchart of a method for calculating the mass of a synthetic CME-based halo CME according to an embodiment of the present invention.

Referring to FIG. 1 , first, a three-dimensional factor of the CME is calculated using a full ice cream cone model assuming that the CME structure is a full ice cream cone that is a flat cone combined with a hemisphere ( S100 ).

Content disclosed in Korean Patent Registration No. 1786361 may be referred to for calculating the three-dimensional factor of CME using a model in which the CME structure is a full ice cream cone in an embodiment of the present invention. Specifically, observed projection speeds are calculated from height and time data in a plurality of azimuthal directions of the halo CME obtained from single satellite observation data, and emission speeds are calculated from the observed projection speeds. (radial velocity), angular width (angular width) and the initial value of the three-dimensional factor including the source location (source location) is obtained. Thereafter, a full ice cream cone model is constructed using the initial values of the obtained three-dimensional factors and projected on the sky plane, and measured at the outer boundary of the full ice cream cone model projected on the sky plane. After calculating the estimated projection speeds, the calculated observed projection speed and the measured projection speed are compared to calculate the 3D factor of the halo CME when the error value is the minimum value.

Referring to the contents disclosed in the aforementioned Korean Patent No. 1786361, the CME structure of the full ice cream cone model assumed in the embodiment of the present invention is as shown in FIG. 2 .

축은 지구와 2차원 평면(

,

) 방면을 나타내며, 콘 좌표계에서

축은 원뿔의 중심축에 놓여 있고, 콘의 단면에 평행한 면(

,

)을 나타낸다. 또한, r은 태양 핵의 중심(410)으로부터 CME의 가장자리(412)까지의 거리, 즉 반경 높이를 나타내고,

및

는 헬리오센트릭 좌표계에서의 여위도(colatitude) 및 경도(longitude)를 의미하며,

는 각도 폭(angular width)을 의미하고,

는 평면과 원뿔의 중심축 사이의 각도(

)를 의미한다.2 is a diagram illustrating the geometric structure of a full ice cream cone model that is a CME structure according to an embodiment of the present invention. Referring to FIG. 2 , a full ice cream cone model that is a CME structure according to an embodiment of the present invention is shown through a relationship between a heliocentric coor-dinate system and a cone coor-dinate system. . in a heliocentric coordinate system

The axis is the Earth and a two-dimensional plane (

,

) represents the direction, in the cone coordinate system

The axis lies on the central axis of the cone and is parallel to the cross section of the cone (

,

) is indicated. In addition, r represents the distance from the center 410 of the solar nucleus to the edge 412 of the CME, that is, the radial height,

and

means latitude and longitude in the heliocentric coordinate system,

means the angular width,

is the angle between the plane and the central axis of the cone (

) means

), 각도 폭(

), 여위도(

) 및 경도(

)를 포함한다.Therefore, in the above-described step S100, the three-dimensional factors of the CME calculated using the full ice cream cone model are the radial height (r), the radial velocity (

), angle width (

), latitude (

) and hardness (

) is included.

)을 측정한다(S110). Next, the observed mass of the CME ((

) is measured (S110).

The CME region projected on the sky plane may be defined as the projection of the full ice cream cone constructed by the three-dimensional factor obtained in step S100 described above. These regions can also be used to generate synthetic CME images corresponding to observations. Since the observed mass in this region can be measured in a general way based on the assumption that all electrons are in the sky plane, a detailed description is omitted here.

Next, the density distribution is calculated from the observation from the portion where the projection speed of the CME is the fastest, that is, the measurement position angle (MPA), which is the position angle of the leading edge (S120).

)를 측정한다. 여기서, MPA에서의 밀도 분포와 부피는 다음의 [수학식 1]에 의해 주어진다.More specifically, from the observation, the density distribution (

) is measured. Here, the density distribution and volume in MPA are given by the following [Equation 1].

는 MPA의 질량이고(여기서

는 미리 주어지는 것인지 아니면 S110에서 산출된 관측 질량(

)으로부터 산출되는 것인지 확인부탁드립니다),

은 MPA의 높이이며,

는 CME 물질이 90% 수소와 10% 헬륨(전자 당 1.97 × 10^-24g)을 포함한다는 가정에서 전자 당 유효 질량이고,

은 MPA의 부피로서,

이며, 여기서

는 가림판인 오컬터(occulter)의 높이이고,

은 MPA의 각도 폭이다. here,

is the mass of MPA (where

Whether is given in advance or the observed mass calculated in S110 (

), please check whether it is calculated from

is the height of MPA,

is the effective mass per electron, assuming that the CME material contains 90% hydrogen and 10% helium (1.97 × 10 ^-24 g per electron),

is the volume of MPA,

and where

is the height of the occulter, the screen,

is the angular width of MPA.

로부터 멱함수 지수인

를 획득할 수 있다.In this case, the power function distribution of the density

from the power exponent

can be obtained.

에 기초하여 합성 CME를 생성한다(S130).Then, the above-mentioned full ice cream cone structure and the above-mentioned power function distribution

Synthetic CME is generated based on (S130).

은 다음의 [수학식 2]에 의해 추정된다.The overall brightness of the pixels in the synthetic CME is

is estimated by the following [Equation 2].

는 Billings(1966)에 의해 사용된 톰슨 단면이고,

은 태양 중심에서 시선(line of sight, LOS) 상의 산란 위치까지의 거리이며,

는 산란 위치에서 관찰자까지의 거리이고,

는 해당 위치의 전자 수 밀도로서

이며,

는 해당 위치에서 단일 전자의 밝기이다.here,

is the Thomson cross-section used by Billings (1966),

is the distance from the center of the sun to the scattering location on the line of sight (LOS),

is the distance from the scattering position to the observer,

is the electron number density at that position

is,

is the brightness of a single electron at that location.

The LOS boundary is determined by the intersection of the full ice cream cone and the line from the pixel to the viewer. Thereafter, the total brightness of all pixels in the CME region defined in step S110 is estimated. A synthetic CME can be generated by changing the power exponent from n _cme = n _os to n _cme = n _os + 0.1×i _max . As an example, i _max is 15, which may be determined by multiple trials.

Next, the density distribution at the position angle of the leading edge having the fastest projection speed is calculated from the generated synthetic CME (S140).

)는 다음의 [수학식 3]에 의해 산출될 수 있다.Specifically, the density distribution at the position angle of the leading edge with the fastest projection velocity from the synthetic CME (

) can be calculated by the following [Equation 3].

는 합성 CME에서의 MPA의 질량으로, 합성 CME에서 측정된 총 밝기와 모든 전자가 스카이 평면에 있다는 가정에서 추정될 수 있으며,

는 합성 CME에서의 MPA의 부피로서, 이는 전술한

과 동일하다.here,

is the mass of MPA in the synthetic CME, which can be estimated from the total brightness measured in the synthetic CME and assuming that all electrons are in the sky plane,

is the volume of MPA in synthetic CME, which is

same as

Subsequently, the CME power function in which the root mean square (RMS) error is minimized by comparing the density distribution in the observation calculated in step S120 with the density distribution calculated in the synthetic CME in step S140 . The density distribution is determined (S150).

)와 상기 단계(S140)에서 산출된 합성 CME에서의 밀도 분포(

) 사이의 RMS 오차가 추정된다. 이러한 RMS 오차는

(여기서

는 태양 반경임)인 밀도부터 계산되며, 그 후, RMS 오차가 최소화되는 n_cme가 결정된다.Specifically, the density distribution in the observation calculated in step S120 (

) and the density distribution in the synthetic CME calculated in step (S140) (

), the RMS error between This RMS error is

(here

is the solar radius) is calculated from the density, and then n _cme for which the RMS error is minimized is determined.

Next, the mass of the CME is calculated using the power density distribution of the CME determined in step S150 and the full ice cream cone structure (S160).

Specifically, the mass (M _cme ) of the CME can be calculated from the following [Equation 4] using n _cme determined in the step (S150).

이다.here,

to be.

As such, according to an embodiment of the present invention, the mass of the halo CME can be more accurately calculated based on only the brightness data measured from a single satellite by generating the synthetic CME.

Using the steps of the above-described synthetic CME-based halo CME mass calculation method, the following two kinds of synthetic CME masses can be estimated. (1) a halo CME mass that considers only the observed CME region (M _cme1 ) and (2) a halo CME mass that includes both the observed and occluded regions (M _cme2 ). The mass of the occluded region (hereafter occluded mass) is estimated by the density distribution determined in the above step, but only regions with a radial distance greater than 4R _⊙ are considered. When the limb CME is observed in the LASCO (Large Angle and Spectrometric Coronagraph) C3 field of view, the CME structure can be identified through the shielding plate of the LASCO C3. The lower part of the CME is not yet known, and from this fact, a region in which r > 4R _⊙ is considered when calculating M _cme2 through Equation 4 above.

On the other hand, the conversion factor (C _o ), which is the ratio of the synthetic CME mass (M _cme ) calculated in the aforementioned step S160 to the observed mass ( _{Mo obs} ) of the CME calculated in the aforementioned step S110, is the following [Equation 5] can be derived as

In the case of such a conversion factor, two conversion factors may be derived as follows in response to the two types of synthetic CME masses described above. (1) C _o1 is the ratio of the mass of the synthetic CME without obscured mass to the observed mass and (2) C _o2 is the ratio of the mass of the synthetic CME with the obscured mass to the observed mass.

Hereinafter, an apparatus for calculating the mass of a synthetic CME-based halo CME according to an embodiment of the present invention will be described with reference to the drawings.

3 is a diagram illustrating a schematic configuration of an apparatus for calculating a mass of a synthetic CME-based halo CME according to an embodiment of the present invention.

As shown in FIG. 3 , the apparatus 100 for calculating the mass of a synthetic CME-based halo CME according to an embodiment of the present invention includes at least one processor 110 , a memory 120 , a communicator 130 , and an input/output device 140 . ) and a communication bus 150 .

The processor 110 may be a general-purpose central processing unit (CPU), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits for controlling program execution in the solution of the present application.

The memory 120 stores information related to mass calculation of a synthetic CME-based halo CME according to an embodiment of the present invention.

Specifically, the memory 120 is further configured to store a set of codes, and the codes are used to control the processor 110 to execute a process as follows. This process is a process of calculating the three-dimensional factor of the CME, which is a full ice cream cone model, based on single satellite observation data input through the communicator 130 or the input/output device 140, and the entire halo CME obtained from the single satellite observation data. The process of measuring the observed mass of the CME using the brightness and the assumption that all electrons of the CME are located on the sky plane, the process of calculating the density distribution in MPA, which is the position angle of the leading edge of the CME from the single satellite observation data, pool, The process of generating a synthetic CME based on the ice cream cone structure and power function distribution; This includes the process of determining the power density distribution of the CME in which the RMS error is minimized by comparing the calculated density distributions, and the process of calculating the mass of the CME using the determined power density distribution of the CME and the full ice cream cone structure.

Here, the mass of CME includes a halo CME mass that considers only the observed CME region, or a halo CME mass that includes both the observed CME region and the obscured region.

Memory 120 can be read-only memory (ROM) or other type of static storage that can store instructions, or random access memory (RAM) or other type of dynamic storage that can store information and instructions; or Electronically Erasable Programmable Read-Only Memory (EEPROM), Compact Disc Read-Only Memory (CD-ROM), or other compact disc storage device or optical disc storage device (compressed optical disc, laser disc, optical disc, digital versatile disc, blue ray disk, etc.), a magnetic disk storage medium or other magnetic storage device, or any other medium that can carry or store expected program code in the form of instructions or data structures and that can be accessed by a computer; This is not limited. The memory 120 may exist independently and is coupled to the processor 110 by a communication bus 150 .

The communicator 130 communicates with other devices or communication networks, and may be implemented using various communication technologies. That is, WCDMA (Wideband CDMA), HSDPA (High Speed Downlink Packet Access), HSUPA (High Speed Uplink Packet Access), HSPA (High Speed Packet Access), Mobile WiMAX, WiBro, LTE (Long Term) Evolution), Bluetooth, infrared data association (IrDA), near field communication (NFC), Zigbee, wireless LAN technology, USB (Universal Serial Bus), etc. may be applied. In addition, when a service is provided by being connected to the Internet, TCP/IP, which is a standard protocol for information transmission on the Internet, may be followed.

The input/output device 140 specifically includes an input device 141 and an output device 142 , and the input device 141 may communicate with the processor 110 and receive user input in a plurality of ways. . For example, the input device 141 may be a mouse, a keyboard, a touch screen, or a sensing device. The output device 142 may communicate with the processor 110 and display information or output voice in a plurality of ways. For example, the output device 142 may be a liquid crystal display (LCD), a light emitting diode (LED) display, an organic light emitting diode (OLED) display, a speaker, or the like.

The communication bus 150 is configured to couple all the components of the mass calculating device 100 of the synthetic CME-based Halo CME, that is, the processor 110 , the memory 120 , the communicator 130 , and the input/output device 140 .

Example

From January 2000 to September 2014, the aforementioned synthetic CME-based halo CME mass calculation method was applied to 56 halo CMEs. A synthetic CME was generated based on the full ice cream cone structure generated using three-dimensional factors, and Fig. 4 shows the basic difference image of the CME observed on September 10, 2014 and the corresponding synthetic CME image. 4, (a) is a basic difference image of CME observed on September 10, 2014 by SOHO (Solar and Heliospheric Observatory)/LASCO C3, (b) is generated according to an embodiment of the present invention The corresponding synthetic CME images are shown. As shown in FIG. 4 , it can be seen that the composite CME image is very consistent with the observed images.

To determine the electron number density distribution of the synthetic CME corresponding to the observation, the RMS error between the density distribution of the MPA from the observation and the synthetic CME is minimized. 5(a) and (b) show the observed and synthesized CMEs with projected heights for events on November 17, 2001 (hereafter event 1) and September 10, 2014 (hereafter event 2), respectively. The density distribution in MPA is shown.

The density distribution of MPAs in the observations follows a power function approximately with an exponent (n _os ) of 2.7 for event 1 and 2.4 for event 2. For all 56 events, the power exponent of density at MPA for the observations, n _os , ranges from 1.5 to 3.6 with an average of about 2.5. It was found that the density distribution in MPA from synthetic CME was very similar to the observation, as can be seen in Fig. 5(a) and (b). This implies that the power function assumption for CME density is sufficient to generate synthetic CMEs corresponding to the observations. From this minimization, it was confirmed that the electron number density of the synthetic CME followed a power function distribution with an exponent (n _cme ) equal to 3.3 for event 1 and 2.6 for event 2, respectively. On average, the density distribution of synthetic CMEs for all 56 events follows a power function with an exponent (n _cme ) of 2.9.

6(a) and (b) show the observed mass (M ^c _obs ) and the two synthetic CME masses (synthetic CME mass without occluded mass (M ^c ) for event 1 and event 2, respectively, as a function of radius height. _cme1 ) and synthetic CME masses with occluded masses (M ^c _cme2 )). For M ^c _cme2 , the masked mass does not include all the masses behind the blanking plate because we consider the region where r > 4R _⊙ . Here, it was confirmed that the synthetic CME mass based on the power density distribution converges to a constant value similar to the observed mass. From the observations, the final total mass M ^c _obs is obtained as 10 ^16.3 g for event 1 and 10 ^16.2 g for event 2. For event 1, the synthetic CME mass without occluded mass (M ^c _cme1 ) and the synthetic CME mass with occluded mass (M ^c _cme2 ) are 10 ^16.5 g and 10 ^16.6 g, respectively. For event 2, M ^c _cme1 and M ^c _cme2 are 10 ^16.6 g and 10 ^16.8 g, respectively.

For the two composite masses, a transformation factor C _o1 , which is the ratio of the mass of the composite CME without obscured mass to the observed mass, and C _o2 , which is the ratio of the mass of the composite CME with an obscured mass to the observed mass, is determined. For event 1, C _o1 is 1.6 and C _o2 is 2.1. For event 2, C _o1 and C _o2 are 2.5 and 3.6, respectively.

For all 56 events, the transformation factor C _o1 ranges from 1.4 to 3.0 with a mean of 2.0. The factor C _o2 ranges from 1.8 to 5.0 with a mean of 3.0. These results imply that the observed halo CME mass can be underestimated by about a factor of 2 when considering the observed CME area and by about 3 times when considering the area including the occluded area.

The embodiment of the present invention described above is not implemented only through the apparatus and method, and may be implemented through a program for realizing a function corresponding to the configuration of the embodiment of the present invention or a recording medium in which the program is recorded.

Although the embodiments of the present invention have been described in detail above, the scope of the present invention is not limited thereto. is within the scope of the right.

Claims (15)

A method for a halo Coronal Mass Ejection (CME) mass calculation device to calculate the mass of a halo CME based on a synthetic CME, comprising:
Calculating three-dimensional factors from observation data of a single satellite using a full ice-cream cone model, in which the CME structure is a flat cone combined with a hemisphere;
calculating an observation density distribution at a position angle of a leading edge having the fastest projection speed in the observation data;
generating a synthetic CME based on a power-law distribution following the full ice cream cone structure to which the three-dimensional factor is applied and the electron density of the CME structure;
calculating a composite CME density distribution at a position angle of a leading edge in the composite CME;
determining a power function density distribution of the synthesized CME when an error value is a minimum value by comparing the observed density distribution with the synthesized CME density distribution, and
calculating the mass of the synthetic CME using a power density distribution of the synthetic CME;
A synthetic CME-based halo CME mass calculation method comprising a. According to claim 1,
The step of calculating the three-dimensional factor,
calculating observed projection speeds from height and time data in a plurality of azimuthal directions of the halo CME obtained from the observation data;
obtaining initial values of three-dimensional parameters including a radial velocity, an angular width and a source location from the observed projection velocity;
constructing the full ice cream cone model using the initial value of the three-dimensional factor and projecting it on a sky plane;
calculating estimated projection speeds at an outer boundary of a full ice cream cone model projected on the sky plane; and
Comparing the observed projection velocity and the measured projection velocity, calculating a three-dimensional factor of the halo CME when an error value is a minimum value;
Including, synthetic CME-based halo CME mass calculation method. According to claim 1,
The three-dimensional factor is the radial height (r), the radial velocity (

), angle width (

), latitude (

) and hardness (

) containing,
Synthetic CME-based Halo CME Mass Calculation Method. 4. The method of claim 3,
In the step of calculating the observed density distribution, the density distribution at a position angle corresponding to the leading edge of the CME (

) is the relation below

is calculated by
here,

is the mass at the position angle, and

is the height at the position angle,

is the effective mass per electron,

is the volume at the position angle,

and where

is the height of the occulter, the screen,

is the angular width at the position angle,
Synthetic CME-based Halo CME Mass Calculation Method. 5. The method of claim 4,
The power function distribution is the following relation

is given by
here,

is the electron number density, n _cme is the power exponent in the synthetic CME,

is a constant,
Synthetic CME-based Halo CME Mass Calculation Method. 6. The method of claim 5,
The total brightness of the pixels of the synthetic CME (

) is the following relation

is estimated by
here,

is the Thomson cross section,

is the distance from the center of the sun to the scattering location on the line of sight (LOS),

is the distance from the scattering position to the observer,

is the brightness of a single electron,
Synthetic CME-based Halo CME Mass Calculation Method. 6. The method of claim 5,
In the step of calculating the synthetic CME density distribution, the density distribution of the synthetic CME (

) is the following relation

is calculated by
here,

is the mass at the position angle of the synthetic CME,

is the volume at the position angle of the synthetic CME,
Synthetic CME-based Halo CME Mass Calculation Method. 8. The method of claim 7,
In the step of determining the power density distribution of the synthetic CME,
The error value is

(here

is the solar radius), calculated from the density that satisfies the
Synthetic CME-based Halo CME Mass Calculation Method. 8. The method of claim 7,
The mass (M _cme ) of the synthetic CME is the following relation

calculated according to
Synthetic CME-based Halo CME Mass Calculation Method. 10. The method according to any one of claims 1 to 9,
The mass of the synthetic CME is,
A halo CME mass that considers only the observed CME region and a halo CME mass that includes both the observed CME region and the region obscured by the blanking plate.
Including, synthetic CME-based halo CME mass calculation method. 11. The method of claim 10,
The halo CME mass including both the observed CME region and the region obscured by the shield is calculated by considering only the region where r > 4R _⊙ ,
Synthetic CME-based Halo CME Mass Calculation Method. A device for calculating the mass of halo CME based on synthetic CME, comprising:
including a processor, memory, communicator and input/output device;
The communicator performs communication with an external device,
The input/output device displays information to the outside, or receives information or commands input from the outside,
the memory is configured to store a set of codes;
The code is
A process of calculating a three-dimensional factor from the observation data of a single satellite input through the communicator or the input/output device using the full ice cream cone model, and the observation density distribution at the position angle of the leading edge with the fastest projection speed in the observation data A process for calculating a synthetic CME based on a power function distribution following the electron density of the full ice cream cone structure to which the three-dimensional factor is applied and the CME structure, the synthetic CME density at the position angle of the leading edge in the synthetic CME a process of calculating a distribution, a process of comparing the observed density distribution with the synthesized CME density distribution to determine a power function density distribution of the synthesized CME when an error value is a minimum value, and using the power function density distribution of the synthesized CME The process of calculating the mass of the synthetic CME by
used to control the processor to execute
Synthetic CME-based Halo CME Mass Calculator. 13. The method of claim 12,
The processor is
A process of measuring the observed mass of the CME assuming that the total brightness of the halo CME obtained from the observation data and all electrons of the CME are located on the sky plane
run more,
wherein the observed density distribution is calculated using the observed mass of the CME;
Synthetic CME-based Halo CME Mass Calculator. 13. The method of claim 12,
The power function distribution is the following relation

is given by
here,

is the electron number density, n _cme is the power exponent in the synthetic CME,

is a constant,
Synthetic CME-based Halo CME Mass Calculator. 13. The method of claim 12,
The mass of the synthetic CME is,
A halo CME mass that considers only the observed CME region and a halo CME mass that includes both the observed CME region and the region obscured by the blanking plate.
Including, synthetic CME-based halo CME mass calculation device.

Images