We study the asymmetric distribution of weak photospheric magnetic field values in the two hemispheres separately using synoptic maps from SDO/HMI, SOLIS/VSM and WSO during solar cycles 21-24. We calculate the weak-field asymmetry (shift) by fitting the distributions of weak-field values to a shifted Gaussian. Hemispheric shifts derived from the three data sets agree very well, and increase systematically when reducing the spatial resolution of the map. Shifts of the northern and southern hemisphere are typically opposite to each other. Shifts follow the evolution of the trailing flux and have a strong solar cycle variation with maxima in the early to mid-declining phase of the solar cycle. The sign of the hemispheric weak-field shift is always the same as the polarity of the polar field in the respective hemisphere and solar cycle. We also find that shifts in the south are systematically larger in absolute value than in the north.

Context. Systematic observations of magnetic field strength and polarity in sunspots began at Mount Wilson Observatory (MWO), USA in early 1917. Except for a few brief interruptions, this historical dataset continues till present. Aims. The sunspot field strength and polarity observations are critical in our project of reconstructing the solar magnetic field over the last hundred years. Here we provide a detailed description of the newly digitized dataset of drawings of sunspot magnetic field observations. Methods. The digitization of MWO drawings is based on a software package develope d by us. It includes a semi-automatic selection of solar limbs and other features of the drawing, and a manual entry of the time of observations, the measured field strength and other notes hand-written on each drawing. The data are preserved in a MySQL database. Results. We provide a brief history of the project and describe the results from digitizing this historical dataset. We also provide a summary of the final dataset, and describe its known limitations. Finally, we compare the sunspot magnetic field measurements with other instruments, and demonstrate that, if needed, the dataset could be continued using modern observations such as, for example, Vector Stokes Magnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) platform.

The line-of-sight (LOS) component of the large-scale photospheric magnetic field has been observed since the 1950s, but the daily full-disk observations of the full vector magnetic field started only in 2010 using the SOLIS Vector Stokes Magnetograph (VSM) and the SDO helioseismic and magnetic imager (HMI). Traditionally, potential field extrapolations are based on the assumption that the magnetic field in the photosphere is approximately radial. The validity of this assumption has not been tested yet. We investigate here the structure and evolution of the three components of the solar large-scale magnetic field in 2010 - 2017, covering the ascending to mid-declining phase of solar cycle 24, using SOLIS/VSM vector synoptic maps of the photospheric magnetic field. We compare the observed VSM vector magnetic field to the potential vector field derived using the VSM LOS magnetic field observations as an input. The new vector field data allow us to derive the meridional inclination and the azimuth angle of the magnetic field and to investigate their solar cycle evolution and latitudinal profile of these quantities. SOLIS/VSM vector data show that the photospheric magnetic field is in general fairly non-radial. In the meridional plane the field is inclined toward the equator, reflecting the dipolar structure of the solar magnetic field. Rotationally averaged meridional inclination does not have significant solar cycle variation. While the vector radial component $B_r$ and the potential radial component $B_r^{PFSS}$ are fairly similar, the meridional and zonal components do not agree very well. We find that SOLIS/VSM vector observations are noisy at high latitudes and suffer from the vantage point effect more than LOS observations. This is due to different noise properties in the LOS and transverse components of the magnetic field, which needs to be addressed in future studies.

We use the synoptic maps of the photospheric magnetic field observed at Wilcox Solar Observatory, Mount Wilson Observatory, Kitt Peak, SOHO/MDI, SOLIS/VSM, and SDO/HMI to study the distribution of weak photospheric magnetic field values in 1974-2018. We fit the histogram distribution of weak field values for each synoptic map of the six data-sets separately with a parametrized Gaussian function in order to calculate the possible shift (to be called here the weak-field asymmetry) of the maximum of the Gaussian distribution from zero. We estimate the statistical significance of the weak-field asymmetry for each rotation. We also calculate several versions of lower-resolution synoptic maps from the high-resolution maps and calculate their rotational weak-field asymmetries. We find that the weak-field asymmetries increase with decreasing map resolution. A very large fraction of weak-field asymmetries are statistically significant, with the fraction of significant weak-field asymmetries increasing with decreasing resolution. Significant weak-field asymmetries of high- and low-resolution maps mainly occur at the same times and have the same sign. Weak-field asymmetries for the different data-sets and resolutions vary quite similarly in time, and their mutual correlations are very high, especially for low-resolution maps. These results give strong evidence for weak-field asymmetries reflecting a real feature of weak field values, which is best seen in medium- and low-resolution synoptic maps and is most likely related to the supergranulation scale of the photospheric field.

We use recently digitized sunspot drawings from Mount Wilson Observatory to investigate the latitudinal dependence of tilt angles of active regions and its change with solar cycle. The drawings cover the period from 1917 to present and contain information about polarity and strength of magnetic field in sunspots. We identify clusters of sunspots of same polarity, and used these clusters to form ``bipole pairs''. The orientation of these bipole pairs was used to measure their tilts. We find that the latitudinal profile of tilts does not monotonically increase with latitude as most previous studies assumed, but instead, it shows a clear maximum at about 25--30 degree latitudes. Functional dependence of tilt ($\gamma$) on latitude ($\varphi$) was found to be $\gamma = (0.20\pm 0.08) \sin (2.80 \varphi) + (-0.00\pm 0.06)$. We also find that latitudinal dependence of tilts varies from one solar cycle to another, but larger tilts do not seem to result in stronger solar cycles. Finally, we find the presence of a systematic offset in tilt of active regions (non-zero tilts at the equator), with odd cycles exhibiting negative offset and even cycles showing the positive offset.

Using a recently developed two-scale formalism to determine the magnetic helicity spectrum (Brandenburg et al. 2017), we analyze synoptic vector magnetograms built with data from the Vector Spectromagnetograph (VSM) instrument on the \emph{Synoptic Optical Long-term Investigations of the Sun} (SOLIS) telescope during January 2010-July 2016. In contrast to an earlier study using only three Carrington rotations, our analysis includes 74 synoptic Carrington rotation maps. We recover here bihelical spectra at different phases of solar cycle~24, where the net magnetic helicity in the majority of the data is consistent with a large-scale dynamo with helical turbulence operating in the Sun. More than $20\%$ of the analyzed maps, however, show violations of the expected sign rule.

Coronal holes (CH) are regions of open magnetic field lines in the solar corona and the source of fast solar wind. Understanding the evolution of coronal holes is critical for solar magnetism as well as for accurate space weather forecasts. We study here the extreme ultraviolet (EUV) synoptic maps at three wavelengths (195A/193A, 171A and 304A) measured by Solar and Heliospheric Observatory/Extreme Ultraviolet Imaging Telescope (SOHO/EIT) and Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) instruments. The two datasets are first homogenized by scaling the SDO/AIA data to the SOHO/EIT level by means of histogram equalization. We then develop a novel automated method to identify CHs from these homogenized maps by determining the intensity threshold of CH regions separately for each synoptic map. This is done by identifying the best location and size of an image segment, which optimally contains portions of coronal holes and the surrounding quiet Sun allowing us to detect the momentary intensity threshold. Our method is thus able to adjust itself to the changing scale size of coronal holes and to temporally varying intensities. To make full use of the information in the three wavelengths we construct, a composite CH distribution, which is more robust than distributions based on one wavelength. Using the composite CH dataset we discuss the temporal evolution of CHs during the solar cycles 23 and 24.

Intensity scintillations of cosmic radio sources are used to study astrophysical plasmas like the ionosphere, the solar wind, and the interstellar medium. Normally these observations are relatively narrow band. With Low Frequency Array (LOFAR) technology at the Kilpisj\"arvi Atmospheric Imaging Receiver Array (KAIRA) station in northern Finland we have observed scintillations over a 3 octave bandwidth. ``Parabolic arcs'', which were discovered in interstellar scintillations of pulsars, can provide precise estimates of the distance and velocity of the scattering plasma. Here we report the first observations of such arcs in the ionosphere and the first broad-band observations of arcs anywhere, raising hopes that study of the phenomenon may similarly improve the analysis of ionospheric scintillations. These observations were made of the strong natural radio source Cygnus-A and covered the entire 30-250\,MHz band of KAIRA. Well-defined parabolic arcs were seen early in the observations, before transit, and disappeared after transit although scintillations continued to be obvious during the entire observation. We show that this can be attributed to the structure of Cygnus-A. Initial results from modeling these scintillation arcs are consistent with simultaneous ionospheric soundings taken with other instruments, and indicate that scattering is most likely to be associated more with the topside ionosphere than the F-region peak altitude. Further modeling and possible extension to interferometric observations, using international LOFAR stations, are discussed.