Temporal Variations in Solar Irradiance Since 1947

The study of variations in total solar irradiance (TSI) and spectral irradiance is important for understanding how the Sun affects the Earth’s climate. A data-driven approach is used in this article to analyze and model the temporal variation of the TSI and Mg?ii index back to 1947. In both cases, observed data in the time interval of the satellite era, 1978?–?2013, were used for neural network (NN) model-design and testing. For this particular purpose, the evolution of the solar magnetic field is assumed to be the main driver for the day-to-day irradiance variability. First, we design a model for the Mg?ii index data from F10.7 cm solar radio-flux using the NN approach in the time span of 1978 through 2013. Results of Mg?ii index model were tested using various numbers of hidden nodes. The predicted values of the hidden layer with five nodes correspond well to the composite Mg?ii values. The model reproduces 94% of the variability in the composite Mg?ii index, including the secular decline between the 1996 and 2008 solar cycle minima. Finally, the extrapolation of the Mg?ii index was performed using the developed model from F10.7 cm back to 1947. Similarly, the NN model was designed for TSI variability study over the time span of the satellite era using data from the Physikalisch-Meteorologisches Observatorium Davos (PMOD) as a target, and solar activity indices as model inputs. This model was able to reproduce the daily irradiance variations with a correlation coefficient of 0.937 from sunspot and facular measurements in the time span of 1978?–?2013. Finally, the temporal variation of the TSI was analyzed using the designed NN model back to 1947 from the Photometric Sunspot Index (PSI) and the extrapolated Mg?ii index. The extrapolated TSI result indicates that the amplitudes of Solar Cycles 19 and 21 are closely comparable to each other, and Solar Cycle 20 appears to be of lower irradiance during its maximum.     

Estimating the Mg <Emphasis Type="SmallCaps">ii</Emphasis> Index from 1961 Through 1981 Using Ca <Emphasis Type="SmallCaps">ii</Emphasis> K Images from the Mt Wilson Observatory

An empirical model of solar UV spectral irradiance has been developed that is based on observed spectral radiance measurements and full disk Ca ii K images. The Mg ii index is then calculated from the estimated spectra in a narrow wavelength range (180 Å) near the Mg ii doublet at 2800 Å. Our long term goal is to expand this wavelength range from 10 to 4000 Å in continuing studies based on spectral data covering this wavelength range (e.g. Skylab, UARS/SUSIM, TIMED/SEE, etc.). Our previous modeling effort produced spectra in this 180 Å range and the resulting Mg ii index values for the period from 1991 through 1995 and we have used observations during this time period to validate the model results. The current paper presents results from this model based on a 21-year portion of the recently digitized Ca ii K images from the Mt Wilson Observatory (MWO) film archive. Here we present details of the model, the required model modifications, and the resulting Mg ii index from 1961 through 1981. Since the NOAA Mg ii index did not begin until 1978, the present model results are compared to a Mg ii index estimated from the F_10.7 radio flux over this 21-year period. The NOAA Mg ii index, which is derived from measured UV spectra, is also included for comparison from late 1978 through 1981.     

Solar Irradiance from 165 to 400 nm in 2008 and UV Variations in Three Spectral Bands During Solar Cycle 24

Accurate measurements of the solar spectral irradiance (SSI) and its temporal variations are of primary interest to better understand solar mechanisms, and the links between solar variability and Earth’s atmosphere and climate. The SOLar SPECtrum (SOLSPEC) instrument of the Solar Monitoring Observatory (SOLAR) payload onboard the International Space Station (ISS) has been built to carry out SSI measurements from 165 to 3088 nm. We focus here on the ultraviolet (UV) part of the measured solar spectrum (wavelengths less than 400 nm) because the UV part is potentially important for understanding the solar forcing of Earth’s atmosphere and climate. We present here SOLAR/SOLSPEC UV data obtained since 2008, and their variations in three spectral bands during Solar Cycle 24. They are compared with previously reported UV measurements and model reconstructions, and differences are discussed.     

Comparison of Solar UV Spectral Irradiance from SUSIM and SORCE

Knowledge of solar spectral irradiance (SSI) is important in determining the impact of solar variability on climate. Observations of UV SSI have been made by the Solar Ultraviolet Spectral Irradiance Monitor (SUSIM) on the Upper Atmosphere Research Satellite (UARS), the Solar-Stellar Irradiance Comparison Experiment (SOLSTICE), and the Solar Irradiance Monitor (SIM), both on the Solar Radiation and Climate Experiment (SORCE) satellite. Measurements by SUSIM and SORCE overlapped from 2003 to 2005. SUSIM and SORCE observations represent ～?20 years of absolute UV SSI. Unfortunately, significant differences exist between these two data sets. In particular, changes in SORCE UV SSI measurements, gathered at moderate and minimum solar activity, are a factor of two greater than the changes in SUSIM observations over the entire solar cycle. In addition, SORCE UV SSI have a substantially different relationship with the Mg ii index than did earlier UV SSI observations. Acceptance of these new SORCE results impose significant changes on our understanding of UV SSI variation. Alternatively, these differences in UV SSI observations indicate that some or all of these instruments have changes in instrument responsivity that are not fully accounted for by the current calibration. In this study, we compare UV SSI changes from SUSIM with those from SIM and SOLSTICE. The primary results are that (1) long-term observations by SUSIM and SORCE generally do not agree during the overlap period (2003?–?2005), (2) SUSIM observations during this overlap period are consistent with an SSI model based on Mg ii and early SUSIM SSI, and (3) when comparing the spectral irradiance for times of similar solar activity on either side of solar minimum, SUSIM observations show slight differences while the SORCE observations show variations that increase with time between spectra. Based on this work, we conclude that the instrument responsivity for SOLSTICE and SIM need to be reevaluated before these results can be used for climate-modeling studies.     

Effects of Hysteresis Between Maximum CME Speed Index and Typical Solar Activity Indicators During Cycle 23

Using the smoothed time series of maximum CME speed index for solar cycle 23, it is found that this index, analyzed jointly with six other solar activity indicators, shows a hysteresis phenomenon. The total solar irradiance, coronal index, solar radio flux (10.7?cm), Mg?ii core-to-wing ratio, sunspot area, and H?? flare index follow different paths for the ascending and the descending phases of solar cycle?23, while a saturation effect exists at the maximum phase of the cycle. However, the separations between the paths are not the same for the different solar activity indicators used: the H?? flare index and total solar irradiance depict broad loops, while the Mg?ii core-to-wing ratio and sunspot area depict narrow hysteresis loops. The lag times of these indices with respect to the maximum CME speed index are discussed, confirming that the hysteresis represents a clue in the search for physical processes responsible for changing solar emission.     

Lifetimes of some excited states in the rare earths: La ii,Ce ii,Pr ii,Nd ii,Sm ii,Yb i,Yb ii,and Lu ii

Radiative lifetimes for excited states in La ii, Ce ii, Pr ii, Nd ii, Sm ii, Yb i, Yb ii, and Lu ii have been determined by means of the beam-foil technique or the zero-field level-crossing method. The lifetimes for La, Ce, Pr, Nd, Sm, and Yb are shorter than those computed by summing the transition probabilities of Corliss and Bozman by a factor of up to ~5. The large discrepancies between the abundance of La, Ce, Pr, Nd, and Sm in the solar photosphere and in meteorites are eliminated or greatly reduced if the abundance determinations of the solar photosphere are based upon the gf values of Corliss and Bozman, corrected for by the present lifetimes.     

A Change in the Solar He ii EUV Global Network Structure as an Indicator of the Geo-Effectiveness of Solar Minima

Solar activity during 2007?–?2009 was very low, causing anomalously low thermospheric density. A comparison of solar extreme ultraviolet (EUV) irradiance in the He?ii spectral band (26 to 34 nm) from the Solar Extreme ultraviolet Monitor (SEM), one of instruments on the Charge Element and Isotope Analysis System (CELIAS) on board the Solar and Heliospheric Observatory (SOHO) for the two latest solar minima showed a decrease of the absolute irradiance of about 15±6 % during the solar minimum between Cycles 23 and 24 compared with the Cycle 22/23 minimum when a yearly running-mean filter was used. We found that some local, shorter-term minima including those with the same absolute EUV flux in the SEM spectral band show a higher concentration of spatial power in the global network structure from the 30.4 nm SOHO/Extreme ultraviolet Imaging Telescope (EIT) images for the local minimum of 1996 compared with the minima of 2008?–?2011. We interpret this higher concentration of spatial power in the transition region’s global network structure as a larger number of larger-area features on the solar disk. These changes in the global network structure during solar minima may characterize, in part, the geo-effectiveness of the solar He?ii EUV irradiance in addition to the estimations based on its absolute levels.     

A Fast Model for the Reconstruction of Spectral Solar Irradiance in the Near- and Mid-Ultraviolet

We present a model for the reconstruction of spectral solar irradiance between 200 and 400?nm. This model is an extension of the total solar irradiance (TSI) model of Crouch et al. (Astrophys.?J. 677, 723, 2008) which is based on a data-driven Monte Carlo simulation of sunspot emergence, fragmentation, and erosion. The resulting time-evolving daily area distribution of magnetic structures of all sizes is used as input to a four-component irradiance model including contributions from the quiet Sun, sunspots, faculae, and network. In extending the model to spectral irradiance in the near- and mid-ultraviolet, the quiet Sun and sunspot emissivities are calculated from synthetic spectra at T _eff=5750?K and 5250?K, respectively. Facular emissivities are calculated using a simple synthesis procedure proposed by Solanki and Unruh (Astron. Astrophys. 329, 747, 1998). The resulting time series of ultraviolet flux is calibrated against the data from the SOLSTICE instrument on the Upper Atmospheric Research Satellite (UARS). Using a genetic algorithm, we invert quiet Sun corrections, profile of facular temperature variations with height, and network model parameters which yield the best fit to these data. The resulting best-fit time series reproduces quite well the solar-cycle timescale variations of UARS ultraviolet observations, as well as the short-timescale fluctuations about the 81 day running mean. We synthesize full spectra between 200 and 400?nm, and validate these against the spectra obtained by the ATLAS-1 and ATLAS-3 missions, finding good agreement, to better than 3?% at most wavelengths. We also compare the UV variability predicted by our reconstructions in the descending phase of sunspot cycle 23 to SORCE/SIM data as well as to other reconstructions. Finally, we use the model to reconstruct the time series of spectral irradiance starting in 1874, and investigate temporal correlations between pairs of wavelengths in the bands of interest for stratospheric chemistry and dynamics.     

Sensitivity of Selected Ba <Emphasis Type="SmallCaps">ii</Emphasis>, Fe <Emphasis Type="SmallCaps">i</Emphasis>, Fe <Emphasis Type="SmallCaps">ii</Emphasis>, and Cr <Emphasis Type="SmallCaps">i</Emphasis> Spectral Lines to Velocity in Quiet Solar Atmosphere

We examine the sensitivity of selected Ba?ii, Fe?i, Fe?ii, and Cr?i spectral lines to changes of the line-of-sight velocity by sharpness of their line profiles and response functions to line-of-sight velocity evaluated by the 1-D model of the quiet solar atmosphere in the LTE approximation. The set of selected lines includes the Ba?ii 4554 Å line, generally considered to be an excellent Doppler mapper. Our findings confirm earlier results showing that the sensitivity increases not only with wavelength, as anticipated from the Doppler relation, but mainly with the sharpness of line profiles given by the ratio of their depths and widths. The line Fe?i 5247 Å is the most sensitive in our set, whereas the Fe?i and Fe?ii infrared lines show very low sensitivity because of their large thermal widths. The line Ba?ii 4554 Å shows only moderate sensitivity due to its large width, given by a broad hyperfine structure and isotopic split. For the first time we identify a very promising and so far unknown Doppler mapper of the solar photosphere and low chromosphere, which is the line Ba?ii 6497 Å. Its sensitivity is comparable with the sensitivity of Fe?i 5247 Å and clearly surpasses the sensitivity of Ba?ii 4554 Å. The line Ba?ii 6497 Å offers many advantages, making it a highly recommendable choice for future studies of line-of-sight velocities in the photosphere and low chromosphere.     

A Possible Cause of the Diminished Solar Wind During the Solar Cycle 23 – 24 Minimum

Interplanetary magnetic field and solar wind plasma density observed at 1 AU during Solar Cycle 23?–?24 (SC-23/24) minimum were significantly smaller than those during its previous solar cycle (SC-22/23) minimum. Because the Earth’s orbit is embedded in the slow wind during solar minimum, changes in the geometry and/or content of the slow wind region (SWR) can have a direct influence on the solar wind parameters near the Earth. In this study, we analyze solar wind plasma and magnetic field data of hourly values acquired by Ulysses. It is found that the solar wind, when averaging over the first (1995.6?–?1995.8) and third (2006.9?–?2008.2) Ulysses’ perihelion (\({\sim}\,1.4~\mbox{AU}\)) crossings, was about the same speed, but significantly less dense (\({\sim}\,34~\%\)) and cooler (\({\sim}\,20~\%\)), and the total magnetic field was \({\sim}\,30~\%\) weaker during the third compared to the first crossing. It is also found that the SWR was \({\sim}\,50~\%\) wider in the third (\({\sim}\,68.5^{\circ}\) in heliographic latitude) than in the first (\({\sim}\,44.8^{\circ}\)) solar orbit. The observed latitudinal increase in the SWR is sufficient to explain the excessive decline in the near-Earth solar wind density during the recent solar minimum without speculating that the total solar output may have been decreasing. The observed SWR inflation is also consistent with a cooler solar wind in the SC-23/24 than in the SC-22/23 minimum. Furthermore, the ratio of the high-to-low latitude photospheric magnetic field (or equatorward magnetic pressure force), as observed by the Mountain Wilson Observatory, is smaller during the third than the first Ulysses’ perihelion orbit. These findings suggest that the smaller equatorward magnetic pressure at the Sun may have led to the latitudinally-wider SRW observed by Ulysses in SC-23/24 minimum.     

IRIS Observations of Spicules and Structures Near the Solar Limb

We have analyzed Interface Region Imaging Spectrograph (IRIS) spectral and slit-jaw observations of a quiet region near the South Pole. In this article we present an overview of the observations, the corrections, and the absolute calibration of the intensity. We focus on the average profiles of strong (Mg?ii h and k, C?ii and Si?iv), as well as of weak spectral lines in the near ultraviolet (NUV) and the far ultraviolet (FUV), including the Mg?ii triplet, thus probing the solar atmosphere from the low chromosphere to the transition region. We give the radial variation of bulk spectral parameters as well as line ratios and turbulent velocities. We present measurements of the formation height in lines and in the NUV continuum from which we find a linear relationship between the position of the limb and the intensity scale height. We also find that low forming lines, such as the Mg?ii triplet, show no temporal variations above the limb associated with spicules, suggesting that such lines are formed in a homogeneous atmospheric layer and, possibly, that spicules are formed above the height of \(2''\). We discuss the spatio-temporal structure of the atmosphere near the limb from images of intensity as a function of position and time. In these images, we identify p-mode oscillations in the cores of lines formed at low heights above the photosphere, slow-moving bright features in O?i and fast-moving bright features in C?ii. Finally, we compare the Mg?ii k and h line profiles, together with intensity values of the Balmer lines from the literature, with computations from the PROM57Mg non-LTE model, developed at the Institut d’ Astrophysique Spatiale, and estimated values of the physical parameters. We obtain electron temperatures in the range of \({\sim}\, 8000~\mbox{K}\) at small heights to \({\sim}\, 20\,000~\mbox{K}\) at large heights, electron densities from \(1.1\times 10^{11}\) to \(4\times 10^{10}~\mbox{cm}^{-3}\) and a turbulent velocity of \({\sim}\, 24~\mbox{km}\,\mbox{s}^{-1}\).     

Progress in ultraviolet stellar spectrophotometry with S59

A progress report is given of investigations related to observations of stellar spectra obtained with the ultraviolet stellar spectrophotometer S59 aboard the ESRO TD-1A satellite. We describe first the processing of the observations: intensity and wavelength calibration, identification of lines, classification of spectra. Thereafter some important groups of lines are dealt with: non-LTE computation of the Mgii lines are presented: apart from the peculiar emission line and shell stars they compare well with the observations; intensity ratios, of Feii and Feiii lines are an important temperature classification criterion. Interstellar UV lines indicate large deficiencies of some metals, as compared to solar values. A new ultraviolet continuous extinction curve is determined. We discuss the strong outstreaming motions observed in α Cygni, and the consequent mass loss of this star (< 3 × 10^?10 \(\mathfrak{M}_ \odot \) yr^?1and finally we describe the composite spectrum of the (WC8+09I) bnary γ² Velorum; the ultraviolet continuous spectrum of the WC star is about one magnitude brighter than any theory predicts.     

The Solar Cycle 23?�C?24 Minimum. A Benchmark in Solar Variability and Effects in the Heliosphere

Given the numerous ground-based and space-based experiments producing the database for the Cycle 23??C?24 Minimum epoch from September 2008 to May 2009, we have an extraordinary opportunity to understand its effects throughout the heliosphere. We use solar radiative output in this period to obtain minimum values for three measures of the Sun??s radiative output: the total solar irradiance, the Mg ii index, and the 10.7 cm solar radio flux. The derived values are included in the research summaries as a means to exchange ideas and data for this long minimum in solar activity.     

Sun–Earth L_{1} and L_{2} to Moon transfers exploiting natural dynamics

This paper examines the design of transfers from the Sun–Earth libration orbits, at the \(L_{1}\) and \(L_{2}\) points, towards the Moon using natural dynamics in order to assess the feasibility of future disposal or lifetime extension operations. With an eye to the probably small quantity of propellant left when its operational life has ended, the spacecraft leaves the libration point orbit on an unstable invariant manifold to bring itself closer to the Earth and Moon. The total trajectory is modeled in the coupled circular restricted three-body problem, and some preliminary study of the use of solar radiation pressure is also provided. The concept of survivability and event maps is introduced to obtain suitable conditions that can be targeted such that the spacecraft impacts, or is weakly captured by, the Moon. Weak capture at the Moon is studied by method of these maps. Some results for planar Lyapunov orbits at \(L_{1}\) and \(L_{2}\) are given, as well as some results for the operational orbit of SOHO.     

Evolution of the \mathcal {L}_1 halo family in the radial solar sail circular restricted three-body problem

We present a detailed investigation of the dramatic changes that occur in the \(\mathcal {L}_1\) halo family when radiation pressure is introduced into the Sun–Earth circular restricted three-body problem (CRTBP). This photo-gravitational CRTBP can be used to model the motion of a solar sail orientated perpendicular to the Sun-line. The problem is then parameterized by the sail lightness number, the ratio of solar radiation pressure acceleration to solar gravitational acceleration. Using boundary-value problem numerical continuation methods and the AUTO software package (Doedel et al. in Int J Bifurc Chaos 1:493–520, 1991) the families can be fully mapped out as the parameter \(\beta \) is increased. Interestingly, the emergence of a branch point in the retrograde satellite family around the Earth at \(\beta \approx 0.0387\) acts to split the halo family into two new families. As radiation pressure is further increased one of these new families subsequently merges with another non-planar family at \(\beta \approx 0.289\) , resulting in a third new family. The linear stability of the families changes rapidly at low values of \(\beta \) , with several small regions of neutral stability appearing and disappearing. By using existing methods within AUTO to continue branch points and period-doubling bifurcations, and deriving a new boundary-value problem formulation to continue the folds and Krein collisions, we track bifurcations and changes in the linear stability of the families in the parameter \(\beta \) and provide a comprehensive overview of the halo family in the presence of radiation pressure. The results demonstrate that even at small values of \(\beta \) there is significant difference to the classical CRTBP, providing opportunity for novel solar sail trajectories. Further, we also find that the branch points between families in the solar sail CRTBP provide a simple means of generating certain families in the classical case.     

Velocities and Linewidths in the Network and Cell Interiors of a Polar Coronal Hole Compared with the Quiet Sun

The relative Doppler velocities and linewidths in a polar coronal hole and the nearby quiet-Sun region have been obtained from the Solar and Heliospheric Observatory (SOHO)/Coronal Diagnostic Spectrometer (CDS) observations using emission lines originating at different heights in the solar atmosphere from the lower transition region (TR) to the low solar corona. The observed region is separated into the network and the cell interior, and the behavior of the above parameters were examined in the different regions. It has been found that the histograms of Doppler velocity and width are generally broader in the cell interior than in the network. The histograms of Doppler velocities of the network and cell interior do not show significant differences in most cases. However, in the case of the quiet Sun, the Doppler velocities of the cell interior are more blueshifted than those of the network for the lowermost line He?ii 304 Å, and an opposite behavior is seen for the uppermost line Mg?ix 368 Å. The linewidth histograms show that the network–cell difference is more prominent in the coronal hole. The network has a significantly larger linewidth than the cell interior for the lowermost TR line He?ii 304 Å for the quiet Sun. For the coronal hole, this is true for the three lower TR lines: He?ii 304 Å, O?iii 599 Å, and O?v 630 Å. We also obtained the correlations between the relative Doppler velocity and the width. A mild positive correlation is found for the lowermost transition-region line He?ii 304 Å, which decreases even more or become insignificant for the intermediate lines. For the low coronal line Mg?ix 368 Å, the correlation becomes strongly negative. This might be caused by standing waves or waves propagating from the lower to the upper solar atmosphere. The results may have implications for the generation of the fast solar wind and coronal heating.     

The solar niobium abundance

The solar Nb abundance is derived from five Nb i and ten Nb ii lines in the photospheric spectrum. Equivalent widths are obtained from measurements on spectra recorded at Kitt Peak National Observatory. Synthetic spectrum calculations gave abundances of 2.23 and 2.08 from neutral and ionized lines respectively in the logarithmic A _H = 12.00 scale. This gives an average abundance value of A _Nb = 2.13 ± 0.10.     

Effects of Geometries and Substructures of ICMEs on Geomagnetic Storms

To better understand geomagnetic storm generations by ICMEs, we consider the effect of substructures (magnetic cloud, MC, and sheath) and geometries (impact location of flux-rope at the Earth) of the ICMEs. We apply the toroidal magnetic flux-rope model to 59 CDAW CME–ICME pairs to identify their substructures and geometries, and select 20 MC-associated and five sheath-associated storm events. We investigate the relationship between the storm strength indicated by minimum Dst index \((\mathrm{Dst}_{\mathrm{min}})\) and solar wind conditions related to a southward magnetic field. We find that all slopes of linear regression lines for sheath-storm events are steeper (\({\geq}\,1.4\)) than those of the MC-storm events in the relationship between \(\mathrm{Dst}_{\mathrm{min}}\) and solar wind conditions, implying that the efficiency of sheath for the process of geomagnetic storm generations is higher than that of MC. These results suggest that different general solar wind conditions (sheaths have a higher density, dynamic and thermal pressures with a higher fluctuation of the parameters and higher magnetic fields than MCs) have different impact on storm generation. Regarding the geometric encounter of ICMEs, 100% (2/2) of major storms (\(\mathrm{Dst}_{\mathrm{min}} \leq -100~\mbox{nT}\)) occur in the regions at negative \(P_{Y}\) (relative position of the Earth trajectory from the ICME axis in the \(Y\) component of the GSE coordinate) when the eastern flanks of ICMEs encounter the Earth. We find similar statistical trends in solar wind conditions, suggesting that the dependence of geomagnetic storms on 3D ICME–Earth impact geometries is caused by asymmetric distributions of the geoeffective solar wind conditions. For western flank events, 80% (4/5) of the major storms occur in positive \(P_{Y}\) regions, while intense geoeffective solar wind conditions are not located in the positive \(P_{Y}\). These results suggest that the strength of geomagnetic storms depends on ICME–Earth impact geometries as they determine the solar wind conditions at Earth.     

Sun-as-a-Star Observation of Flares in Lyman α by the PROBA2/LYRA Radiometer

There are very few reports of flare signatures in the solar irradiance at H i Lyman α at 121.5 nm, i.e. the strongest line of the solar spectrum. The LYRA radiometer onboard PROBA2 has observed several flares for which unambiguous signatures have been found in its Lyman-α channel. Here we present a brief overview of these observations followed by a detailed study of one of them: the M2 flare that occurred on 8 February 2010. For this flare, the flux in the LYRA Lyman-α channel increased by 0.6 %, which represents about twice the energy radiated in the GOES soft X-ray channel and is comparable with the energy radiated in the He ii line at 30.4 nm. The Lyman-α emission represents only a minor part of the total radiated energy of this flare, for which a white-light continuum was detected. Additionally, we found that the Lyman-α flare profile follows the gradual phase but peaks before other wavelengths. This M2 flare was very localized and had a very brief impulsive phase, but more statistics are needed to determine if these factors influence the presence of a Lyman-α flare signal strong enough to appear in the solar irradiance.