Long-term variations in the magnetic fields of the Sun and the heliosphere: their origin, effects, and implications

Recent studies of the variation of geomagnetic activity over the past 140 years have quantified the "coronal source" magnetic flux F-s that leaves the solar atmosphere and enters the heliosphere and have shown that it has risen, on average, by an estimated 34% since 1963 and by 140% since 1900. This variation of open solar flux has been reproduced by Solanki et al. [2000] using a model which demonstrates how the open flux accumulates and decays, depending on the rate of flux emergence in active regions and on the length of the solar cycle. We here use a new technique to evaluate solar cycle length and find that it does vary in association with the rate of change of F-s in the way predicted. The long-term variation of the rate of flux emergence is found to be very similar in form to that in F-s, which may offer a potential explanation of why F-s appears to be a useful proxy for extrapolating solar total irradiance back in time. We also find that most of the variation of cosmic ray fluxes incident on Earth is explained by the strength of the heliospheric field (quantified by F-s) and use observations of the abundance of the isotope Be-10 (produced by cosmic rays and deposited in ice sheets) to study the decrease in F-s during the Maunder minimum. The interior motions at the base of the convection zone, where the solar dynamo is probably located, have recently been revealed using the helioseismology technique and found to exhibit a 1.3-year oscillation. This periodicity is here reported in observations of the interplanetary magnetic field and geomagnetic activity but is only present after 1940, When present, it shows a strong 22-year variation, peaking near the maximum of even-numbered sunspot cycles and showing minima at the peaks of odd-numbered cycles. We discuss the implications of these long-term solar and heliospheric variations for Earth's environment.

Long-term changes in the solar photosphere associated with changes in the coronal source flux

Using sunspot observations from Greenwich and Mount Wilson, we show that the latitudinal spread of sunspot groups has increased since 1874, in a manner that closely mirrors the long-term (similar to 100 year) changes in the coronal source flux, F-s, as inferred from geomagnetic activity. This latitude spread is shown to be well correlated with the flux emergence rate required by the model of the coronal source flux variation by Solanki er al. [2000]. The time constant for the decay of this open flux is found to be 3.6 +/-0.8 years. Using this value, and quantifying the photospheric flux emergence rate using the latitudinal spread of sunspot groups, the model reproduces the observed coronal source flux variation. The ratio of the 100-year drift to the solar cycle amplitude for the flux emergence rate is found to be half of the same ratio for F-s.

Long term changes in EUV and X-ray emissions from the solar corona and chromosphere as measured by the response of the Earth’s ionosphere during total solar eclipses from 1932 to 1999

Measurements of the ionospheric E region during total solar eclipses in the period 1932-1999 have been used to investigate the fraction of Extreme Ultra Violet and soft X-ray radiation, phi, that is emitted from the limb corona and chromosphere. The relative apparent sizes of the Moon and the Sun are different for each eclipse, and techniques are presented which correct the measurements and, therefore, allow direct comparisons between different eclipses. The results show that the fraction of ionising radiation emitted by the limb corona has a clear solar cycle variation and that the underlying trend shows this fraction has been increasing since 1932. Data from the SOHO spacecraft are used to study the effects of short-term variability and it is shown that the observed long-term rise in phi has a negligible probability of being a chance occurrence.

Long-term variations in the magnetic field of the Sun and possible implications for terrestrial climate

Recent studies of the variation of geomagnetic activity over the past 140 years have quantified the "coronal source" or "open" magnetic flux F-s that leaves the solar atmosphere and enters the heliosphere and have shown that it has risen, on average, by 34% since 1963 and by 140% since 1900. This variation is reflected in studies of the heliospheric field using isotopes deposited in ice sheets and meteorites by the action of galactic comic rays. The variation has also been reproduced using a model that demonstrates how the open flux accumulates and decays, depending on the rate of flux emergence in active regions and on the length of the solar cycle. The cosmic ray flux at energies > 3 GeV is found to have decayed by about 15% during the 20(th) century (and by about 4% at > 13 GeV). We show that the changes in the open flux do reflect changes in the photospheric and sub-surface field which offers an explanation of why open flux appears to be a good proxy for solar irradiance extrapolation. Correlations between F-s, solar cycle length, L, and 11-year smoothed sunspot number, R-11, explain why the various irradiance reconstructions for the last 150 years are similar in form. Possible implications of the inferred changes in cosmic ray flux and irradiance for global temperatures on Earth are discussed.

Solar causes of the long-term increase in geomagnetic activity

We analyze the causes of the century-long increase in geomagnetic activity, quantified by annual means of the aa index, using observations of interplanetary space, galactic cosmic rays, the ionosphere, and the auroral electrojet, made during the last three solar cycles. The effects of changes in ionospheric conductivity, the Earth's dipole tilt, and magnetic moment are shown to be small; only changes in near-Earth interplanetary space make a significant contribution to the long-term increase in activity. We study the effects of the interplanetary medium by applying dimensional analysis to generate the optimum solar wind-magnetosphere energy coupling function, having an unprecedentedly high correlation coefficient of 0.97. Analysis of the terms of the coupling function shows that the largest contributions to the drift in activity over solar cycles 20-22 originate from rises in the average interplanetary magnetic field (IMF) strength, solar wind concentration, and speed; average IMF orientation has grown somewhat less propitious for causing geomagnetic activity. The combination of these factors explains almost all of the 39% rise in aa observed over the last three solar cycles. Whereas the IMF strength varies approximately in phase with sunspot numbers, neither its orientation nor the solar wind density shows any coherent solar cycle variation. The solar wind speed peaks strongly in the declining phase of even-numbered cycles and can be identified as the chief cause of the phase shift between the sunspot numbers and the aa index. The rise in the IMF magnitude, the largest single contributor to the drift in geomagnetic activity, is shown to be caused by a rise in the solar coronal magnetic field, consistent with a rise in the coronal source field, modeled from photospheric observations, and an observed decay in cosmic ray fluxes.

Long-term drift of the coronal source magnetic flux and the total solar irradiance

We test the method of Lockwood et al. [1999] for deriving the coronal source flux from the geomagnetic aa index and show it to be accurate to within 12% for annual means and 4.5% for averages over a sunspot cycle. Using data from four solar constant monitors during 1981-1995, we find a linear relationship between this magnetic flux and the total solar irradiance. From this correlation, we show that the 131% rise in the mean coronal source field over the interval 1901-1995 corresponds to a rise in the average total solar irradiance of {\Delta}I = 1.65 +/- 0.23 Wm^{-2}.

Long term individual load forecast under different electrical vehicles uptake scenarios

More and more households are purchasing electric vehicles (EVs), and this will continue as we move towards a low carbon future. There are various projections as to the rate of EV uptake, but all predict an increase over the next ten years. Charging these EVs will produce one of the biggest loads on the low voltage network. To manage the network, we must not only take into account the number of EVs taken up, but where on the network they are charging, and at what time. To simulate the impact on the network from high, medium and low EV uptake (as outlined by the UK government), we present an agent-based model. We initialise the model to assign an EV to a household based on either random distribution or social influences - that is, a neighbour of an EV owner is more likely to also purchase an EV. Additionally, we examine the effect of peak behaviour on the network when charging is at day-time, night-time, or a mix of both. The model is implemented on a neighbourhood in south-east England using smart meter data (half hourly electricity readings) and real life charging patterns from an EV trial. Our results indicate that social influence can increase the peak demand on a local level (street or feeder), meaning that medium EV uptake can create higher peak demand than currently expected.

Design of a computerised treatment for short-term memory deficits in aphasia

The treatment of auditory-verbal short-term memory (STM) deficits in aphasia is a growing avenue of research (Martin & Reilly, 2012; Murray, 2012). STM treatment requires time precision, which is suited to computerised delivery. We have designed software, which provides STM treatment for aphasia. The treatment is based on matching listening span tasks (Howard & Franklin, 1990), aiming to improve the temporal maintenance of multi-word sequences (Salis, 2012). The person listens to pairs of word-lists that differ in word-order and decides if the pairs are the same or different. This approach does not require speech output and is suitable for persons with aphasia who have limited or no output. We describe the software and how its review from clinicians shaped its design.

Memo test: standardisation of computerised assessment of auditory verbal short-term memory

Short-term memory (STM) impairments are prevalent in adults with acquired brain injuries. While there are several published tests to assess these impairments, the majority require speech production, e.g. digit span (Wechsler, 1987). This feature may make them unsuitable for people with aphasia and motor speech disorders because of word finding difficulties and speech demands respectively. If patients perceive the speech demands of the test to be high, the may not engage with testing. Furthermore, existing STM tests are mainly ‘pen-and-paper’ tests, which can jeopardise accuracy. To address these shortcomings, we designed and standardised a novel computerised test that does not require speech output and because of the computerised delivery it would enable clinicians identify STM impairments with greater precision than current tests. The matching listening span tasks, similar to the non-normed PALPA 13 (Kay, Lesser & Coltheart, 1992) is used to test short-term memory for serial order of spoken items. Sequences of digits are presented in pairs. The person hears the first sequence, followed by the second sequence and s/he decides whether the two sequences are the same or different. In the computerised test, the sequences are presented in live voice recordings on a portable computer through a software application (Molero Martin, Laird, Hwang & Salis 2013). We collected normative data from healthy older adults (N=22-24) using digits, real words (one- and two-syllables) and non-words (one- and two- syllables). Their performance was scored following two systems. The Highest Span system was the highest span length (e.g. 2-8) at which a participant correctly responded to over 7 out of 10 trials at the highest sequence length. Test re-test reliability was also tested in a subgroup of participants. The test will be available as free of charge for clinicians and researchers to use.

Nodal variations and long-term changes in the main tides on the coasts of China

The long-term changes in the main tidal constituents (O1, K1, M2, N2, and S2) along the coasts of China and in adjacent seas are investigated based on 17 tide-gauge records covering the period 1954–2012. The observed 18.61 year nodal modulations of the diurnal constituents O1 and K1 are in agreement with the equilibrium tidal theory, except in the South China Sea. The observed modulations of the M2 and N2 amplitudes are smaller than theoretically predicted at the northern stations and larger at the southern stations. The discrepancies between the theoretically predicted nodal variations and the observations are discussed. The 8.85 year perigean cycle is identifiable in the N2 parameters at most stations, except those in the South China Sea. The radiational component of S2 contributes on average 16% of the observed S2 except in the Gulf of Tonkin, on the south coast, where it accounts for up to 65%. We confirmed the existence of nodal modulation in S2, which is stronger on the north coast. The semidiurnal tidal parameters show significant secular trends in the Bohai and Yellow Seas, on the north coast, and in the Taiwan Strait. The largest increase is found for M2 for which the amplitude increases by 4–7 mm/yr in the Yellow Sea. The potential causes for the linear trends in tidal constants are discussed.

Simulation of long-term thermal characteristics of three Estonian lakes

A one-dimensional surface energy-balance lake model, coupled to a thermodynamic model of lake ice, is used to simulate variations in the temperature of and evaporation from three Estonian lakes: Karujärv, Viljandi and Kirjaku. The model is driven by daily climate data, derived by cubic-spline interpolation from monthly mean data, and was run for periods of 8 years (Kirjaku) up to 30 years (Viljandi). Simulated surface water temperature is in good agreement with observations: mean differences between simulated and observed temperatures are from −0.8°C to +0.1°C. The simulated duration of snow and ice cover is comparable with observed. However, the model generally underpredicts ice thickness and overpredicts snow depth. Sensitivity analyses suggest that the model results are robust across a wide range (0.1–2.0 m−1) of lake extinction coefficient: surface temperature differs by less than 0.5°C between extreme values of the extinction coefficient. The model results are more sensitive to snow and ice albedos. However, changing the snow (0.2–0.9) and ice (0.15–0.55) albedos within realistic ranges does not improve the simulations of snow depth and ice thickness. The underestimation of ice thickness is correlated with the overestimation of snow cover, since a thick snow layer insulates the ice and limits ice formation. The overestimation of snow cover results from the assumption that all the simulated winter precipitation occurs as snow, a direct consequence of using daily climate data derived by interpolation from mean monthly data.