The temperature map of the loop and coronal sources of an occulted solar flare

An occulted solar flare occurred at about 06:07 UT on 2002, November 2. The RHESSI X-ray images show two separate parts. The lower part consists of a complete loop and the upper part a coronal source which well extends above the solar limb. The loop source shrank for about 3 min with a speed of ∼24 km s⁻¹ during the early impulsive phase and then expanded at ∼7 km s⁻¹, while the coronal source presented an upward motion at about 6 km s⁻¹. We obtained the temperature map of the loop source from RHESSI image spectrum. The temperature of the loop increases with altitude, indicating that the reconnection X-point of this flare is located above the loop source. However, the apparent coronal source is the top of another independent large-scale loop.     

The evolution of a complex solar radio burst corresponding to special configuration of microwave sources

A complex radio burst associated with periodic (∼1 and 6 min) pulsations and several kinds fine structures, e.g., normal- and reverse-drifting type III bursts, zebra patterns, and slowly drifting structure was observed with the radio spectrometers (1.0–2.0, 2.6–3.8, 5.2–7.6, and 0.65–1.5 GHz) at the National Astronomical Observatories of China (NAOC) in Beijing and Yunnan on 19 October 2001. In combination with the images of 17 and 34 GHz from NoRH and the magnetograms from MDI we reveal the existence and evolution of preexisting and new emerging sources, and find the horseshoe-shaped structure of microwave sources intensity during the late phase of the burst. Through the detailed comparison of the evolution of each source with the time profiles of radio bursts corresponding to these sources we indicate that the intimate correlation between the microwave sources evolution and the generation of the radio burst associated fine structures. Some fine structures can be considered as the MHD turbulence and plasma emission mechanism, based on the anisotropic beam instability and hybrid waves generations. From the characteristics of observations we may presume that the coronal magnetic structures should contain an extended coronal loop system and multiple discrete electrons acceleration/injection sites. The mechanisms of this complex radio burst are deal with the incoherent gyrosynchrotron emission from the trapped electrons and the coherent plasma emission from the non trapped electrons.     

A spectroscopic measurement of high velocity spray plasma from an M-class flare and coronal mass ejection

The M1.5-class flare and associated coronal mass ejection (CME) of 16 February 2011 was observed with the Extreme ultraviolet Imaging Spectrometer on board the Hinode spacecraft. Spray plasma associated with the CME is found to exhibit a Doppler blue-shift of 850 km s^?1 – one of the largest values reported from spectroscopy of the solar disk and inner corona. The observation is unusual in that the emission line (Fe xii 193.51 Å) is not observed directly, but the Doppler shift is so large that the blue-shifted component appears in a wavelength window at 192.82 Å, intended to observe lines of O v, Fe xi and Ca xvii. The Fe xii 195.12 Å emission line is used as a proxy for the rest component of 193.51 Å. The observation highlights the risks of using narrow wavelength windows for spectrometer observations when observing highly-dynamic solar phenomena. The consequences of large Doppler shifts for ultraviolet solar spectrometers, including the upcoming Multi-slit Solar Explorer (MUSE) mission, are discussed.     

The coronal loop controversy: TRACE analysis

The temperature distribution along coronal loops provides an important clue for solving the coronal heating problem. Recent analysis, however, has produced conflicting results. Here, we analyze in detail one component of this analysis – the effect of background subtraction on the temperature of loops observed with the Transition Region and Coronal Explorer (TRACE). Specifically, we selected 10 coronal loops that were visible in the TRACE 171 Å and 195 Å passbands. We chose between 20 and 30 pixel along each loop and background pixels to correspond with the loop pixels. Temperature analysis was done in three different ways: (1) standard TRACE analysis of the loop pixels with no background subtraction; (2) constant background subtraction for each TRACE image; (3) pixel pair background subtraction. Each method produced a temperature estimate for the selected pixels. We find that a flat line is an excellent fit to the temperature results – the analysis indicates that the temperature of the loop is uniform along the length visible by TRACE. However, if we select random pixels and plot the temperature results in the same way, these pixels indicate that the temperature of this “structure” is also uniform. We conclude therefore, that in the cases considered here, the image ratio analysis does not produce a physically meaningful value of plasma temperature; in addition, background subtraction makes no significant difference to the temperatures results.     

Gamma-ray observations of explosive nucleosynthesis products

In this review I discuss the various γ-ray emission lines that can be expected and, in some cases have been observed, from radioactive explosive nucleosynthesis products. The most important γ-ray lines result from the decay chains of ⁵⁶Ni, ⁵⁷Ni, and ⁴⁴Ti. ⁵⁶Ni is the prime explosive nucleosynthesis product of Type Ia supernovae, and its decay determines to a large extent the Type Ia light curves. ⁵⁶Ni is also a product of core-collapse supernovae, and in fact, γ-ray line emission from its daughter product, ⁵⁶Co, has been detected from SN1987A by several instruments. The early occurrence of this emission was surprising and indicates that some fraction of ⁵⁶Ni, which is synthesized in the innermost supernova layers, must have mixed with the outermost supernova ejecta.Special attention is given to the γ-ray line emission of the decay chain of ⁴⁴Ti (⁴⁴Ti → ⁴⁴Sc → ⁴⁴Ca), which is accompanied by line emission at 68, 78, and 1157 keV. As the decay time of ⁴⁴Ti is ∼86 yr, one expects this line emission from young supernova remnants. Although the ⁴⁴Ti yield (typically 10⁻⁵–10⁻⁴M_⊙) is not very high, its production is very sensitive to the energetics and asymmetries of the supernova explosion, and to the mass cut, which defines the mass of the stellar remnant. This makes ⁴⁴Ti an ideal tool to study the inner layers of the supernova explosion. This is of particular interest in light of observational evidence for asymmetric supernova explosions.The γ-ray line emission from ⁴⁴Ti has so far only been detected from the supernova remnant Cas A. I discuss these detections, which were made by COMPTEL (the 1157 keV line) and BeppoSAX (the 68 and 78 keV lines), which, combined, give a flux of (2.6 ± 0.4 ± 0.5) × 10⁻⁵ ph cm⁻² s⁻¹ per line, suggesting a ⁴⁴Ti yield of (1.5 ± 1.0) × 10⁻⁴M_⊙. Moreover, I present some preliminary results of Cas A observations by INTEGRAL, which so far has yielded a 3σ detection of the 68 keV line with the ISGRI instrument with a flux that is consistent with the BeppoSAX detections. Future observations by INTEGRAL-ISGRI will be able to constrain the continuum flux above 90 keV, as the uncertainty about the continuum shape, is the main source of systematic error for the 68 and 78 keV line flux measurements. Moreover, with the INTEGRAL-SPI instrument it will be possible to measure or constrain the line broadening of the 1157 keV line. A preliminary analysis of the available data indicates that narrow line emission (i.e., Δv < 1000 km s⁻¹) can be almost excluded at the 2σ level, for an assumed line flux of 1.9 × 10⁻⁵ ph cm⁻² s⁻¹.     

The thermal X-ray spectrum of the 2003 April 26 solar flare

Observations and their analysis of the thermal X-ray spectrum of the M2 flare on 2003 April 26 are described. The spectrum observed by the RHESSI mission cover the energy range from ∼5 to ∼50 keV. With its ∼1-keV spectral resolution, intensities and equivalent widths of two line complexes, the Fe line group at 6.7 keV (mostly due to Fe xxv lines and Fe xxiv satellites) and the Fe/Ni line group at 8 keV (mostly due to higher-excitation Fe xxv lines and Ni xxvii lines) were obtained as a function of time through a number of flares. The abundance of Fe can also be determined from RHESSI spectra; it appears to be consistent with a coronal value for at least some times during the flare. Comparisons of RHESSI spectra with those from the RESIK crystal spectrometer on CORONAS-F show very satisfactory agreement, giving much confidence in the intensity calibration of both instruments.     

Excitation of flare-induced waves in coronal loops and the effects of radiative cooling

EUV imaging observations from several space missions (SOHO/EIT, TRACE, and SDO/AIA) have revealed a presence of propagating intensity disturbances in solar coronal loops. These disturbances are typically interpreted as slow magnetoacoustic waves. However, recent spectroscopic observations with Hinode/EIS of active region loops revealed that the propagating intensity disturbances are associated with intermittent plasma upflows (or jets) at the footpoints which are presumably generated by magnetic reconnection. For this reason, whether these disturbances are waves or periodic flows is still being studied. This study is aimed at understanding the physical properties of observed disturbances by investigating the excitation of waves by hot plasma injections from below and the evolution of flows and wave propagation along the loop. We expand our previous studies based on isothermal 3D MHD models of an active region to a more realistic model that includes full energy equation accounting for the effects of radiative losses. Computations are initialized with an equilibrium state of a model active region using potential (dipole) magnetic field, gravitationally stratified density and temperature obtained from the polytropic equation of state. We model an impulsive injection of hot plasma into the steady plasma outflow along the loops of different temperatures, warm (～1 MK) and hot (～6 MK). The simulations show that hot jets launched at the coronal base excite slow magnetoacoustic waves that propagate to high altitudes along the loops, while the injected hot flows decelerate rapidly with heights. Our results support that propagating disturbances observed in EUV are mainly the wave features. We also find that the effect of radiative cooling on the damping of slow-mode waves in 1–6 MK coronal loops is small, in agreement with the previous conclusion based on 1D MHD models.     

Spectral analysis of a two-ribbon microflare

Using high-resolution Hα, CaII 8542 Å and FeI 6302.5 Å Stokes spectral data obtained simultaneously with THEMIS in 2002 September, we have analyzed the spectra and the characteristics of a two-ribbon microflare (MF). The hard X-ray emission provides evidence of non-thermal particle acceleration in the microflare. The two-ribbons are located on either sides of the magnetic polarity inversion line. The non-thermal characteristics mainly appeared at the outer edges of the flare ribbons. It indicates that the instantaneous magnetic reconnection and the particle acceleration mainly took place at the outer edges of the flare ribbons. Using the Hα and CaII 8542 Å line profiles and the non-LTE calculation, we obtain the semi-empirical atmospheric model for the bright kernel of the MF. The result indicates that the temperature enhancement in the chromosphere is about 2000–2500 K.     

Signatures of magnetic reconnection in solar radio observations?

Magnetic reconnection occurs during eruptive processes (flares, CMEs) in the solar corona. This leads to a change of magnetic connectivity. Nonthermal electrons propagate along the coronal magnetic field thereby exciting dm- and m-wave radio burst emission after acceleration during reconnection or other energy release processes in heights of some Mm to ⩾700 Mm. We summarize the results of some case studies which can be interpreted as radio evidence of magnetic reconnection: under certain conditions, simple spectral structures (pulsation pulses, reverse drift bursts) are formed by simultaneously acting but widely spaced radio sources. Narrowband spikes are emitted as a side-effect during large-scale coronal loop collisions. In dynamic radio spectra, the lower fast mode shock formed in the reconnection outflow appears as type II burst-like but nondrifting emission lane. It has been several times observed at the harmonic mode of the local plasma frequency between 250 and 500 MHz and at heights of ≈200 Mm.     

The magnetic field topology associated with two M flares

On 27 October, 2003, two GOES M-class flares occurred in an interval of 3 h in active region NOAA 10486. The two flares were confined and their associated brightenings appeared at the same location, displaying a very similar shape both at the chromospheric and coronal levels. We focus on the analysis of magnetic field (SOHO/MDI), chromospheric (HASTA, Kanzelhöhe Solar Observatory, TRACE) and coronal (TRACE) observations. By combining our data analysis with a model of the coronal magnetic field, we compute the magnetic field topology associated with the two M flares. We find that both events can be explained in terms of a localized magnetic reconnection process occurring at a coronal magnetic null point. This null point is also present at the same location one day later, on 28 October, 2003. Magnetic energy release at this null point was proposed as the origin of a localized event that occurred independently with a large X17 flare on 28 October, 2003 [Mandrini, C.H., Démoulin, P., Schmieder, B., Deluca, E., Pariat, E., Uddin, W. Companion event and precursor of the X17 flare on 28 October, 2003. Solar Physics, 238, 293–312, 2006], at 11:01 UT. The three events, those on 27 October and the one on 28 October, are homologous. Our results show that coronal null points can be stable topological structures where energy release via magnetic reconnection can happen, as proposed by classical magnetic reconnection models.     

Low-frequency electric field and density fluctuation measurements on Solar Orbiter

Solar Orbiter will orbit the Sun down to a distance of 0.22 AU allowing detailed in situ studies of important but unexplored regions of the solar wind in combination with coordinated remote sensing of the Sun. In-situ measurements require high quality measurements of particle distributions and electric and magnetic fields. We show that such important scientific topics as the identification of coronal heating remnants, solar wind turbulence, magnetic reconnection and shock formation within coronal mass ejections all require electric field and plasma density measurements in the frequency range from DC up to about 100 Hz. We discuss how such measurements can be achieved using the double-probe technique. We sketch a few possible antenna design solutions.     

Evidence for a solar coronal thick-target hard X-ray source observed by RHESSI

We study a solar flare hard X-ray (HXR) source observed by the Reuven Ramaty high energy solar spectroscopic imager (RHESSI) in which the HXR emission is almost entirely in a coronal loop so dense as to be collisionally thick at electron energies up to ∼45−60 keV. This contrasts with most events previously reported in which the HXR emission is primarily from the loop footpoints in the collisionally dense chromosphere. In particular, we show that the high loop column densities inferred from the GOES and RHESSI soft X-ray emission measure and the volume of the flare loop are consistent with the coronal thick-target interpretation of the HXR images and spectra. The high column densities observed already at the very beginning of the impulsive phase are explained by chromospheric evaporation during a preflare which, as Nobeyama 17 GHz radio images reveal, took place in the same set of nested loops as the main flare.     

Hard X-ray footpoint motions in solar flares: Comparing magnetic reconnection models with observations

Hard X-ray observations from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) of the October 29, 2003 GOES X10 two-ribbon flare are used together with magnetic field observations from the Michelson Doppler Imager (MDI) onboard SoHO to compare footpoint motions with predictions from magnetic reconnection models. The temporal variations of the velocity v of the hard X-ray footpoint motions and the photospheric magnetic field strength B in footpoints are investigated. The underlying photospheric magnetic field strength is generally higher (B ∼ 700–1200 G) in the slower moving (v ∼ 20–50 km s⁻¹) western footpoint than in the faster (v ∼ 20–100 km s⁻¹) moving eastern source (∼100–600 G). Furthermore, a rough temporal correlation between the HXR flux and the product vB² is observed.     

On some large solar energetic particle events: Direct impulsive component and broad associated coronal mass ejections

Using the proton intensity and X-ray flux data from the GOES, combined with the observations of the associated solar eruptions by the Large Angle and Spectrometric Coronagraph Experiment (LASCO) on board the Solar and Heliospheric Observatory (SOHO), 14 large SEP events occurring in the period 2000 January–2002 April have been studied. It is found that: (1) events with the SEPs increasing shortly after the maximum of their parent flares (<1 h; hereafter prompt events) have rapid and great (up to four orders of magnitude) SEP increments in high-energy channels (> ∼100 MeV); however, for events whose onset of the SEP injection lags the flare maximum for a long time (>3 h; hereafter delayed events), the high-energy SEPs show no obvious enhancements (within one order of magnitude); (2) peak intensity of the prompt events is distinctly larger than that of the delayed events; (3) CMEs associated with the poorly magnetically connected events (source region <W30°) in our survey are all halo CMEs. From these observational differences, we propose a special scenario of the production of the largest SEP events: both CMEs and flares are induced in the same coronal process; high-energy particles accelerated in the reconnection region can escape easily from the open field lines and/or be transported by fast CMEs into interplanetary space, indicating a direct impulsive component in large gradual SEP events. Meanwhile, the broad width of the associated CMEs implies that the CME width is more important in SEP events production than previously considered.     

Comparison of vertical E × B drift velocities and ground-based magnetometer observations of DELTA H in the low latitude under geomagnetically disturbed conditions

In the present work, we analyzed the daytime vertical E × B drift velocities obtained from Jicamarca Unattended Long-term Ionosphere Atmosphere (JULIA) radar and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca, and Piura in Peru for 22 geomagnetically disturbed events in which either SC has occurred or Dst_max < ?50 nT during the period 2006–2011. The ΔH component of geomagnetic field is measured as the differences in the magnitudes of horizontal H component between magnetometer placed directly on the magnetic equator and one displaced 6–9° away. It will provide a direct measure of the daytime electrojet current, due to the eastward electric field. This will in turn gives the magnitude of vertical E × B drift velocity in the F region. A positive correlation exists between peak values of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of the events. It was observed that 45% of the events have daytime vertical E × B drift velocity peak in the magnitude range 10–20 m/s and 20–30 m/s and 20% have peak ΔH in the magnitude range 50–60 nT and 80–90 nT. It was observed that the time of occurrence of the peak value of both the vertical E × B drift velocity and the ΔH have a maximum (40%) probability in the same time range 11:00–13:00 LT. We also investigated the correlation between E × B drift velocity and D_st index and the correlation between delta H and D_st index. A strong positive correlation is found between E × B drift and D_st index as well as between delta H and D_st Index. Three different techniques of data analysis – linear, polynomial (order 2), and polynomial (order 3) regression analysis were considered. The regression parameters in all the three cases were calculated using the Least Square Method (LSM), using the daytime vertical E × B drift velocity and ΔH. A formula was developed which indicates the relationship between daytime vertical E × B drift velocity and ΔH, for the disturbed periods. The E × B drift velocity was then evaluated using the formulae thus found for the three regression analysis and validated for the ‘disturbed periods’ of 3 selected events. The E × B drift velocities estimated by the three regression analysis have a fairly good agreement with JULIA radar observed values under different seasons and solar activity conditions. Root Mean Square (RMS) errors calculated for each case suggest that polynomial (order 3) regression analysis provides a better agreement with the observations from among the three.     

Centre-to-limb variation of photospheric facular radiance and image resolution

We study the effect of the angular resolution on the determination of the angular properties of the facular radiance. We analyze photospheric intensity in the continuum, around the Ni 676.8 nm line, and longitudinal magnetic field along the line of sight, measured by the MDI instrument aboard SOHO with two spatial resolutions, 4″ and 1.2″ (2″ and 0.6″ pixels, respectively). The effect of the limited photometric sensitivity of the instrument and the limited information on the angular structure of the magnetic field tubes are considered. Our study of the high-resolution data shows that intensity contrast of magnetic features between 80 and 600 Gauss increases from centre to limb up to a maximum that occurs at higher heliocentric angles (θ) when obtained with higher resolution data than for lower resolution data. There is a suggestion that at heliocentric angles below about 75° there is only a monotonic increase in the contrast as one goes from cos (θ) = 1 to cos (θ) = 0.2.     

The plerionic supernova remnant G21.5-0.9: In and out

The absence of a supernova remnant (SNR) shell surrounding the Crab and other plerions (pulsar wind nebulae) has been a mystery for three decades. G21.5-0.9 is a particularly intriguing plerionic SNR in which the central powering engine is not yet detected. Early CHANDRA observations revealed a faint extended X-ray halo which was suggested to be associated with the SNR shell; however its spectrum was non-thermal, unlike what is expected from an SNR shell. On the other hand, a plerionic origin to the halo is problematic since the X-ray plerion would be larger than the radio plerion. We present here our analysis of an integrated 245 ks of archival CHANDRA data acquired with the High-Resolution Camera (HRC) and 520 ks acquired with the Advanced CCD Imaging Spectrometer (ACIS). This study provides the deepest and highest resolution images obtained to date. The resulting images reveal for the first time: (1) a limb-brightened morphology in the eastern section of the halo, and (2) a rich structure in the inner (40″-radius) bright plerion including wisps and a double-lobed morphology with an axis of symmetry running in the northwest–southeast direction. Our spatially resolved spectroscopic study of the ACIS-I data indicates that the photon index steepens with increasing distance from the central point source out to a radius of 40″ then becomes constant at ∼2.4 in the X-ray halo (for a column density N_H = 2.2 × 10²² cm⁻²). No line emission was found from the eastern limb; however marginal evidence for line emission in the halo’s northern knots was found. This study illustrates the need for deep CHANDRA observations to reveal the missing SNR material in Crab-like plerions.     

The 16-day waves in the mesosphere and lower thermosphere over Wuhan (30.6°N, 114.5°E) and Adelaide (35°S, 138°E)

Winds from a meteor radar at Wuhan (30.6°N, 114.5°E) and a MF radar at Adelaide (35°S, 138°E) are used to study the 16-day waves in the mesosphere and lower thermosphere (MLT). The height range is 78–98 km at Wuhan and 70–98 km at Adelaide. By comparison, it is found that the zonal components at both sites are generally larger than the meridional ones, and eastward motion of the zonal background winds is favorable for the 16-day waves penetration to the MLT region. The zonal maximum amplitude appears in the autumn (September–October) around 86–98 km at Wuhan and in the winter months and early spring (July–October) around 72–82 km at Adelaide. Differences are found in wave amplitudes and time of appearance between the two years of 2002 and 2003. In 2003, the intensity of the wave amplitudes is relatively smaller than that for 2002 at both sites. The summer 16-day waves are comparatively weaker at Adelaide in both years, but stronger in 2002 at Wuhan near the mesopause and the lower thermosphere (86–98 km). The strong summer waves at Wuhan may come from the winter southern hemisphere.     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.     

GNSS global real-time augmentation positioning: Real-time precise satellite clock estimation,prototype system construction and performance analysis

Lots of ambiguities in un-differenced (UD) model lead to lower calculation efficiency, which isn’t appropriate for the high-frequency real-time GNSS clock estimation, like 1 Hz. Mixed differenced model fusing UD pseudo-range and epoch-differenced (ED) phase observations has been introduced into real-time clock estimation. In this contribution, we extend the mixed differenced model for realizing multi-GNSS real-time clock high-frequency updating and a rigorous comparison and analysis on same conditions are performed to achieve the best real-time clock estimation performance taking the efficiency, accuracy, consistency and reliability into consideration. Based on the multi-GNSS real-time data streams provided by multi-GNSS Experiment (MGEX) and Wuhan University, GPS + BeiDou + Galileo global real-time augmentation positioning prototype system is designed and constructed, including real-time precise orbit determination, real-time precise clock estimation, real-time Precise Point Positioning (RT-PPP) and real-time Standard Point Positioning (RT-SPP). The statistical analysis of the 6 h-predicted real-time orbits shows that the root mean square (RMS) in radial direction is about 1–5 cm for GPS, Beidou MEO and Galileo satellites and about 10 cm for Beidou GEO and IGSO satellites. Using the mixed differenced estimation model, the prototype system can realize high-efficient real-time satellite absolute clock estimation with no constant clock-bias and can be used for high-frequency augmentation message updating (such as 1 Hz). The real-time augmentation message signal-in-space ranging error (SISRE), a comprehensive accuracy of orbit and clock and effecting the users’ actual positioning performance, is introduced to evaluate and analyze the performance of GPS + BeiDou + Galileo global real-time augmentation positioning system. The statistical analysis of real-time augmentation message SISRE is about 4–7 cm for GPS, whlile 10 cm for Beidou IGSO/MEO, Galileo and about 30 cm for BeiDou GEO satellites. The real-time positioning results prove that the GPS + BeiDou + Galileo RT-PPP comparing to GPS-only can effectively accelerate convergence time by about 60%, improve the positioning accuracy by about 30% and obtain averaged RMS 4 cm in horizontal and 6 cm in vertical; additionally RT-SPP accuracy in the prototype system can realize positioning accuracy with about averaged RMS 1 m in horizontal and 1.5–2 m in vertical, which are improved by 60% and 70% to SPP based on broadcast ephemeris, respectively.