A New Prediction Method for the Arrival Time of Interplanetary Shocks

Solar transient activities such as solar flares, disappearing filaments, and coronal mass ejections (CMEs) are solar manifestations of interplanetary (IP) disturbances. Forecasting the arrival time at the near Earth space of the associated interplanetary shocks following these solar disturbances is an important aspect in space weather forecasting because the shock arrival usually marks the geomagnetic storm sudden commencement (SSC) when the IMF Bz component is appropriately southward and/or the solar wind dynamic pressure behind the shock is sufficiently large. Combining the analytical study for the propagation of the blast wave from a point source in a moving, steady-state, medium with variable density (wei, 1982; wei and dryer 1991) with the energy estimation method in the ISPM model (smith and dryer 1990, 1995), we present a new shock propagation model (called SPM below) for predicting the arrival time of interplanetary shocks at Earth. The duration of the X-ray flare, the initial shock speed and the total energy of the transient event are used for predicting the arrival of the associated shocks in our model. Especially, the background speed, i.e., the convection effect of the solar wind is considered in this model. Applying this model to 165 solar events during the periods of January 1979 to October 1989 and February 1997 to August 2002, we found that our model could be practically equivalent to the prevalent models of STOA, ISPM and HAFv.2 in forecasting the shock arrival time. The absolute error in the transit time in our model is not larger than those of the other three models for the same sample events. Also, the prediction test shows that the relative error of our model is ≤10% for 27.88% of all events, ≤30% for 71.52%, and ≤50% for 85.46%, which is comparable to the relative errors of the other models. These results might demonstrate a potential capability of our model in terms of real-time forecasting.     

Implications of Proton Anisotropy Development Observed by the Erne Instrument During the 9 July 1996 Solar Particle Event

We analysed the solar particle event following the 9 July 1996 solar flare. High-energy protons were detected by the ERNE instrument on board SOHO. Anisotropy of arriving protons revealed very peculiar non-monotonic development. A short period of almost isotropic distribution was imbedded into the prolonged period of beam-like distribution of 14–17 MeV protons. This implies the existence of a narrow magnetic channel with a much smaller mean free path than in the surrounding quiet solar wind plasma. We used Monte Carlo simulations of interplanetary transport to fit the observed anisotropies and intensity–time profiles. Proton injection and transport parameters are estimated. The injection scenario is found to be very close to the scenario of the 24 May 1990 event, but the intensity and the interplanetary transport parameters are different. The extreme anisotropy observed implies prolonged injection of high-energy protons at the Sun and at the interplanetary shock front, and either a very large mean free path (≥ 5 AU) outside the slow transport channel, or alternatively, a somewhat smaller mean free path (≈2 AU) and enhanced focusing between the Sun and the Earth.     

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.     

Predicting the Arrival Time of Shock Passages at Earth

The purpose of this parametric study is to predict the arrival time at Earth of shocks due to disturbances observed on the Sun. A 3D magnetohydrodynamic (MHD) simulation code is used to simulate the evolution of these disturbances as they propagate out to 1 AU. The model in Han, Wu and Dryer (1988) uses solar data for input at 0.08 AU (18 solar radii). The initial shock speed (ISS) is assumed to be constant from the corona to 0.08 AU. We investigate how variations of this ISS affect the arrival times of the shock at Earth. This basic parametric study, however, does not consider inhomogeneous background solar wind structures such as corotating interaction regions and their precursor stream–stream interactions, nor interplanetary manifestations of complex coronal mass ejecta such as magnetic clouds. In the latter case, only their associated shocks are considered. Because the ambient (pre-existing background) solar wind speed is known to affect the shock arrival time at 1 AU, we also simulated events with various background solar wind speeds (BSWS) to investigate this effect. The results show that the shock arrival time at Earth depends on the BSWS, the speed of solar disturbances, their size, and their source location at the Sun. However, it is found that for a sufficiently large momentum input, the shock arrival time at Earth is not significantly affected by the pre-existing solar wind speed.     

THE BASTILLE DAY (14 JULY 2000) EVENT IN HISTORICAL LARGE SUN–EARTH CONNECTION EVENTS

One of the large Sun–Earth connection events in solar cycle 23 occurred between 14 and 16 July 2000. Anomalies occurring on several satellites were reported in association with this event. Statistical study of extreme events is important not only for a view of space weather but for seeking ways to predict such kinds of large events. The Bastille Day event was characterized by a large flux (24 000 p.f.u. at its maximum) of solar energetic protons and a fast average transit speed of approximately 1500 km s⁻¹ of the interplanetary disturbance. A geomagnetic Kp index of more than 9 was observed after an interval of approximately eleven years since 1989. We found that return periods of extreme space weather (e.g., large flares, solar energetic proton events, and large geomagnetic storms) satisfied the Weibull distribution.     

Propagation of the Bastille Day 2000 CME Shock in the Outer Heliosphere

The Bastille Day (14 July) 2000 CME is a fast, halo coronal mass ejection event headed earthward. The ejection reached Earth on 15 July 2000 and produced a very significant magnetic storm and widespread aurora. At 1 AU the Wind spacecraft recorded a strong forward shock with a speed jump from ∼ 600 to over 1000 km s⁻¹. About 6 months later, this CME-driven shock arrived at Voyager 2 (∼ 63 AU) on 12 January 2001 with a speed jump of ∼ 60 km s⁻¹. This provides a good opportunity to study the shock propagation in the outer heliosphere. In this study, we employ a 2.5-D MHD numerical model, which takes the interaction of solar wind protons and interstellar neutrals into account, to investigate the shock propagation in detail and compare the model predictions with the Voyager 2 observations. The Bastille Day CME shock undergoes a dramatic change in character from the inner to outer heliosphere. Its strength and propagation speed decay significantly with distance. The model results at the location of Voyager 2 are in good agreement with in-situ observations. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1014293527951     

The Bastille day Shocks and Merged Interaction Region

Near 1 AU the solar wind structure associated with the solar flare of 14 July 2000 (Bastille Day) consisted of a large high-speed stream of 15 July and five nearby small streams during a 10-day period. At the leading edge of the large high-speed stream, in less than 6 hours, the flow speed increased from 600 km s⁻¹ to 1100 km s⁻¹, the magnetic field intensity increased from 10 nT to 60 nT, and an interaction region was identified. The interaction region was bounded between the pair of a forward shock F and a reverse shock R. Additional forward shocks were also identified at the leading edge of each of the five smaller streams. This paper presents a magnetohydrodynamics (MHD) simulation using ACE plasma and magnetic field data near 1 AU as input to study the radial evolution of the Bastille Day solar wind event. The two shocks, F and R, propagated in opposite directions away from each other in the solar wind frame and interacted with neighboring shocks and streams; the spatial and temporal extent of the interaction region continued to increase with the heliocentric distance. The solar wind was restructured from a series of streams at 1 AU to a huge merged interaction region (MIR) extending over a period of 12 days at 5.5 AU. Throughout the interior of the MIR bounded by the shock pair F and R the magnetic field intensity was a few times stronger than that outside the MIR. The simulation shows how merging of shocks, collision of shocks, and formation of new shocks contributed to the evolution process.     

Evolution of Coronal Mass Ejections in the Inner Heliosphere: A Study Using White-Light and Scintillation Images

Knowledge of the radial evolution of the coronal mass ejection (CME) is important for the understanding of its arrival at the near-Earth space and of its interaction with the disturbed/ambient solar wind in the course of its travel to 1 AU and further. In this paper, the radial evolution of 30 large CMEs (angular width > 150^∘, i.e., halo and partial halo CMEs) has been investigated between the Sun and the Earth using (i) the white-light images of the near-Sun region from the Large Angle Spectroscopic Coronagraph (LASCO) onboard SOHO mission and (ii) the interplanetary scintillation (IPS) images of the inner heliosphere obtained from the Ooty Radio Telescope (ORT). In the LASCO field of view at heliocentric distances R≤30 solar radii (R_⊙), these CMEs cover an order of magnitude range of initial speeds, V_CME≈260–2600 km s⁻¹. Following results have been obtained from the speed evolution of these CMEs in the Sun–Earth distance range: (1) the speed profile of the CME shows dependence on its initial speed; (2) the propagation of the CME goes through continuous changes, which depend on the interaction of the CME with the surrounding solar wind encountered on the way; (3) the radial-speed profiles obtained by combining the LASCO and IPS images yield the factual view of the propagation of CMEs in the inner heliosphere and transit times and speeds at 1 AU computed from these profiles are in good agreement with the actual measurements; (4) the mean travel time curve for different initial speeds and the shape of the radial-speed profiles suggest that up to a distance of ∼80 R_⊙, the internal energy of the CME (or the expansion of the CME) dominates and however, at larger distances, the CME's interaction with the solar wind controls the propagation; (5) most of the CMEs tend to attain the speed of the ambient flow at 1 AU or further out of the Earth's orbit. The results of this study are useful to quantify the drag force imposed on a CME by the interaction with the ambient solar wind and it is essential in modeling the CME propagation. This study also has a great importance in understanding the prediction of CME-associated space weather at the near-Earth environment.     

Prediction Test for the Two Extremely Strong Solar Storms in October 2003

In late October and early November 2003, a series of space weather hazard events erupted in solar-terrestrial space. Aiming at two intense storm (shock) events on 28 and 29 October, this paper presents a Two-Step method, which combines synoptic analysis of space weather–`observing’ and quantitative prediction – ‘palpating’, and uses it to test predictions. In the first step, ‘observing’, on the basis of observations of the source surface magnetic field, interplanetary scintillation (IPS) and ACE spacecraft, we find that the propagation of the shock waves is asymmetric and northward relative to the normal direction of their solar sources due to the large-scale configuration of the coronal magnetic fields, and the Earth is located near the direction of the fastest speed and greatest energy of the shocks. Being two fast ejection shock events, the fast explosion of extremely high temperature and strong magnetic field, and background solar wind velocity as high as 600 and 1000 km s⁻¹, are also helpful to their rapid propagation. According to the synoptic analysis, the shock travel times can be estimated as 21 and 20 h, which are close to the observational results of 19.97 and 19.63 h, respectively. In the second step, ‘palpating’, we adopt a new membership function of the fast shock events for the ISF method. The predicted results here show that for the onset time of the geomagnetic disturbance, the relative errors between the observational and the predicted results are 1.8 and 6.7%, which are consistent with the estimated results of the first step; and for the magnetic disturbance magnitude, the relative errors between the observational and the predicted results are 4.1 and 3.1%, respectively. Furthermore, the comparison among the predicted results of our Two-Step method with those of five other prevailing methods shows that the Two-Step method is advantageous in predicting such strong shock event. It can predict not only shock arrival time, but also the magnitude of magnetic disturbance. The results of the present paper tell us that understanding the physical features of shock propagation thoroughly is of great importance in improving the prediction efficiency.     

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.     

Statistical Models for the Solar Flare Interval Distribution in Individual Active Regions

This article discusses statistical models for the solar flare interval distribution in individual active regions. We analyzed solar flare data in 55 active regions that are listed in the Geosynchronous Operational Environmental Satellite (GOES) soft X-ray flare catalog for the years from 1981 to 2005. We discuss some problems with a conventional procedure to derive probability density functions from any data set and propose a new procedure, which uses the maximum likelihood method and Akaike Information Criterion (AIC) to objectively compare some competing probability density functions. Previous studies of the solar flare interval distribution in individual active regions only dealt with constant or time-dependent Poisson process models, and no other models were discussed. We examine three models – exponential, lognormal, and inverse Gaussian – as competing models for probability density functions in this study. We found that lognormal and inverse Gaussian models are more likely models than the exponential model for the solar flare interval distribution in individual active regions. The possible solar flare mechanisms for the distribution models are briefly mentioned. We also briefly investigated the time dependence of probability density functions of the solar flare interval distribution and found that some active regions show time dependence for lognormal and inverse Gaussian distribution functions. The results suggest that solar flares do not occur randomly in time; rather, solar flare intervals appear to be regulated by solar flare mechanisms. Determining a solar flare interval distribution is an essential step in probabilistic solar flare forecasting methods in space weather research. We briefly mention a probabilistic solar flare forecasting method as an application of a solar flare interval distribution analysis. The application of our distribution analysis to a probabilistic solar flare forecasting method is one of the main objectives of this study.     

The role of solar-wind velocity-waves in comet outburst activity

Comet outburst activity and the structure of solar wind streams were compared on the basis of Pioneer 10, 11, Vela 3 and IMP 7, 8 measurements at the heliocentric distance r ≈ 1–6 AU. It is shown that the solar wind velocity waves which are evolving into corotating shock waves beyond the Earth orbit may be responsible for comet outburst activity. The correlation between variations of comet outburst activity with heliocentric distance and the behavior of the solar wind velocity waves is established. The closeness of the characteristic times for the velocity waves and comet outburst activity (7–8 days at r = 1 AU) as well as the simultaneous growth of both the characteristic times with r are noted. The observed distribution of the comet outburst activity parameters during the 11-year cycle is also in good agreement with the phase distributions during the 11-year cycle of variations of the coronal hole areas and the rate of change of the sunspot area δS _p.     

Solar Cycle Variation of the Interplanetary Forward Shock Drivers Observed at 1 AU

Forecasting space weather more accurately from solar observations requires an understanding of the variations in physical properties of interplanetary (IP) shocks as solar activity changes. We examined the characteristics (occurrence rate, physical parameters, and types of shock driver) of IP shocks. During the period of 1995 – 2001, a total of 249 forward IP shocks were observed. In calculating the shock parameters, we used the solar wind data from Wind at the solar minimum period (1995 – 1997) and from ACE since 1998 including the solar maximum period (1999 – 2001). Most of IP shocks (68%) are concentrated in the solar maximum period. The values of physical quantities of IP shocks, such as the shock speed, the sonic Mach number, and the ratio of plasma density compression, are larger at solar maximum than at solar minimum. However, the ratio of IMF compression is larger at solar minimum. The IP shock drivers are classified into four groups: magnetic clouds (MCs), ejecta, high speed streams (HSSs), and unidentified drivers. The MC is the most dominant and strong shock driver and 150 out of total 249 IP shocks are driven by MCs. The MC is a principal and very effective shock driver not only at solar maximum but also at solar minimum, in contrast to results from previous studies, where the HSS is considered as the dominant IP shock driver.     

Evolution and flare productivity of active regions in October–November 2003

In the current solar cycle, the concentration of flare activity peaked during the period from October 19 to November 5, 2003, 3.5 years after the maximum point of the current solar-activity cycle. During this time, 56 high-(16) and medium-class flares occurred on the Sun, including 11 X flares. The flux of every such flare exceeded by a factor of 30 to 600 the 1–8 Å soft X-ray background flux of the entire Sun during flare-free periods. The disturbances caused by these flares produced six major S2-to S4-level proton events and four G1-to G5-class magnetic storms in the Earth’s space environment. Among the solar events observed were the most powerful X-ray flare of the current solar cycle, the eighth solar proton event in terms of particle flux during the entire history of observations, and the seventh magnetic storm in terms of Ap index. The most powerful flare resulted in the fastest coronal mass ejection during the current solar cycle with the solar plasma moving through interplanetary space at a velocity of 10⁶ km/s, which is about four times higher than the average velocity. Severe magnetic storms during the period from September 29 through October 3 posed a lot of problems for research and technological satellites (Advanced Composition Explorer (ACE), Aqua, Chandra, Chips, Cluster, Geostationary Operational Environmental Satellites (GOES) 9, 10, and 12, etc.) and spacecraft in interplanetary space (Mars Explorer Rover and Microwave Anisotropy Probe). The Advanced Earth Observing Satellite 2 (ADEOS 2), a Japanese satellite for monitoring the Earth’s environment, was disabled at the time of the arrival of the powerful interplanetary shock from the superflare of October 28, 2003. During this period, the ISS astronauts were forced to escape into the aft part of the station five times, which ensured the strongest protection against radiation. This paper is dedicated to the study of the solar situation and individual flare events.     

ACE Observations of the Bastille Day 2000 Interplanetary Disturbances

We present ACE observations for the six-day period encompassing the Bastille Day 2000 solar activity. A high level of transient activity at 1 AU, including ICME-driven shocks, magnetic clouds, shock-accelerated energetic particle populations, and solar energetic ions and electrons, are described. We present thermal ion composition signatures for ICMEs and magnetic clouds from which we derive electron temperatures at the source of the disturbances and we describe additional enhancements in some ion species that are clearly related to the transient source. We describe shock acceleration of 0.3–2.0 MeV nucl⁻¹ protons and minor ions and the relative inability of some of the shocks to accelerate significant energetic ion populations near 1 AU. We report the characteristics of < 20 MeV nucl⁻¹ solar energetic ions and < 0.32 MeV electrons and attempt to relate the release of energetic electrons to particular source regions.     

Post-Impulsive-Phase Acceleration in a Wide Range of Solar Longitudes

We have analyzed five solar energetic particle (SEP) events observed aboard the SOHO spacecraft during 1996–1997. All events were associated with impulsive soft X-ray flares, Type II radio bursts and coronal mass ejections (CMEs). Most attention is concentrated on the SEP acceleration during the first 100 minutes after the flare impulsive phase, post-impulsive-phase acceleration, being observed in eruptions centered at different solar longitudes. As a representative pattern of a (nearly) well-connected event, we consider the west flare and CME of 9 July 1996 (S10 W30). Similarities and dissimilarities of the post-impulsive-phase acceleration at large heliocentric-angle distance from the eruption center are illustrated with the 24 September 1997 event (S31 E19). We conclude that the proton acceleration at intermediate scales, between flare acceleration and interplanetary CME-driven shock acceleration, significantly contributes to the production of ≳10 MeV protons. This post-impulsive-phase acceleration seems to be caused by the CME lift-off.     

Responses of the polar ionosphere to the Bastille Day solar event

Simultaneous observations at Zhongshan Station, Antarctica, are presented for the interval of 13–17 July 2000 to show responses of the polar ionosphere to the Bastille Day (14 July 2000) solar event. The polar ionosphere was highly disturbed, as shown by frequently large deviations of the geomagnetic H-component, large riometer absorption events and strong ULF waves. Associated with the huge solar proton event produced by the X5/3B flare, a polar cap absorption (PCA) was observed. It began at ∼ 10:40 UT on 14 July and ended at ∼ 19:40 UT on 17 July. Superposed on it, there was a large absorption event with a peak of 26 dB, starting at ∼ 03:00 UT and ending at ∼ 11:10 UT on 15 July. This kind of absorption was probably produced by an intense `cloud of energetic electrons' during an auroral substorm. The ULF waves were very intense during the main phase and the recovery phase of the severe magnetic storm on 15 and 16 July. The ionospheric absorption was so strong that the digisonde signal was blacked out most of the time. The ionosphere returned to normal in the afternoon on 17 July.     

Lasco and eit Observations of the Bastille day 2000 Solar Storm

The period of 10–14 July 2000 saw a large number of energetic solar events ending with a very energetic flare that was associated with a large solar energetic particle event and a fast halo coronal mass ejection (CME) that produced the largest geomagnetic disturbance since 1989. This paper tries to summarize the complex coronal activity observed during this period, in order to establish a background for a number of papers in this topical issue. The GOES X-ray data are presented. Data animations of observations from EIT and LASCO C2 and C3 are presented on the accompanying CD-ROM. The observations around the time of the three X-class flares are considered. EIT observations of the Bastille Day flare show coronal brightening followed by dimming. LASCO had good data coverage for all three events. For one of the flares, no coronal response was seen. The other two flares are associated with halo CMEs. The timing suggests that the start of the flares and CMEs are simultaneous to approximately 30 min. Analysis of the LASCO and EIT images following the Bastille Day flare show the arrival of energetic particles at SOHO at approximately 10:41 UT on 14 July. Individual features of these CMEs have been tracked and the height–time plots used to estimate the dynamics of the CMEs. The initial speed and deceleration of the halo CMEs estimated from the fitting of height–time plots are compared with the in-situ observations at L1. The three flares are identified as the solar sources of three shocks observed at 1 AU. Finally, it is stressed that global heliospheric effects during periods of exceptional activity should consider a cumulative scenario rather than events in isolation.     

Correction of SOHO CELIAS/SEM EUV measurements saturated by extreme solar flare events

The solar irradiance in the Extreme Ultraviolet (EUV) spectral bands has been observed with a 15 s cadence by the SOHO Solar EUV Monitor (SEM) since 1995. During remarkably intense solar flares the SEM EUV measurements are saturated in the central (zero) order channel (0.1–50.0 nm) by the flare soft X‐ray and EUV flux. The first order EUV channel (26–34 nm) is not saturated by the flare flux because of its limited bandwidth, but it is sensitive to the arrival of Solar Energetic Particles (SEP). While both channels detect nearly equal SEP fluxes, their contributions to the count rate is sensibly negligible in the zero order channel but must be accounted for and removed from the first channel count rate. SEP contribution to the measured SEM signals usually follows the EUV peak for the gradual solar flare events. Correcting the extreme solar flare SEMEUV measurements may reveal currently unclear relations between the flare magnitude, dynamics observed in different EUV spectral bands, and the measured Earth atmosphere response. A simple and effective correction technique based on analysis of SEM count‐rate profiles, GOES X‐ray, and GOES proton data has been developed and used for correcting EUV measurements for the five extreme solar flare events of July 14, 2000, October 28, November 2, November 4, 2003, and January 20, 2005. Although none of the 2000 and 2003 flare peaks were contaminated by the presence of SEPs, the January 20, 2005 SEPs were unusually prompt and contaminated the peak. The estimated accuracy of the correction is about ±7.5% for large X‐class events. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A study on the dynamic spectral indices for SEP events on 2000 July 14 and 2005 January 20

We have studied the dynamic proton spectra for the two solar energetic particle(SEP) events on2000 July 14(hereafter GLE59) and 2005 January 20(hereafter GLE69). The source locations of GLE59 and GLE69 are N22 W07 and N12 W58 respectively. Proton fluxes 30 Me V have been used to compute the dynamic spectral indices of the two SEP events. Our results show that spectral indices of the two SEP events increased more swiftly at early times, suggesting that the proton fluxes 30 Me V might be accelerated particularly by the concurrent flares at early times for the two SEP events. For the GLE69 with source location at N12 W58, both flare site and shock nose are well connected with the Earth at the earliest time. However, only the particles accelerated by the shock driven by eastern flank of the CME can propagate along the interplanetary magnetic field line to the Earth after the flare. For the GLE59 with source location at N22 W07, only the particles accelerated by the shock driven by western flank of the associated CME can reach the Earth after the flare. Our results also show that there was slightly more than one hour during which the proton spectra for GLE69 are softer than that for GLE59 after the flares, suggesting that the shock driven by eastern flank of the CME associated with GLE69 is weaker than the shock driven by the western flank of the CME associated with GLE59. The results support that quasi-perpendicular shock has stronger potential in accelerating particles than the quasi-parallel shock. These results also suggest that only a small part of the shock driven by western flank of the CME associated with the GLE59 is quasi-perpendicular.