We investigate the rising flux tube and the formation of sunspots in an unprecedentedly deep computational domain that covers the whole convection zone with a radiative magnetohydrodynamics simulation. Previous calculations had shallow computational boxes (< 30 Mm) and convection zones at a depth of 200 Mm. By using our new numerical code R2D2, we succeed in covering the whole convection zone and reproduce the formation of the sunspot from a simple horizontal flux tube because of the turbulent thermal convection. The main findings are (1) The rising speed of the flux tube is larger than the upward convection velocity because of the low density caused by the magnetic pressure and the suppression of the mixing. (2) The rising speed of the flux tube exceeds 250 m/s at a depth of 18 Mm, while we do not see any clear evidence of the divergent flow 3 hr before the emergence at the solar surface. (3) Initially, the root of the flux tube is filled with the downflows and then the upflow fills the center of the flux tube during the formation of the sunspot. (4) The essential mechanisms for the formation of the sunspot are the coherent inflow and the turbulent transport. (5) The low-temperature region is extended to a depth of at least 40 Mm in the matured sunspot, with the high-temperature region in the center of the flux tube. Some of the findings indicate the importance of the deep computational domain for the flux emergence simulations.

The morphological asymmetry of leading and following sunspots is a well-known characteristic of the solar surface. In the context of large-scale evolution of the surface magnetic field, the asymmetry has been assumed to have only a negligible effect. Using the surface flux transport model, we show that the morphological asymmetry of leading and following sunspots has a significant impact on the evolution of the large-scale magnetic field on the solar surface. By evaluating the effect of the morphological asymmetry of each bipolar magnetic region (BMR), we observe that the introduction of the asymmetry in the BMR model significantly reduces its contribution to the polar magnetic field, especially for large and high-latitude BMRs. Strongly asymmetric BMRs can even reverse the regular polar field formation. The surface flux transport simulations based on the observed sunspot record shows that the introduction of the morphological asymmetry reduces the root-mean-square difference from the observed axial dipole strength by 30--40 percent. These results indicate that the morphological asymmetry of leading and following sunspots has a significant effect on the solar cycle prediction.

We investigate the diagnostic potential of the spectral lines at 850 nm for understanding the magnetism of the lower atmosphere. For that purpose, we use a newly developed 3D simulation of a chromospheric jet to check the sensitivity of the spectral lines to this phenomenon as well as our ability to infer the atmospheric information through spectropolarimetric inversions of noisy synthetic data. We start comparing the benefits of inverting the entire spectrum at 850 nm versus only the Ca II 8542 A spectral line. We found a better match of the input atmosphere for the former case, mainly at lower heights. However, the results at higher layers were not accurate. After several tests, we determined that we need to weight more the chromospheric lines than the photospheric ones in the computation of the goodness of the fit. The new inversion configuration allows us to obtain better fits and consequently more accurate physical parameters. Therefore, to extract the most from multi-line inversions, a proper set of weights needs to be estimated. Besides that, we conclude again that the lines at 850 nm, or a similar arrangement with Ca II 8542 A plus Zeeman sensitive photospheric lines, poses the best observing configuration for examining the thermal and magnetic properties of the lower solar atmosphere.

The reduced speed of sound technique (RSST) has been used for efficient simulation of low Mach number flows in solar and stellar convection zones. The basic RSST equations are hyperbolic, and are suitable for parallel computation by domain decomposition. The application of RSST is limited to cases where density perturbations are much smaller than the background density. In addition, non-conservative variables are required to be evolved using this method, which is not suitable in cases where discontinuities like shock waves co-exist in a single numerical domain. In this study, we suggest a new semi-conservative formulation of the RSST that can be applied to low Mach number flows with large density variations. We derive the wave speed of the original and newly suggested methods to clarify that these methods can reduce the speed of sound without affecting the entropy wave. The equations are implemented using the finite volume method. Several numerical tests are carried out to verify the suggested methods. The analysis and numerical results show that the original RSST is not applicable when mass density variations are large. In contrast, the newly suggested methods are found to be efficient in such cases. We also suggest variants of the RSST that conserve momentum in the machine precision. The newly suggested variants are formulated as semi-conservative equations, which reduce to the conservative form of the Euler equations when the speed of sound is not reduced. This property is advantageous when both high and low Mach number regions are included in the numerical domain. The newly suggested forms of RSST can be applied to a wider range of low Mach number flows.

We report the results of the Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the solar chromosphere on the southern polar limb. Coordinated observations with the Interface Region Imaging Spectrograph (IRIS) are also conducted. ALMA provided unprecedented high spatial resolution in the millimeter band ($\approx$ 2.0 arcsec) at 100 GHz frequency with a moderate cadence (20 s). The results are as follows: (1) The ALMA 100 GHz images show saw-tooth patterns on the limb, and a comparison with SDO/AIA 171\AA\ images shows a good correspondence of the limbs with each other. (2) The ALMA 100 GHz movie shows a dynamic thorn-like structure elongating from the saw-tooth patterns on the limb, with lengths reaching at least 8 arcsec, thus suggesting jet-like activity in the ALMA microwave range. These ALMA jets are in good correspondence with IRIS jet clusters. (3) A blob ejection event is observed. By comparing with the IRIS Mg II slit-jaw images, the trajectory of the blob is located along the spicular patterns.

Sunspot observations with telescopes in 18th century were carried out in Japan as well. One of these sunspot observations is recorded in an account called Sansaizusetsu narabini Kansei irai Jissoku Zusetsu (Charts of Three Worlds and Diagrams of Actual Observations since Kansei Era). We analyze manuscripts of this account to show in total 15 sunspot drawings in 1749-1750. These observations were carried out by contemporary official astronomers in Japan, with telescopes covered by zongurasus (< zonglas in Dutch, corresponding to "sunglass" in English). We count their group number of sunspots to locate them in long-term solar activity and show that their observations were situated around the solar maximum in 1749 or 1750. We also computed their locations and areas, while we have to admit the difference of variant manuscripts with one another. These observational records show the spread of sunspot observations not only in Europe but also in Japan and hence may contribute to crosscheck or possibly improve the known sunspot indices.

Aims. We report the small temporal variation of the axial dipole moment near the solar minimum and its application to the solar cycle prediction by the surface flux transport (SFT) model. Methods. We measure the axial dipole moment using the photospheric synoptic magnetogram observed by the Wilcox Solar Observatory (WSO), the ESA/NASA Solar and Heliospheric Observatory Michelson Doppler Imager (MDI), and the NASA Solar Dynamics Observatory Helioseismic and Magnetic Imager (HMI). We also use the surface flux transport model for the interpretation and prediction of the observed axial dipole moment. Results. We find that the observed axial dipole moment becomes approximately constant during the period of several years before each cycle minimum, which we call the axial dipole moment plateau. The cross-equatorial magnetic flux transport is found to be small during the period, although the significant number of sunspots are still emerging. The results indicates that the newly emerged magnetic flux does not contributes to the build up of the axial dipole moment near the end of each cycle. This is confirmed by showing that the time variation of the observed axial dipole moment agrees well with that predicted by the SFT model without introducing new emergence of magnetic flux. These results allows us to predict the axial dipole moment in Cycle 24/25 minimum using the SFT model without introducing new flux emergence. The predicted axial dipole moment of Cycle 24/25 minimum is 60--80 percent of Cycle 23/24 minimum, which suggests the amplitude of Cycle 25 even weaker than the current Cycle 24. Conclusions. The plateau of the solar axial dipole moment is an important feature for the longer prediction of the solar cycle based on the SFT model.

This paper presents a three-dimensional simulation of chromospheric jets with twisted magnetic field lines. Detailed treatments of the photospheric radiative transfer and the equation of states allow us to model realistic thermal convection near the solar surface, which excites various MHD waves and produces chromospheric jets in the simulation. A tall chromospheric jet with a maximum height of 10--11 Mm and lifetime of 8--10 min is formed above a strong magnetic field concentration. The magnetic field lines are strongly entangled in the chromosphere, which helps the chromospheric jet to be driven by the Lorentz force. The jet exhibits oscillatory motion as a natural consequence of its generation mechanism. We also find that the produced chromospheric jet forms a cluster with a diameter of several Mm with finer strands. These results imply a close relationship between the simulated jet and solar spicules.

We investigate the effect of coronal temperature on the formation process of solar chromospheric jets using two-dimensional magnetohydrodynamic simulations of the region from the upper convection zone to the lower corona. We develop a new radiative magnetohydrodynamic code for the dynamic modeling of the solar atmosphere, employing a LTE equation of state, optically thick radiative loss in the photosphere, optically thin radiative loss in the chromosphere and the corona, and thermal conduction along the magnetic field lines. Many chromospheric jets are produced in the simulations by shock waves passing through the transition region. We find that these jets are projected farther outward when the coronal temperature is lower (similar to that in coronal holes) and shorter when the coronal temperature is higher (similar to that in active regions). When the coronal temperature is high, the deceleration of the chromospheric jets is consistent with the model in which deceleration is determined by the periodic chromospheric shock waves. However, when the coronal temperature is low, the gravitational deceleration becomes more important and the chromospheric jets approach ballistic motion.