The sunspot area is a critical physical quantity for assessing the solar activity level; forecasts of the sunspot area are of great importance for studies of the solar activity and space weather. We developed an innovative hybrid model prediction method by integrating the complementary ensemble empirical mode decomposition (CEEMD) and extreme learning machine (ELM). The time series is first decomposed into intrinsic mode functions (IMFs) with different frequencies by CEEMD; these IMFs can be divided into three groups, a high-frequency group, a low-frequency group, and a trend group. The ELM forecasting models are established to forecast the three groups separately. The final forecast results are obtained by summing up the forecast values of each group. The proposed hybrid model is applied to the smoothed monthly mean sunspot area archived at NASA's Marshall Space Flight Center (MSFC). We find a mean absolute percentage error (MAPE) and a root mean square error (RMSE) of 1.80% and 9.75, respectively, which indicates that: (1) for the CEEMD-ELM model, the predicted sunspot area is in good agreement with the observed one; (2) the proposed model outperforms previous approaches in terms of prediction accuracy and operational efficiency.

The autoregressive method provides a univariate procedure to predict the future sunspot number (SSN) based on past record. The strength of this method lies in the possibility that from past data it yields the SSN in the future as a function of time. On the other hand, its major limitation comes from the intrinsic complexity of solar magnetic activity that may deviate from the linear stationary process assumption that is the basis of the autoregressive model. By analyzing the residual errors produced by the method, we have obtained the following conclusions: (1) the optimal duration of the past time for the forecast is found to be 8.5 years; (2) the standard error increases with prediction horizon and the errors are mostly systematic ones resulting from the incompleteness of the autoregressive model; (3) there is a tendency that the predicted value is underestimated in the activity rising phase, while it is overestimated in the declining phase; (5) the model prediction of a new Solar Cycle is fairly good when it is similar to the previous one, but is bad when the new cycle is much different from the previous one; (6) a reasonably good prediction of a new cycle can be made using the AR model 1.5 years after the start of the cycle. In addition, we predict the next cycle (Solar Cycle 25) will reach the peak in 2024 at the activity level similar to the current cycle.

High resolution, multi-wavelength images from the Dutch Open Telescope were used to study the detailed mechanisms that might be involved in the multiple layer solar atmosphere observed in high cadence multi-wavelength observations. With the exceptional data observed for active region NOAA 10789 on 2005 July 13th, we study the changing pattern of the fibril using multi-wavelength tomography of the Hα line center and blue wing, Ca II H, and the G Band. It is believed that a long fibril that is rooted in the umbra, with longer apparent periodicity, may be due to morphological changes. To determine this, we conduct phase difference and coherency analysis between points along the fibril to understand how the wave propagates.

Using sunspot number data from 270 historical stations for the period 1981-2013, we investigate their personal reduction coefficients (k) statistically. Chang & Oh (2012) perform a simulation showing that the k varies with the solar cycle. We try to verify their results using observational data. For this, a weighted mean and weighted standard deviation of monthly sunspot number are used to estimate the error from observed data. We find that the observed error (noise) is much smaller than that used in the simulation. Thus no distinct k-variation with the solar cycle is observed contrary to the simulation. In addition, the probability distribution of k is determined to be non-Gaussian with a fat-tail on the right side. This result implies that the relative sunspot number after 1981 might be overestimated since the mean value of k is less than that of the Gaussian distribution.