Centennial variations in sunspot number, open solar flux and streamer belt width: 1. correction of the sunspot number record since 1874

We analyse the widely-used international/ Zürich sunspot number record, R, with a view to quantifying a suspected calibration discontinuity around 1945 (which has been termed the “Waldmeier discontinuity” [Svalgaard, 2011]). We compare R against the composite sunspot group data from the Royal Greenwich Observatory (RGO) network and the Solar Optical Observing Network (SOON), using both the number of sunspot groups, N{sub}G{\sub}, and the total area of the sunspots, A{sub}G{\sub}. In addition, we compare R with the recently developed interdiurnal variability geomagnetic indices IDV and IDV(1d). In all four cases, linearity of the relationship with R is not assumed and care is taken to ensure that the relationship of each with R is the same before and after the putative calibration change. It is shown the probability that a correction is not needed is of order 10{sup}−8{\sup} and that R is indeed too low before 1945. The optimum correction to R for values before 1945 is found to be 11.6%, 11.7%, 10.3% and 7.9% using A{sub}G{\sub}, N{sub)G{\sub}, IDV, and IDV(1d), respectively. The optimum value obtained by combining the sunspot group data is 11.6% with an uncertainty range 8.1-14.8% at the 2σ level. The geomagnetic indices provide an independent yet less stringent test but do give values that fall within the 2σ uncertainty band with optimum values are slightly lower than from the sunspot group data. The probability of the correction needed being as large as 20%, as advocated by Svalgaard [2011], is shown to be 1.6 × 10{sup}−5{\sup}.

Centennial variations in sunspot number, open solar flux and streamer belt width: 2. Comparison with the geomagnetic data

We investigate the relationship between interdiurnal variation geomagnetic activity indices, IDV and IDV(1d), corrected sunspot number, R{sub}C{\sub}, and the group sunspot number R{sub}G{\sub}. R{sub}C{\sub} uses corrections for both the “Waldmeier discontinuity”, as derived in Paper 1 [Lockwood et al., 2014c], and the “Wolf discontinuity” revealed by Leussu et al. [2013]. We show that the simple correlation of the geomagnetic indices with R{sub}C{\sub}{sup}n{\sup} or R{sub}G{\sub}{sup}n{\sup} masks a considerable solar cycle variation. Using IDV(1d) or IDV to predict or evaluate the sunspot numbers, the errors are almost halved by allowing for the fact that the relationship varies over the solar cycle. The results indicate that differences between R{sub}C{\sub} and R{sub}G{\sub} have a variety of causes and are highly unlikely to be attributable to errors in either R{sub}G{\sub} alone, as has recently been assumed. Because it is not known if R{sub}C{\sub} or R{sub}G{\sub} is a better predictor of open flux emergence before 1874, a simple sunspot number composite is suggested which, like R{sub}G{\sub}, enables modelling of the open solar flux for 1610 onwards in Paper 3, but maintains the characteristics of R{sub}C{\sub}.

Centennial variations in sunspot number, open solar flux and streamer belt width: 3. Modelling

From the variation of near-Earth interplanetary conditions, reconstructed for the mid-19th century to the present day using historic geomagnetic activity observations, Lockwood and Owens [2014] have suggested that Earth remains within a broadened streamer belt during solar cycles when the Open Solar Flux (OSF) is low. From this they propose that the Earth was immersed in almost constant slow solar wind during the Maunder minimum (c. 1650-1710). In this paper, we extend continuity modelling of the OSF to predict the streamer belt width using both group sunspot numbers and corrected international sunspot numbers to quantify the emergence rate of new OSF. The results support the idea that the solar wind at Earth was persistently slow during the Maunder minimum because the streamer belt was broad.

Long-term variations in the open solar flux and possible links to Earth’s climate

Recent paleoclimate studies provide strong evidence for an association between cosmogenic isotope production and Earth’s climate throughout the holecene. These isotopes are generated by the bombardment of Earth’s atmosphere by galactic cosmic rays, the fluxes of which vary in approximately inverse proportion to the total open magnetic flux of the Sun. This paper discusses how results from the Ulysses spacecraft allow us to quantify the open solar flux from observations of near-Earth interplanetary space and to study its long-term variations using the homogeneous record of geomagnetic activity. A study of the results and of their accuracy is presented. The two proposed mechanisms that could lead to the open solar flux being a good proxy for solar-induced climate change are discussed: the first is the modulation of the production of some types of cloud by the air ions produced by cosmic rays; the second is a variation in the total or spectral solar irradiance, in association with changes in the open flux. Some implications for our understanding of anthropogenic climate change are discussed.

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 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.

Ionospheric convection response to slow, strong variations in a northward interplanetary magnetic field: A case study for January 14, 1988

We analyze ionospheric convection patterns over the polar regions during the passage of an interplanetary magnetic cloud on January 14, 1988, when the interplanetary magnetic field (IMF) rotated slowly in direction and had a large amplitude. Using the assimilative mapping of ionospheric electrodynamics (AMIE) procedure, we combine simultaneous observations of ionospheric drifts and magnetic perturbations from many different instruments into consistent patterns of high-latitude electrodynamics, focusing on the period of northward IMF. By combining satellite data with ground-based observations, we have generated one of the most comprehensive data sets yet assembled and used it to produce convection maps for both hemispheres. We present evidence that a lobe convection cell was embedded within normal merging convection during a period when the IMF By and Bz components were large and positive. As the IMF became predominantly northward, a strong reversed convection pattern (afternoon-to-morning potential drop of around 100 kV) appeared in the southern (summer) polar cap, while convection in the northern (winter) hemisphere became weak and disordered with a dawn-to-dusk potential drop of the order of 30 kV. These patterns persisted for about 3 hours, until the IMF rotated significantly toward the west. We interpret this behavior in terms of a recently proposed merging model for northward IMF under solstice conditions, for which lobe field lines from the hemisphere tilted toward the Sun (summer hemisphere) drape over the dayside magnetosphere, producing reverse convection in the summer hemisphere and impeding direct contact between the solar wind and field lines connected to the winter polar cap. The positive IMF Bx component present at this time could have contributed to the observed hemispheric asymmetry. Reverse convection in the summer hemisphere broke down rapidly after the ratio |By/Bz| exceeded unity, while convection in the winter hemisphere strengthened. A dominant dawn-to-dusk potential drop was established in both hemispheres when the magnitude of By exceeded that of Bz, with potential drops of the order of 100 kV, even while Bz remained northward. The later transition to southward Bz produced a gradual intensification of the convection, but a greater qualitative change occurred at the transition through |By/Bz| = 1 than at the transition through Bz = 0. The various convection patterns we derive under northward IMF conditions illustrate all possibilities previously discussed in the literature: nearly single-cell and multicell, distorted and symmetric, ordered and unordered, and sunward and antisunward.

Holocene variations in peatland methane cycling associated with the Asian summer monsoon system

Atmospheric methane concentrations decreased during the early to middle Holocene; however, the governing mechanisms remain controversial. Although it has been suggested that the mid-Holocene minimum methane emissions are associated with hydrological change, direct evidence is lacking. Here we report a new independent approach, linking hydrological change in peat sediments from the Tibetan Plateau to changes in archaeal diether concentrations and diploptene delta C-13 values as tracers for methanogenesis and methanotrophy, respectively. A minimum in inferred methanogenesis occurred during the mid-Holocene, which, locally, corresponds with the driest conditions of the Holocene, reflecting a minimum in Asian monsoon precipitation. The close coupling between precipitation and methanogenesis is validated by climate simulations, which also suggest a regionally widespread impact. Importantly, the minimum in methanogenesis is associated with a maximum in methanotrophy. Therefore, methane emissions in the Tibetan Plateau region were apparently lower during the mid-Holocene and partially controlled by interactions of large-scale atmospheric circulation.

Palaeoenvironmental changes during the Valanginian: new insights from variations in phosphorus contents and bulk- and clay mineralogies in the western Tethys

Paleoenvironmental and paleoclimatic changes during the Valanginian carbon isotopic excursion (CIE) have been investigated in the western Tethys. For this purpose, bulk-rock and clay mineralogies, as well as phosphorus (P) contents were evaluated in a selection of five sections located in the Vocontian Basin (Angles, SE France; Alvier, E Switzerland; Malleval, E France), and the Lombardian Basin (Capriolo, N Italy; Breggia, S Switzerland). Within the CIE interval, bulk-rock and clay mineralogies are inferred to reflect mostly climate change. The onset of the CIE (Busnardoites campylotoxus ammonite Zone) is characterized by higher detrital index (DI: sum of the detrital minerals divided by calcite contents) values and the presence of kaolinite in their clay-mineral assemblages. In the late Valanginian (from the Saynoceras verrucosum Zone up to the end of the Valanginian), the samples show relatively variable DI and lower values or the absence of kaolinite. The variation in the mineralogical composition is interpreted as reflecting a change from a climate characterized by optimal weathering conditions associated with an increase in terrigenous input on the southern European margin during the CIE towards an overall unstable climate associated with drier conditions in the late Valanginian. This is contrasted by a dissymmetry (proximal vs distal) along the studied transect, the northern Tethyan margin being more sensitive to changes in continental input compared to the distal environments. P accumulation rates (PAR) present similar features. In the Vocontian basin, P content variations are associated with changes in terrigenous influx, whereas in the Lombardian basin (i.e. Capriolo and Breggia), PAR values are less well correlated. This is mainly because the deeper part of the Tethys was less sensitive to changes in continental inputs. The onset of the CIE (top of the B. campylotoxus Zone) records a general increase in PAR suggesting an increase in marine nutrient levels. This is linked to higher continental weathering rates and the enhanced influx of nutrients into the ocean. In the period corresponding to the shift itself, P contents show a dissymmetry between the Vocontian and Lombardian basins (proximal vs distal). For the sections of Malleval, Alvier and Angles, a decrease in P concentrations associated to a decrease in detrital input is observed. In Capriolo and Breggia, PAR show maximum values during the plateau, indicating a more complex interaction between different P sources. The time interval including the top of S. verrucosum Zone up to the end of the Valanginian is characterized by variable PAR values, suggesting variable nutrient influxes. These changes are in agreement with an evolution towards seasonally contrasted conditions in the late Valanginian.

Time-variations in commodity price jumps

In this paper, we study jumps in commodity prices. Unlike assumed in existing models of commodity price dynamics, a simple analysis of the data reveals that the probability of tail events is not constant but depends on the time of the year, i.e. exhibits seasonality. We propose a stochastic volatility jump–diffusion model to capture this seasonal variation. Applying the Markov Chain Monte Carlo (MCMC) methodology, we estimate our model using 20 years of futures data from four different commodity markets. We find strong statistical evidence to suggest that our model with seasonal jump intensity outperforms models featuring a constant jump intensity. To demonstrate the practical relevance of our findings, we show that our model typically improves Value-at-Risk (VaR) forecasts.

Exploring the effects of genotypical and phenotypical variations in bitter taste sensitivity on perception, liking and intake of brassica vegetables in the UK

Brassicaceous vegetables (BV) have chemoprotective effects and yet consumption of BV in the UK is low. Previous studies suggest perception, liking and intake of BV are influenced by bitter taste sensitivity which this study further explores. Phenotypical taste sensitivity of 136 subjects was classified using propythiouracil (PROP) and sodium chloride and fungiform papillae density (FPD) was measured from tongue images. Polymorphisms of TAS2R38 and gustin (CA6) genes were analysed. Liking and bitterness of four raw vegetables (two BV (broccoli and white cabbage) and two non-BV (spinach and courgette)), as well as habitual consumption, were evaluated. There was a significant association between TAS2R38 genotype and PROP taster status (p<0.0001) and between FPD and PROP taster status (p=0.029). Individuals with greater sensitivity for PROP predominantly had TAS2R38 PAV/PAV genotype and greater FPD. BV were perceived as more bitter than non-BV (p<0.0001) with PAV/PAV subjects perceiving significantly stronger bitter intensity. There was a significant difference in liking for the four vegetables (p=0.002), and between consumers of different TAS2R38 genotype (p=0.0024). Individuals with TAS2R38 AVI/AVI genotype liked BV more. Regarding intake, both PAV/PAV and AVI/AVI individuals consumed more total vegetables and BV than PAV/AVI. Although PROP nontasters tended to consume more vegetables and BV than the other two phenotype groups, liking and vegetable intake were not significantly affected by taste phenotype. Although there was not a significant effect of CA6 genotype on bitterness ratings, there was a significant interaction between CA6 and TAS2R38, and in addition CA6 genotype was significantly associated with BV intake. However, these effects require validation as the proportions of the population with the CA6 G/G genotype was extremely small (7%). Our results confirmed that bitter taste perception in vegetables was influenced by both genotype and phenotype of bitter taste sensitivity. Moreover, our findings demonstrated that neither genotype nor phenotype of taste sensitivity alone accurately predict vegetable liking and intake as demographic factors were found to have a substantial influence.

Urban warming in villages

Long term meteorological records (> 100 years) from stations associated with villages are generally classified as rural and assumed to have no urban influence. Using networks installed in two European villages, the local and microclimatic variations around two of these rural-village sites are examined. An annual average temperature difference ($\Delta{T}$) of 0.6 and 0.4 K was observed between the built-up village area and the current meteorological station in Geisenheim (Germany) and Haparanda (Sweden), respectively. Considerably larger values were recorded for the minimum temperatures and during summer. The spatial variations in temperature within the villages are of the same order as recorded over the past 100+ years in these villages (0.06 to 0.17 K/10 years). This suggests that the potential biases in the long records of rural-villages also warrant careful consideration like those of the more commonly studied large urban areas effects.