Abstract In recent years, ensemble modeling has been widely employed in space weather to estimate uncertainties in forecasts. We here focus on the ensemble modeling of Coronal Mass Ejections (CME) arrival times and arrival velocities using a drag‐based model, which is well‐suited for this purpose due to its simplicity and low computational cost. Although ensemble techniques have previously been applied to the drag‐based model, it is still not clear how to best determine distributions for its input parameters, namely the drag parameter and the solar wind speed. The aim of this work is to evaluate statistical distributions for these model parameters starting from a list of past CME‐ICME events. We employ LASCO coronagraph observations to measure initial CME position and speed, and in situ data to associate them with an arrival date and arrival speed. For each event we ran a statistical procedure to invert the model equations, producing parameters distributions as output. Our results indicate that the distributions employed in previous works were appropriately selected, even though they were based on restricted samples and heuristic considerations. On the other hand, possible refinements to the current method are also identified, such as the dependence of the drag parameter distribution on the CME being accelerated or decelerated by the solar wind, which deserve further investigation.
We present a Python tool to generate a standard dataset from solar images that allows for user-defined selection criteria and a range of pre-processing steps. Our Python tool works with all image products from both the Solar and Heliospheric Observatory (SoHO) and Solar Dynamics Observatory (SDO) missions. We discuss a dataset produced from the SoHO mission's multi-spectral images which is free of missing or corrupt data as well as planetary transits in coronagraph images, and is temporally synced making it ready for input to a machine learning system. Machine-learning-ready images are a valuable resource for the community because they can be used, for example, for forecasting space weather parameters. We illustrate the use of this data with a 3-5 day-ahead forecast of the north-south component of the interplanetary magnetic field (IMF) observed at Lagrange point one (L1). For this use case, we apply a deep convolutional neural network (CNN) to a subset of the full SoHO dataset and compare with baseline results from a Gaussian Naive Bayes classifier.
Propagating transverse waves are thought to be a key transporter of Poynting flux throughout the Sun’s atmosphere. Recent studies have shown that these transverse motions, interpreted as the magnetohydrodynamic kink mode, are prevalent throughout the corona. The associated energy estimates suggest the waves carry enough energy to meet the demands of coronal radiative losses in the quiescent Sun. However, it is still unclear how the waves deposit their energy into the coronal plasma. We present the results from a large-scale study of propagating kink waves in the quiescent corona using data from the Coronal Multi-channel Polarimeter (CoMP). The analysis reveals that the kink waves appear to be weakly damped, which would imply low rates of energy transfer from the large-scale transverse motions to smaller scales via either uniturbulence or resonant absorption. This raises questions about how the observed kink modes would deposit their energy into the coronal plasma. Moreover, these observations, combined with the results of Monte Carlo simulations, lead us to infer that the solar corona displays a spectrum of density ratios, with a smaller density ratio (relative to the ambient corona) in quiescent coronal loops and a higher density ratio in active-region coronal loops.
The Coronal Multi-channel Polarimeter (CoMP) has opened up exciting opportunities to probe transverse MHD waves in the Sun’s corona. The archive of CoMP data is utilized to generate a catalog of quiescent coronal loops that can be used for studying propagating kink waves. The catalog contains 120 loops observed between 2012 and 2014. This catalog is further used to undertake a statistical study of propagating kink waves in the quiet regions of the solar corona, investigating phase speeds, loop lengths, footpoint power ratio (a measure of wave power entering the corona through each footpoint of a loop) and equilibrium parameter (which provides a measure of the change in wave amplitude) values. The statistical study enables us to establish the presence of a relationship between the rate of damping and the length of the coronal loop, with longer coronal loops displaying weaker wave damping. We suggest the reason for this behavior is related to a decreasing average density contrast between the loop and ambient plasma as loop length increases. The catalog presented here will provide the community with the foundation for the further study of propagating kink waves in the quiet solar corona.
Alfvenic waves have gained renewed interest since the existence of ubiquitous propagating kink waves were discovered in the corona. It has long been suggested that Alfvenic waves play an important role in coronal heating and the acceleration of the solar wind. To this effect, it is imperative to understand the mechanisms that enable their energy to be transferred to the plasma. Mode conversion via resonant absorption is believed to be one of the main mechanisms for kink wave damping and it is considered to play a key role in the process of energy transfer. This study examines the damping of propagating kink waves in quiescent coronal loops using the Coronal Multi-channel Polarimeter. A coherence-based method is used to track the Doppler velocity signal of the waves, which enables us to investigate the spatial evolution of velocity perturbations. The power ratio of outward to inward propagating waves is used to estimate the associated damping lengths and quality factors. To enable accurate estimates of these quantities, we provide the first derivation of a likelihood function suitable for fitting models to the ratio of two power spectra obtained from discrete Fourier transforms. Maximum likelihood estimation is used to fit an exponential damping model to the observed variation in power ratio as a function of frequency. We confirm earlier indications that propagating kink waves are undergoing frequency-dependent damping. Additionally, we find that the rate of damping decreases, or equivalently the damping length increases, for. longer coronal loops that reach higher in the corona.