Cone model-based SEP event calculations for applications to multipoint observations

The problem of modeling solar energetic particle (SEP) events is important to both space weather research and forecasting, and yet it has seen relatively little progress. Most important SEP events are associated with coronal mass ejections (CMEs) that drive coronal and interplanetary shocks. These shocks can continuously produce accelerated particles from the ambient medium to well beyond 1 AU. This paper describes an effort to model real SEP events using a Center for Integrated Space weather Modeling (CISM) MHD solar wind simulation including a cone model of CMEs to initiate the related shocks. In addition to providing observation-inspired shock geometry and characteristics, this MHD simulation describes the time-dependent observer field line connections to the shock source. As a first approximation, we assume a shock jump-parameterized source strength and spectrum, and that scatter-free transport occurs outside of the shock source, thus emphasizing the role the shock evolution plays in determining the modeled SEP event profile. Three halo CME events on May 12, 1997, November 4, 1997 and December 13, 2006 are used to test the modeling approach. While challenges arise in the identification and characterization of the shocks in the MHD model results, this approach illustrates the importance to SEP event modeling of globally simulating the underlying heliospheric event. The results also suggest the potential utility of such a model for forcasting and for interpretation of separated multipoint measurements such as those expected from the STEREO mission.

Ambient solar wind's effect on ICME transit times

Most empirical and numerical models of Interplanetary Coronal Mass Ejection (ICME) propagation use the initial CME velocity as their primary, if not only, observational input. These models generally predict a wide spread of 1 AU transit times for ICMEs with the same initial velocity. We use a 3D coupled MHD model of the corona and heliosphere to determine the ambient solar wind's effect on the propagation of ICMEs from 30 solar radii to 1 AU. We quantitatively characterize this deceleration by the velocity of the upstream ambient solar wind. The effects of varying solar wind parameters on the ICME transit time are quantified and can explain the observed spread in transit times for ICMEs of the same initial velocity. We develop an adjustment formula that can be used in conjunction with other models to reduce the spread in predicted transit times of Earth-directed ICMEs.

A kinematically distorted flux rope model for magnetic clouds

Constant-α force-free magnetic flux rope models have proven to be a valuable first step toward understanding the global context of in situ observations of magnetic clouds. However, cylindrical symmetry is necessarily assumed when using such models, and it is apparent from both observations and modeling that magnetic clouds have highly noncircular cross sections. A number of approaches have been adopted to relax the circular cross section approximation: frequently, the cross-sectional shape is allowed to take an arbitrarily chosen shape (usually elliptical), increasing the number of free parameters that are fit between data and model. While a better “fit” may be achieved in terms of reducing the mean square error between the model and observed magnetic field time series, it is not always clear that this translates to a more accurate reconstruction of the global structure of the magnetic cloud. We develop a new, noncircular cross section flux rope model that is constrained by observations of CMEs/ICMEs and knowledge of the physical processes acting on the magnetic cloud: The magnetic cloud is assumed to initially take the form of a force-free flux rope in the low corona but to be subsequently deformed by a combination of axis-centered self-expansion and heliocentric radial expansion. The resulting analytical solution is validated by fitting to artificial time series produced by numerical MHD simulations of magnetic clouds and shown to accurately reproduce the global structure.

Numerical simulation of the 12 May 1997 CME Event: The role of magnetic reconnection

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the 12 May 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance (≈1R ) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

The distribution of solar wind speeds during solar minimum: calibration for numerical solar wind modeling constraints on the source of the slow solar wind

It took the solar polar passage of Ulysses in the early 1990s to establish the global structure of the solar wind speed during solar minimum. However, it remains unclear if the solar wind is composed of two distinct populations of solar wind from different sources (e.g., closed loops which open up to produce the slow solar wind) or if the fast and slow solar wind rely on the superradial expansion of the magnetic field to account for the observed solar wind speed variation. We investigate the solar wind in the inner corona using the Wang-Sheeley-Arge (WSA) coronal model incorporating a new empirical magnetic topology–velocity relationship calibrated for use at 0.1 AU. In this study the empirical solar wind speed relationship was determined by using Helios perihelion observations, along with results from Riley et al. (2003) and Schwadron et al. (2005) as constraints. The new relationship was tested by using it to drive the ENLIL 3-D MHD solar wind model and obtain solar wind parameters at Earth (1.0 AU) and Ulysses (1.4 AU). The improvements in speed, its variability, and the occurrence of high-speed enhancements provide confidence that the new velocity relationship better determines the solar wind speed in the outer corona (0.1 AU). An analysis of this improved velocity field within the WSA model suggests the existence of two distinct mechanisms of the solar wind generation, one for fast and one for slow solar wind, implying that a combination of present theories may be necessary to explain solar wind observations.

Computer simulation study of equilibrium properties of magnetic dipolar fluids

Langevin dynamics simulations are used to investigate the equilibrium magnetization properties and structure of magnetic dipolar fluids. The influence of using different boundary conditions are systematically studied. Simulation results on the initial susceptibility and magnetization curves are compared with theoretical predictions. The effect of particle aggregation is discussed in detail by performing a cluster analysis of the microstructure.

Ground-based and satellite observations of high-latitude auroral activity in the dusk sector of the auroral oval

On 7 December 2000, during 13:30-15:30 UT the MIRACLE all-sky camera at Ny Alesund observed auroras at high-latitudes (MLAT similar to 76) simultaneously when the Cluster spacecraft were skimming the magnetopause in the same MLT sector (at similar to 16:00-18:00 MLT). The location of the auroras (near the ionospheric convection reversal boundary) and the clear correlation between their dynamics and IMF variations suggests their close relationship with R1 currents. Consequently, we can assume that the Cluster spacecraft were making observations in the magnetospheric region associated with the auroras, although exact magnetic conjugacy between the ground-based and satellite observations did not exist. The solar wind variations appeared to control both the behaviour of the auroras and the magnetopause dynamics. Auroral structures were observed at Ny Alesund especially during periods of negative IMF B-Z. In addition, the Cluster spacecraft experienced periodic (T similar to 4 - 6 min) encounters between magnetospheric and magnetosheath plasmas. These undulations of the boundary can be interpreted as a consequence of tailward propagating magnetopause surface waves. Simultaneous dusk sector ground-based observations show weak, but discernible magnetic pulsations (Pc 5) and occasionally periodic variations (T - 2 - 3 min) in the high-latitude auroras. In the dusk sector, Pc 5 activity was stronger and had characteristics that were consistent with a field line resonance type of activity. When IMF BZ stayed positive for a longer period, the auroras were dimmer and the spacecraft stayed at the outer edge of the magnetopause where they observed electromagnetic pulsations with T similar to 1 min. We find these observations interesting especially from the viewpoint of previously presented studies relating poleward-moving high-latitude auroras with pulsation activity and MHD waves propagating at the magnetospheric boundary layers.

On the cause of a magnetospheric flux transfer event

We present a detailed investigation of a magnetospheric flux transfer event (FTE) seen by the Active Magnetospheric Tracer Explorer (AMPTE) UKS and IRM satellites around 1046 UT on October 28, 1984. This event has been discussed many times previously in the literature and has been cited as support for a variety of theories of FTE formation. We make use of a model developed to reproduce ion precipitations seen in the cusp ionosphere. The analysis confirms that the FTE is well explained as a brief excursion into an open low-latitude boundary layer (LLBL), as predicted by two theories of magnetospheric FTEs: namely, that they are bulges in the open LLBL due to reconnection rate enhancements or that they are indentations of the magnetopause by magnetosheath pressure increases (but in the presence of ongoing steady reconnection). The indentation of the inner edge of the open LLBL that these two models seek to explain is found to be shallow for this event. The ion model reproduces the continuous evolution of the ion distribution function between the sheath-like population at the event center and the surrounding magnetospheric populations; it also provides an explanation of the high-pressure core of the event as comprising field lines that were reconnected considerably earlier than those that are draped over it to give the event boundary layer. The magnetopause transition parameter is used to isolate a field rotation on the boundaries of the core, which is subjected to the tangential stress balance test. The test identifies this to be a convecting structure, which is neither a rotational discontinuity (RD) nor a contact discontinuity, but could possibly be a slow shock. In addition, evidence for ion reflection off a weak RD on the magnetospheric side of this structure is found. The event structure is consistent in many ways with features predicted for the open LLBL by analytic MHD theories and by MHD and hybrid simulations. The de Hoffman-Teller velocity of the structure is significantly different from that of the magnetosheath flow, indicating that it is not an indentation caused by a high-pressure pulse in the sheath but is consistent with the motion of newly opened field lines (different from the sheath flow because of the magnetic tension force) deduced from the best fit to the ion data. However, we cannot here rule out the possibility that the sheath flow pattern has changed in the long interval between the two satellites observing the FTE and subsequently emerging into the magnetosheath; thus this test is not conclusive in this particular case. Analysis of the fitted elapsed time since reconnection shows that the core of the event was reconnected in one pulse and the event boundary layer was reconnected in a subsequent pulse. Between these two pulses is a period of very low (but nonzero) reconnection rate, which lasts about 14 mins. Thus the analysis supports, but does not definitively verify, the concept that the FTE is a partial passage into an open LLBL caused by a traveling bulge in that layer produced by a pulse in reconnection rate.

Inferring the structure of the solar corona and inner heliosphere during the Maunder minimum using global thermodynamic magnetohydrodynamic simulations

Observations of the Sun’s corona during the space era have led to a picture of relatively constant, but cyclically varying solar output and structure. Longer-term, more indirect measurements, such as from 10Be, coupled by other albeit less reliable contemporaneous reports, however, suggest periods of significant departure from this standard. The Maunder Minimum was one such epoch where: (1) sunspots effectively disappeared for long intervals during a 70 yr period; (2) eclipse observations suggested the distinct lack of a visible K-corona but possible appearance of the F-corona; (3) reports of aurora were notably reduced; and (4) cosmic ray intensities at Earth were inferred to be substantially higher. Using a global thermodynamic MHD model, we have constructed a range of possible coronal configurations for the Maunder Minimum period and compared their predictions with these limited observational constraints. We conclude that the most likely state of the corona during—at least—the later portion of the Maunder Minimum was not merely that of the 2008/2009 solar minimum, as has been suggested recently, but rather a state devoid of any large-scale structure, driven by a photospheric field composed of only ephemeral regions, and likely substantially reduced in strength. Moreover, we suggest that the Sun evolved from a 2008/2009-like configuration at the start of the Maunder Minimum toward an ephemeral-only configuration by the end of it, supporting a prediction that we may be on the cusp of a new grand solar minimum.