With recent advances in the field of machine learning, the use of deep neural networks for time series forecasting has become more prevalent. The quasi-periodic nature of the solar cycle makes it a good candidate for applying time series forecasting methods. We employ a combination of WaveNet and LSTM neural networks to forecast the sunspot number using the years 1749 to 2019 and total sunspot area using the years 1874 to 2019 time series data for the upcoming Solar Cycle 25. Three other models involving the use of LSTMs and 1D ConvNets are also compared with our best model. Our analysis shows that the WaveNet and LSTM model is able to better capture the overall trend and learn the inherent long and short term dependencies in time series data. Using this method we forecast 11 years of monthly averaged data for Solar Cycle 25. Our forecasts show that the upcoming Solar Cycle 25 will have a maximum sunspot number around 106 $\pm$ 19.75 and maximum total sunspot area around 1771 $\pm$ 381.17. This indicates that the cycle would be slightly weaker than Solar Cycle 24.

Magnetic field extrapolation is an important tool to study the three-dimensional (3D) solar coronal magnetic field which is difficult to directly measure. Various analytic models and numerical codes exist but their results often drastically differ. Thus a critical comparison of the modeled magnetic field lines with the observed coronal loops is strongly required to establish the credibility of the model. Here we compare two different non-potential extrapolation codes, a non-linear force-free field code (CESE-MHD-NLFFF) and a non-force-free field (NFFF) code in modeling a solar active region (AR) that has a sigmoidal configuration just before a major flare erupted from the region. A 2D coronal-loop tracing and fitting method is employed to study the 3D misalignment angles between the extrapolated magnetic field lines and the EUV loops as imaged by SDO/AIA. It is found that the CESE-MHD-NLFFF code with preprocessed magnetogram performs the best, outputting a field which matches the coronal loops in the AR core imaged in AIA 94 {\AA} with a misalignment angle of ~10 degree. This suggests that the CESE-MHD-NLFFF code, even without using the information of coronal loops in constraining the magnetic field, performs as good as some coronal-loop forward-fitting models. For the loops as imaged by AIA 171 {\AA} in the outskirts of the AR, all the codes including the potential-field give comparable results of mean misalignment angle (~30 degree). Thus further improvement of the codes is needed for a better reconstruction of the long loops enveloping the core region.

We address points recently discussed in Georgoulis (2011) in reference to Leka et al. (2009b). Most importantly, we find that the results of Georgoulis (2011) support a conclusion of Leka et al. (2009b): that limited spatial resolution and the presence of unresolved magnetic structures can challenge ambiguity- resolution algorithms. Moreover, the findings of both Metcalf et al. (2006) and Leka et al. (2009b) are confirmed in Georgoulis (2011): a method's performance can be diminished when the observed field fails to conform to that method's assumptions. The implication of boundaries in models of solar magnetic structures is discussed; we confirm that the distribution of the field components in the model used in Leka et al. (2009b) is closer to what is observed on the Sun than what is proposed in Georgoulis (2011). It is also shown that method does matter with regards to simulating limited spatial resolution and avoiding an inadvertent introduction of bias. Finally, the assignment of categories to data- analysis algorithms is revisited; we argue that assignments are only useful and elucidating when used appropriately.