# Patent

## Title
System and method for reducing trapped energetic proton flux at low earth orbits

## Abstract

A system and method for reducing energetic proton flux trapped in the inner radiation belt by injecting Ultra Low Frequency (ULF) electromagnetic waves is disclosed. The ULF electromagnetic waves is generated by space or ground based transmitters and the frequency range is selected such that the injected waves are in gyrofrequency resonance with trapped 10 to 100 Mev protons. Pitch angle scattering of the trapped protons in gyro-resonance with the injected waves increases their precipitation rate by forcing their orbits into pitch angles inside the atmospheric loss-cone where they are lost by intaracting with the dense neutral atmosphere at altitudes below 100 km. The reduction of energetic proton flux trapped in the inner radiation belt allows use of commercial electronics with submicron feature size on Low Earth Orbit satellites and microsatellites without the operational constraints imposed by the presence of energetic proton fluxes trapped at the inner radiation belts.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This Application claims rights under 35 USC §119(e) from U.S. Application Ser. No. 61/448,480 filed Mar. 2, 2011, the contents of which are incorporated herein by reference.


TECHNICAL FIELD
Embodiments are generally related to radiation shielding. Embodiments also relate to shielding an electronic component against space environment radiations. Embodiments additionally relate to a system and method for reducing energetic proton flux trapped in the inner radiation belt by injecting Ultra Low Frequency (ULF) electromagnetic waves.
BACKGROUND OF THE INVENTION
The structure and behavior of the energetic electrons and protons trapped in Earth's Radiation Belt (RB) has been the subject of numerous experimental and theoretical studies. Morphologically, two regions are distinguished in an ionosphere such as an inner RB for L shells lower than two and an outer RB for L shells higher than two. The inner RB is dominated by protons with energy in excess of 10 MeV and lifetimes from a few years at low altitudes of 400 to 500 km to many tens of years at higher altitudes. Overall the inner belt energetic protons are relatively stable with a typical lifetime of ten years. Contrary to this, the outer RB is very dynamic and dominated by energetic electron fluxes associated with solar events and space weather process.
Earth's inner radiation belt located inside L=2 is dominated by a relatively stable flux of trapped protons with energy from a few to over 100 MeV. Radiation effects in spacecraft electronics caused by the inner radiation belt protons are the major cause of performance anomalies and lifetime of Low Earth Orbit satellites. For electronic components with large feature size, of the order of a micron, anomalies occur mainly when crossing the South Atlantic Anomaly (SAA). However, current and future commercial electronic systems are incorporating components with submicron size features. Such systems cannot function in the presence of the trapped 30 to 100 MeV protons, as hardening against such high-energy protons is essentially impractical.
Low Earth Orbiting (LEO) satellites spend a significant part of their orbit in the inner RB that is populated by energetic protons with energy, from one to more than one hundred MeV. The interaction of energetic protons with electronic devices of modern spacecraft results in high rates of anomalies due to Single Event Effects (SEE). Such anomalies range from nuisance effects that require operator intervention to debilitating effects leading to functional or total loss of the spacecraft. A set of operational problems occur when protons deposit enough charge in a small volume of silicon to change the state of memory cell, so that a one becomes zero and vice versa. The memories can become corrupted and lead to erroneous commands. Such soft errors are referred to as Single Event Upsets (SEU) and often generate high background counts to render the sensor unusable. Sometimes a single proton can upset more than one bit giving rise to Multiple Bit Upsets (MBU). Some devices can be triggered into a high current drain, leading to burn-out and hardware failure, known as single event latch-up or burn-out. Other devices suffer dielectric breakdown and rupture.
For LEO satellites, the dominant source of proton influence is the South Atlantic Anomaly (SAA). The SAA is a localized region at a fixed altitude, where protons in the inner RB reach their maximum intensity as a result of the asymmetry of the earth's magnetic field that can be approximated by a tilted, offset dipole in the inner magnetosphere. At present, satellites with micron size Commercial-Off-The Shelf (COTS) electronics experience serious effects mainly when transiting the SAA. For example, the intolerable frequency of SEU of the IBM 603 microprocessors (5 micron CMOS) in Iridium forced Motorola to disable the cache while transiting the SAA. Similar anomalies were experienced by the Hubble Space Telescope and numerous other satellites. To mitigate such effects, spacecraft utilize shielded electronic components that can reduce the flux of protons with energy below few MeV. However, it is very hard to shield against proton fluxes with energy in excess of 20 to 30 MeV.
The severity of the environment is usually expressed as an integral linear energy transfer spectrum, that represents the flux of particles depositing more than a certain amount of energy and charge per unit length of the material. This is referred as Linear Energy Transfer (LET), and given in units of MeV per g/cm2 or per mg/cm2. The effect on devices is characterized as a cross section (effective area presented to a beam), that is a function of the LET. The frequency of SEU caused by energetic protons is a non-linear function of the feature size. For large feature sizes, SEU are due to charge deposition caused by secondary particles with higher LET. For feature sizes smaller than 90 nm, direct proton ionization can cause SEU, resulting in an increase of the frequency of proton SEU by two or more orders of magnitude for deep submicron devices. This could preclude their use even for orbit latitudes different than the SAA. Further hardening the microelectronic components, besides the added weight, is very ineffective for proton energies higher than 20 to 30 MeV. For example, even one inch of Al reduces the 60-80 MeV flux by less than a factor of three. The recent tests have shown that the SEU cross section for energies between 1 to 10 MeV for bulk 65 nm Complementary Metal Oxide Semiconductor (CMOS) technology is by two orders of magnitude higher than for micron size devices, rendering current shielding level inadequate even at low proton energies.
Use of COTS in space applications is dictated by their high volume production and wide-spread use. The high volume production drives down their recurring component costs because of high yields and economies of scale. The wide-spread use of COTS reduces the system cost. Furthermore open standards drive down development and life-time support costs reduce the time to market for new products. The SEE issue for submicron CMOS or other electronic components presents a major dilemma, since it will prohibit use of COTS circuits with sub-micron size features and will limit the use of micro-satellites at LEO orbits.
Thus it is difficult to shield against 30 to 100 Mev protons to the level required by sub-micron features of current and future commercial electronic components. Heavy weight penalty must be paid to effect such shielding. Therefore, it is believed that a need exists for an improved system and method for reducing the energetic proton flux trapped in the inner radiation belt. Such system and method should allow the use of commercial electronics with submicron feature size on Low Earth Orbit (LEO) satellites and microsatellites without the operational constraints imposed by the presence of energetic proton fluxes trapped at the inner radiation belts.
BRIEF SUMMARY
The following summary is provided to facilitate an understanding of some of the innovative features unique to the disclosed embodiment and is not intended to be a full description. A full appreciation of the various aspects of the embodiments disclosed herein can be gained by taking the entire specification, claims, drawings, and abstract as a whole.
It is, therefore, one aspect of the disclosed embodiments to provide for radiation shielding.
It is another aspect of the disclosed embodiments to provide for shielding an electronic component against space environment radiations.
It is a further aspect of the present invention to provide for a system and method for reducing energetic proton flux trapped in the inner radiation belt by injecting Ultra Low Frequency (ULF) electromagnetic waves.
It is a another aspect of the present invention to provide for a system and method that allow the use of commercial electronics with submicron feature size on Low Earth Orbit (LEO) satellites and microsatellites without the operational constraints imposed by the presence of energetic proton fluxes trapped at the inner radiation belts.
It is a yet another aspect of the present invention to provide for a system and method for reducing energetic proton flux trapped in the inner radiation belt by injecting ULF electromagnetic waves into LEO and selecting ULF frequency range by ensuring that the injected waves are in gyrofrequency resonance with trapped 10 to 100 Mev protons. The ULF electromagnetic waves can be generated by space or ground based transmitters.
The aforementioned aspects and other objectives and advantages can now be achieved as described herein. According to the present invention, the energetic proton flux trapped in the inner radiation belt may be reduced by injecting Ultra-Low Frequency (ULF) electromagnetic waves, generated by space or ground based transmitters. The transmitted ULF frequency range is selected by the requirement that the injected waves are in gyrofrequency resonance with trapped 10 to 100 Mev protons.
Pitch angle scattering of the trapped protons in gyro-resonance with the injected waves increases their precipitation rate by forcing their orbits into pitch angles inside the atmospheric loss-cone and are lost by intaracting with the dense neutral atmosphere at altitudes below 100 km. Efficient techniques for generating and injecting the required ULF power include Horizontal Electric Dipole (HED) transmitters, Rotating Magnetic Fields (RMF) using arrays of permanent or supreconducting magnets and Transient Horizontal Electric Dipole Transmitters (THED). The proton flux reduction can be efficiently accomplished by using ground based arrays of permanent or superconducting magnets rotating at the selected ULF frequencies.
The present invention allow the use of COTS micro-electronic circuits with sub-micron features aboard LEO satellites and micro-satellites, reduce the current shielding weight and increase the useful lifetime of LEO satellites. The invention is based on the recognition that, the rate of SEU and other anomalies of electronic circuits aboard LEO satellites as well as the lifetime limitations are predominantly a function of the trapped proton flux in the 30 to 100 MeV energy range. The SEU and electronic circuit anomaly issue will be resolved by providing techniques that will reduce the trapped energetic proton flux encountered by LEO satellites.



BRIEF DESCRIPTION OF THE DRAWINGS
The accompanying figures, in which like reference numerals refer to identical or functionally-similar elements throughout the separate views and which are incorporated in and form a part of the specification, further illustrate the disclosed embodiments and, together with the detailed description of the invention, serve to explain the principles of the disclosed embodiments.
FIG. 1 illustrates a schematic diagram of the proton and electron RB structure as a function of L shell;
FIG. 2a illustrates a graph showing variation of omni-directional proton flux above 10 MeV and 50 MeV energies as a function of the L value;
FIG. 2b illustrates a graph showing differential spectrum of inner RB protons;
FIG. 3a illustrates a map showing a radiation flux at the SAA at altitude 400 km;
FIG. 3b illustrates a map showing geographical distribution of SEU in nMOS DRAms on UoSAT-2 showing clustering of proton events in the SAA;
FIG. 4 illustrates a graph showing B-L plot of the mean atmospheric density encountered by the trapped RB particles (atomic electrons/cm3), in accordance with the disclosed embodiments;
FIGS. 5a and 5b illustrate graphs showing time variation of 55 MeV proton flux and their redistribution caused by the Starfish nuclear test;
FIG. 6 illustrates a schematic drawing of the Proton Radiation Belt Remediation (PRBR) system, in accordance with the disclosed embodiments;
FIG. 7 illustrates a graph showing variation of minimum resonant SAW frequency and proton energy, under conditions typical of the magnetic equator of L=1.5, for two bounding pitch angles, in accordance with the disclosed embodiments;
FIG. 8 illustrates a graph showing the dependence of injected wave power on ionospheric reflection coefficient R to maintain wave energy W=75 kJ inside the shell volume at L=1.5 with δL=0, in accordance with the disclosed embodiments;
FIG. 9a illustrates an apparatus for injecting Shear Alfven Waves (SAW) to the ionosphere using Horizontal Electric Dipole (HED), in accordance with the disclosed embodiments;
FIG. 9b illustrates an apparatus for injecting SAW to the ionosphere using Rotating Magnetic Field (RMF) antennas, in accordance with the disclosed embodiments;
FIG. 10 illustrates a graph of amplitude of the electric field in the ionosphere versus time, showing a sneak-through concept of Transient Horizontal Electric Dipole (THED) Transmitters, in accordance with the disclosed embodiments;
FIG. 11a illustrates a graph showing ground response time as a function of antenna length and ground conductivity, in accordance with the disclosed embodiments;
FIG. 11b and FIG. 11c illustrate a graph showing relative efficiency as compared to an HED without reduction by ground return effects for rise times of T/3 and T respectively, in accordance with the disclosed embodiments; and
FIG. 12 illustrates a flow chart showing a method of reducing trapped energetic proton flux at LEO orbits, in accordance with the disclosed embodiments.



DETAILED DESCRIPTION
The particular values and configurations discussed in these non-limiting examples can be varied and are cited merely to illustrate at least one embodiment and are not intended to limit the scope thereof.
1. The Proton Radiation Belts
Radiation Belts (RB), also known as Van-Allen belts, are the locus of energetic electron and protons trapped by the earth's magnetic field. It is customary and convenient to describe a given magnetic field line by its L value. The L value corresponds to the radial location of its intersection with the magnetic equator in units of earth radius (RE). FIG. 1 illustrates a schematic diagram of the proton and electron RB structure 100 as a function of L shell around the earth 106. The intensity of a proton radiation belt 104 is confined inside the L=2 value. Contrary to this, the electron radiation belts 102 have two peaks one inside L=2 and another between L=4 and L=5, separated by a gap between L=2 and L=3. It is customary to label regions with L<2 as the inner RB and regions with L>2 as the outer RB.
The inner RB dominated by the presence of trapped energetic protons is shown in FIG. 1. A more detailed structure of the inner RB proton belt 104 of FIG. 1 is shown in FIGS. 2a and 2b. FIG. 2a illustrates a graph 200 showing variation of omni-directional proton flux above 10 MeV and 50 MeV energies as a function of the L value and FIG. 2b illustrates a graph 250 showing differential spectrum of inner RB protons. Note that for high proton energies, in excess of 20 to 30 MeV, the flux peaks at values of L=1.8 and is marked by an extremely sharp decrease for L>2. Furthermore, the energy spectrum is dominated by energies from 10 to 100 MeV and drops as a power law with energy.
FIG. 3a illustrates a map 300 showing a radiation flux in the SAA at altitude 400 km and FIG. 3b illustrates a map 350 showing a geographical distribution of SEU in nMOS DRAms on UoSAT-2 showing clustering of proton events in the SAA. Note that SAA is a localized region at a fixed altitude where protons in the inner RB protons reach their maximum intensity.
2. Lifetime of Energetic Protons Affecting Satellites at LEO Orbits
The large gradient of proton flux at the boundary between the inner and outer RB can be explained using the simplified “leaky-bucket” model. According to this model, the average proton flux at a particular L shell and energy is given by balancing the source of the energetic protons, such as Cosmic-Ray Neutron Albedo Decay (CRAND), to their loss to the atmosphere by processes such as inelastic nuclear collisions and slowing down by collisions with atomic oxygen at low altitude. The loss time T of a proton with energy E is controlled by the rate of energy degradation by collisional interactions with atomic oxygen and is given by the approximate formula,

T≈2×104(E/MeV)1.3(#/cm3/<p>)years  Equation 1

where <p> is the atomic electron density averaged over a proton orbit. The mean atmospheric density encountered by RB protons (atomic electrons per cm3) averaged over solar cycle in a B-L map 400 as computed by Cornwall et al. (1965) is shown in FIG. 4, where B is the magnetic field value at the particle mirror point. These results are consistent with particle lifetimes measured following the Starfish nuclear test for fifty five MeV protons trapped at L=1.4 with mirror points between two hundred and twenty and seven thirty kilometer measured by Filz and Coleman [1965] and shown in FIGS. 5a and 5b. Referring to L=1.8 and to the energy range thirty to hundred Mega electron Volt (MeV), the lifetime is from few years to few thousand years. FIGS. 5a and 5b illustrate graphs 500 and 550 showing time variation of fifty five MeV proton flux their redistribution caused by the Starfish test.
The dilemma occurs when similar considerations are applied to the outer RBs and compared with the results as shown in FIG. 2a, based on data provided by the National Space Science Data Center (NSSDC). As seen in the FIG. 2a, the flux maximizes in the inner RB, inside the L=2 shell, and drops precipitously at larger values of L. Furthermore, an examination of their temporal behavior indicates that the energetic proton flux in the inner RB is constant over times longer than tens to a hundred years. This is not the case for the outer RB, where the dynamic behavior is of the order of few hours. The dominant source of the trapped protons in the inner RB is due to what is known as Cosmic Ray Albedo Neutron Decay (GRAND). CRAND refers to trapping by the Earth's magnetic field of the protons produced by the decay of albedo neutrons, created by collisions of cosmic rays with atmospheric nuclei. The loss time is consistent with that given by Equation 1 and the data shown in FIGS. 5a and 5b. This loss time scale gives lifetimes in excess of thirty to fifty years for particles mirroring above thousand kilometer (km). The long lifetimes of the particles mirroring at high altitudes is related to the low density of the atomic nuclei and the low collisional pitch angle scattering into the loss cone. It is clear that increasing the rate of pitch angle of the protons trapped at high altitudes by artificial means will result in faster loss rate and lower trapped flux.
The concept resembles the work on controlled precipitation of relativistic electrons following an accidental or deliberate High Altitude Nuclear Explosions (HANEs), known as Radiation Belt Remediation (RBR). The main difference is that RBR is accomplished by injecting whistler waves with frequency in the kiloHertz (kHz) range and the remediation requires timescales of few days, while the Proton Radiation Belt Remediation (PRBR) requires SAW in the frequency range of one to ten Hertz (Hz) and remediation timescales of one to three year.
3. Proton Radiation Belt Remediation (PRBR)
A schematic drawing of the Proton Radiation Belt Remediation (PRBR) system 600 is illustrated in FIG. 6. ULF waves are injected from the ground to fill and maintain Shear Alfven Waves (SAW) with the proper frequency ω and amplitude B in the region to be remediated. The energetic protons interacting with the SAW will pitch angle scatter and forced to precipitate in the atmosphere. Designing such a PRBR system 600 requires the physics input such as frequency selection, pitch angle scattering rate and proton lifetime and energy-power requirements.
3.1 Frequency Selection
Neglecting relativistic effects and concentrating on the primary resonance, energetic protons interact with SAW when the Doppler shifted wave frequency ω seen in the reference frame of the energetic proton is equal to its gyro-frequency Q namely,

ω−kzvz=−Ω  Equation 2

In Equation 2, kz is the wave-number in the magnetic field direction. Assuming ω<<Q, and using the dispersion relation of SAW, with VA as the Alfven speed, as

ω=kzVA  Equation 3

The protons velocity v and pitch angle α, resonate with SAW when




  
    
      
        
          ω
          ⁡
          
            (
            
              v
              ,
              α
            
            )
          
        
        =
        
          
            Ω
            
              cos
              ⁢
              
                
              
              ⁢
              α
            
          
          ⁢
          
            
              V
              A
            
            v
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        4
      
    
  




Equation 4 can be re-defined to obtain the minimum frequency required to interact with protons outside the loss cone angle αL of energy E, as




  
    
      
        ω
        ≥
        
          
            Ω
            
              cos
              ⁢
              
                
              
              ⁢
              
                α
                L
              
            
          
          ⁢
          
            
              
                MV
                A
                2
              
              
                2
                ⁢
                
                  
                
                ⁢
                E
              
            
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        5
      
    
  




where M is the proton mass. FIG. 7 shows a graph 700 of the minimum frequency as a function of energy for typical equatorial plasma conditions at L=1.5. The pitch angle loss cone at L=1.5 is around 28 degree. One can see for example that SAW with frequency thirteen Hertz and bandwidth 2δω/ωo1 will resonate with thirty MeV at all angles inside the loss cone, as well as with higher energy trapped protons.
3.2. Pitch Angle Scattering Rate and Proton Lifetime
Computation of the scattering rate and proton lifetime in the presence of a given SAW spectrum requires a couple of relatively complex but otherwise standard computations. The first is to follow textbook procedure (Lyons and Williams, 1984) to calculate the pitch angle diffusion coefficient as a function of the SAW spectrum and amplitude for the energy range of interest. The second is to determine the effective diffusion coefficient by averaging over the bounce and azimuthal orbit. This gives the effective diffusion coefficient as a function of the SAW amplitude <δB>. It is important to emphasize that, this is the average amplitude that a proton sees when it completes its entire orbit. The results are the same, if the waves concentrated in a small azimuthal shell with higher amplitude or if the waves uniformly distribute over the azimuth. Finally the lifetime is computed as discussed in Lyons and Williams (1984), by solving the bounced averaged pitch angle diffusion equation as an eigen-value problem. The details of this analysis can be found in Shao et al. (2009).












TABLE 1







E in MeV
f1 = 6.5 Hz
f2 = 10 Hz
f3 = 13 Hz





















30
1688 days 
880 days
595 days



50
900 days
586 days
920 days



100
580 days
1032 days 
1600 days 










Table 1 shows the proton lifetimes in the presence of SAW with average amplitude 25 pT for selected injection frequencies and proton energies. Notice that the diffusion rate as well as the lifetime scale as the square of the SAW amplitude.
3.3. Energy-Power Requirements
Two factors affect the energy-power required to accomplish a desirable remediation. The first is obviously the size of the region in units of δL. The second is the SAW confinement time that in its turn depends on the reflection coefficient R of the SAW from the ionosphere. The results per δL=0.1 is expressed as a function of the reflection coefficient R. Considering the region of L=1.5, the volume is given approximately by 3×1020(δL/0.1) m3. Therefore, to achieve the lifetimes referred to in Table 1 the volume should contain a total energy of 75 kJ in SAW.
When SAW trapping region is treated as a leaky cavity, then




  
    
      
        
          
            d
            W
          
          
            d
            t
          
        
        =
        
          P
          -
          vW
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        6
      
    
  




Where W is the SAW energy, P the injected power and v is the energy loss rate. The energy loss rate due to transmission at the ionospheric boundary is v=−lnR/To, where R is the reflection coefficient at the ionospheric boundary and To is the Alfven wave transit time along the magnetic field line. FIG. 8 illustrates a graph 750 showing the power required to maintain SAW energy of 75 kJ as a function of the reflection coefficient. The power varies from a few hundred of watts to several kW.
4. Apparatus for Accomplishing PRBR
The basic PRBR system 600 concept illustrated in FIG. 6. SAW can be generated either from ground transmitters and injected upwards along magnetic field lines or injected directly from space based platforms. The SAW in the chosen frequency range interact resonantly with 10 to 100 MeV protons causing to precipitate at a rate that depends on the amplitude of the SAW. For amplitude of the order of 25 to 30 pT, the proton flux can be reduced by one order of magnitude, or more on times, one to three years.
4.1. Steady State Horizontal Electric Dipole (HED) Transmitters
HED system is similar to traditional ELF transmitters that are used for submarine communications such as the FELF system located in Michigan. Greifinger examined a similar system for lateral injection of ULF signals in the earth-ionosphere or the Alfvenic waveguide. FIG. 9a illustrates an apparatus 800 for injecting SAW to the ionosphere using Horizontal Electric Dipole (HED) and FIG. 9b illustrates an apparatus 850 for injecting SAW to the ionosphere using Rotating Magnetic Field (RMF) antennas.
In HED system, antenna can be utilized to inject SAW upwards through the lower ionosphere along the magnetic field lines as illustrated in FIG. 9a. The electric and magnetic fields of such an antenna will depend on the current I, the horizontal length L and the skin depth δ of the ground through which the return current flows. It is easy to calculate the magnetic field at altitude 75 km, the bottom of the magnetized ionosphere and given by




  
    
      
        
          
            H
            y
          
          ≈
          
            
              IL
              
                4
                ⁢
                π
                ⁢
                
                  
                
                ⁢
                
                  h
                  2
                
              
            
            ⁢
            
              (
              
                δ
                /
                h
              
              )
            
          
        
        ,
        
          
            E
            x
          
          ≈
          
            
              Z
              ⁡
              
                (
                h
                )
              
            
            ⁢
            
              IL
              
                4
                ⁢
                π
                ⁢
                
                  
                
                ⁢
                
                  h
                  2
                
              
            
            ⁢
            
              (
              
                δ
                /
                h
              
              )
            
          
        
        ,
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        7
      
    
  
  
    
      
        
          Z
          ⁡
          
            (
            h
            )
          
        
        =
        
          
            
              iωμ
              
                
                  
                    σ
                    P
                  
                  ⁡
                  
                    (
                    h
                    )
                  
                
                +
                iωɛ
              
            
          
          .
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        8
      
    
  




Z(h) is the impedance at the bottom of the ionosphere and σP(h) is the corresponding Pedersen conductivity. Note that the electric and magnetic fields are driven by two anti-parallel currents such as antenna current and image current, separated by the skin depth distance δ, assuming δ<<L. From Equations 7 and 8, the power density injected in the ionosphere by a HED with dipole moment IL at a frequency f is thus given by

S=Zo(1/√{square root over (1+σP(h)/iωε)})(IL/4πh2)2(δ/h)2  Equation 9

Where Z0=120π is the impedance of free space. Taking the approximate area at an altitude h as h2, the injected power in the RB in practical units will be given by

P(z=h)≈α4(IL/3×104 A-km)2(75 km/h)4(δ/7 km)2kW  Equation 10

In Equation 10 α≈cos2 θ√{square root over (εω/σP(h))} is the efficiency with which the power at the bottom of the ionosphere will couple to the SAW, if the angle that the earth's magnetic field makes to the ground at the transmitter location is θ. Based on the fact that the ionospheric attenuation at few Hz frequencies is negligible and using nighttime conditions, the factor α is of order unity. As a zero order estimate, a HED with L≈10 15 km and I≈1 to 3 kA located on ground with conductivity approximately 10−4 S/m could in principle inject a few kW of power into the SAW mode required to achieve lifetime of the order of 2 to 3 years for 30 to 100 MeV trapped protons. In such a system the main loss is ohmic heating of the ground and overall efficiencies of the order or better than 10−3 can be achieved. The total ground power required is of the order of few MW.
4.2. Rotating Magnetic Field (RMF)
An alternative system 850 that can inject SAW efficiently in the radiation belts is illustrated in FIG. 9b. The system 850 utilizes an array of superconducting coils (permanent magnets) with magnetic moment parallel to the ground and rotating at the relevant frequency ω=2πf. Note that in contrary to the HED, the image source of a Horizontal Magnetic Dipole is in the same direction, thereby doubling the strength of the source. The magnetic field at the bottom of the ionosphere of such a Rotating Magnetic Field (RMF) source will be given approximately by




  
    
      
        
          
            
              H
              →
            
            ⁡
            
              (
              
                z
                ≈
                h
              
              )
            
          
          ≈
          
            
              M
              ⁡
              
                (
                
                  
                    
                      
                        e
                        ^
                      
                      x
                    
                    ⁢
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                  +
                  
                    
                      
                        e
                        ^
                      
                      y
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                
                )
              
            
            
              π
              ⁢
              
                
              
              ⁢
              
                h
                3
              
            
          
        
        ,
        
          

        
        ⁢
        
          
            
              E
              →
            
            ⁡
            
              (
              
                z
                =
                h
              
              )
            
          
          ≈
          
            
              Z
              ⁡
              
                (
                h
                )
              
            
            ⁢
            
              
                
                  M
                  ⁡
                  
                    (
                    
                      
                        
                          
                            e
                            ^
                          
                          y
                        
                        ⁢
                        cos
                        ⁢
                        
                          
                        
                        ⁢
                        ω
                        ⁢
                        
                          
                        
                        ⁢
                        t
                      
                      +
                      
                        
                          
                            e
                            ^
                          
                          x
                        
                        ⁢
                        sin
                        ⁢
                        
                          
                        
                        ⁢
                        ω
                        ⁢
                        
                          
                        
                        ⁢
                        t
                      
                    
                    )
                  
                
                
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  
                    h
                    3
                  
                
              
              .
            
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        11
      
    
  




In practical units the power injected in SAW will be approximately

P≈α64(75 km/h)2(M/2×104 A-km2)2kW  Equation 12

An advantage of the RMF system is its compactness and portability. For example, a superconducting magnet with 25 m2 area, Four hundred Ampere DC current and 105 turns has an approximate magnetic moment of 109 A-m2. Approximately twenty coils will be needed to get inject kilowatt level power. A further advantage of such a system is that it does not require low conductivity ground and can thus be located in any desirable location as well as it can be portable. For example it can be located in a barge or any platform such as oil rig platforms.
4.3. Transient Horizontal Electric Dipole (THED) Transmitters
THED systems which are similar to HED systems operate in transient mode with pulse length of the order of 0.1 to 1 seconds and can inject broadband waves in the desired frequency band. A significant advantage of such a system is that it can increase the injection efficiency of the steady state HED by as much as 20 dB by avoiding the effect of the magnetic field generated by the ground return current at the bottom of the ionosphere. This is accomplished by an innovative “sneak-though” operation part of the present invention. FIG. 10 illustrates, a graph 900 of amplitude of the electric field in the ionosphere versus time, showing a sneak-through concept of Transient Horizontal Electric Dipole (THED) Transmitters. The sneak-though concept relies on the fact that, in order to reach the bottom of the ionosphere, the magnetic field generated by the return current has to travel through the conducting ground first, before moving at the speed of light to the bottom of the ionosphere, while the magnetic field of the antenna travels only through vacuum. Since the wave propagation through the ground is controlled by diffusive processes and by more than four orders of magnitude slower than the speed of light. As a result, operating on a transient mode allows the antenna signal to inject SAW before the ground signal reduces its amplitude.
FIG. 11a illustrates a graph 910 showing the ground response time as a function of antenna length and ground conductivity, similar to the UR time of a circuit. Timescales of 0.1 to 1 seconds can be accomplished by a combination of antenna lengths in the few hundred meters for a range of ground conductivities. FIG. 11b and FIG. 11c illustrate graphs 920 and 930, showing relative efficiency as compared to a HED without reduction by ground return effects for rise times of T/3 and T respectively. The relative efficiency of pulses with rise time T/3 and T as compared to the value expected by an antenna without effects due to the ground return are shown. In this case the power density is given by

S=Zo(1/√{square root over (1+σP(h)/iωε)})(IL/4πh2)2  Equation 13

This is similar to Equation 9 but with the value of δ/h=1. This factor reduces the efficiency of the steady state HED by more than 15 to 20 dB. As seen in from FIGS. 11b and 11c, the reduction due to the ground return in the THED case is less than one dB.
FIG. 12 illustrates a flow chart 950 showing a method of reducing trapped energetic proton flux at LEO orbits, in accordance with the disclosed embodiments. As illustrated at block 952, the electromagnetic waves are generated using at least one transmitter. Then as said at block 954, the generated electromagnetic waves are injected into earth ionosphere. The frequency range of injected electromagnetic waves are then selected either from ground or space in order to ensure electromagnetic waves are in gyro frequency resonance with trapped protons of energy between 10 and 100 Mev as depicted at blocks 956 and 958. Pitch angle scattering of the trapped protons in gyro-resonance with the injected waves increases their precipitation rate by forcing their orbits into pitch angles inside the atmospheric loss-cone and are lost by intaracting with the dense neutral atmosphere at altitudes below 100 km. Thus, the energetic proton flux trapped in the inner radiation belt are reduced by injecting Ultra-Low Frequency (ULF) electromagnetic waves.
It will be appreciated that variations of the above disclosed and other features and functions, or alternatives thereof, may be desirably combined into many other different systems or applications. Also that various presently unforeseen or unanticipated alternatives, modifications, variations or improvements therein may be subsequently made by those skilled in the art which are also intended to be encompassed by the following claims.

## Claims

1. A method for reducing the energetic proton flux trapped in the inner radiation belt of earth comprising:
generating said electromagnetic waves using at least one transmitter;
injecting said electromagnetic waves to earth ionosphere; and
selecting frequency range of electromagnetic waves, wherein the frequency range is selected such that said electromagnetic waves are gyro frequency resonance with trapped protons of energy between 10 and 100 Mev, and wherein the gyro frequency resonance of the trapped protons with the injected electromagnetic waves increases the proton precipitation rate by forcing the orbits of the trapped protons into pitch angles inside the atmospheric loss-cone and the frequency range of said electromagnetic waves is selected such that the minimum frequency is either equal to or greater than





  
    Ω
    
      cos
      ⁢
      
        
      
      ⁢
      α
      ⁢
      
        
      
      ⁢
      L
    
  
  ⁢
  
    
      
        
          M
          ·
          
            VA
            2
          
        
        
          2
          ⁢
          
            
          
          ⁢
          E
        
      
    
    .

2. The method of claim 1 wherein said at least one transmitter comprising a ground transmitter for injecting said electromagnetic waves upwards along magnetic field lines.

3. The method of claim 1 wherein said at least one transmitter comprising a space transmitter for injecting said electromagnetic waves directly from space based platforms.

4. The method of claim 1, wherein the step of generating said electromagnetic waves comprises of rotating one or more ground based arrays of permanent or superconducting magnets at the selected frequency range thereby reducing proton flux.

5. The method of claim 1 wherein said electromagnetic waves are generated and injected into earth's ionosphere by utilizing a transmitter comprising horizontal electric dipole transmitter, rotating magnetic fields using arrays of permanent or superconducting magnets and transient horizontal electric dipole transmitter.

6. A method for reducing the energetic proton flux trapped n the inner radiation belt of earth comprising:
generating Ultra Low Frequency electromagnetic waves using at least one transmitter;
injecting said Ultra Low Frequency electromagnetic waves to earth ionosphere; and
selecting frequency range of Ultra Low Frequency electromagnetic waves, wherein the frequency range is selected such that said Ultra Low Frequency electromagnetic waves are in gyro frequency resonance with trapped protons of energy between 10 and 100 Mev, and wherein the gyro frequency resonance of the trapped protons with the injected electromagnetic waves increases the proton precipitation rate by forcing the orbits of the trapped protons into pitch angles inside the atmospheric loss-cone and the frequency range of said Ultra Low Frequency electromagnetic waves is selected such that the minimum frequency is either equal to or greater than





  
    Ω
    
      cos
      ⁢
      
        
      
      ⁢
      α
      ⁢
      
        
      
      ⁢
      L
    
  
  ⁢
  
    
      
        
          M
          ·
          
            VA
            2
          
        
        
          2
          ⁢
          
            
          
          ⁢
          E
        
      
    
    .

7. The method of claim 6 wherein said at least one transmitter comprising a ground transmitter for injecting said electromagnetic waves upwards along magnetic field lines.

8. The method of claim 6 wherein said at least one transmitter comprising a space transmitter for injecting said electromagnetic waves directly from space based platforms.

9. The method of claim 6 wherein the step of generating said Ultra Low Frequency electromagnetic waves comprises of rotating one or more ground based arrays of permanent or superconducting magnets at the selected frequency range thereby reducing proton flux.

10. The method of claim 6 wherein said electromagnetic waves are generated and injected into the earth's ionosphere by utilizing a transmitter comprising horizontal electric dipole transmitter, rotating magnetic fields using arrays of permanent or superconducting magnets and transient horizontal electric dipole transmitter.

11. A method for reducing the energetic proton flux trapped in the inner radiation belt of earth comprising:
generating said Ultra Low Frequency electromagnetic waves using at least one transmitter;
injecting said Ultra Low Frequency electromagnetic waves to earth ionosphere where they are lost by interacting with a dense neutral atmosphere at altitudes below 100 km;
selecting frequency range of Ultra Low Frequency electromagnetic waves, wherein the frequency range is selected such that said Ultra Low Frequency electromagnetic waves are in gyro frequency resonance with trapped protons of energy between 10 and 100 Mev and wherein the gyro frequency resonance of the trapped protons with the injected electromagnetic waves increases the proton precipitation rate by for the orbits of the trapped protons into pitch angles inside the atmospheric loss-cone where the protons are lost by interacting with a dense neutral atmosphere at altitudes below 100 km and the frequency range of said Ultra Low Frequency electromagnetic waves is selected such that the minimum frequency is either equal to or greater than





  
    Ω
    
      cos
      ⁢
      
        
      
      ⁢
      α
      ⁢
      
        
      
      ⁢
      L
    
  
  ⁢
  
    
      
        
          M
          ·
          
            VA
            2
          
        
        
          2
          ⁢
          
            
          
          ⁢
          E
        
      
    
    .

12. The method of claim 11 wherein said at least one transmitter comprising a ground transmitter for injecting said electromagnetic waves upwards along magnetic field lines.

13. The method of claim 11 wherein said at least one transmitter comprising a space transmitter for injecting said electromagnetic waves directly from space based platforms.

14. The method of claim 11, wherein the step of generating said Ultra Low Frequency electromagnetic waves comprises of rotating one or more ground based arrays of permanent or superconducting magnets at the selected frequency range thereby reducing proton flux.

15. The method of claim 1, wherein the frequency range of said electromagnetic waves is more than 10 Hz.

16. The method of claim 6, wherein the frequency of said Ultra Low Frequency electromagnetic waves is more than 10 Hz.

17. The method of claim 11, wherein the frequency range of said Ultra Low Frequency electromagnetic waves is more than 10 Hz.