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.

In the recent years there has been a resurgence of the study of Extended Solar Cycles (ESCs) through observational proxies mainly in Extreme Ultraviolet. But most of them are limited only to space-based era covering only about two solar cycles. Long-term historical data-sets are worth in examining the consistency of ESCs. Kodaikanal Solar Observatory (KSO) and Mount Wilson Observatory (MWO) are the two major sources of long-term Ca II K digitised spectroheliograms covering the temporal spans 1907-2007 and 1915-1985 respectively. In this study, we detected supergranule boundaries, commonly known as networks, using the Carrington maps from both KSO and MWO datasets. Subsequently we excluded the plage areas to consider only quiet sun (QS) and detected small scale bright features through intensity thresholding over the QS network. Latitudinal density of those features, which we named as `Network Bright Elements' (NBEs), could clearly depict the existence of overlapping cycles with equator-ward branches starting at latitude $\approx 55^{\circ}$ and taking about $15\pm1$ years to reach the equator. We performed superposed epoch analysis to depict the similarity of those extended cycles. Knowledge of such equator-ward band interaction, for several cycles, may provide critical constraints on solar dynamo models.

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.

We study the variability of solar activity using new photospheric proxies originally developed for the analysis of stellar magnetism with the CoRoT and Kepler photometric observations. These proxies are obtained by tracking the temporal modulations in the observations associated to the spots and magnetic features as the Sun rotates. We analyze here 21 years, spanning solar cycles 23 and 24, of the observations collected by the space-based photometric VIRGO and radial velocity GOLF instruments on board the SoHO satellite. The photospheric activity proxy Sph is then calculated for each of the three VIRGO photometers as well as the associated Svel proxy from the radial velocity GOLF observations. Comparisons with several standard solar activity proxies sensitive to different layers of the Sun demonstrate that these new activity proxies, Sph and Svel, provide a new manner to monitor solar activity. We show that both the long- and short-term magnetic variabilities respectively associated to the 11-year cycle and the quasi-biennial oscillation are well monitored, and that the magnetic field interaction between the subsurface, the photosphere, and the chromosphere of the Sun, was modified between Cycle 24 and Cycle 23. Furthermore, the photometric proxies show a wavelength dependence of the response function of the solar photosphere among the three channels of the VIRGO photometers, providing inputs for the study of the stellar magnetism of Sun-like stars.

We describe here a procedure to combine measurements in the 393.37 nm Ca II K spectral line taken at different observatories. Measurements from the National Solar Observatory (NSO) Integrated Sunlight Spectrometer (ISS) on the Synoptic Optical Long-term Investigations of the Sun (SOLIS) telescope, the NSO/Sac Peak Ca II K-Line Monitoring Program, and Ca II K filtergrams from Kodaikanal Solar Observatory (KKL) are merged together to create a pair of composites of the Ca II K 1-A emission index. These composites are publicly available from the SOLIS website at http://solis.nso.edu/0/iss/.

The most commonly used index of stellar magnetic activity is the instrumental flux scale of singly-ionized calcium H & K line core emission, S, developed by the Mount Wilson Observatory (MWO) HK Project, or the derivative index R'_HK. Accurately placing the Sun on the S scale is important for comparing solar activity to that of the Sun-like stars. We present previously unpublished measurements of the reflected sunlight from the Moon using the second-generation MWO HK photometer during solar cycle 23 and determine cycle minimum S_min,23 = 0.1634 +/- 0.0008, amplitude Delta S_23 = 0.0143 +/- 0.0012, and mean <S_23> = 0.1701 +/- 0.0005. By establishing a proxy relationship with the closely related National Solar Observatory Sacramento Peak calcium K emission index, itself well-correlated with the Kodaikanal Observatory plage index, we extend the MWO S time series to cover cycles 15-24 and find on average <S_min> = 0.1621 +/- 0.0008, <Delta S_cyc> = 0.0145 +/- 0.0012, <S_cyc> = 0.1694 +/- 0.0005. Our measurements represent an improvement over previous estimates which relied on stellar measurements or solar proxies with non-overlapping time series. We find good agreement from these results with measurements by the Solar-Stellar Spectrograph at Lowell Observatory, an independently calibrated instrument, which gives us additional confidence that we have accurately placed the Sun on the S-index flux scale.

We have detected several periodicities in the solar equatorial rotation rate of sunspot groups in the Greenwich Photoheliographic Results (GPR) during the period 1931-1976, the Solar Optical Observing Network (SOON) during the period 1977-2014, and the Debrecen Photoheliographic Data (DPD) during the period 1974-2014. Our results suggest a ~250-day period in the equatorial rotation rate determined from both the Mt. Wilson Doppler-velocity data and the sunspot-group data during 1986-2007. However, a wavelet analysis reveals that this periodicity appears mostly around 1991 in the velocity data, while it is present in most of the solar cycles covered by the sunspot-group data, mainly near the minimum epochs of the solar cycles. We also found the signature of a period of ~1.4 years period in the velocity data during 1990-1995, and in the equatorial rotation rate of sunspot groups mostly around the year 1956. The equatorial rotation rate of sunspot groups reveals a strong ~1.6-year periodicity around 1933 and 1955 and a weaker one around 1976, and a strong ~1.8-year periodicity around 1943. Our analysis also suggests periodicities of ~5 years, ~7 years, and ~17 years as well as some other short-term periodicities. However, short-term periodicities are mostly present at the time of solar minima. Hence, short-term periodicities cannot be confirmed because of the larger uncertainty in the data.

Long-term synoptic observations in the resonance line of Ca II K constitute a fundamental database for a variety of retrospective analyses of the state of the solar magnetism. Synoptic Ca II K observations began in late 1904 at the Kodaikanal Observatory, in India. In early 1970s, the National Solar Observatory (NSO) at Sacramento Peak (USA) started a new program of daily Sun-as-a-star observations in the Ca II K line. Today the NSO is continuing these observations through its Synoptic Optical Long-term Investigations of the Sun (SOLIS) facility. These different data sets can be combined into a single disk-integrated Ca II K index time series that describes the average properties of the chromospheric emission over several solar cycles. We present such a Ca II K composite and discuss its correlation with the new entirely revised sunspot number data series. For this preliminary investigation, the scaling factor between pairs of time series was determined assuming a simple linear model for the relationship between the monthly mean values during the duration of overlapping observations.

Recently, the National Solar Observatory (NSO) Solar-atmosphere Pipeline Working Group has undertaken the production of synoptic maps from Helioseismic and Magnetic Imager (HMI) magnetograms. A set of maps has been processed spanning the data available for 2010-2015 using twice daily images (taken at UT midnight and noon) and running them through the same algorithms used to produce SOLIS/VSM 6302l mean-magnetic and spatial-variance maps. The contents of this document provide an overview of what these maps look like, and the processing steps used to generate them from the original HMI input data.

Analysis of synoptic data from the Vector Stokes Magnetograph (VSM) of the Synoptic Optical Long-term Investigations of the Sun (SOLIS) and the NASA/NSO Spectromagnetograph (SPM) at the NSO/Kitt Peak Vacuum Telescope facility shows that the reversals of solar polar magnetic fields exhibit elements of a stochastic process, which may include the development of specific patterns of emerging magnetic flux, and the asymmetry in activity between northern and southern hemispheres. The presence of such irregularities makes the modeling and prediction of polar field reversals extremely hard if possible. In a classical model of solar activity cycle, the unipolar magnetic regions (UMRs) of predominantly following polarity fields are transported polewards due to meridional flows and diffusion. The UMRs gradually cancel out the polar magnetic field of the previous cycle, and re-build the polar field of opposite polarity setting the stage for the next cycle. We show, however, that this deterministic picture can be easily altered by the developing of a strong center of activity, or by the emergence of an extremely large active region, or by a "strategically placed" coronal hole. We demonstrate that the activity occurring during the current cycle 24 may be the result of this randomness in the evolution of the solar surface magnetic field.

The Vector Spectromagnetograph (VSM) instrument on the Synoptic Optical Long-term Investigations of the Sun (SOLIS) telescope is designed to obtain high-quality magnetic field observations in both the photosphere and chromosphere by measuring the Zeeman-induced polarization of spectral lines. With 1$^{\prime \prime}$ spatial resolution (1.14$^{\prime \prime}$ before 2010) and 0.05\AA\ spectral resolution, the VSM provides, among other products, chromospheric full-disk magnetograms using the CaII 854.2 nm spectral line and both photospheric full-disk vector and longitudinal magnetograms using the FeI 630.15 nm line. Here we describe the procedure used to compute daily weighted averages of the photospheric radial polar magnetic field at different latitude bands from SOLIS/VSM longitudinal full-disk observations. Time series of these measurements are publicly available from the SOLIS website at http://solis.nso.edu/0/vsm/vsm\_plrfield.html. Future plans include the calculation of the mean polar field strength from SOLIS/VSM chromospheric observations and the determination of the {\it true} radial polar field from SOLIS/VSM full-Stokes measurements.

Unlike other stars, the surface of the Sun can be spatially resolved to a high degree of detail. But the Sun can also be observed as if it was a distant star. The availability of solar disk-resolved and disk-integrated spectra offers an opportunity to devise methods to derive information about the spatial distribution of solar features from Sun-as-a-star measurements. Here, we present an update on work done at the National Solar Observatory to reconcile disk-integrated and disk-resolved solar spectra from the Synoptic Optical Long-term Investigation of the Sun (SOLIS) station. The results of this work will lead to a new approach to infer the information about the spatial distribution of features on other stars, from the overall filling factor of active regions to, possibly, the latitude/longitude distribution of features.

Magnetic flux synoptic charts are critical for reliable modeling of the corona and heliosphere. Until now, however, these charts were provided without any estimate of uncertainties. The uncertainties are due to instrumental noise in the measurements and to the spatial variance of the magnetic flux distribution that contributes to each bin in the synoptic chart. We describe here a simple method to compute synoptic magnetic flux maps and their corresponding magnetic flux spatial variance charts that can be used to estimate the uncertainty in the results of coronal models. We have tested this approach by computing a potential-field-source-surface model of the coronal field for a Monte Carlo simulation of Carrington synoptic magnetic flux maps generated from the variance map. We show that these uncertainties affect both the locations of source-surface neutral lines and the distributions of coronal holes in the models.

Measurements from the Mount Wilson Observatory (MWO) are used to study the long-term variations of sunspot field strengths from 1920 to 1958. Following a modified approach similar to that in Pevtsov et al. (2011), for each observing week we select a single sunspot with the strongest field strength measured that week and then compute monthly averages of these weekly maximum field strengths. The data show the solar cycle variation of the peak field strengths with an amplitude of about 500-700 gauss (G), but no statistically significant long-term trends. Next, we use the sunspot observations from the Royal Greenwich Observatory (RGO) to establish a relationship between the sunspot areas and the sunspot field strengths for Cycles 15-19. This relationship is then used to create a proxy of peak magnetic field strength based on sunspot areas from the RGO and the USAF/NOAA network for the period from 1874 to early 2012. Over this interval, the magnetic field proxy shows a clear solar cycle variation with an amplitude of 500-700 G and a weaker long-term trend. From 1874 to around 1920, the mean value of magnetic field proxy increases by about 300-350 G, and, following a broad maximum in 1920-1960, it decreases by about 300 G. Using the proxy for the magnetic field strength as the reference, we scale the MWO field measurements to the measurements of the magnetic fields in Pevtsov et al. (2011) to construct a combined data set of maximum sunspot field strengths extending from 1920 to early 2012. This combined data set shows strong solar cycle variations and no significant long-term trend (linear fit to the data yields a slope of $-0.2\pm$0.8 G year$^{-1}$). On the other hand, the peak sunspot field strengths observed at the minimum of the solar cycle show a gradual decline over the last three minima (corresponding to cycles 21-23) with a mean downward trend of $\approx$ 15 G year$^{-1}$.

This article presents a comparative analysis of solar activity data, Mt Wilson diameter data, Super-Kamiokande solar neutrino data, and nuclear decay data acquired at the Lomonosov Moscow State University (LMSU). We propose that salient periodicities in all of these datasets may be attributed to r-mode oscillations. Periodicities in the solar activity data and in Super-Kamiokande solar neutrino data may be attributed to r-mode oscillations in the known tachocline, with normalized radius in the range 0.66 to 0.74, where the sidereal rotation rate is in the range 13.7 to 14.6 year-1. We propose that periodicities in the Mt Wilson and LMSU data may be attributed to similar r-mode oscillations where the sidereal rotation rate is approximately which we attribute to a hypothetical inner tachocline separating a slowly rotating core from the radiative zone. We also discuss the possible role of the RSFP (Resonant Spin Flavor Precession) process, which leads to estimates of the neutrino magnetic moment and of the magnetic field strength in or near the solar core.

The characterization of solar surface differential rotation (SDR) from disk-integrated chromospheric measurements has important implications for the study of differential rotation and dynamo processes in other stars. Some chromospheric lines, such as Ca II K, are very sensitive to the presence of activity on the disk and are an ideal choice for investigating SDR in Sun-as-a star observations. Past studies indicate that when the activity is low, the determination of Sun's differential rotation from integrated-sunlight measurements becomes uncertain. However, our study shows that using the proper technique, SDR can be detected from these type of measurements even during periods of extended solar minima. This paper describes results from the analysis of the temporal variations of Ca II K line profiles observed by the Integrated Sunlight Spectrometer (ISS) during the declining phase of Cycle 23 and the rising phase of Cycle 24, and discusses the signature of SDR in the power spectra computed from time series of parameters derived from these profiles. The described methodology is quite general, and could be applied to photometric time series of other Main-Sequence stars for detecting differential rotation.

We compare photospheric line-of-sight magnetograms from the Synoptic Long-term Investigations of the Sun (SOLIS) vector spectromagnetograph (VSM) instrument with observations from the 150-foot Solar Tower at Mt. Wilson (MWO), Helioseismic and Magnetic Imager (HMI) on Solar Dynamics Observatory (SDO), and Michelson Doppler Imager (MDI) on Solar and Heliospheric Observatory (SOHO). We find very good agreement between VSM and the other data sources for both disk-averaged flux densities and pixel-by-pixel measurements. We show that the VSM mean flux density time series is of consistently high signal-to-noise with no significant zero-offsets. We discuss in detail some of the factors -spatial resolution, flux dependence and position on the solar disk- affecting the determination of scaling between VSM and SOHO/MDI or SDO/HMI magnetograms. The VSM flux densities agree well with spatially smoothed data from MDI and HMI, although the scaling factors show clear dependence on flux density. The factor to convert VSM to HMI increases with increasing flux density (from $\approx$1 to $\approx$1.5). The nonlinearity is smaller for the VSM vs. ~SOHO/MDI scaling factor (from $\approx$1 to $\approx$1.2).

Measurements of the ionized Ca II K line are one of the major resources for long-term studies of solar and stellar activity. They also play a critical role in many studies related to solar irradiance variability, particularly as a ground-based proxy to model the solar ultraviolet flux variation that may influence the Earth's climate. Full disk images of the Sun in Ca II K have been available from various observatories for more than 100 years and latter synoptic Sun-as-a-star observations in Ca II K began in the early 1970s. One of these instruments, the Integrated Sunlight Spectrometer (ISS) has been in operation at Kitt Peak (Arizona) since late 2006. The ISS takes daily observations of solar spectra in nine spectra bands, including the Ca II K and H line s. We describe recent improvements in data reduction of Ca II K observations, and present time variations of nine parameters derived from the profile of this spectral line.

This article presents a power-spectrum analysis of 39,024 measurements of the solar diameter made at the Mount Wilson Observatory from 1968.670 to 1997.965. This power spectrum contains a number of very strong peaks. We find that eight of these peaks agree closely with the frequencies of r-mode oscillations for a region of the Sun where the sidereal rotation frequency is 12.08 year$^{-1}$. We estimate that there is less than one chance in ten to the sixth power of finding this pattern by chance.

It is well established that both total and spectral solar irradiance are modulated by variable magnetic activity on the solar surface. However, there is still disagreement about the contribution of individual solar features for changes in the solar output, in particular over decadal time scales. Ionized Ca II K line spectroheliograms are one of the major resources for these long-term trend studies, mainly because such measurements have been available now for more than 100 years. In this paper we introduce a new Ca II K plage and active network index time series derived from the digitization of almost 40,000 photographic solar images that were obtained at the 60-foot solar tower, between 1915 and 1985, as a part of the monitoring program of the Mount Wilson Observatory. We describe here the procedure we applied to calibrate the images and the properties of our new defined index, which is strongly correlated to the average fractional area of the visible solar disk occupied by plages and active network. We show that the long-term variation of this index is in an excellent agreement with the 11-year solar cycle trend determined from the annual international sunspot numbers series. Our time series agrees also very well with similar indicators derived from a different reduction of the same data base and other \caii K spectroheliograms long-term synoptic programs, such as those at Kodaikanal Observatory (India), and at the National Solar Observatory at Sacramento Peak (USA). Finally, we show that using appropriate proxies it is possible to extend this time series up to date, making this data set one of the longest Ca II K index series currently available.

We present results using the AutoClass analysis application available at NASA/Ames Intelligent Systems Div. (2002) which is a Bayesian, finite mixture model classification system developed by Cheeseman and Stutz (1996). We apply this system to Mount Wilson Solar Observatory (MWO) intensity and magnetogram images and classify individual pixels on the solar surface to calculate daily indices that are then correlated with total solar irradiance (TSI) to yield a set of regression coefficients. This approach allows us to model the TSI with a correlation of better than 0.96 for the period 1996 to 2007. These regression coefficients applied to classified pixels on the observed solar surface allow the construction of images of the Sun as it would be seen by TSI measuring instruments like the Solar Bolometric Imager recently flown by Foukal et al., (2004). As a consequence of the very high correlation we achieve in reproducing the TSI record, our approach holds out the possibility of creating an on-going, accurate, independent estimate of TSI variations fromground-based observations which could be used to compare, and identify the sources of disagreement among, TSI observations from the various satellite instruments and to fill in gaps in the satellite record. Further, our spatially-resolved images should assist in characterizing the particular solar surface regions associated with TSI variations. Also, since the particular set of MWO data on which this analysis is based is available on a daily basis back to at least 1985, and on an intermittent basis before then, it will be possible to estimate the TSI emission due to identified solar surface features at several solar minima to constrain the role surface magnetic effects have on long-term trends in solar energy output.

The power spectral analyses of the Sun's surface equatorial rotation rate determined from the Mt. Wilson daily Doppler velocity measurements during the period 3 December 1985 to 5 March 2007 suggests the existence of 7.6 year, 2.8 year, 1.47 year, 245 day, 182 day and 158 day periodicities in the surface equatorial rotation rate during the period before 1996. However, there is no variation of any kind in the more accurately measured data during the period after 1995. That is, the aforementioned periodicities in the data during the period before the year 1996 may be artifacts of the uncertainties of those data due to the frequent changes in the instrumentation of the Mt. Wilson spectrograph. On the other hand, the temporal behavior of most of the activity phenomena during cycles 22 (1986-1996) and 23 (after 1997) is considerably different. Therefore, the presence of the aforementioned short-term periodicities during the last cycle and absence of them in the current cycle may, in principle, be real temporal behavior of the solar rotation during these cycles.

This study based on longitudinal Zeeman effect magnetograms and spectral line scans investigates the dependence of solar surface magnetic fields on the spectral line used and the way the line is sampled in order to estimate the magnetic flux emerging above the solar atmosphere and penetrating to the corona from magnetograms of the Mt. Wilson 150-foot tower synoptic program (MWO). We have compared the synoptic program \lambda5250\AA line of Fe I to the line of Fe I at \lambda5233\AA since this latter line has a broad shape with a profile that is nearly linear over a large portion of its wings. The present study uses five pairs of sampling points on the $\lambda5233$\AA line. We recommend adoption of the field determined with a line bisector method with a sampling point as close as possible to the line core as the best estimate of the emergent photospheric flux. The combination of the line profile measurements and the cross-correlation of fields measured simultaneously with \lambda5250\AA and \lambda5233\AA yields a formula for the scale factor 1/\delta that multiplies the MWO synoptic magnetic fields. The new calibration shows that magnetic fields measured by the MDI system on the SOHO spacecraft are equal to 0.619+/-0.018 times the true value at a center-to-limb position 30 deg. Berger and Lites (2003) found this factor to be 0.64+/-0.013 based on a comparison the the Advanced Stokes Polarimeter.

In a previous paper (Varadi et al., 1999), Random Lag Singular Spectrum Analysis was offered as a tool to find oscillations in very noisy and long time series. This work presents a generalization of the technique to search for common oscillations in two or more time series.