Accurate forecasting of the arrival time and subsequent geomagnetic impacts of Coronal Mass Ejections (CMEs) at Earth is an important objective for space weather forecasting agencies. Recently, the CME Arrival and Impact working team has made significant progress towards defining community-agreed metrics and validation methods to assess the current state of CME modeling capabilities. This will allow the community to quantify our current capabilities and track progress in models over time. Firstly, it is crucial that the community focuses on the collection of the necessary metadata for transparency and reproducibility of results. Concerning CME arrival and impact we have identified 6 different metadata types: 3D CME measurement, model description, model input, CME (non-)arrival observation, model output data and metrics and validation methods. Secondly, the working team has also identified a validation time period, where all events within the following two periods will be considered: 1 January 2011-31 December 2012 and January 2015-31 December 2015. Those two periods amount to a total of about 100 hit events at Earth and a large amount of misses. Considering a time period will remove any bias in selecting events and the event set will represent a sample set that will not be biased by user selection. Lastly, we have defined the basic metrics and skill scores that the CME Arrival and Impact working team will focus on.

The Wang-Sheeley-Arge (WSA)-ENLIL+Cone model is used extensively in space weather operations world-wide to model CME propagation. As such, it is important to assess its performance. We present validation results of the WSA-ENLIL+Cone model installed at the Community Coordinated Modeling Center (CCMC) and executed in real-time by the CCMC space weather team. CCMC uses the WSA-ENLIL+Cone model to predict CME arrivals at NASA missions throughout the inner heliosphere. In this work we compare model predicted CME arrival-times to in-situ ICME leading edge measurements at STEREO-A, STEREO-B, and Earth (Wind and ACE) for simulations completed between March 2010-December 2016 (over 1,800 CMEs). We report hit, miss, false alarm, and correct rejection statistics for all three locations. For all predicted CME arrivals, the hit rate is 0.5, and the false alarm rate is 0.1. For the 273 events where the CME was predicted to arrive at Earth, STEREO-A, or STEREO-B, and was actually observed (hit event), the mean absolute arrival-time prediction error was 10.4 +/- 0.9 hours, with a tendency to early prediction error of -4.0 hours. We show the dependence of the arrival-time error on CME input parameters. We also explore the impact of the multi-spacecraft observations used to initialize the model CME inputs by comparing model verification results before and after the STEREO-B communication loss (since September 2014) and STEREO-A sidelobe operations (August 2014-December 2015). There is an increase of 1.7 hours in the CME arrival time error during single, or limited two-viewpoint periods, compared to the three-spacecraft viewpoint period. This trend would apply to a future space weather mission at L5 or L4 as another coronagraph viewpoint to reduce CME arrival time errors compared to a single L1 viewpoint.

We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfv\'en wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed field regions the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes we account for a reflection of the outward propagating waves, which is accompanied by generation of weaker counter-propagating waves. The non-linear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes. The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, $R=R_{b}\sim1.1\ R_\odot$ well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfv\'en wave propagation, which adequately describe physics within the heliocentric distances range, $R_{\odot}<R<R_{b}$, in the low solar corona. By interfacing this threaded-field-lines model with the full MHD global corona model at $r=R_{b}$, we find the global solution and achieve a faster-than-real-time performance of the model on $\sim200$ cores.

We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GONG, SOHO/LASCO, and STEREO/COR. We automate this process so that minimum manual intervention is required in specifying the CME initial state. With the newly developed data-driven Eruptive Event Generator Gibson-Low (EEGGL), we present a method to derive Gibson-Low (GL) flux rope parameters through a handful of observational quantities so that the modeled CMEs can propagate with the desired CME speeds near the Sun. A test result with CMEs launched with different Carrington rotation magnetograms are shown. Our study shows a promising result for using the first-principles-based MHD global model as a forecasting tool, which is capable of predicting the CME direction of propagation, arrival time, and ICME magnetic field at 1 AU (see companion paper by Jin et al. 2016b).

On 7 January 2014 an X1.2 flare and CME with a radial speed $\approx$2500 km s$^{-1}$ was observed from near an active region close to disk center. This led many forecasters to estimate a rapid arrival at Earth ($\approx$36 hours) and predict a strong geomagnetic storm. However, only a glancing CME arrival was observed at Earth with a transit time of $\approx$49 hours and a $K_{\rm P}$ geomagnetic index of only $3-$. We study the interplanetary propagation of this CME using the ensemble Wang-Sheeley-Arge (WSA)-ENLIL+Cone model, that allows a sampling of CME parameter uncertainties. We explore a series of simulations to isolate the effects of the background solar wind solution, CME shape, tilt, location, size, and speed, and the results are compared with observed in-situ arrivals at Venus, Earth, and Mars. Our results show that a tilted ellipsoid CME shape improves the initial real-time prediction to better reflect the observed in-situ signatures and the geomagnetic storm strength. CME parameters from the Graduated Cylindrical Shell model used as input to WSA--ENLIL+Cone, along with a tilted ellipsoid cloud shape, improve the arrival-time error by 14.5, 18.7, 23.4 hours for Venus, Earth, and Mars respectively. These results highlight that CME orientation and directionality with respect to observatories play an important role in understanding the propagation of this CME, and for forecasting other glancing CME arrivals. This study also demonstrates the importance of three-dimensional CME fitting made possible by multiple viewpoint imaging.

Ensemble modeling of CMEs provides a probabilistic forecast of CME arrival time which includes an estimation of arrival time uncertainty from the spread and distribution of predictions and forecast confidence in the likelihood of CME arrival. The real-time ensemble modeling of CME propagation uses the WSA-ENLIL+Cone model installed at the CCMC and executed in real-time. The current implementation evaluates the sensitivity of WSA-ENLIL+Cone model simulations of CME propagation to initial CME parameters. We discuss the results of real-time ensemble simulations for a total of 35 CME events between January 2013 - July 2014. For the 17 events where the CME was predicted to arrive at Earth, the mean absolute arrival time prediction error was 12.3 hours, which is comparable to the errors reported in other studies. For predictions of CME arrival at Earth the correct rejection rate is 62% and the false-alarm rate is 38%. The arrival time was within the range of the ensemble arrival predictions for 8 out of 17 events. The Brier Score for CME arrival predictions is 0.15 (where 1 is a perfect forecast), indicating that on average, the predicted likelihood of CME arrival is fairly accurate. The reliability of ensemble CME arrival predictions is heavily dependent on the initial distribution of CME input parameters, particularly the median and spread. Preliminary analysis of the probabilistic forecasts suggests undervariability, indicating that these ensembles do not sample a wide enough spread in CME input parameters. Prediction errors can also arise from ambient model parameters, the accuracy of the solar wind background derived from coronal maps, or other model limitations. Finally, predictions of the Kp geomagnetic index differ from observed values by less than one for 11 out of 17 of the ensembles and Kp prediction errors computed from the mean predicted Kp show a mean absolute error of 1.3.