Spaceflight

Three clocks and nothing to fall onto

A charged particle in a dipole field gyrates, bounces and drifts, on timescales a millisecond, a second and a quarter of an hour. The three periods are three decades apart, and that separation is not a curiosity — it is the reason each motion has a conserved quantity, and the reason a magnetic storm can accelerate particles rather than merely stir them.

Assumes Magnetosphere and Orbital debris.

A particle in a gravitational field falls. A charged particle in a magnetic field cannot: the force is always perpendicular to its motion, so it does no work, and the particle has nowhere to go and nothing to lose. It simply keeps moving, and the geometry of the field decides where.

Around a planet with a dipole field the answer is that it stays. Not in an orbit — the field does not attract it — but in a region, executing three superposed motions for as long as nothing disturbs them. The belts of trapped radiation that result were the first discovery of the space age, and they are the environment every satellite is designed against.

The discovery is worth a sentence because of how it was made. Explorer 1 carried a Geiger counter to measure cosmic rays, and above about a thousand kilometres it reported zero counts — not a high rate, none at all. The instrument had not failed; it had saturated, and a saturated counter reads empty. What looked like the absence of radiation was a flux thousands of times higher than anything expected, and the correct reading of a null result took a second flight with a shielded detector to establish.

Three motions, three decades apart, and that is the point. The three periods of a trapped 1-MeV electron's motion against the shell it is trapped on, on a logarithmic time axis. It gyrates about a field line in 0.22 milliseconds, bounces between its mirror points in 0.33 seconds, and drifts right round the planet in 16 minutes — ratios of 1498 and 3015 at L = 4. Each motion carries a conserved quantity: the magnetic moment, the longitudinal invariant, and the magnetic flux the drift shell encloses. The separation is what makes them conserved. An invariant survives anything that changes slowly compared with its own period, so a disturbance lasting minutes destroys the third and leaves the first two untouched — and a particle that keeps its magnetic moment while being moved inward to a stronger field must gain energy. That is not a loophole; it is how the belts are filled.
Fig. 1 The three periods of a trapped particle’s motion against the shell it is trapped on. It gyrates in a fifth of a millisecond, bounces between its mirror points in a third of a second, and drifts right round the planet in sixteen minutes. The ratios are about fifteen hundred and three thousand, and the whole of the argument below depends on their being large.

The first motion, and what it conserves

The fastest motion is the gyration: a helix about a field line, at a radius set by the particle’s momentum and the field strength, and at a frequency set by the field strength alone.

That frequency is the thing to hold on to. It does not depend on the particle’s energy at all in the non-relativistic case, and only through the Lorentz factor otherwise. A field of a milligauss gyrates every electron in it at the same rate, whatever they are doing, which is why a plasma has a characteristic frequency rather than a distribution of them. It is also why a splitting proportional to the field is the standard magnetic diagnostic in an atom: the same gyration, at the same frequency, on a bound electron rather than a free one.

Around that gyration there is a conserved quantity: the magnetic moment, the perpendicular kinetic energy divided by the field strength. It is conserved not exactly but adiabatically — meaning that it is conserved to all orders in a small parameter, and the small parameter is the ratio of the gyration period to the timescale over which the field the particle sees is changing.

The proof of that is a classical result about slowly varying oscillators, and its content is unusual: the quantity conserved is not an energy or a momentum but an action, the area enclosed in phase space by one cycle of the motion. Actions are the quantities that survive slow change, which is the same reason they are the quantities that get quantised.

“Adiabatic” is doing real work in that sentence and is not a hedge. An exactly conserved quantity is protected by a symmetry and cannot change at all; an adiabatic invariant has no symmetry behind it and changes by an amount that is exponentially small in the ratio of timescales rather than merely small. That is a much stronger statement than “approximately constant”, and it is why the invariants can be treated as coordinates — a trapped particle is labelled by three numbers that do not change, and its position and velocity are derived quantities.

For a particle moving along a field line into a stronger field, conservation of the moment means the perpendicular energy must rise. Total energy is fixed, so the parallel energy falls, and at some point it reaches zero. The particle stops moving along the line and turns round.

A mirror, made of nothing

That turning point is a magnetic mirror, and it is worth insisting on how peculiar it is. Nothing is there. The field has no gradient sharp enough to be called a surface, no material, no potential barrier in the ordinary sense. What reflects the particle is the geometry of its own motion.

Where a particle turns round, and why it is not the same for all of them. A dipole field line at L = 4, with the path of a particle whose equatorial pitch angle is 45°. The first adiabatic invariant fixes the ratio of the perpendicular energy to the field strength, so as the particle moves toward the pole and the field strengthens, its motion tips further into the perpendicular — and where the field has grown by 1/sin²α, which here is a factor of 2.00, the parallel motion has gone entirely and it turns round. That happens at 23.1° of latitude. A particle with a smaller pitch angle mirrors deeper, and below 5.5° it would mirror inside the atmosphere and be absorbed instead — the loss cone, which at this shell is a narrow 5.5-degree hole in an otherwise full distribution. The gyration drawn along the path is the real one: the radius falls as the field rises, which is why the helix tightens toward each end.
Fig. 2 The reflection, drawn. A dipole field line at four planetary radii, with a particle whose motion at the equator makes forty-five degrees to the field. Conservation of the magnetic moment requires it to turn where the field has grown by a factor of two, which happens at twenty-three degrees of latitude — and the drawn helix tightens toward each end because the gyroradius falls as the field rises.

Where the mirror sits depends only on the angle between the particle’s velocity and the field at the equator, and on nothing else — not its energy, not its mass, not its charge. A particle moving almost perpendicular to the field mirrors almost immediately; one moving almost along it penetrates deep toward the pole.

That gives a criterion for whether a particle is trapped at all. If its equatorial pitch angle is small enough that its mirror point would lie inside the atmosphere, it does not mirror; it collides, and it is gone. The set of angles for which that happens is the loss cone, and at four planetary radii it is between five and six degrees wide.

Five degrees out of ninety is a small hole, and the smallness is the reason the belts exist. A particle scattered by even a few degrees can fall into it, but scattering has to happen, and in a nearly collisionless plasma it happens slowly.

The loss cone’s width is set by the ratio of the equatorial field to the field at the top of the atmosphere, and it therefore shrinks rapidly with distance: at two planetary radii it is about seventeen degrees and at six about three. So the outer belt has a smaller hole to lose particles through than the inner one, which is part of why it responds faster to being filled and holds less well once quiet returns.

There is a directly observable consequence. Particles inside the loss cone are the ones that reach the atmosphere, and what they do when they get there is produce aurora and X-rays. So a measurement of precipitation at the top of the atmosphere is a measurement of how full the loss cone is, and therefore of how hard the scattering is working — an observable for a process happening thousands of kilometres above.

The second and third motions

Bouncing between mirror points is the second motion, with its own conserved quantity — the integral of the parallel momentum along the bounce path, which is again an action.

The third motion is a drift in longitude. It arises because the field is not uniform across a gyration: the particle’s orbit is slightly tighter on the side nearer the planet, where the field is stronger, so the circle does not close and the guiding centre walks sideways. Positive and negative charges drift in opposite directions, which makes a current — the ring current, whose field at the planet’s surface is what a magnetic storm’s magnetograms record.

The conserved quantity here is the magnetic flux enclosed by the drift path, and the surface traced out by a particle’s bounce path as it drifts is called a drift shell.

Drift shells are the reason the belts can be described by one coordinate at all. A particle’s position varies over sixty thousand kilometres in the course of a drift, and its energy varies as it moves through fields of different strength — but the shell it is on does not change, so a flux quoted “at L equals four” is a statement about a population rather than about a place. The coordinate was invented by Carl McIlwain in 1961 for exactly this purpose, and it collapses a messy three-dimensional field into a single number that particle data organise themselves by.

The organisation only works while the third invariant holds. A storm scrambles it, and the same particle data then have to be re-sorted, which is the practical face of the acceleration mechanism described below.

Two belts and the gap a mission is designed around. Trapped particle flux against shell parameter, on a logarithmic scale. Two populations, and they are different particles: the inner belt is protons of tens of megaelectronvolts, produced when cosmic rays strike the upper atmosphere and the resulting neutrons decay in flight; the outer belt is electrons, injected and accelerated by magnetic storms. Between them at L = 2.33 is the slot, 49 times emptier than either peak, and it is emptier because a plasma wave at just the right frequency scatters particles into the loss cone there. Every one of those numbers is a design input. A transfer orbit crosses both belts twice per revolution; a geostationary satellite lives at L = 6.6, on the outer edge of the outer belt, where the flux varies by three orders of magnitude with the state of the solar wind. The dose is not an environment to be endured but a geometry to be flown through.
Fig. 3 The result, as a flux against shell. Two populations, of different particles: an inner belt of energetic protons made when cosmic rays strike the upper atmosphere and the resulting neutrons decay in flight, and an outer belt of electrons injected and accelerated by storms. Between them is a slot two orders of magnitude emptier, and the slot is a design input rather than a curiosity.

Three motions, three actions, three timescales. The separation between the timescales is what makes each action conserved, because each is conserved against disturbances slow compared with its own period — and a disturbance slow enough to break the third is far too slow to touch the first.

Breaking the slowest one is an accelerator

This is where the structure earns its keep, and it is not obvious.

A magnetic storm disturbs the field on a timescale of minutes to hours. That is slow compared with a gyration and a bounce, and fast compared with nothing — it is comparable to the drift period. So it breaks the third invariant and leaves the first two intact.

A particle whose enclosed flux is no longer conserved wanders in shell: it can be moved inward, toward stronger field. But its magnetic moment is still conserved, and the magnetic moment is perpendicular energy over field strength. Move a particle to a field ten times stronger at fixed moment and its perpendicular energy has gone up by ten.

That is radial diffusion, and it is an accelerator built out of a conservation law. A hundred-kiloelectronvolt electron in the outer magnetosphere, transported inward by successive storms, arrives at a few planetary radii with megaelectronvolts. Nothing did work on it in the ordinary sense; the field’s slow change did, and the bookkeeping is the invariant.

Four torques, and the altitude where the answer changes. The four environmental torques acting on one spacecraft — 10 square metres of area, principal moments differing by 200 kilogram metres squared, a centre of pressure 10 centimetres from the centre of mass and a residual magnetic dipole of 1 ampere metre squared — plotted against altitude, both axes logarithmic. The aerodynamic torque falls by six orders of magnitude across the range because the density does; the gravity gradient falls only as the cube of the distance from the Earth's centre, which over this range is barely a factor of three; the magnetic torque falls as the same cube; and sunlight does not change at all, which the figure checks rather than states. The leader changes at 438 kilometres, from the atmosphere to the gravity gradient, and that crossing is measured off the drawn curves. Nothing here is a disturbance in the ordinary sense: each of these acts in the same direction for a large part of every orbit, so what matters is not the size of the torque but the part of its integral that survives an orbit — and the wheel that stores that integral is emptied on the schedule these curves set.
Fig. 4 The same accounting logic in a different setting. A quantity that is conserved on one timescale and not another can be moved between reservoirs by choosing the timescale of the disturbance — which is how a spacecraft’s momentum wheels absorb a secular torque and then dump it, and how a storm turns a slow rearrangement of a field into a fast particle.

The mechanism was proposed before it was seen and has since been watched happening: spacecraft measuring particle fluxes as functions of the three invariants have followed populations moving inward at fixed moment and gaining exactly the energy the conservation law demands. It is one of the cleaner confirmations in the subject, because the prediction is not a rate but an exact relation between two measured quantities.

The energy available this way has a ceiling, and the ceiling is where the particle runs out of shells to be moved through. Inward transport stops when the particle reaches the plasmasphere, where the waves that scatter it into the loss cone are strong; outward transport takes it to the magnetopause, where it is simply lost to the wind. So the belts are bounded on both sides by loss processes rather than by any limit on the acceleration, in much the way that a trapped population of stars in a cluster is bounded by evaporation rather than by a limit on how hard it can be stirred, which is why their outer edge tracks the magnetopause’s position during a storm and why a compressed magnetosphere empties the outer belt rather than filling it.

What empties the slot

The gap between the belts is as informative as the belts, and it took some time to explain.

The slot is not a place where particles are absent because none arrive. Radial diffusion delivers them there as readily as anywhere. It is a place where they are removed, and the removal is done by waves.

A plasma supports waves of many kinds, and some of them have frequencies matching the gyration frequency of the trapped particles. A particle in resonance with a wave has its pitch angle scattered — the wave changes the direction of the motion without much changing its energy — and a scattered particle can end up in the loss cone.

The waves responsible in the slot are called plasmaspheric hiss, and they exist there because the slot is inside the plasmasphere, the cold dense plasma that corotates with the planet. So the slot’s location is set by the plasmasphere’s outer boundary, which moves in and out with geomagnetic activity — and after a large storm the slot fills, then empties again over weeks.

That is a satisfying closure. The belts are populated by the breaking of one invariant and depopulated by the breaking of another, and the second breaking is done by waves whose existence depends on a third population of particles entirely.

It is also a warning about how the belts should be described. A picture of two static doughnuts is wrong in a way that matters: the outer belt can vary in flux by four orders of magnitude in a day, empty in hours and refill in a day, and occasionally split into two. In 2012 a third belt was observed to form and persist for four weeks, a structure that had not been predicted and that required the spacecraft measuring it to have the energy resolution to see it as separate rather than as a broadening.

The lesson for anyone designing to a radiation environment is that the belts are a statistical model rather than a map. The engineering standards quote fluxes averaged over years, and the worst case a given mission meets can exceed the average by a large factor for a short time.

What a spacecraft actually integrates

The engineering question is not what the flux is but what dose a given mission accumulates, and the two are related by geometry and by shielding.

The first millimetre buys most of what there is to buy. Dose behind aluminium shielding, relative to the unshielded value, against thickness. The belt electron spectrum falls roughly exponentially with energy, and an electron's range is roughly proportional to its energy, so the dose falls exponentially with thickness — until it does not. The flat part is bremsstrahlung: an electron stopped in the shield radiates, and the resulting photons are far more penetrating than the electron was, so past about 10 millimetres each extra kilogram removes almost nothing — the fourteenth millimetre buys a factor of 1.005 where the first bought 10.3. The asymptote sits at 1.2 per cent of the unshielded dose. This is why spacecraft shielding is millimetres rather than centimetres, why the electronics rather than the structure carries the local shielding, and why the honest way to survive the belts is to spend as little time in them as the orbit allows.
Fig. 5 Dose behind aluminium against thickness, relative to unshielded. The first millimetre buys a factor of ten; the fourteenth buys almost nothing, because an electron stopped in the shield radiates a photon that the rest of the shield does not stop. The floor is bremsstrahlung, and it is why spacecraft shielding is measured in millimetres rather than centimetres.

The floor is the important feature. Stopping an electron converts part of its energy into penetrating photons, so the returns from adding mass diminish sharply and then stop. Beyond a few millimetres the sensible engineering move is not more shielding but less exposure — a different orbit, a faster transfer, or a component that tolerates the dose.

That is why the belts shape mission design so strongly. A transfer to geostationary orbit crosses both belts twice per revolution and takes weeks if flown with electric propulsion, which is why a low-thrust spiral trades propellant against radiation. Medium Earth orbit, where the navigation constellations live, sits in the outer belt permanently. Low orbit is largely below the belts, except where the field is weak — and low orbit has a different hazard of its own whose statistics work in a completely different way.

Where a particle turns round, and why it is not the same for all of them. A dipole field line at L = 1.5, with the path of a particle whose equatorial pitch angle is 30°. The first adiabatic invariant fixes the ratio of the perpendicular energy to the field strength, so as the particle moves toward the pole and the field strengthens, its motion tips further into the perpendicular — and where the field has grown by 1/sin²α, which here is a factor of 4.00, the parallel motion has gone entirely and it turns round. That happens at 33.2° of latitude. A particle with a smaller pitch angle mirrors deeper, and below 28.1° it would mirror inside the atmosphere and be absorbed instead — the loss cone, which at this shell is a narrow 28.1-degree hole in an otherwise full distribution. The gyration drawn along the path is the real one: the radius falls as the field rises, which is why the helix tightens toward each end.
Fig. 6 The inner belt’s geometry, drawn on a much smaller shell. The field line is shorter, the mirror points are further from the equator for the same pitch angle, and the loss cone is several times wider — so the inner belt loses particles more readily per bounce and holds them anyway, because nothing scatters them into it.

That exception has a name. The Earth’s dipole is offset from the planet’s centre, so the field at a given altitude is weakest over the South Atlantic, and the belts reach lower there. Spacecraft in low orbit take most of their dose in a few minutes per day crossing the South Atlantic Anomaly, and many instruments are simply switched off for the crossing.

What the belts are made of, and where it came from

The two belts are different particles from different sources, and the difference is worth setting out because it explains their very different behaviour.

The inner belt is mostly protons of tens of megaelectronvolts. They are made by a process with no magnetospheric content at all: a cosmic ray strikes a nucleus in the upper atmosphere, the collision produces neutrons, some neutrons travel upward, and a neutron decays in about fifteen minutes into a proton, an electron and an antineutrino. The proton appears in the middle of the trapping region already energetic and already charged, and is trapped from that instant.

Ten decades of one power law, with two places where it bends. The cosmic-ray spectrum, multiplied by energy to the 2.7 so that its features can be seen at all — undivided, it falls by thirty orders of magnitude across this plot and every bend in it is invisible. The knee at 3·10¹⁵ electronvolts is where the spectrum steepens from E^−2.7 to E^−3.1, and it sits within a factor of a few of the energy at which a proton's gyroradius in a microgauss field becomes comparable to the thickness of the galactic disc: above it, confinement begins to leak. The ankle at 3·10¹⁸ is where it flattens again, which is read as a galactic population running out and an extragalactic one taking over. Above 5·10¹⁹ the spectrum is cut off, because a proton that energetic loses energy to the microwave background within about fifty megaparsecs and cannot have come from further. Three features, each of them a statement about a magnetic field: one about the galaxy's, one about where the galaxy's ends, and one about the fact that empty space is not empty.
Fig. 7 The parent population of the inner belt, measured on the ground. Cosmic rays striking the atmosphere are what make the neutrons whose decay stocks the inner belt, so a proton trapped at one and a half planetary radii began as a particle accelerated somewhere in the galaxy — and its confinement there was also magnetic.

That source is steady, because the cosmic-ray flux is steady, and the losses are slow, because at low altitude the field is strong and the loss cone narrow. So the inner belt is the stable one: its flux varies by a factor of two over a solar cycle and essentially not at all from day to day.

The outer belt is electrons, and its source is the magnetosphere’s own plasma sheet — particles of a few kiloelectronvolts, injected during substorms and then accelerated by the two mechanisms above. That source is episodic and the losses are fast, which is why the outer belt is the volatile one.

The asymmetry is fortunate for spacecraft design in one respect and unfortunate in another. The inner belt can be modelled once and designed against; the outer belt cannot, and a mission that spends time in it has to be designed against a distribution rather than a number.

What the invariants do not survive

The adiabatic picture is very good and it is not universal, and the exceptions are where the current research is.

The first invariant fails when the field changes on the gyration timescale, which happens in the sharp field reversal at the centre of the magnetotail: the field there is weak, the gyroradius is large, and particles execute chaotic figure-eight orbits rather than helices. The tail’s plasma sheet is therefore a place where particles are energised non-adiabatically, and it is the source population for much of what the belts later accelerate.

It also fails for waves at the gyrofrequency itself, which is the case of chorus — an intense, structured emission generated near the equator outside the plasmasphere. Chorus can accelerate electrons locally, in place, without any radial transport, and the two mechanisms compete. Chorus is worth one more sentence because of what it sounds like. Converted to audio, the emissions are a rising tone repeated every fraction of a second, which is exactly what early radio listeners on the ground described as a dawn chorus of birds — heard, and named, decades before anything was known about where it came from. It is generated by the same electrons it later accelerates, which makes the belts one of the rare astrophysical systems that both drives and is driven by its own waves.

Distinguishing local acceleration from radial transport requires measuring the phase-space density as a function of the invariants: if the density peaks at an intermediate shell, something is making particles there; if it rises monotonically outward, they are being transported inward. That measurement has been made, the peak is there, and the current view is that both processes operate with chorus dominating during the main phase of a storm.

The belts and the dose they deliver are the two things a mission actually has to plan around, and both are worth reading at a second setting.

Two belts and the gap a mission is designed around. Trapped particle flux against shell parameter, on a logarithmic scale. Two populations, and they are different particles: the inner belt is protons of tens of megaelectronvolts, produced when cosmic rays strike the upper atmosphere and the resulting neutrons decay in flight; the outer belt is electrons, injected and accelerated by magnetic storms. Between them at L = 2.33 is the slot, 49 times emptier than either peak, and it is emptier because a plasma wave at just the right frequency scatters particles into the loss cone there. Every one of those numbers is a design input. A transfer orbit crosses both belts twice per revolution; a geostationary satellite lives at L = 6.6, on the outer edge of the outer belt, where the flux varies by three orders of magnitude with the state of the solar wind. The dose is not an environment to be endured but a geometry to be flown through.
Fig. 8 The two belts with the marker at a shell of six Earth radii. The outer belt’s electrons dominate there and its boundary moves with magnetic activity, so a spacecraft in a high orbit crosses a region whose intensity varies by orders of magnitude on a timescale of days.
The first millimetre buys most of what there is to buy. Dose behind aluminium shielding, relative to the unshielded value, against thickness. The belt electron spectrum falls roughly exponentially with energy, and an electron's range is roughly proportional to its energy, so the dose falls exponentially with thickness — until it does not. The flat part is bremsstrahlung: an electron stopped in the shield radiates, and the resulting photons are far more penetrating than the electron was, so past about 20 millimetres each extra kilogram removes almost nothing — the fourteenth millimetre buys a factor of 1.005 where the first bought 10.3. The asymptote sits at 1.2 per cent of the unshielded dose. This is why spacecraft shielding is millimetres rather than centimetres, why the electronics rather than the structure carries the local shielding, and why the honest way to survive the belts is to spend as little time in them as the orbit allows.
Fig. 9 And the dose behind shields from half a millimetre to two centimetres. The first millimetre stops most of the electrons and almost none of the protons, so beyond a few millimetres the returns are very small — which is why radiation-hard electronics exist rather than thicker boxes.

Why the separation of timescales is the whole thing

Strip everything else away and one structural fact is left, and it is worth stating on its own.

A particle in a dipole has three periodic motions whose periods stand in ratios of about a thousand each. Because they are separated, each has an adiabatic invariant. Because the invariants are separated, a disturbance breaks them one at a time, in order, from the slowest. And because breaking the slowest while keeping the fastest forces an energy change, the system has a built-in accelerator that runs on nothing but slow rearrangement.

None of that would work if the three periods were comparable. A system with one timescale has one conserved quantity or none, and disturbing it stirs rather than sorts.

Three motions, three decades apart, and that is the point. The three periods of a trapped 2-MeV electron's motion against the shell it is trapped on, on a logarithmic time axis. It gyrates about a field line in 1.63 milliseconds, bounces between its mirror points in 0.47 seconds, and drifts right round the planet in 6 minutes — ratios of 287 and 710 at L = 6.6. Each motion carries a conserved quantity: the magnetic moment, the longitudinal invariant, and the magnetic flux the drift shell encloses. The separation is what makes them conserved. An invariant survives anything that changes slowly compared with its own period, so a disturbance lasting minutes destroys the third and leaves the first two untouched — and a particle that keeps its magnetic moment while being moved inward to a stronger field must gain energy. That is not a loophole; it is how the belts are filled.
Fig. 10 The same three periods at geostationary distance for a more energetic particle. The separation is preserved — it comes from the geometry rather than from the numbers — which is why the argument transfers unchanged to Jupiter, where the same three motions run a thousand times slower and produce a radiation environment that has destroyed spacecraft.

The transfer to other planets is the test. Jupiter’s field is twenty thousand times the Earth’s in moment, its rotation is faster, and its magnetosphere is loaded with plasma from Io — and the same three-motion structure describes it, with the same invariants, giving belts whose electron fluxes are high enough that a spacecraft passing through them accumulates in hours what would take years at Earth. The mission design response has been the same one an Earth-orbiting engineer would recognise: fly a polar orbit that ducks under the worst of it, and keep the perijove passes short.

There is one further echo worth noting. The same three-timescale structure — a fast motion with an action, a slower one, and a slowest that a disturbance can break — is the reason a spacecraft’s attitude control fills a momentum store on a schedule and has to empty it against something outside the vehicle. In both cases the useful move is to disturb the system at the timescale of the invariant one wants to break, and to leave the others alone.

What this makes readable

Essays that name this one as a prerequisite.

What links here

Essays that link to this one from their own argument.

The objects this essay names

Each one links to every other essay that touches it.

Adiabatic invariantBremsstrahlungDrift shellGeomagnetic stormLoss coneMagnetic mirrorMagnetic momentMagnetosphereRadial diffusionRadiation beltsSpacecraft shielding