To simulate solar Coronal Mass Ejections (CMEs), predict their time of arrival and geomagnetic impact, it is important to accurately model the background solar wind conditions in which CMEs propagate. We use the Alfvén Wave Solar-atmosphere Model (AWSoM) within the the Space Weather Modeling Framework (SWMF) to simulate solar maximum conditions during two Carrington rotations and produce solar wind background conditions comparable to the observations. We describe the inner boundary conditions for AWSoM using the ADAPT global magnetic maps and validate the simulated results with EUV observations in the low corona and measured plasma parameters at L1 as well as at the position of the STEREO spacecraft. This work complements our prior AWSoM validation study for solar minimum conditions Sachdeva et al. (2019), and shows that during periods of higher magnetic activity, AWSoM can reproduce the solar plasma conditions (using properly adjusted photospheric Poynting flux) suitable for providing proper initial conditions for launching CMEs.

The magnetohydrodynamics with embedded particle-in-cell (MHD-EPIC) model has been successfully applied to global magnetospheric simulations in recent years. However, the PIC region was restricted to be one or more static boxes, which is not always sufficient to cover the whole physical structure of interest efficiently. The FLexible Exascale Kinetic Simulator (FLEKS), which is a new PIC code and allows a dynamic PIC region of any shape, is designed to break this restriction. FLEKS is usually used as the PIC component of the MHD with adaptively embedded particle-in-cell (MHD-AEPIC) model. FLEKS supports dynamically activating or deactivating cells to fit the regions of interest during a simulation. An adaptive time-stepping scheme is also introduced to improve the accuracy and efficiency of a long simulation. The particle number per cell may increase or decrease significantly and lead to load imbalance and large statistical noise in the cells with fewer particles. A particle splitting scheme and a particle merging algorithm are designed to limit the change of the particle number and hence improve the accuracy of the simulation as well as load balancing. Both particle splitting and particle merging conserve the total mass, momentum, and energy. FLEKS also contains a test-particle module to enable tracking particle trajectories due to the time-dependent electromagnetic field that is obtained from a global simulation.

We perform a geomagnetic event simulation using a newly developed magnetohydrodynamic with adaptively embedded particle-in-cell (MHD-AEPIC) model. We have developed effective criteria to identify reconnection sites in the magnetotail and cover them with the PIC model. The MHD-AEPIC simulation results are compared with Hall MHD and ideal MHD simulations to study the impacts of kinetic reconnection at multiple physical scales. At the global scale, the three models produce very similar SYM-H and SuperMag Electrojet (SME) indexes, which indicates that the global magnetic field configurations from the three models are very close to each other. We also compare the ionospheric solver results and all three models generate similar polar cap potentials and field aligned currents. At the mesoscale we compare the simulations with in situ Geotail observations in the tail. All three models produce reasonable agreement with the Geotail observations. At the kinetic scales, the MHD-AEPIC simulation can produce a crescent shape distribution of the electron velocity space at the electron diffusion region which agrees very well with MMS observations near a tail reconnection site. These electron scale kinetic features are not available in either the Hall MHD or ideal MHD models. Overall, the MHD-AEPIC model compares well with observations at all scales, it works robustly, and the computational cost is acceptable due to the adaptive adjustment of the PIC domain. It remains to be determined whether kinetic physics can play a more significant role in other types of events, including but not limited to substorms.

The largest moon in the solar system, Ganymede, is the only moon known to possess a strong intrinsic magnetic field and a corresponding magnetosphere. Using the latest version of Space Weather Modeling Framework (SWMF), we study the upstream plasma interactions and dynamics in this sub-Alfvenic system. Results from the Hall MHD and the coupled MHD with embedded Particle-in-Cell (MHD-EPIC) models are compared. We find that under steady upstream conditions, magnetopause reconnection occurs in a non-steady manner. Flux ropes of Ganymede’s radius in length form on the magnetopause at a rate about 2/minute and create spatiotemporal variations in plasma and field properties. Upon reaching proper grid resolutions, the MHD-EPIC model can resolve both electron and ion kinetics at the magnetopause and show localized non-gyrotropic behavior inside the diffusion region. The estimated global reconnection rate from the models is about 80 kV with 60% efficiency, and there is weak evidence of about 1 minute periodicity in the temporal variations due to the dynamic reconnection process.

5 For the first time, we simulate the detailed spectral line emission from a solar active 6 region (AR) with the Alfvén Wave Solar Model (AWSoM). We select an AR appearing 7 near disk center on 2018 July 13 and use an NSO/HMI synoptic magnetogram to specify 8 the magnetic field at the model’s inner boundary. To resolve smaller-scale magnetic 9 features, we apply adaptive mesh refinement to resolve the AR with a horizontal spatial 10 resolution of 0.35 • (4.5 Mm), four times higher than the background corona. We then 11 apply the SPECTRUM code informed with CHIANTI spectral emissivities to calculate 12 spectral lines forming at temperatures ranging from 0.5 to 3 MK. Comparisons are 13 made between the simulated line intensities and those observed by the Hinode/EIS 14 instrument where we find close agreement (about 20% relative error for both loop top 15 and footpoints at a temperature of about 1.5 MK) across a wide range of loop sizes and 16 temperatures. We also simulate and compare Doppler velocities and find that simulated 17 flow patterns are of comparable magnitude to what is observed. Our results demonstrate 18 the broad applicability of the low-frequency Alfvén wave balanced turbulence theory 19 for explaining the heating of coronal loops. 20

Flux Transfer Events (FTEs) are transient magnetic flux ropes typically found at the Earth’s magnetopause on the dayside. While it is known that FTEs are generated by magnetic reconnection, it remains unclear how the details of magnetic reconnection controls their properties. A recent study showed that the helicity sign of FTEs positively correlates with the east-west (By) component of the Interplanetary Magnetic Field (IMF). With data from the Cluster and Magnetospheric Multiscale missions, we performed a statistical study of 166 quasi force-free FTEs. We focus on their helicity sign and possible association with upstream solar wind conditions and local magnetic reconnection properties. Using both in situ data and magnetic shear modeling, we find that FTEs whose helicity sign corresponds to the IMF By are associated with moderate magnetic shears while those that does not correspond to the IMF By are associated with higher magnetic shears. While uncertainty in IMF propagation to the magnetopause may lead to randomness in the determination of the flux rope core field and helicity, we rather propose that for small IMF By, which corresponds to high shear and low guide field, the Hall pattern of magnetic reconnection determines the FTE core field and helicity sign. In that context we explain how the temporal sequence of multiple X-line formation and the reconnection rate are important in determining the flux rope helicity sign. This work highlights a fundamental connection between kinetic processes at work in magnetic reconnection and the macroscale structure of FTEs.

Magnetospheric sawtooth oscillations are observed during strong and steady solar wind driving conditions. The simulation results of our global MHD model with embedded kinetic physics show that when the total magnetic flux carried by constant solar wind exceeds a threshold, sawtooth-like magnetospheric oscillations are generated. Different from previous works, this result is obtained without involving time-varying ionospheric outflow in the model. The oscillation period and amplitude agree well with observations. The simulated oscillations cover a wide range of local times, although the distribution of magnitude as a function of longitude is different from observations. Our comparative simulations using ideal or Hall MHD models do not produce global time-varying features, which suggests that kinetic reconnection physics in the magnetotail is a major contributing factor to sawtooth oscillations.

The MHD with embedded PIC (MHD-EPIC) model makes it feasible to incorporate kinetic physics into a global simulation. Still, this requires a large enough box-shaped PIC domain to accommodate the movement and changes of the magnetic reconnection regions over time. This wastes computational resources on simulating regions with the expensive PIC model where MHD would be sufficient to describe the physics. We have developed a new MHD with Adaptively Embedded PIC (MHD-AEPIC) algorithm that couples the BATS-R-US MHD model with the new FLexible Exascale Kinetic Simulator (FLEKS) PIC code. In the new coupled model the PIC domains can move with the magnetic reconnection regions and adapt to them with an arbitrary shape. In this work, we will first introduce the algorithms for selecting the reconnection regions in the MHD model that need to be resolved with the kinetic PIC model. Then we will compare simulations obtained with MHD-EPIC using fixed PIC regions versus MHD-AEPIC employing adaptive PIC regions to verify that the new model generates reliable results. Finally, we will apply the MHD-AEPIC model to a global magnetic storm simulation and demonstrate the improved efficiency.