Abstract
Introduction: On the propagation path to the satellite, the ionosphere will distort the nuclear electromagnetic pulse (NEMP) and change its physical properties.
Methods: This paper proposes a method for calculating the propagation of NEMP to the satellite. The method decomposes NEMP into the superposition of simple harmonic waves, and each simple harmonic wave is calculated separately in the ionosphere. With the consideration of different time of arrival and critical frequency of the ionosphere, the NEMP after propagating in the ionosphere is obtained by superposition of simple harmonic waves in time domain rather than the inverse Fourier transform which will erase the time domain information.
Results: The results show that NEMP is dispersive in ionosphere with the pulse broadened, the speeds changed and the bandwidth narrowed. The time-frequency spectrum can provide the frequency band where the signal energy is located.
Discussion: Our proposed method provides a simple and effective way to calculate the NEMP propagation in the ionosphere, which should afford help to the design of NEMP receivers and the selection of satellite orbit altitude.
1 Introduction
Nuclear electromagnetic pulse (NEMP) is one of the important effects of nuclear explosions, with fast propagation speed, high peak value, wide frequency band, wide coverage and other characteristics. NEMP is appropriate for use as a detection target for space-based nuclear explosion monitoring systems. As a consequence of the Comprehensive Nuclear Test Ban Treaty, non-nuclear experimental and theoretical efforts are expanded to study and observe the NEMP phenomenology and to develop appropriate descriptive models. For the non-nuclear experimental efforts, especially in space-based NEMP explosion monitoring, in 1997, the United States launched a near Earth satellite, FORTE (Fast On orbit Recording of Transient Events), with an optical and radio frequency detection system. The satellite is equipped with a 90 MHz receiver and a pair of 20 MHz receivers, both of which can be tuned in the range of 30–300 MHz. The FORTE received the pulse signal sent by the ground electromagnetic pulse generator (Massey et al., 1998). After band-pass filtering and converting it into a digital signal, it was transmitted back to the ground for analysis, so as to obtain the ionospheric information of the propagation path (Roussel-Dupre et al., 2001; Huang and Roussel-Dupre, 2005; Minter et al., 2007) and carry out the ground signal source positioning algorithm (Jacobson and Shao, 2001). With the above mentioned techniques, FORTE has the ability to detected the NEMP signal.
Theoretical methods for studying NEMP propagation can be classified into two categories: time-domain methods and frequency-domain methods. For the time-domain methods, the Maxwell equations are the most essential parts. The electric field components and magnetic field components are calculated with appropriate current sources, medium parameters and boundary conditions. The widely used numerical calculation method is Finite Difference Time Domain (FDTD) method. With this method, the physical process of electromagnetic pulse excited by Compton electron currents is simulated (Meng et al., 2003; Meng, 2013). The discrete iterative equations for the electromagnetic field components are presented in the CPML (Convolutional Perfect Matched Layer) media under the two-dimension and three-dimension prolate-spheroidal coordinate system (Gao et al., 2005a; Gao et al., 2005b). Chen et al. derived the time-domain difference equations for the interaction between ionosphere and electromagnetic pulse (Chen et al., 2011). Those results show that the combination of ionospheric time-domain calculation method and NEMP calculation method is reasonable and effective. However, FDTD method is suitable for the electromagnetic field calculation of subtle structures, and it takes a very long time for the calculation of long distance electromagnetic pulse propagation.
In frequency domain methods, the ionosphere is considered as the transfer function, and the Fourier transform is used to obtain the spectrum of NEMP signal. The amplitude and phase of each frequency component after propagation in the ionosphere are calculated respectively. Then, the inverse Fourier transform is used to obtain the NEMP signal received at the satellite orbit (Meng et al., 2004; Yao et al., 2019). Nevertheless, the Fourier transform will erase the time domain information of NEMP signal, and the inverse Fourier transform ignores that the different frequency components have different speeds in dispersive material.
This paper proposes a method for NEMP propagation in the ionosphere. The method decomposes the NEMP signal into the superposition of simple harmonic waves. The propagation of each simple harmonic wave in the ionosphere is calculated separately. With the consideration of the time differences among the simple harmonic waves, the NEMP signal is obtained by superposition of simple harmonic waves after propagation in the ionosphere. This paper is organized as follows: In Section 2, the model of ionosphere, decomposition of NEMP and superposition of simple harmonic waves are presented. In Section 3, the waveforms chosen for simulation are introduced and the results and discussion are presented. In the last section, we conclude this paper and present the applications.
2 Methods
2.1 Model of ionosphere
2.1.1 Multilayer model
The multilayer model is used for calculation, as a result of the non-uniform distribution of the plasma density in the ionosphere. Each layer is considered to be uniform and has a thickness of , as is shown in Figure 1. It is assumed that the layer is vacuum due to the negligible plasma density.
FIGURE 1
The plasma frequency of the ionosphere is , where is the charge of electron, is the mass of electron, is the vacuum permittivity, and is the electron density.
2.1.2 Critical frequency
increases with the increase of altitude, and after reaching the peak , starts to decrease. The plasma frequency corresponding to is the critical frequency, namely, . The frequency components below the critical frequency cannot continue to propagate upward.
2.1.3 Collision frequency
Particles in the ionosphere are in thermal motion and inevitably collide. Frequencies of electron-neutral and electron-ion collisions affect the ionosphere propagation of the signals (He and Zhao, 2009).
The collision frequencies of electron-neutral particle include.
(1) Collision frequency of electrons with Nitrogen molecules
(2) Collision frequency of electrons with Oxygen molecules
(3) Collision frequency of electrons with Oxygen atoms
(4) Collision frequency of electrons with Helium atoms
From Eqs 1–4, is the electron temperature (in eV, 1 eV = 11,600 K), , , , are the density of neutral particles , , , , respectively.
The collision frequency of electron-ion iswhere is the collision frequency of electrons with the ion i, which represents , , in this paper. represents the density of the ion i.
The total electron collision frequency can be expressed as
2.1.4 Refraction index
The propagation of electromagnetic waves in the ionosphere follows the Appleton-Lassen equation. In the absence of geomagnetic field (Huang, 2000), the refraction index can be expressed bywhere and are real and imaginary parts, respectively; ; ; is the frequency of the simple harmonic wave.
In vacuum, the speed of a simple harmonic wave is equal to the speed of light , namely, . In the ionosphere considering collisions, the phase speed can be expressed by (Liu et al., 2010) where, is the conductivity. is the real part of where The item is always great than 1, resulting in . According to Eq. 8, as the frequency increases, also increases, shown in Figure 2.
FIGURE 2
The propagation time in each layer of ionosphere is given by
In vacuum, the amplitude is a function of distance from the source pointwhere is the initial amplitude, is the amplitude at distance from the source point.
In the ionosphere, the attenuation of amplitude with distance can be expressed aswhere is the amplitude at distance from the source point; is the attenuation coefficient with frequency .
2.2 Dispersive propagation in the ionosphere
2.2.1 Decomposition of NEMP
In the light of Eqs 7–11, the influence of the ionosphere differ with the frequency . Therefore, we decompose NEMP into simple harmonic waves, and study the propagation of each simple harmonic wave in the ionosphere.
The Fourier transform is applied to and the signal is converted from the time domain to the frequency domain . The complex value corresponding to frequency is . The simple harmonic wave decomposed from can be described aswhere is the frequency, and () are the initial amplitude and phase, respectively. The superposition of all simple harmonics can be expressed aswhere is the simple harmonic wave with highest frequency , and is the simple harmonic wave with the lowest frequency .
2.2.2 A simple harmonic wave propagation in the ionosphere
2.2.2.1 Propagation from to
For a simple harmonic wave with frequency , when propagating vertically upward from to the bottom layer of the ionosphere, where there is no plasma, the propagation time is calculated byand the amplitude is calculated by
The expression is derived from Eq. 12 as
2.2.2.2 Propagation from to
When propagating from to , i.e., in layer, the refraction index of layer is expressed by
The propagation time is calculated by
The absorption loss can be derived from Eqs 10, 11where is the attenuation coefficient in layer, Eq. 16 is then expressed as
2.2.2.3 Propagation from to
By analogy, when the simple harmonic wave propagates to the height , the expression iswhere is the attenuation coefficient in layer, is the propagation time in layer, and is the propagation time from to .
2.2.2.4 Definition domain
Each simple harmonic wave has the same duration . The time difference between and , propagating from to , is given by
Equation 22 represents that when reaches , needs to reach it. Simple harmonic waves of different frequencies propagation in the ionosphere are shown in Figure 3.
FIGURE 3
Taking the moment when reaches as the zero moment, then the time domain of is defined as
In particular, the time domain of is defined as
2.2.3 Superposition of simple harmonic waves
The ionosphere will reflect the simple harmonic waves with frequencies below the critical frequency . The simple harmonic waves that can penetrate the ionosphere are , , …, , …, . is the last simple harmonic wave that can penetrate the ionosphere. , …, will be reflected by the ionospheric electron frequency peak (foF2).
With the consideration of the critical frequency and the different definition domain, the simple harmonic waves , , …, are superimposed, so NEMP signal that propagates from vertically upward to is finally described as
3 Simulation
The NEMP is emitted at the ground and is received at an altitude of 800 km and let km. The data of electron density and ions density are downloaded from the International Reference Ionosphere (International Reference Ionosphere - IRI (2016) | CCMC (nasa.gov)). The data of neutral particle density are downloaded from the NRLMSISE-00 Atm Model (NRLMSISE-00 Atm Model | CCMC (nasa.gov)). The electron density of the ionosphere as a function of altitude is plotted in Figure 1. The parameters of ionosphere are shown in Table 1.
TABLE 1
| Date and time | Latitude | Longitude | Top height | Critical frequency |
|---|---|---|---|---|
| 12:00 21 March 2020 | 40.2°N | 115.4°E | 800km | 6.24MHz |
Parameters of ionosphere.
3.1 Simulation waveform
We choose Bell Laboratory pulse as the simulation waveform. The first reason is that the method proposed in this paper is to study the dispersion of NEMP in the ionosphere, rather than the generation of NEMP. Therefore, the simulation waveform can be adjusted and replaced to meet the needs of the simulation. The second reason is that it is easier to be simulated, which facilitates post-calculation analysis.
Bell Laboratory pulse is represented by the double exponential functionwhere is the maximum amplitude, is the modifying factor, and are the parameters characterizing both the rise time and fall time simultaneously, is the peak time.
Parameters of Bell Laboratory pulse are shown in Table 2. The original waveform and superposition waveform which is obtained from Eq. 13, are shown in Figure 4. The relative error is less than 1 × 10−11%. Compared to the order of original waveform, the error is negligible.
TABLE 2
| α/s-1 | β/s-1 | tr/ns | tw/ns | |
|---|---|---|---|---|
| Bell Laboratory | 4,000,000 | 476,000,000 | 4.1 | 184 |
Parameters of Bell Laboratory pulse.
FIGURE 4
3.2 Results and discussion
Let represents the Bell Laboratory pulse after propagation in the ionosphere. becomes an oscillatory waveform with a long tail after propagation in the ionosphere, and is stretched from 1 to 30 μs. The rising edge of the envelope becomes slower, and the duration of the rising edge and the falling edge are almost the same. In addition, there is a peak time delay of 10 μs for , shown in Figure 5.
FIGURE 5
The delay can be explained with the simple harmonic waves’ arrival time differences and amplitudes, shown in Table 3. From 0 to 6 μs, the frequencies of the simple harmonic waves are all above 50 MHz, but their amplitudes are small, so the envelope of is flat. The envelope begins to rise at 6 μs, corresponding to the arrival time of a 50 MHz harmonic wave. Simple harmonic waves with lower frequencies but larger amplitudes then arrive and superimposed one after another. The energy of the pulse accumulates and reaches the peak at 10 μs. After 10 μs, although the amplitude of the arriving simple harmonic waves is larger than those that have arrived before 10 μs, the time difference of arrival is also increasing, so that the overlap between adjacent simple harmonic waves is constantly reduced. The energy of the pulse fails to accumulate, and the envelope begins to fall.
TABLE 3
| Frequency/MHz | Bell laboratory Amplitude/μV | Time difference/s |
|---|---|---|
| 500 | 0.86 | 0 |
| 450 | 0.88 | 1.5361 × 10−8 |
| 400 | 0.95 | 3.6837 × 10−8 |
| 350 | 1.09 | 6.8161 × 10−8 |
| 300 | 1.31 | 1.1642 × 10−7 |
| 250 | 1.74 | 1.9660 × 10−7 |
| 200 | 2.46 | 3.4409 × 10−7 |
| 150 | 3.92 | 6.6290 × 10−7 |
| 100 | 7.39 | 1.5750 × 10−6 |
| 50 | 17.48 | 6.5323 × 10−6 |
| 45 | 19.17 | 8.0962 × 10−6 |
| 40 | 20.89 | 1.0292 × 10−5 |
| 35 | 22.47 | 1.3517 × 10−5 |
| 30 | 23.47 | 1.8534 × 10−5 |
| 25 | 22.97 | 2.6988 × 10−5 |
| 10 | 0.69 | 2.0598 × 10−4 |
Amplitudes and time differences of simple harmonic waves.
Because of the critical frequency, the simple harmonic waves below it will be reflected. The critical frequency in the ionosphere is 6.24 MHz, and there are no amplitudes below the critical frequency, which are shown in Figure 6. It should be noted that Figure 6 is the one-to-one correspondence between simple harmonic waves and their amplitudes after propagation in the ionosphere. The frequency spectrum has the same shape as the original waveform.
FIGURE 6
The spectrum of the time-frequency of shows that the simple harmonic waves with higher frequencies arrive earlier, which is consistent with the time-frequency spectrum of the signal observed by the satellite (Massey et al., 1998). The brightest area in Figure 7 represents the frequency band with the highest concentration of energy and corresponds to the peak of , too.
FIGURE 7
Comparing the above results with other scholars’ great work (Liang and Zheng, 2009; Chen et al., 2011), our proposed method can explain the changes of waveforms after propagation in the ionosphere, such as the delay of the peak time, in light of the arrival time differences and amplitudes of simple harmonic waves. Besides, the application of inverse Fourier transform requires that the frequency components have the same wave front, in other words, the frequency components have the same propagation speed. Therefore, inverse Fourier transform is not suitable for the dispersive materials (Yao et al., 2019). Our proposed method is to superimpose the simple harmonic waves directly, which coincides with realistic physical process.
We have chosen different values of to calculate the peak value of . The results are shown in Table 4. The origin peak value of the pulse is 50,000 V/m, and peak value of the pulse after propagation in the ionosphere is in the order of mV/m. The order of amplitude of the calculation is what we care about most.
TABLE 4
| /km | 10 | 7 | 4 | 1 | 0.5 |
|---|---|---|---|---|---|
| /(mV/m) | 3.2 | 4.2 | 2.9 | 3.0 | 2.8 |
Peak value and thickness of each layer.
4 Conclusion
This paper proposes a method for calculating the propagation of NEMP in the ionosphere. Our method surpasses the weakness of the inverse Fourier transform that erases the propagation time difference for different frequencies. The results show that the NEMP signal is broadened as 30 times as its original length in the ionosphere. The spectrum of NEMP signal has the shape of the original signal envelope. The obtained dispersed spectrum is consistent with the satellite observations. The method can be used for the design of the space-based NEMP detectors and satellite-ground docking experiments. The propagation model may be optimized upon receipt of downlink data from the satellite.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
YLW proposed the model and the method, and wrote the paper. DHZ and ZXL ran the simulation. LHW and YW analyzed the results. TCZ and BX wrote the paper. BFC and PL did the literature review and checked the paper. All authors contributed to the article and approved the submitted version.
Funding
This paper is supported by the Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDA17040503).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
ChenY.MaL.ZhouH.WuW.LiJ.LiB. (2011). Calculation of high altitude nuclear electromagnetic pulse propagation in ionosphere. HIGH POWERLASER Part. BEAMS23 (02), 441–444. 10.3788/hplpb20112302.0441
2
GaoC.ChenY.WangL. (2005a). Application of CPML to two-dimension numerical simulation of nuclear electromagnetic pulse from air explosions. HIGH POWERLASER Part. BEAMS07, 1111–1116.
3
GaoC.ChenY.WangL. (2005b). The application of CPML in three-dimensional numerical simulation of low-altitude electromagnetic pulse. Chin. J. Comput. Phys.22 (4), 292–298. 10.19596/j.cnki.1001-246x.2005.04.002
4
HeF.ZhaoZ. (2009). Ionospheric loss of high frequency radio wave propagated in the ionospheric regions. Chin. J. Radio Sci.24 (04), 720–723+747.
5
HuangJ. (2000). Correction for atmospheric refractive error of radio wave. Beijing: National Defense Industry Press.
6
HuangZ.Roussel-DupreR. (2005). Total electron content (TEC) variability at Los Alamos, New Mexico: a comparative study: FORTE-derived TEC analysis. Radio Sci.40 (6), RS6007. 10.1029/2004rs003202
7
JacobsonA. R.ShaoX. (2001). Using geomagnetic birefringence to locate sources of impulsive, terrestrial VHF signals detected by satellites on orbit. Radio Sci.36 (04), 671–680. 10.1029/2000rs002555
8
LiangR.ZhengY. (2009). Calculation of nuclear electromagnetic pulse propagation into space. Nucl. Electronics& Detect. Technol.29 (06), 1393–1396.
9
LiuS.LiuS.HongW. (2010). Finite difference time domain method for dispersive media. Beijing: Science Press.
10
MasseyR. S.KnoxS. O.FranzR. C.HoldenD. N.RhodesC. T. (1998). Measurements of transionospheric radio propagation parameters using the FORTE satellite. Radio Sci.33 (6), 1739–1753. 10.1029/98rs02032
11
MengC. (2013). Numerical simulation of the HEMP environment. IEEE Trans. Electromagn. Compat.55(3), 440–445. 10.1109/temc.2013.2258024
12
MengC.ChenY.ZhouH. (2003). Effects of the HOB and the burst yield on the properties of NEMP. Chin. J. Comput. Phys.20(2), 173–177. 10.19596/j.cnki.1001-246x.2003.02.016
13
MengC.ChenY.ZhouH.GongJ. (2004). The analysis of the propagation of transient electromagnetic pulse through ionosphere. Nucl. Electron. Detect. Technol.24 (04), 369–372+406.
14
MinterC. F.RobertsonD. S.SpencerP. S. J.JacobsonA. R.Full-RowellT. J.Araujo-PradereE. A.et al (2007). A comparison of Magic and FORTE ionosphere measurements. Radio Sci.42 (3), RS3026. 10.1029/2006rs003460
15
Roussel-DupreR. A.JacobsonA. R.TriplettL. A. (2001). Analysis of FORTE data to extract ionospheric parameters. Radio Sci.36 (6), 1615–1630. 10.1029/2000rs002587
16
YaoJ.ShanH.CaoB.JiangT. (2019), “Waveform prediction of V/UHF EMP signal propagation in the ionosphere”, in Proceedings of the 2019 IEEE Asia-Pacific Microwave Conference (APMC), Singapore, December 2019. 1062–1064.
Summary
Keywords
nuclear electromagnetic pulse, ionosphere, closed-form expression, electromagnetic pulse propagation, dispersive propagation
Citation
Wei Y, Zhu D, Li Z, Wang L, Wang Y, Zhang T, Xing B, Cao B and Li P (2023) Dispersive propagation of nuclear electromagnetic pulse in the ionosphere. Front. Astron. Space Sci. 10:1201921. doi: 10.3389/fspas.2023.1201921
Received
07 April 2023
Accepted
23 October 2023
Published
09 November 2023
Volume
10 - 2023
Edited by
Chen Zhou, Wuhan University, China
Reviewed by
Sampad Kumar Panda, K L University, India
Xiang Xu, Wuhan University, China
Updates
Copyright
© 2023 Wei, Zhu, Li, Wang, Wang, Zhang, Xing, Cao and Li.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peng Li, lipeng@sklnbcpc.cn; Baofeng Cao, caobaofeng@sklnbcpc.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.