A large proton event associated with solar filament activity

We report observations made from several interplanetary spacecraft, of the large low-energy particle event of 23–27 April, 1979 associated with solar filament activity. We discuss the intensity, spectral and directional evolution of the event as observed in the energy range 35–1600 keV on ISEE-3, located ～ 0.99 AU from the Sun upstream of the Earth. We demonstrate that the shock disturbance propagating through the interplanetary medium and observed at ISEE-3 on 24/25 April strongly controls the particle event. From a comparison of the ISEE-3 observations with those on other spacecraft, in particular on Helios-2, located at 0.41 AU heliocentric distance near the Sun-Earth line, we identify the solar filament erupting on late 22 April near central meridian as the trigger for the propagating shock disturbance. This disturbance which comprises a forward shock and a reverse shock at the orbit of ISEE-3 is found to be the main source of the energetic proton population observed.     

Predicted interplanetary shocks/particles at Mars compared with in-situ observations: An overview

The biological and technological consequences of long-duration, solar-related, energetic particle radiation for manned/unmanned spacecraft warrant that consideration be given to providing reliable space weather predictions for future space missions to planet Mars. An account is, herein, provided of how the HAFv.2 numerical model was applied to predict the arrivals of four, flare-related, shocks at Mars generated during a >20-day active period on the Sun in March 1989, and of the arrival of another composite shock produced in association with a 10-day period of solar activity in December 2006. These predictions are compared with in-situ measurements of shock signatures at Mars recorded, in the former case, by the solar-low-energy-detector (SLED) and by the low-energy-telescope (LET) aboard the Phobos-2 spacecraft and, in the latter case, in data recorded by the ASPERA-3/IMA instrument aboard Mars Express. The success of the predictions is discussed and the requirement for further validation of the modeling technique using a large statistical sample pointed out. In-situ measurements made aboard Mars Express by the ASPERA-3/IMA experiment during the rising phase of Solar Cycle 24 can provide data relevant to such validation. The successful application of a SOLar Particle ENgineering COde (SOLPENCO), that estimates solar energetic particle (SEP) fluxes and fluences at the Earth, to the case of an energetic particle event at Mars (6 March 1989) is discussed. Measurements of SEP events recorded by the Solar TErrestrial RElations Observatory (STEREO) supplemented by Mars Express measurements can potentially allow the predictions of SOLPENCO to be further studied downstream using a large statistical sample. However, we are presently only at the beginning of our understanding of the complex Sun-Earth-Mars scenarios that give rise to shock/particle events in the close Martian environment.     

Energetic protons associated with a forward-reverse interplanetary shock pair at 1 A.U.

A forward-reverse interplanetary shock was observed on 25 March 1969 by the magnetometer and plasma detector on the HEOS-1 satellite. This relatively rare event was described by Chao et al (1972) who concluded that the shock pair was formed at a distance 0.10–0.13 A.U. upstream of the Earth as a result of the interaction between a fast and a slow solar wind streams. Simultaneous observations of 1 MeV solar proton fluxes were also performed on HEOS-1. A characteristic intensity peak was observed as the forward shock passed by the spacecraft. The evolution of the proton intensity, together with a detailed analysis of anisotropies and pitch angle distributions show a complex dynamic picture of the effect of the forward shock on the ambient proton population. Significant changes in particle fluxes are seen to be correlated with fluctuations in the magnetic field. It is suggested that simple geometrical models of shock-associated acceleration should be expanded to include the effect of magnetic fluctuations on particle fluxes. The interaction region limited by the forward and reverse shocks contained a large variety of magnetic fluctuations. Following the tangential discontinuity separating the fast solar wind stream from the preceding slow stream, a sunward flow was observed in the proton data, followed by a small but significant drop in intensity prior to the reverse shock.     

Multi-spacecraft Study of the 8 November 2000 SEP Event: Electron Injection Histories 100° Apart

We present the analysis of a large solar near-relativistic electron event observed by the Ulysses and the ACE spacecraft on 8 November 2000, when Ulysses was located at a heliocentric distance of 2.4 AU and at a heliographic latitude of ??80° S. We use a particle propagation model to infer the local interplanetary transport conditions and the injection histories of the near-relativistic electrons observed by both spacecraft. We find different local transport conditions for each set of observations. The inferred injection profiles for both spacecraft extend for several hours; but the injection at Ulysses was smaller and started later. The association with type II radio emission suggests that the heliospheric electrons were provided by coronal shock acceleration. An analysis of the in situ magnetic field and plasma measurements indicates that the global configuration of the heliosphere (disturbed by transient structures) could play a role in shaping the characteristics of solar energetic particle events observed from different locations.     

Dominant acceleration processes of ambient energetic protons (E ⩾ 50 keV) at the Bow Shock: Conditions and limitations

Energetic proton (E_p ? 50 keV) and magnetic field observations during crossings of the Earth's Bow Shock by the IMP-7 and 8 spacecraft are incorporated in this work in order to examine the effect of the Bow Shock on a pre-existing proton population under different “interplanetary magnetic field-Bow Shock” configurations, as well as the conditions for the presence of the Bow Shock associated energetic proton intensity enhancements. The presented observations indicate that the dominant process for the efficient acceleration of ambient energetic particles to energies exceeding ~ 50 keV is by “gradient-B” drifting parallel to the induced electric field at quasi-perpendicular Bow Shocks under certain well defined limitations deriving from the finite and curved Bow Shock surface. It is shown that the proton acceleration at the Bow Shock is most efficient for high values of the upstream magnetic field (in general B₁ > 8γ), high upstream plasma speed and expanded Bow Shock fronts, as well as for directions of the induced electric field oriented almost parallel to the flanks of the Bow Shock, i.e. when the drift distance of protons parallel to the electric field at the shock front is considerably smaller than the local radius of curvature of the Bow Shock. The implications of the presented observations of Bow Shock crossings as to the source of the energetic proton intensity enhancements are discussed.     

Energetic Particle Fluxes during the Bastille Day Solar Eruption

We report on our observations of solar energetic particle fluxes of p, He, C, O, Ne, Mg, Si, and Fe ions measured by the Energetic and Relativistic Nucleon and Electron (ERNE) experiment associated with the Bastille Day solar flare and coronal mass ejection (CME) on 14 July 2000. We observed two clear maxima of the Fe/O ratio at the energies 8.5–15 MeV nucl⁻¹. The first Fe/O maximum occurred ∼ 3 hours after the beginning of the particle event, and the second maximum ∼ 22 hours after the first one at the arrival of the shock associated with the Bastille Day eruption. We also observed a change in the energy spectrum of oxygen concurrent with a change in the direction of the interplanetary magnetic field at the start of the second enhancement of the Fe/O ratio. We propose an interpretation of the particle event where observed interplanetary particle fluxes are associated with two different particle sources near the Sun and in interplanetary space. We suggest that heavy ions observed during the first period of the Fe/O enhancement were released when a coronal shock reached a magnetic foot point connected to 1 AU. The second maximum of Fe/O occurred when spacecraft encountered Fe-rich material stored in magnetic field flux tubes early in the event and was possibly reaccelerated by the interplanetary shock.     

Cosmic-ray diffusive acceleration at shock waves with finite upstream and downstream escape boundaries

In the present paper we discuss the modifications introduced into the first-order Fermi shock acceleration process due to a finite extent of diffusive regions near the shock or due to boundary conditions leading to an increased particle escape upstream and/or downstream of the shock. In the simple example of the planar shock wave considered we idealize the escape phenomenon by imposing a particle escape boundary at some distance from the shock. The presence of such a boundary (or boundaries) leads to coupled steepening of the accelerated particle spectrum and decreasing of the acceleration time scale. It allows for a semi-quantitative evaluation and, in some specific cases, also for modelling of the observed steep particle spectra as a result of the first-order Fermi shock acceleration. We also note that the particles close to the upper energy cut-off are younger than the estimate based on the respective acceleration time scale. In Appendix A we present a new time-dependent solution for infinite diffusive regions near the shock allowing for different constant diffusion coefficients upstream and downstream of the shock.     

Multi-spacecraft observations of particle events and interplanetary shocks during November/December 1982

We present a sample of solar energetic particle events observed between November 18 and December 31, 1982 by the HELIOS 1, the VENERA 13, and IMP 8 spacecraft. During the entire time period all three spacecraft were magnetically connected to the western hemisphere of the Sun with varying radial and angular distances from the flares. Eleven proton events, all of them associated with interplanetary shocks, were observed by the three spacecraft. These events are visible in the low-energy (about 4 MeV) as well as the high-energy (30 MeV) protons. In the largest events protons were observed up to energies of about 100 MeV. The shocks were rather fast and in some cases extended to more than 90% east of the flare site. Assuming a symmetrical configuration, this would correspond to a total angular extent of some interplanetary shocks of about 180%. In addition, due to the use of three spacecraft at different locations we find some indication for the shape of the shock front: the shocks are fastest close to the flare normal and are slower at the eastern flank. For particle acceleration we find that close to the flare normal the shock is most effective in accelerating energetic particles. This efficiency decreases for observers connected to the eastern flank of the shock. In this case, the efficiency of shock acceleration for high-energy protons decreases faster than for low-energy protons. Observation of the time-intensity profiles combined with variations of the anisotropy and of the steepness of the proton spectrum allows one in general to define two components of an event which we term solar and interplanetary. We attempt to describe the results in terms of a radially variable efficiency of shock acceleration. Under the assumption that the shock is responsible not only for the interplanetary, but also for the solar component, we find evidence for a very efficient particle acceleration while the shock is still close to the Sun, e.g., in the corona. In addition, we discuss this series of strong flares and interplanetary shocks as a possible source for the formation of a superevent.     

Recurrence of Energetic Particle Flux Anisotropy as Observed by Erne on 9 July 1996

A solar energetic particle event was observed on 9 July 1996, by the ERNE sensors LED and HED on board the SOHO spacecraft. The arrival of the first protons in the energy range >20 MeV took place at 09:55 UT, 43 min after the maximum in the X-ray and Hα radiation of a flare located at S10 W30. The rise phase of the particle intensities at all energies was exceptionally rapid. At 12:50 UT, the intensities dropped in all energy channels. Simultaneously, the magnetic field instrument MFI on board WIND, not far from SOHO, detected a sharp and large change in the magnetic field direction. The analysis of the directional measurements of ERNE in the energy range 14–17 MeV shows the presence of a strong flux anisotropy during the whole period 10:10–12:50 UT. From 12:50 UT until about 16:00 UT, the directional analysis of the proton fluxes gives only a weak anisotropy at the limit of the sensor resolution. Later on, the flux anisotropy was found to recur, indicating a continuous injection of particles into the flux tubes connected to the SOHO spacecraft. These experimental results lead to strict limits on particle injection and transport models. The first period of the anisotropy and its recurrent phase cover 24 hours. This suggests an extended injection of particles. The strength and stability of the anisotropy indicate that, during these periods, SOHO was in an interplanetary sector where the particle transport was almost scatter-free. On the other hand, during the intermediate 3-hr period, we observed particles which traveled in a sector of diffusive transport or which were retarded by magnetic field disturbances not far from the observation site.     

Topology and Dynamics of the 12 December 1980 Magnetic Cloud Event

In this work we have performed an analysis of the interplanetary magnetic field and plasma parameters associated with the 12th December 1980 shock–magnetic cloud event. The study was complemented by including a detailed analysis of the energetic particles fluxes, spectral indexes and directional information. Locally, the magnetic cloud had a latitude of 48°, longitude of 220° and a radius of 0.15 AU. The maximum approach of the ISEE-3 spacecraft to the magnetic cloud axis was 0.047 AU. Moreover, our results suggest that, when encountered by the spacecraft, the magnetic cloud was expanding at a rate of roughly r ₁/r ₀=1.5 (r ₁ is the final radius and r ₀ the initial one) and that it had been expanding during 59 hr. We have also found evidence which indicates that the particle injection inside the cloud depended on the particle energy and pitch angle. These features also suggest that the energetic particle bi-directional fluxes could not be produced by mirroring in the magnetic cloud feet. Instead we think that these fluxes could be generated by the intrinsic properties of the injection mechanism.     

Radar absorption due to a corotating interaction region encounter with Mars detected by MARSIS

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) is a subsurface and topside ionosphere radar sounder aboard the European Space Agency spacecraft Mars Express, in orbit at Mars since 25 December 2003, and in operation since 17 June 2005. The ionospheric sounding mode of MARSIS is capable of detecting the reflection of the sounding wave from the martian surface. This ability has been used in previous work to show that the surface reflection is absorbed and disappears during periods when high fluxes of energetic particles are incident on the ionosphere of Mars. These absorption events are believed to be the result of increased collisional damping of the sounding wave, caused by increased electron density below the spacecraft, in turn caused by impact ionization from the impinging particles. In this work we identify two absorption events that were isolated during periods when the surface reflection is consistently visible and when Mars is nearly at opposition. The visibility of the surface reflection is viewed in conjunction with particle and photon measurements taken at both Mars and Earth. Both absorption events are found to coincide with Earth passing through solar wind speed and ion flux signatures indicative of a corotating interaction region (CIR). The two events are separated by an interval of approximately 27 days, corresponding to one solar rotation. The first of the two events coincides with abruptly enhanced particle fluxes seen in situ at Mars. Simultaneous with the particle enhancement there are an abrupt decrease in the intensity of electron oscillations, typically seen by the Mars Express particle instrument ASPERA-3 between the magnetic pileup boundary and the martian bow shock, and a sharp drop in the solar wind pressure, seen in the proxy quantity based on MGS magnetometer observations. The decrease in oscillation intensity is therefore the probable effect of a relaxation of the martian bow shock. The second absorption event does not show a particle enhancement and complete ASPERA-3 data during that time are unavailable. Other absorption events are the apparent result of solar X-ray and XUV enhancements. We conclude that surface reflection absorption events are sometimes caused by enhanced ionospheric ionization from high energy particles accelerated by the shocks associated with a CIR. A full statistical analysis of CIRs in relation to observed absorption events in conjunction with a quantitative analysis of the deposition of ionization during space weather events is needed for a complete understanding of this phenomenon. If such analyses can be carried out, radar sensing of the martian ionosphere might be useful as a space weather probe.     

Pick-up Ion Measurements in the Heliosphere – A Review

Measurements of the composition and spatial distribution of pick-up ions inside the heliosphere are reviewed. The first interstellar ⁴He⁺pick-up ions were detected with the SULEICA instrument on the AMPTE spacecraft near Earth's orbit. Most data on pick-up ions were taken in the solar-wind and suprathermal energy range of SWICS on Ulysses while the spacecraft cruised from 1.4 to 5.4 AU and explored the high-latitude heliosphere and solar wind from the ecliptic to ± 80^° heliolatitude. This includes the discovery of H⁺, ⁴He⁺⁺, ³He⁺, N⁺,O⁺, and Ne⁺ pick-up ions that originate from the interstellar neutralgas penetrating the heliosphere. From their fluxes properties of the interaction region between the heliosphere and the Local Interstellar Cloud such as the limits on filtration and the strength of the interstellar magnetic field have been revealed. Detailed analysis of the velocity distributions of pick-up ions led to 1) the discovery of a new distinct source, the so-called Inner Source, consisting of atoms released from interstellar and interplanetary dust inside the heliosphere, 2) the determination of pick-up ion transport parameters such as the long mean free path for pitch-angle scattering of order1 AU, and 3) detailed knowledge on the very preferential injection and acceleration of pick-up ions during interplanetary energetic particle events such as Co-rotating Interaction Regions and Coronal Mass Ejections. SWICS measurements have fully confirmed the theory of Fisk, Koslovsky, and Ramaty that pick-up ions derived from the interstellar gas are the dominant source of the Anomalous Cosmic Rays; they are pre-accelerated inside the heliosphere and re-accelerated at the solar-wind Termination Shock according to Pesses, Eichler, and Jokipii. The data indicate that the Inner Source of pick-up ionsis largely responsible for the occurence of C⁺ in the Anomalous Cosmic Rays. The abundances of recently discovered Inner-Source Mg⁺ and Si⁺ are solar-wind like and consistent with their abundances in the energetic particles associated with Co-rotating Interaction Regions. Knowledge on the injection and acceleration processes in Co-rotating Interaction Regions is applied to discuss the current observational evidence for the Interplanetary Focusing Cone of the interstellar neutral gas due to the Sun's gravitational force. The 25–150 keV/amu suprathermal ⁴He⁺ pick-up ion fluxes measured by CELIAS/STOF on board SOHO over 360^° of ecliptic longitude represent a `local' ionization and acceleration of interstellar atoms at 1 AU or smaller heliocentric distances. Completing the first limited data set of SULEICA/AMPTE on ⁴He⁺ pick-up ions they indicate a density enhancement in the Interplanetary Focusing Cone which is confirmed by recent SWICS/ACE data. Clear evidence for signatures in ecliptic longitude are found in the data on energetic neutral H fluxes observed with the CELIAS/HSTOF sensor on board SOHO. These fluxes are enhanced in the upstream and downstream directions of the interstellar wind. Detection of energetic H atoms, which propagate unaffected by the Heliospheric Magnetic Field, provided for the first time a diagnostic tool for observations near Earth to analyze the structure in ecliptic longitude of the interface region between the heliosphere and the Local Interstellar Cloud. This revised version was published online in July 2006 with corrections to the Cover Date.     

Proton activity of the Sun in current solar cycle 24

We present a study of seven large solar proton events in the current solar cycle 24(from 2009 January up to the current date). They were recorded by the GOES spacecraft with the highest proton fluxes being over 200 pfu for energies 10 Me V. In situ particle measurements show that:(1) The profiles of the proton fluxes are highly dependent on the locations of their solar sources, namely flares or coronal mass ejections(CMEs), which confirms the "heliolongitude rules" associated with solar energetic particle fluxes;(2) The solar particle release(SPR) times fall in the decay phase of the flare emission, and are in accordance with the times when the CMEs travel to an average height of 7.9 solar radii; and(3) The time differences between the SPR and the flare peak are also dependent on the locations of the solar active regions. The results tend to support the scenario of proton acceleration by the CME-driven shock,even though there exists a possibility of particle acceleration at the flare site, with subsequent perpendicular diffusion of accelerated particles in the interplanetary magnetic field. We derive the integral time-of-maximum spectra of solar protons in two forms: a single power-law distribution and a power law roll-over with an exponential tail. It is found that the unique ground level enhancement that occurred in the event on 2012 May 17 displays the hardest spectrum and the largest roll-over energy which may explain why this event could extend to relativistic energies.     

Slow-mode shock candidate in the Jovian magnetosheath

We discuss some interesting plasma observations in the Jovian magnetosheath by the onboard plasma instruments of the Cassini spacecraft during the 2000-2001 Jupiter flyby. We propose that the observations are consistent with a slow-mode shock transition. In the terrestrial magnetosheath, a number of observations have been made that are consistent with slow-mode waves or shocks. In addition, a number of observations have established that, at least occasionally, slow-mode structures form at the plasma sheet-lobe boundary in the terrestrial magnetotail, related to X lines associated with reconnection. There has been only one previously reported observation of a slow-mode shock-like transition in the Jovian plasma environment. This observation was made in the dayside magnetosheath. The observation we report here was made well downstream of the magnetosphere in Jupiter’s magnetosheath, at local time ∼19:10. For our analysis we have used the data from the Cassini Plasma Spectrometer (CAPS) the Magnetospheric Imaging Instrument (MIMI) and the Magnetometer (MAG). The bow shock crossings observed by Cassini ranged downstream to −600 R_J from the planet     

OBLIQUE EMISSION OF WHISTLER-MODE WAVES AROUND INTERPLANETARY SHOCKS

We present a study of whistler-mode wave generation and wave particle interaction in the vicinity of interplanetary shocks as observed by the Ulysses spacecraft. Generally the whistler-mode waves (measured in the frequency range 0.22–448 Hz) are observed downstream of the shocks where they persist for some hours. From the electron distribution functions (EDF) in the energy range 1.6 to 862 eV measured by the spacecraft, we compute the wave growth rate of the electromagnetic electron cyclotron and Landau instabilities for the case of oblique propagation of the waves with respect to the interplanetary magnetic field (IMF) B. In general, in agreement with the wave measurements, the instability grows mostly downstream of the shock fronts. Following the wave activity, velocity space diffusion of the electrons results in a marginally stable state with some sporadic fluctuations.     

Shock surfing acceleration

Analytical and numerical analysis identify shock surfing acceleration as an ideal pre-energization mechanism for the slow pick-up ions at quasiperpendicular shocks. After gaining sufficient energy by shock surfing, pick-up ions undergo diffusive acceleration to reach their observed energies. Energetic ions upstream of the cometary bow shock, acceleration of solar energetic particles by magnetosonic waves in corona, ion enhancement in interplanetary shocks, generation of anomalous cosmic rays from interstellar pick-up ions at the termination shock are some of the cases where shock surfing acceleration apply. Inclusion of the lower-hybrid wave turbulence into the laminar model of shock surfing can explain the preferential acceleration of heavier particles as observed by Voyager at the termination shock. At relativistic energies, unlimited acceleration of ions is theoretically possible; because for sufficiently strong shocks main limitation of the mechanism, caused by the escape of accelerated particles downstream of the shock during acceleration no longer exists.     

The Bastille day Magnetic Clouds and Upstream Shocks: Near-Earth Interplanetary Observations

The energetic charged particle, interplanetary magnetic field, and plasma characteristics of the `Bastille Day' shock and ejecta/magnetic cloud events at 1 AU occurring over the days 14–16 July 2000 are described. Profiles of MeV (WIND/LEMT) energetic ions help to organize the overall sequence of events from the solar source to 1 AU. Stressed are analyses of an outstanding magnetic cloud (MC2) starting late on 15 July and its upstream shock about 4 hours earlier in WIND magnetic field and plasma data. Also analyzed is a less certain, but likely, magnetic cloud (MC1) occurring early on 15 July; this was separated from MC2 by its upstream shock and many heliospheric current sheet (HCS) crossings. Other HCS crossings occurred throughout the 3-day period. Overall this dramatic series of interplanetary events caused a large multi-phase magnetic storm with min Dst lower than −300 nT. The very fast solar wind speed (≥ 1100 km s⁻¹) in and around the front of MC2 (for near average densities) was responsible for a very high solar wind ram pressure driving in the front of the magnetosphere to geocentric distances estimated to be as low as ≈ 5 R _E, much lower than the geosynchronous orbit radius. This was consistent with magnetic field observations from two GOES satellites which indicated they were in the magnetosheath for extended times. A static force-free field model is used to fit the two magnetic cloud profiles providing estimates of the clouds' physical and geometrical properties. MC2 was much larger than MC1, but their axes were nearly antiparallel, and their magnetic fields had the same left-handed helicity. MC2's axis and its upstream shock normal were very close to being perpendicular to each other, as might be expected if the cloud were driving the shock at the time of observation. The estimated axial magnetic flux carried by MC2 was 52×10²⁰ Mx, which is about 5 times the typical magnetic flux estimated for other magnetic clouds in the WIND data over its first 4 years and is 17 times the flux of MC1. This large flux is due to both the strong axially-directed field of MC2 (46.8 nT on the axis) and the large radius (R ₀=0.189 AU) of the flux tube. MC2's average speed is consistent with the expected transit time from a halo-CME to which it is apparently related.     

The Large Longitudinal Spread of Solar Energetic Particles During the 17 January 2010 Solar Event

We investigate multi-spacecraft observations of the 17 January 2010 solar energetic particle event. Energetic electrons and protons have been observed over a remarkable large longitudinal range at the two STEREO spacecraft and SOHO, suggesting a longitudinal spread of nearly 360 degrees at 1?AU. The flaring active region, which was on the backside of the Sun as seen from Earth, was separated by more than 100 degrees in longitude from the magnetic footpoints of each of the three spacecraft. The event is characterized by strongly delayed energetic particle onsets with respect to the flare and only small or no anisotropies in the intensity measurements at all three locations. The presence of a coronal shock is evidenced by the observation of a type II radio burst from the Earth and STEREO-B. In order to describe the observations in terms of particle transport in the interplanetary medium, including perpendicular diffusion, a 1D model describing the propagation along a magnetic field line (model 1) (Dr?ge, Astrophys. J. 589, 1027??C?1039, 2003) and the 3D propagation model (model 2) by Dr?ge et?al. (Astrophys. J. 709, 912??C?919, 2010) including perpendicular diffusion in the interplanetary medium have been applied. While both models are capable of reproducing the observations, model 1 requires injection functions at the Sun of several hours. Model 2, which includes lateral transport in the solar wind, reveals high values for the ratio of perpendicular to parallel diffusion. Because we do not find evidence for unusual long injection functions at the Sun, we favor a scenario with strong perpendicular transport in the interplanetary medium as an explanation for the observations.     

Fundamentals of collisionless shocks for astrophysical application, 2. Relativistic shocks

In this concise review of the recent developments in relativistic shock theory in the Universe we restrict ourselves to shocks that do not exhibit quantum effects. On the other hand, emphasis is given to the formation of shocks under both non-magnetised and magnetised conditions. We only briefly discuss particle acceleration in relativistic shocks where much of the results are still preliminary. Analytical theory is rather limited in predicting the real shock structure. Kinetic instability theory is briefed including its predictions and limitations. A recent self-similar relativistic shock theory is described which predicts the average long-term shock behaviour to be magnetised and to cause reasonable power-law distributions for energetic particles. The main focus in this review is on numerical experiments on highly relativistic shocks in (i) pair and (ii) electron-nucleon plasmas and their limitations. These simulations do not validate all predictions of analytic and self-similar theory and so far they do not solve the injection problem and the self-modification by self-generated cosmic rays. The main results of the numerical experiments discussed in this review are: (i) a confirmation of shock evolution in non-magnetised relativistic plasma in 3D due to either the lepton-Weibel instability (in pair plasmas) or to the ion-Weibel instability; (ii) the sensitive dependence of shock formation on upstream magnetisation which causes suppression of Weibel modes for large upstream magnetisation ratios σ>10⁻³; (iii) the sensitive dependence of particle dynamics on the upstream magnetic inclination angle θ _Bn , where particles of θ _Bn >34° cannot escape upstream, leading to the distinction between ‘subluminal’ and ‘superluminal’ shocks; (iv) particles in ultra-relativistic shocks can hardly overturn the shock and escape to upstream; they may oscillate around the shock ramp for a long time, so to speak ‘surfing it’ and thereby becoming accelerated by a kind of SDA; (v) these particles form a power-law tail on the downstream distribution; their limitations are pointed out; (vi) recently developed methods permit the calculation of the radiation spectra emitted by the downstream high-energy particles; (vii) the Weibel-generated downstream magnetic fields form large-amplitude vortices which could be advected by the downstream flow to large distances from the shock and possibly contribute to an extended strong field region; (viii) if cosmic rays are included, Bell-like modes can generate upstream magnetic turbulence at short and, by diffusive re-coupling, also long wavelengths in nearly parallel magnetic field shocks; (ix) advection of such large-amplitude waves should cause periodic reformation of the quasi-parallel shock and eject large-amplitude magnetic field vortices downstream where they contribute to turbulence and to maintaining an extended region of large magnetic fields.