skip to main content


Title: Lagrangian Stochastic Model for the Motions of Magnetic Footpoints on the Solar Wind Source Surface and the Path Lengths of Boundary-driven Interplanetary Magnetic Field Lines
Abstract

In this work, we extend Leighton’s diffusion model describing the turbulent mixing of magnetic footpoints on the solar wind source surface. The present Lagrangian stochastic model is based on the spherical Ornstein–Uhlenbeck process with drift that is controlled by the rotation frequency Ω of the Sun, the Lagrangian integral timescaleτL, and the root-mean-square footpoint velocityVrms. The Lagrangian velocity and the positions of magnetic footpoints on the solar wind source surface are obtained from the solutions of a set of stochastic differential equations, which are solved numerically. The spherical diffusion model of Leighton is recovered in the singular Markov limit when the Lagrangian integral timescale tends to zero while keeping the footpoint diffusivity finite. In contrast to the magnetic field lines driven by standard Brownian processes on the solar wind source surface, the interplanetary magnetic field lines are smooth differentiable functions with finite path lengths in our model. The path lengths of the boundary-driven interplanetary magnetic field lines and their probability distributions at 1 au are computed numerically, and their dependency with respect to the controlling parameters is investigated. The path-length distributions are shown to develop a significant skewness as the width of the distributions increases.

 
more » « less
Award ID(s):
2149771
NSF-PAR ID:
10493800
Author(s) / Creator(s):
;
Publisher / Repository:
ApJ
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
945
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
150
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

     We present a stochastic field line mapping model where the interplanetary magnetic field lines are described by a density distribution function satisfying a Fokker–Planck equation that is solved numerically. Due to the spiral geometry of the nominal Parker field and to the evolving nature of solar wind turbulence, the heliospheric diffusion of the magnetic field lines is both heterogeneous and anisotropic, including a radial component. The longitudinal distributions of the magnetic field lines are shown to be close to circular Gaussian distributions, although they develop a noticeable skewness. The magnetic field lines emanating from the Sun are found to differ, on average, from the spirals predicted by Parker. Although the spirals remain close to Archimedean, they are here underwound, on average. Our model predicts a spiral angle that is smaller by ∼5° than the Parker spiral angle at Earth’s orbit for the same solar wind speed ofVsw= 400 km s−1. It also predicts an angular position on the solar disk of the best magnetically connected footpoint to an observer at 1 au that is shifted westward by ∼10° with respect to the Parker’s field model. This significantly changes the angle of the most probable magnetic connection between possible sources on the Sun and observers in the inner heliosphere. The results have direct implications for the heliospheric transport of “scatter-free” electrons accelerated in the aftermath of solar eruptions.

     
    more » « less
  2. Abstract

    The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement ofψcritπ. In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲1043erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.

     
    more » « less
  3. Over three decades of in-situ observations illustrate that the Kelvin–Helmholtz (KH) instability driven by the sheared flow between the magnetosheath and magnetospheric plasma often occurs on the magnetopause of Earth and other planets under various interplanetary magnetic field (IMF) conditions. It has been well demonstrated that the KH instability plays an important role for energy, momentum, and mass transport during the solar-wind-magnetosphere coupling process. Particularly, the KH instability is an important mechanism to trigger secondary small scale (i.e., often kinetic-scale) physical processes, such as magnetic reconnection, kinetic Alfvén waves, ion-acoustic waves, and turbulence, providing the bridge for the coupling of cross scale physical processes. From the simulation perspective, to fully investigate the role of the KH instability on the cross-scale process requires a numerical modeling that can describe the physical scales from a few Earth radii to a few ion (even electron) inertial lengths in three dimensions, which is often computationally expensive. Thus, different simulation methods are required to explore physical processes on different length scales, and cross validate the physical processes which occur on the overlapping length scales. Test particle simulation provides such a bridge to connect the MHD scale to the kinetic scale. This study applies different test particle approaches and cross validates the different results against one another to investigate the behavior of different ion species (i.e., H+ and O+), which include particle distributions, mixing and heating. It shows that the ion transport rate is about 10 25  particles/s, and mixing diffusion coefficient is about 10 10  m 2  s −1 regardless of the ion species. Magnetic field lines change their topology via the magnetic reconnection process driven by the three-dimensional KH instability, connecting two flux tubes with different temperature, which eventually causes anisotropic temperature in the newly reconnected flux. 
    more » « less
  4. Abstract

    We model lower band chorus observations from the DEMETER satellite using daily and hourly autoregressive‐moving average transfer function (ARMAX) equations. ARMAX models can account for serial autocorrelation between observations that are measured close together in time and can be used to predict a response variable based on its past behavior without the need for recent data. Unstable distributions of radiation belt source electrons (tens of keV) and the substorm activity (SMEd from the SuperMAG array) that is thought to inject these electrons were both statistically significant explanatory variables in a daily ARMAX model describing chorus. Predictions from this model correlated well with observations in a hold‐out test data set (validation correlation of 0.675). Source electron flux was most influential when observations came from the same day or the day before the chorus measurement, with effects decaying rapidly over time. Substorms were more influential when they occurred on previous days, presumably due to their injecting source electrons from the plasma sheet. A daily ARMAX model with interplanetary magnetic field (IMF)|B|, IMFBz, and solar wind pressure as inputs instead of those given above was somewhat less predictive of chorus (r=0.611). An hourly ARMAX model with only solar wind and IMF inputs was even less successful, with a validation correlation of 0.502.

     
    more » « less
  5. Abstract

    We studied the spatiotemporal structure of ground magnetic pulsations on the dayside by displaying magnetic field perturbations detected by the European quasi‐Meridional Magnetometer Array (EMMA) as 2‐D images in the magneticLvalue versus time space, called EMMAgrams. We generated EMMAgrams from observations made on 15 August 2015, including a previously studied pulsation event associated with an interplanetary shock. In addition to signatures of field line resonance driven by a cavity mode oscillation, we found poleward propagating structures withL‐independent periods in the Pc2 band. The Pc2 structures are attributed to periodic magnetohydrodynamic pulses (upstream waves) originating from the ion foreshock and propagating in the magnetosphere along the path proposed by Tamao. Ringing of local field lines atL‐dependent periods (transient pulsations) is also clearly detected as dispersive poleward propagating structures not only immediately after the shock impact but also during time periods of less obvious external disturbances. A transient pulsation decays after a few wave periods, and cross‐spectral analysis of transient pulsations detected at two stations with a small latitudinal separation indicates elevation of the cross phase in a band delimited by the field line resonance frequencies at the stations. Successive excitation of transient pulsations by variations of the solar wind dynamic pressure appears to contribute significantly to formation of similar cross‐phase peaks that are widely used in magnetoseismic studies.

     
    more » « less