Stable Sunspot Area Level of Debrecen Photoheliographic Data and Multivariate Correction Factor of SOON Data

We investigate the spatial and temporal variation of sunspot group areas reported by the Greenwich Photoheliographic Results (GPR), the Solar Optical Observing Network (SOON), the Kislovodsk Mountain Astronomical Station (KMAS), and the Debrecen Photoheliographic Data (DPD) databases. We identify improved correction factors for reconciling these individual records to a common scale. Our results show that the DPD sunspot group areas are stable over the studied interval (1974?–?2014). We find an improved fit between GPR and DPD sunspot group areas when using a correction factor such that \(\mathrm{GPR} = 0.975(\pm 0.006) \times \mathrm{DPD}\), independent of the position of the sunspot group on the solar disk. We also find that the scale of KMAS sunspot group areas fits that of DPD well, but has a small position-dependent trend near the limb. However, in order to set SOON sunspot group area records onto the scale of DPD, we find that there is a need for a multivariate correction factor. This multivariate correction factor has a value ranging between 1.1 and 1.9 and is dependent upon the time of the SOON observation, the distance of the group from disk center, and the observatory within the SOON network. Finally, we provide further context to the systematic bias in SOON sunspot group area observations toward lower values relative to those recorded in the GPR and DPD databases that has previously been reported in the literature. We have identified the two main contributors to the SOON area deficit; some penumbral parts are unobserved, and the spot areas are underestimated. Our analysis is vital for studies that require stable, long-term solar activity records such as solar irradiance models that estimate irradiance reduction from records of sunspot group numbers, areas, and locations.     

INTER-CYCLE VARIATIONS OF SOLAR IRRADIANCE: SUNSPOT AREAS AS A POINTER

Most of the present models and reconstructions of solar irradiance use the concept of Photometric Sunspot Index (PSI) to account for the influence of sunspots on solar brightness. Since PSI is based on measured sunspot areas a firm database of such areas is essential. We show, however, that a significant disagreement exists between the data provided by the Royal Greenwich Observatory (from 1874 to 1976) and newer measurements provided by the observatories of Rome, Yunnan, Catania, and the US Air Force. The overlap of the time intervals over which sunspot areas were measured at Greenwich and Rome allows us to quantify the difference between the Greenwich and other data sets. We find that the various data sets differ, at least in a statistical sense, mainly by a correction factor of between 1.15 and 1.25.The revised time series of sunspot areas correlates well with the Zürich sunspot relative numbers over the last 120 years, with the relationship between sunspot areas and sunspot numbers changing only slightly from one cycle to the next. In particular, no indication exists for any extraordinary magnetic behavior of the Sun during the last 2 decades, as might falsely be concluded if the various sunspot area data sets are uncritically combined. There are, however, some indications that cycles 15 and 16 deviate from the rest. We expect that our results should have a significant influence on the reconstruction of the historical solar irradiance.     

大黑子群对太阳总辐照度的影响

本文用云南天文台在第２２周太阳活动峰年期间拍摄到的大太阳黑子群照相资料，太阳黑子目视描述资料，以及Ｎｉｍｂｕｓ—７卫星上辐射计测量的太阳总辐照度，分别计算了太阳总辐射照度与大黑子群的本影视面积，大黑子群全群视面积和日面上全部黑子的总视面积的相关系数。结果表明，太阳总辐射照度与这三种视面积均存在强的负相关。其中与大黑子群本影视面积的相关最强，其次是与全群视面积的相关，最后是与日面上全部黑子的总视面积的相关。并对以上结果和其它有关结果进行了分析和讨论。     

Study of Differences Between Sunspot and White Light Facular Area Data Determined from SDO/HMI and SOHO/MDI Observations

Sunspot and white light facular areas are important data for solar activity and are used, for example, in the study of the evolution of sunspots and their effect on solar irradiance. Solar Dynamic Observatory??s Helioseismic and Magnetic Imager (SDO/HMI) solar images have much higher resolution (??0.5????pixel^?1) than Solar and Heliospheric Observatory??s Michelson Doppler Imager (SOHO/MDI) solar images (??2????pixel^?1). This difference in image resolution has a significant impact on the sunspot and white light facular areas measured in the two datasets. We compare the area of sunspots and white light faculae derived from SDO/HMI and SOHO/MDI observations. This comparison helps the calibration of the SOHO sunspot and facular area to those in SDO observations. We also find a 0.22 degree difference between the North direction in SDO/HMI and SOHO/MDI images.     

Periodicities in Irradiance and in other Solar Activity Indices During Cycle 23

Magnetic fields give rise to distinctive features in different solar atmospheric regimes. To study this, time variations of the flare index, sunspot number and sunspot area, each index arising from different physical conditions, were compared with the solar composite irradiance throughout cycle 23. Rieger-type periodicities in these time series were calculated using Fourier and wavelet transforms (WTs). The peaks of the wavelet power of these periodicities appeared between the years 1999 and 2002. We found that the solar irradiance oscillations are less significant than those in the other indices during this cycle. The irradiance shows non-periodic fluctuations during this time interval. The peaks of the flare index, sunspot number and sunspot total area were seen around 2000.4, 1999.9 and 2001.0, respectively. These periodicities appeared intermittently and were not simultaneous in different solar activity indices during the three years of the maximum phase of solar cycle 23.     

Magnetic fields of sunspots based on combined optical and radio observations

In the present paper we present the results of measurement of magnetic fields in some sunspots at different heights in the solar atmosphere, based on simultaneous optical and radio measurements. The optical measurements were made by traditional photographic spectral observations of Zeeman splitting in a number of spectral lines originating at different heights in the solar photosphere and chromosphere. Radio observations of the spectra and polarization of the sunspot - associated sources were made in the wavelength range of 2–4 cm using large reflector-type radio telescope RATAN-600. The magnetic field penetrating the hot regions of the solar atmosphere were found from the shortest wavelength of generation of thermal cyclotron emission (presumably in the third harmonic of electron gyrofrequency). For all the eight cases under consideration we have found that magnetic field first drops with height, increases from the photosphere to lower chromosphere, and then decreases again as we proceed to higher chromosphere and chromosphere-corona transition region. Radio measurements were found to be well correlated with optical measurements of magnetic fields for the same sunspot. An alternative interpretation implies that different lines used for magnetic field measurements refer to different locations on the solar surface. If this is the case, then the inversion in vertical gradients of magnetic fields may not exist above the sunspots. Possible sources of systematic and random errors are also discussed.     

Spatial and temporal variations in the solar brightness

We have investigated long-term variations of solar brightness as a function of both time and solar latitude using eight years of ground-based photometric data in conjunction with space-based irradiance data. In particular, we have examined whether the combination of sunspot brightness deficits and facular brightness excesses is sufficient to explain the solar cycle irradiance variations. After correcting for the contribution from sunspots, we find that the irradiance data can be adequately explained by a model in which the remaining brightness variations are due entirely to facular contributions confined to the magnetically active latitudes. Thus we find no support for the hypothesis that there are convectively driven hot bands in the active latitudes, and our data show brightness variations that are well described by a facular contrast function.     

Total Solar Irradiance Variation During Rapid Sunspot Growth

Large sunspot areas correspond to dips in the total solar irradiance (TSI), a phenomenon associated with the local suppression of convective energy transport in the spot region. This results in a strong correlation between sunspot area and TSI. During the growth phase of a sunspot other physics may affect this correlation; if the physical growth of the sunspot resulted in surface flows affecting the temperature, for example, we might expect to see an anomalous variation in TSI. In this paper we study NOAA active region 8179, in which large sunspots suddenly appeared near disk center, at a time (March 1998) when few competing sunspots or plage regions were present on the visible hemisphere. We find that the area/TSI correlation does not significantly differ from the expected pattern of correlation, a result consistent with a large thermal conductivity in solar convection zone. In addition we have searched for a smaller-scale effect by analyzing white-light images from MDI (the Michelson Doppler Imager) on SOHO. A representative upper-limit energy consistent with the images is on the order of 3×10³¹ ergs, assuming the time scale of the actual spot area growth. This is of the same order of magnitude as the buoyant energy of the spot emergence even if it is shallow. We suggest that detailed image analyses of sunspot growth may therefore show `transient bright rings' at a detectable level.     

A Comparison of Feature Classification Methods for Modeling Solar Irradiance Variation

Physical understanding of total and spectral solar irradiance variation depends upon establishing a connection between the temporal variability of spatially resolved solar structures and spacecraft observations of irradiance. One difficulty in comparing models derived from different data sets is that the many ways for identifying solar features such as faculae, sunspots, quiet Sun, and various types of “network” are not necessarily consistent. To learn more about classification differences and how they affect irradiance models, feature “masks” are compared as derived from five current methods: multidimensional histogram analysis of NASA/National Solar Observatory/Kitt Peak spectromagnetograph data, statistical pattern recognition applied to SOHO/Michelson Doppler Imager photograms and magnetograms, threshold masks allowing for influence of spatial surroundings applied to NSO magnetograms, and “one-trigger” and “three-trigger” algorithms applied to California State University at Northridge Cartesian Full Disk Telescope intensity observations. In general all of the methods point to the same areas of the Sun for labeling sunspots and active-region faculae, and available time series of area measurements from the methods correlate well with each other and with solar irradiance. However, some methods include larger label sets, and there are important differences in detail, with measurements of sunspot area differing by as much as a factor of two. The methods differ substantially regarding inclusion of fine spatial scale in the feature definitions. The implications of these differences for modeling solar irradiance variation are discussed. K.L. Harvey and S.R. Walton are deseased, to whom this paper is dedicated.     

A Curious History of Sunspot Penumbrae

Daily records of sunspot group areas compiled by the Royal Observatory, Greenwich, from May of 1874 through 1976 indicate a curious history for the penumbral areas of the smaller sunspot groups. On average, the ratio of penumbral area to umbral area in a sunspot group increases from 5 to 6 as the total sunspot group area increases from 100 to 2000 μHem (a μHem is 10^?6 the area of a solar hemisphere). This relationship does not vary substantially with sunspot group latitude or with the phase of the sunspot cycle. However, for the sunspot groups with total areas <?100 μHem, this ratio changes dramatically and systematically through this historical record. The ratio for these smallest sunspots is near 5.5 from 1874 to 1900. After a rapid rise to more than 7 in 1905, it drops smoothly to less than 3 by 1930 and then rises smoothly back to more than 7 in 1961. It then returns to near 5.5 from 1965 to 1976. The smooth variation from 1905 to 1961 shows no indication of any step-like changes that might be attributed to changes in equipment or personnel. The overall level of solar activity was increasing monotonically during this time period when the penumbra-to-umbra area ratio dropped to less than half its peak value and then returned. If this history can be confirmed by other observations (e.g. Mt. Wilson or Kodaikanal), it may impact our understanding of penumbra formation, our dynamo models, and our estimates of historical changes in the solar irradiance.     

Automatic Detection of Sunspots and Extractionof Their Feature Parameters

Sunspots are solar features located in active regions of the Sun, whose number is an indicator of the Sun's magnetic activity. With a substantial increase in the quantity of solar image data, the automated detection and verification of various solar features have become increasingly important for the accurate and timely forecasts of solar activity and space weather. In order to use the high time-cadence SDO/HMI data to extract the main sunspot features for forecasting solar activities, we have established an automatic detection method of sunspots based on mathematical morphology, and calculated the sunspot group area and sunspot number. By comparing our results with those obtained from the Solar Region Summary compiled by NOAA/SWPC, it is found that the sunspot group areas and sunspot numbers computed with our algorithm are in good agreement with the active region values released by SWPC, and the corresponding correlation coefficients for the sunspot group area and sunspot number are 0.77 and 0.79, respectively. By using the method of this paper, the high time-cadence feature parameters can be obtained from the HMI data to provide the timely and accurate inputs for the solar activity forecast.     

The Greenwich Photo-heliographic Results (1874 – 1976): Summary of the Observations, Applications, Datasets, Definitions and Errors

The measurements of sunspot positions and areas that were published initially by the Royal Observatory, Greenwich, and subsequently by the Royal Greenwich Observatory (RGO), as the Greenwich Photo-heliographic Results (GPR), 1874?–?1976, exist in both printed and digital forms. These printed and digital sunspot datasets have been archived in various libraries and data centres. Unfortunately, however, typographic, systematic and isolated errors can be found in the various datasets. The purpose of the present paper is to begin the task of identifying and correcting these errors. In particular, the intention is to provide in one foundational paper all the necessary background information on the original solar observations, their various applications in scientific research, the format of the different digital datasets, the necessary definitions of the quantities measured, and the initial identification of errors in both the printed publications and the digital datasets. Two companion papers address the question of specific identifiable errors; namely, typographic errors in the printed publications, and both isolated and systematic errors in the digital datasets. The existence of two independently prepared digital datasets, which both contain information on sunspot positions and areas, makes it possible to outline a preliminary strategy for the development of an even more accurate digital dataset. Further work is in progress to generate an extremely reliable sunspot digital dataset, based on the programme of solar observations supported for more than a century by the Royal Observatory, Greenwich, and the Royal Greenwich Observatory. This improved dataset should be of value in many future scientific investigations.     

Variation of the solar constant during 1983 and 1984

Measurements of the Nimbus-7/ERB and SMM/ACRIM radiometers indicated several dips in the total solar irradiance in 1983 and in the first part of 1984. The dips in 1983, which should have a real solar origin, were selected according to the peaks of the projected areas of the active sunspot groups above the 2 error limit of their data set. In the first part of 1984 the sunspot activity was strong and few irradiance dips with relatively large amplitudes were observed. In the second part of 1984 the sunspot activity disappeared and at that time the solar constant only fluctuated around its mean.     

From Sunspot Area to Solar Variability: A Linear Transformation

The relationship between sunspot area and other observable solar parameters, such as spectral solar irradiance or total magnetic flux, is frequently sought by examining scatterplots of daily data, which generally show a non-linear distribution of points. We show that the scatterplots are consistent with our published result that these observable solar parameters are related to sunspot area by a transformation that is both linear and time invariant, namely by convolution with a finite impulse response function. Most solar parameters are affected by extended active regions, not just by sunspots. The fact that a complex active region evolves much more slowly than its associated sunspots provides a simple physical explanation of the observed non-linearities in scatterplots.     

Modeling Solar Spectral Irradiance and Total Magnetic Flux Using Sunspot Areas

We show that daily sunspot areas can be used in a simple, single parameter model to reconstruct daily variations in several other solar parameters, including solar spectral irradiance and total magnetic flux. The model assumes that changes in any given parameter can be treated mathematically as the response of the system to the emergence of a sunspot. Using cotemporal observational data, we compute the finite impulse response (FIR) function that describes that response in detail, and show that the response function has been approximately stationary over the time period for which data exist. For each parameter, the impulse response function describes the physical evolution of that part of a solar active region that is the source of the measured variability. We show that the impulse response functions are relatively narrow functions, no more than 3 years wide overall. Each exhibits a pre-active, active, and post-active region component; the active region component dominates the variability of most of the parameters studied.     

Systematic Time Delay of Hemispheric Solar Activity

Five solar-activity indices – the monthly-mean sunspot numbers from January 1945 to March 2008, the monthly-mean sunspot areas during the period of May 1874 to March 2008, the monthly numbers of sunspot groups from May 1874 to May 2008, the monthly-mean flare indices from January 1966 to December 2006, and the numbers of solar filaments per Carrington rotation in the time interval of solar rotations 876 to 1823 – have been used to show a systematic time delay between northern and southern hemispheric solar activities in a cycle. It is found that solar activity does not occur synchronously in the northern and southern hemispheres, and there is a systematic time lag or lead (phase shift) between northern and southern hemispheric solar activity in a cycle. About an eight-cycle period is inferred to exist in such phase shifts. The activity on the Sun may be governed by two different and coupled processes, not by a single process.     

Sunspot cycle variations of ensemble-averaged active regions

We examine published sunspot and calcium plage areas for 1620 solar active regions between 1974 and 1985. With these data we study the properties of ensemble-averaged active regions. The average sunspot area per region, the average plage to sunspot area ratio, and the average plage intensity of regions all vary significantly with the sunspot cycle and in correlation with one another. The average plage area per region varies significantly but is uncorrelated with the sunspot cycle and with the other quantities. While the plage and sunspot areas and the plage intensities of individual active regions observed over a two-year period are strongly correlated, the relationship among these quantities appears to change over an 11-yr period. These results suggest the existence of some energetic connection between active region sunspot areas and plage intensities. Further, if energy balance between sunspot luminosity deficits and facular luminosity excesses holds, then standard models relating these quantities to sunspot and plage areas will have to be modified. Overall energy balance can neither be established nor ruled out.Solar Cycle Workshop Paper.     

On the Relation Between Solar Activity and Clear-Sky Terrestrial Irradiance

The Mauna Loa Observatory record of direct-beam solar irradiance measurements for the years 1958?–?2010 is analysed to investigate the variation of clear-sky terrestrial insolation with solar activity over more than four solar cycles. The raw irradiance data exhibit a marked seasonal cycle, extended periods of lower irradiance due to emissions of volcanic aerosols, and a long-term decrease in atmospheric transmission independent of solar activity. After correcting for these effects, it is found that clear-sky terrestrial irradiance typically varies by ≈?0.2±0.1 % over the course of the solar cycle, a change of the same order of magnitude as the variations of the total solar irradiance above the atmosphere. An investigation of changes in the clear-sky atmospheric transmission fails to find a significant trend with sunspot number. Hence there is no evidence for a yet unknown effect amplifying variations of clear-sky irradiance with solar activity.     

Variation of the solar constant connected with the coronal activity

A strong correlation was found between the dips in the total solar irradiance and the peaks in the active sunspot areas as well as in the 260 MHz coronal radio flux. This connection might indicate that Alfvén-waves, generated during the interaction of the magnetic fields of the active sunspot groups with the convection, are able to transport away part of the missing energy in the solar constant decreases. These waves can heat the solar corona above the sunspot groups. Another part of the missing energy could be re-radiated later, for example during the decay of the active regions.