We derive the plage area evolution over the last 12 solar cycles employing data from all Ca II K archives available publicly in digital form known to us, including several as yet unexplored Ca II K archives. We analyse more than 290,000 full-disc Ca II K observations from 43 datasets spanning the period 1892-2019. All images were consistently processed with an automatic procedure that performs the photometric calibration (if needed) and the limb-darkening compensation. The processing also accounts for artefacts plaguing many of the images, including some very specific artefacts such as bright arcs found in Kyoto and Yerkes data. We have produced a plage area time-series from each analysed dataset. We found that the differences between the plage areas derived from individual archives are mainly due to the differences in the central wavelength and the bandpass used to acquire the data at the various sites. We have empirically cross-calibrated and combined the results obtained from each dataset to produce a composite series of plage areas. "Backbone" series are used to bridge all the series together. We have also shown that the selection of the backbone series has little effect on the final plage area composite. We have quantified the uncertainty of determining the plage areas with our processing due to shifts in the central wavelength and found it to be less than 0.01 in fraction of the solar disc for the average conditions found on historical data. We also found the variable seeing conditions during the observations to slightly increase the plage areas during activity maxima. We provide the so far most complete time series of plage areas based on corrected and calibrated historical and modern Ca II K images. Consistent plage areas are now available on 88% of all days from 1892 onwards and on 98% from 1907 onwards.

Sunspots are the most notable structure on the solar surface with strong magnetic fields. The field is generally strongest in a dark area (umbra), but sometimes stronger fields are found in non-dark regions such as a penumbra and a light bridge. The formation mechanism of such strong fields outside umbrae is still puzzling. Here we report clear evidence of the magnetic field of 6,250 G, which is the strongest field among Stokes I profiles with clear Zeeman splitting ever observed on the Sun. The field was almost parallel to the solar surface and located in a bright region sandwiched by two opposite-polarity umbrae. Using a time series of spectral datasets, we discussed the formation process of the super-strong field and suggested that this strong field region was generated as a result of compression of one umbra pushed by the horizontal flow from the other umbra, like the subduction of the Earth's crust in plate tectonics.

We have carried out a statistical study of the average orientation of the magnetic field in solar filaments with respect to their axes for more than 400 samples, based on data taken with daily full-Sun, full-Stokes spectropolarimetric observations using the He I 1083.0 nm line. The major part of the samples are the filaments in the quiet areas, but those in the active areas are included as well. The average orientation of the magnetic field in filaments shows a systematic property depending on the hemisphere; the direction of the magnetic field in filaments in the northern (southern) hemisphere mostly deviates clockwise (counterclockwise) from their axes, which run along the magnetic polarity inversion line. The deviation angles of the magnetic field from the axes, are concentrated between 10--30 deg. This hemispheric pattern is consistent with that revealed for chirality of filament barbs, filament channels and for other solar features found to possess chirality. For some filaments it was confirmed that their magnetic field direction is locally parallel to their structure seen in Halpha images. Our results for the first time confirmed this hemispheric pattern with the direct observation of the magnetic field in filaments. Interestingly, the filaments which show the opposite magnetic field deviation to the hemispheric pattern, are in many cases found above the polarity inversion line whose ambient photospheric magnetic field has the polarity alignment being opposite to that of active regions following the Hale-Nicholson law.

Current helicity and twist of solar magnetic fields are important quantities to characterize the dynamo mechanism working in the convection zone of the Sun. We have carried out a statistical study on the current helicity of solar active regions observed with the Spectro-Polarimeter (SP) of Hinode Solar Optical Telescope (SOT). We used SOT-SP data of 558 vector magnetograms of a total of 80 active regions obtained from 2006 to 2012. We have applied spatial smoothing and division of data points into weak and strong field ranges to compare the contributions from different scales and field strengths. We found that the current helicity follows the so-called hemispheric sign rule when the weak magnetic fields (absolute field strength $< 300$ gauss) are considered and no smoothing is applied. On the other hand, the pattern of current helicity fluctuates and violates the hemispheric sign rule when stronger magnetic fields are considered and the smoothing of 2.0 arcsec (mimicking ground-based observations) is applied. Furthermore, we found a tendency that the weak and inclined fields better conform to and the strong and vertical fields tend to violate the hemispheric sign rule. These different properties of helicity through the strong and weak magnetic field components give important clues to understanding the solar dynamo as well as the mechanism of formation and evolution of solar active regions.

The current helicity in solar active regions derived from vector magnetograph observations for more than 20 years indicates the so-called hemispheric sign rule; the helicity is predominantly negative in the northern hemisphere and positive in the southern hemisphere. In this paper we revisit this property and compare the statistical distribution of current helicity with Gaussian distribution using the method of normal probability paper. The data sample comprises 6630 independent magnetograms obtained at Huairou Solar Observing Station, China, over 1988-2005 which correspond to 983 solar active regions. We found the following. (1) For the most of cases in time-hemisphere domains the distribution of helicity is close to Gaussian. (2) At some domains (some years and hemispheres) we can clearly observe significant departure of the distribution from a single Gaussian, in the form of two- or multi-component distribution. (3) For the most non-single-Gaussian parts of the dataset we see co-existence of two or more components, one of which (often predominant) has a mean value very close to zero, which does not contribute much to the hemispheric sign rule. The other component has relatively large value of helicity that often determines agreement or disagreement with the hemispheric sign rule in accord with the global structure of helicity reported by Zhang et al. (2010).

The structure and dynamics of the solar corona is dominated by the magnetic field. In most areas in the corona magnetic forces are so dominant that all non-magnetic forces like plasma pressure gradient and gravity can be neglected in the lowest order. This model assumption is called the force-free field assumption, as the Lorentz force vanishes. This can be obtained by either vanishing electric currents (leading to potential fields) or the currents are co-aligned with the magnetic field lines. First we discuss a mathematically simpler approach that the magnetic field and currents are proportional with one global constant, the so-called linear force-free field approximation. In the generic case, however, the relation between magnetic fields and electric currents is nonlinear and analytic solutions have been only found for special cases, like 1D or 2D configurations. For constructing realistic nonlinear force-free coronal magnetic field models in 3D, sophisticated numerical computations are required and boundary conditions must be obtained from measurements of the magnetic field vector in the solar photosphere. This approach is currently of large interests, as accurate measurements of the photospheric field become available from ground-based (for example SOLIS) and space-born (for example Hinode and SDO) instruments. If we can obtain accurate force-free coronal magnetic field models we can calculate the free magnetic energy in the corona, a quantity which is important for the prediction of flares and coronal mass ejections. Knowledge of the 3D structure of magnetic field lines also help us to interpret other coronal observations, e.g., EUV-images of the radiating coronal plasma.

Activity studies of solar-type stars, especially with reference to the status of our current Sun among them, have exposed the importance of (1) homogeneously selecting the sample stars and (2) reliably evaluating their activities down to a considerably low level. Motivated by these requirements, we conducted an extensive study on the activities of 118 solar-analog stars (of sufficiently similar properties to each other) by measuring the emission strength at the core of Ca II 3933.663 line (K line) on the high-dispersion spectrogram obtained by Subaru/HDS, where special attention was paid to correctly detecting the chromospheric emission by removing the wing-fitted photospheric profile calculated from the classical solar model atmosphere. This enabled us to detect low-level activities down to log R' ~ -5.4 (R' is the ratio of the chromospheric core emission flux to the total bolometric flux), by which we could detect subtle activity differences which were indiscernible in previous studies. Regarding the Sun, we found log R'sun = -5.33 near to the low end of the distribution, which means that it belongs to the distinctly low activity group among solar analogs. This excludes the once-suggested possibility for the high frequency of Maunder-minimum stars showing appreciably lower activities than the minimum-Sun.

An extensive beryllium abundance analysis was conducted for 118 solar analogs (along with 87 FGK standard stars) by applying the spectrum synthesis technique to the near-UV region comprising the Be II line at 3131.066 A, in an attempt to investigate whether Be suffers any depletion such as the case of Li showing a large diversity. We found that, while most of these Sun-like stars are superficially similar in terms of their A(Be) (Be abundances) around the solar value within ~ +/- 0.2dex, 4 out of 118 samples turned out strikingly Be-deficient (by more than ~2 dex) and these 4 stars belong to the group of lowest v_e sin i (projected rotation velocity). Moreover, even for the other majority showing an apparent similarity in Be, we can recognize a tendency that A(Be) gradually increases with an increase in v_e sin i. These observational facts suggest that any solar analog star (including the Sun) generally suffers some kind of Be depletion during their lives, where the rotational velocity (or the angular momentum) plays an important role in the sense that depletion tends to be enhanced by slower rotation. Hence, our findings require that the occasionally stated view "G-type dwarfs with T_eff ~< 6000 K are essentially homogeneous in Be with their original composition retained" should be revised. Also, our analysis indicates that the difference of ~0.2 dex in A(Be) between the solar photosphere and the meteorite really exists, implying that "UV missing opacity" is irrelevant at least for this Be II line.

We previously attempted to ascertain why the Li I 6708 line-strengths of Sun-like stars differ so significantly despite the superficial similarities of stellar parameters. We carried out a comprehensive analysis of 118 solar analogs and reported that a close connection exists between the Li abundance A_Li and the line-broadening width (v_r+m; mainly contributed by rotational effect), which led us to conclude that stellar rotation may be the primary control of the surface Li content. To examine our claim in more detail, we study whether the degree of stellar activity exhibits a similar correlation with the Li abundance, which is expected because of the widely believed close connection between rotation and activity. We measured the residual flux at the line center of the strong Ca II 8542 line, r_0(8542), known to be a useful index of stellar activity, for all sample stars using newly acquired spectra in this near-IR region. The projected rotational velocity (v_e sin i) was estimated by subtracting the macroturbulence contribution from v_r+m that we had already established. A remarkable (positive) correlation was found in the A_Li versus (vs.) r_0(8542) diagram as well as in both the r_0(8542) vs. v_e sin i and A_Li vs. v_e sin i diagrams, as had been expected. With the confirmation of rotation-dependent stellar activity, this clearly shows that the surface Li abundances of these solar analogs progressively decrease as the rotation rate decreases. Given this observational evidence, we conclude that the depletion of surface Li in solar-type stars, probably caused by effective envelope mixing, operates more efficiently as stellar rotation decelerates. It may be promising to attribute the low-Li tendency of planet-host G dwarfs to their different nature in the stellar angular momentum.

Previously unobservable mirror asymmetry of the solar magnetic field -- a key ingredient of the dynamo mechanism which is believed to drive the 11-year activity cycle -- has now been measured. This was achieved through systematic monitoring of solar active regions carried out for more than 20 years at observatories in Mees, Huairou, and Mitaka. In this paper we report on detailed analysis of vector magnetic field data, obtained at Huairou Solar Observing Station in China. Electric current helicity (the product of current and magnetic field component in the same direction) was estimated from the data and a latitude-time plot of solar helicity during the last two solar cycles has been produced. We find that like sunspots helicity patterns propagate equatorwards but unlike sunspot polarity helicity in each solar hemisphere does not change sign from cycle to cycle - confirming the theory. There are, however, two significant time-latitudinal domains in each cycle when the sign does briefly invert. Our findings shed new light on stellar and planetary dynamos and has yet to be included in the theory.

We investigate coronal transients associated with a GOES M6.7 class flare and a coronal mass ejection (CME) on 13 July 2004. During the rising phase of the flare, a filament eruption, loop expansion, a Moreton wave, and an ejecta were observed. An EIT wave was detected later on. The main features in the radio dynamic spectrum were a frequency-drifting continuum and two type II bursts. Our analysis shows that if the first type II burst was formed in the low corona, the burst heights and speed are close to the projected distances and speed of the Moreton wave (a chromospheric shock wave signature). The frequency-drifting radio continuum, starting above 1 GHz, was formed almost two minutes prior to any shock features becoming visible, and a fast-expanding piston (visible as the continuum) could have launched another shock wave. A possible scenario is that a flare blast overtook the earlier transient, and ignited the first type II burst. The second type II burst may have been formed by the same shock, but only if the shock was propagating at a constant speed. This interpretation also requires that the shock-producing regions were located at different parts of the propagating structure, or that the shock was passing through regions with highly different atmospheric densities. This complex event, with a multitude of radio features and transients at other wavelengths, presents evidence for both blast-wave-related and CME-related radio emissions.

Knowledge regarding the coronal magnetic field is important for the understanding of many phenomena, like flares and coronal mass ejections. Because of the low plasma beta in the solar corona the coronal magnetic field is often assumed to be force-free and we use photospheric vector magnetograph data to extrapolate the magnetic field into the corona with the help of a non-linear force-free optimization code. Unfortunately the measurements of the photospheric magnetic field contain inconsistencies and noise. In particular the transversal components (say Bx and By) of current vector magnetographs have their uncertainties. Furthermore the magnetic field in the photosphere is not necessary force-free and often not consistent with the assumption of a force-free field above. We develop a preprocessing procedure to drive the observed non force-free data towards suitable boundary conditions for a force-free extrapolation. As a result we get a data set which is as close as possible to the measured data and consistent with the force-free assumption.