While graphical records of astronomical/meteorological events before telescopic observations are of particular interest, they have frequently undergone multiple copying and may have been modified from the original. Here, we analyze a graphical record of the cross-sign of 806 CE in the Anglo-Saxon Chronicle, which has been considered one of the earliest datable halo drawings in British records, whereas another cross-sign in 776 CE has been associated with the aurora. However, philological studies have revealed the later 806 event is derived from Continental annals. Here, records and drawings for the 806 event have been philologically traced back to mid-9th Century Continental manuscripts and the probable observational site identified as the area of Sens in northern France. The possible lunar halos at that time have been comprehensively examined by numerical ray tracing. Combined with calculations of twilight sky brightness, they identify a visibility window supporting monastic observation. Cruciform halos are shown to be fainter and rarer than brighter and more commonplace lunar halos. Physically credible cloud ice crystal variations can reproduce all the manuscript renditions. The manuscript records prove less than desirable detail but what is presented is fully consistent with a lunar halo interpretation. Finally, the possible societal impacts of such celestial events have been mentioned in the context of contemporary coins in Anglo-Saxon England and the Carolingian Empire. These analyses show that we need to trace their provenance back as far as possible, to best reconstruct the original event, even if graphical records are available for given astronomical/meteorological events.

The Carrington event is considered to be one of the most extreme space weather events in observational history within a series of magnetic storms caused by extreme interplanetary coronal mass ejections (ICMEs) from a large and complex active region (AR) emerged on the solar disk. In this article, we study the temporal and spatial evolutions of the source sunspot active region and visual aurorae, and compare this storm with other extreme space weather events on the basis of their spatial evolution. Sunspot drawings by Schwabe, Secchi, and Carrington describe the position and morphology of the source AR at that time. Visual auroral reports from the Russian Empire, Iberia, Ireland, Oceania, and Japan fill the spatial gap of auroral visibility and revise the time series of auroral visibility in mid to low magnetic latitudes (MLATs). The reconstructed time series is compared with magnetic measurements and shows the correspondence between low to mid latitude aurorae and the phase of magnetic storms. The spatial evolution of the auroral oval is compared with those of other extreme space weather events in 1872, 1909, 1921, and 1989 as well as their storm intensity, and contextualizes the Carrington event, as one of the most extreme space weather events, but likely not unique.

Large sunspots can be observed with the unaided eye under suitable atmospheric seeing conditions. Such observations are of particular value because the frequency of their appearance provides an approximate indication of the prevailing level of solar activity. Unaided-eye sunspot (UES) observations can be traced back well before the start of telescopic observations of the Sun, especially in the East Asian historical records. It is therefore important to compare more modern, UES observations with the results of telescopic sunspot observations, to gain a better understanding of the nature of the UES records. A previous comparison of Chinese UES records and Greenwich photo-heliographic results between 1874 and 1918 indicated that a few of the UES were apparently not supported by direct photographic evidence of at least one sunspot with a large area. This article reveals that one of the Chinese unaided-eye observations had possibly captured the transit of Venus on 9 December 1874. The Chinese sunspot records on this date are compared with Western sunspot observations on the same day. It is concluded that sunspots on the solar disk were quite small and the transit of Venus was probably misinterpreted as a sunspot by the Chinese local scholars. This case indicates that sunspots or comparable "obscuring" objects with an area as large as 1000 millionths of the solar disk could easily have been seen with the unaided eye under suitable seeing conditions. It also confirms the visibility of sunspots near the solar limb with the unaided eye. This study provides an explanation of the apparent discrepancy between the Chinese UES observation on 9 December 1874 and the Western sunspot observations using telescopes, as well as a basis for further discussion on the negative pairs in 1900 and 1911, apparently without sufficiently large area.

The enhancement of carbon-14 in tree rings around AD 774/775 has generated wide interest in solar activity at that time. The historical auroral records have been examined critically. Of particular interest was the "white vapour" observed in China on AD 776 January 12/13. Both Usoskin et al. (2013, Astron. Astrophys. 55, L3; U13) and Stephenson (2015, Adv. Sp. Res. 55, 1537; S15) interpreted this record as an auroral display. Subsequently, Neuh\"auser and Neuh\"auser (2015, Astron. Nachr. 336, 225; NN15) proposed five "criteria" for the likeliness of aurorae and on this basis rejected an auroral interpretation. Instead, they interpreted it as a lunar halo, and suggested there were no auroral records as a proxy of solar activity in the interval AD 774-785. We consider if their "lunar halo hypothesis" and their auroral criteria could be of use in future researches on historical auroral candidates. We first show a counter-example for the lunar halo hypothesis from a parallel record on 1882 November 17, which was seen as a whitish colour, in the southerly direction, and near the Moon. We then consider NN15's criteria on colour, direction, and sky brightness and investigate other counter-examples from early-modern auroral observations. We also consider the extension of the white vapour in AD 776 according to the distribution of Chinese asterisms, and show that its large extension was inconsistent with the lunar halo hypothesis. Conversely, the streaks of white vapour penetrating the eight Chinese asterisms can be reproduced if we consider auroral-ray structures at altitudes between 97 km and 170 km, along geomagnetic field lines between the L-shells L=1.55 and 1.64. Our investigations show that we should consider candidate auroral records in historical documents not on the basis of the newly suggested a priori criteria by NN15 but on all the available observational evidence.

The anomalous concentration of radiocarbon in 774/775 attracted intense discussion on its origin, including the possible extreme solar event(s) exceeding any events in observational history. Anticipating such extreme solar events, auroral records were also surveyed in historical documents and those including the red celestial sign after sunset in the Anglo-Saxon Chronicle (ASC) were subjected to consideration. Usoskin et al. (2013: U13) interpreted this record as an aurora and suggested enhanced solar activity around 774/775. Conversely, Neuhauser and Neuhauser (2015a, 2015b: N15a and N15b) interpreted "after sunset" as during sunset or twilight; they considered this sign as a halo display and suggested a solar minimum around 774. However, so far these records have not been discussed in comparison with eyewitness auroral records during the known extreme space-weather events, although they were discussed in relationship with potential extreme events in 774/775. Therefore, we reconstruct the observational details based on the original records in the ASC and philological references, compare them with eyewitness auroral observations during known extreme space-weather events, and consider contemporary solar activity. We clarify the observation was indeed "after sunset", reject the solar halo hypothesis, define the observational time span between 25 Mar. 775 and 25 Dec. 777, and note the parallel halo drawing in 806 in the ASC shown in N15b was not based on the original observation in England. We show examples of eyewitness auroral observations during twilight in known space-weather events, and this celestial sign does not contradict the observational evidence. Accordingly, we consider this event happened after the onset of the event in 774/775, but shows relatively enhanced solar activity, with other historical auroral records in the mid-770s, as also confirmed by the Be data from ice cores.

The study of historical great geomagnetic storms is crucial for assessing the possible risks to the technological infrastructure of a modern society, caused by extreme space-weather events. The normal benchmark has been the great geomagnetic storm of September 1859, the so-called "Carrington Event". However, there are numerous records of another great geomagnetic storm in February 1872. This storm, about 12 years after the Carrington Event, resulted in comparable magnetic disturbances and auroral displays over large areas of the Earth. We have revisited this great geomagnetic storm in terms of the auroral and sunspot records in the historical documents from East Asia. In particular, we have surveyed the auroral records from East Asia and estimated the equatorward boundary of the auroral oval to be near 24.3 deg invariant latitude (ILAT), on the basis that the aurora was seen near the zenith at Shanghai (20 deg magnetic latitude, MLAT). These results confirm that this geomagnetic storm of February 1872 was as extreme as the Carrington Event, at least in terms of the equatorward motion of the auroral oval. Indeed, our results support the interpretation of the simultaneous auroral observations made at Bombay (10 deg MLAT). The East Asian auroral records have indicated extreme brightness, suggesting unusual precipitation of high-intensity, low-energy electrons during this geomagnetic storm. We have compared the duration of the East Asian auroral displays with magnetic observations in Bombay and found that the auroral displays occurred in the initial phase, main phase, and early recovery phase of the magnetic storm.

More than 70 years ago it was recognised that ionospheric F2-layer critical frequencies $foF2$ had a strong relationship to sunspot number. Using historic datasets from the Slough and Washington ionosondes, we evaluate the best statistical fits of $foF2$ to sunspot numbers (at each Universal Time [UT] separately) in order to search for drifts and abrupt changes in the fit residuals over Solar Cycles 17 - 21. Polynomial fits are made both with and without allowance for the white-light facular area, which has been reported as being associated with cycle-to-cycle changes in the sunspot number - $foF2$ relationship. Over the interval studied here, the ISN, $R$, the backbone group number $Rbb$, and the corrected number $Rc$ largely differ in their allowance for the 'Waldmeier discontinuity' around 1945 (the correction factor for which for $R$, $Rbb$ and $Rc$ is, respectively, zero, effectively over 20%, and explicitly 11.6%). It is shown that for Solar Cycles 18 - 21, all three sunspot data sequences perform well, but that the fit residuals are lowest and most uniform for $Rbb$. We here use $foF2$ for those UTs for which $R$, $Rbb$, and $Rc$ all give correlations exceeding 0.99 for intervals both before and after the Waldmeier discontinuity. The error introduced by the Waldmeier discontinuity causes $R$ to underestimate the fitted values based on the $foF2$ data for 1932 - 1945 but $Rbb$ overestimates them by almost the same factor, implying that the correction for the Waldmeier discontinuity inherent in $Rbb$ is too large by a factor of two. Fit residuals are smallest and most uniform for $Rc$ and the ionospheric data support the optimum discontinuity multiplicative correction factor derived from the independent Royal Greenwich Observatory (RGO) sunspot group data for the same interval.

Long-lived (>20 days) sunspot groups extracted from the Greenwich Photo- heliographic Results (GPR) are examined for evidence of decadal change. The problem of identifying sunspot groups which are observed on consecutive solar rotations (recurrent sunspot groups) is tackled by first constructing manually an example dataset of recurrent sunspot groups and then using machine learning to generalise this subset to the whole GPR. The resulting dataset of recurrent sunspot groups is verified against previous work by A. Maunder and other Royal Greenwich Observatory (RGO) compilers. Recurrent groups are found to exhibit a slightly larger value for the Gnevyshev-Waldmeier Relationship than the value found by Petrovay and van Driel-Gesztelyi (Solar Phys. 51, 25, 1997), who used recurrence data from the Debrecen Photoheliographic Results. Evidence for sunspot group lifetime change over the previous century is observed within recurrent groups. A lifetime increase of 1.4 between 1915 and 1940 is found, which closely agrees with results from Blanter et al. (Solar Phys. 237, 329, 2006). Furthermore, this increase is found to exist over a longer period (1915 to 1950) than previously thought and provisional evidence is found for a decline between 1950 and 1965. Possible applications of machine-learning procedures to the analysis of historical sunspot observations, the determination of the magnetic topology of the solar corona and the incidence of severe space-weather events are outlined briefly.

We use Swift/BAT Earth occultation data at different geomagnetic latitudes to derive a sensitive measurement of the Cosmic X-ray background (CXB) and of the Earth albedo emission in the 15--200 keV band. We compare our CXB spectrum with recent (INTEGRAL, BeppoSAX) and past results (HEAO-1) and find good agreement. Using an independent measurement of the CXB spectrum we are able to confirm our results. This study shows that the BAT CXB spectrum has a normalization ~8(+/-3)% larger than the HEAO-1 measurement. The BAT accurate Earth albedo spectrum can be used to predict the level of photon background for satellites in low Earth and mid inclination orbits.

The spectrometer aboard INTEGRAL, SPI, has the capability to detect the signature of polarised emission from a bright gamma-ray source. GRB 041219a is the most intense burst localised by INTEGRAL and is an ideal candidate for such a study. Polarisation can be measured using multiple events scattered into adjacent detectors because the Compton scatter angle depends on the polarisation of the incoming photon. A search for linear polarisation in the most intense pulse of duration 66 seconds and in the brightest 12 seconds of GRB 041219a was performed in the 100-350keV, 100-500keV and 100keV-1MeV energy ranges. The multiple event data from the spectrometer was analysed and compared with the predicted instrument response obtained from Monte-Carlo simulations using the GEANT 4 INTEGRAL mass model. The chi^2 distribution between the real and simulated data as a function of the percentage polarisation and polarisation angle was calculated for all three energy ranges. The degree of linear polarisation in the brightest pulse of duration 66s was found to be 63+/-31% at an angle of 70+/-14 degrees in the 100-350keV energy range. The degree of polarisation was also constrained in the brightest 12s of the GRB and a polarisation fraction of 96+/-40% at an angle of 60+/-14 degrees was determined over the same energy range. However, despite extensive analysis and simulations, a systematic effect that could mimic the weak polarisation signal could not be definitively excluded. Our results over several energy ranges and time intervals are consistent with a polarisation signal of about 60% at a low level of significance (2 sigma). We conclude that the procedure described here demonstrates the effectiveness of using SPI as a polarimeter, and is a viable method of measuring polarisation levels in intense gamma--ray bursts.

Parameters from the literature, such as positions, photoelectric absorption (nH), spin and orbital periods, and distances or redshifts, were collected for all ~500 sources detected by INTEGRAL-ISGRI so far. We investigate where new and previously-known sources detected by ISGRI fit in the parameter space of high-energy objects, and we use the parameters to test correlations expected from theoretical predictions. For example, the influence of the local absorbing matter on periodic modulations is studied for Galactic High-Mass X-ray Binaries (HMXBs) with OB supergiant and Be companions. We examine the spatial distribution of different types of sources in the Milky Way using various projections of the Galactic plane, in order to highlight signatures of stellar evolution and to speculate on the origin of the group of sources whose classifications are still uncertain. ISGRI has detected similar numbers of X-ray Binaries and Active Galactic Nuclei (AGN). The former group contains new members of the class of HMXBs with supergiant stellar companions. Thanks to these additional systems, we are able to show that HMXBs are generally segregated in plots of intrinsic nH versus the orbital period of the system and versus the spin period of the pulsar, based on whether the companion is a Be or an OB supergiant star. We also find a tentative but expected anti-correlation between nH and the orbital period, and a possible and unexpected correlation between the nH and the spin period. While only a handful of new Low-Mass X-ray Binaries (LMXBs) have been discovered, there are many sources that remain unclassified and they appear to follow a spatial distribution typical of Galactic sources (especially LMXBs) rather than extragalactic sources.

Fast X-ray Transients (FXTs) are believed to be non-recurrent bright X-ray sources lasting less than a day and occuring at serendipitous positions, they can best detected and discovered by instruments having a sufficiently wide field of view and high sensitivity. The IBIS/ISGRI instrument onboard INTEGRAL is particularly suited to detect new or already known fast X-ray transient sources. We report on IBIS/ISGRI detection of newly discovered outbursts of three fast transient sources located at low Galactic latitude: SAX J1818.6-1703; IGR J16479-4514; IGR J17391-302/XTE J1739-302. The reported results confirm and strengthen the very fast transient nature of these sources, given that all their newly detected outbursts have a duration less than about 3 hours. Additionally, they provide the first evidence for a possible recurrent fast transient behaviour as all three sources were detected in outburst by ISGRI more than once during the last 2 years.

INTEGRAL is the first gamma-ray astronomy mission with a sufficient sensitivity and angular resolution combination appropriate to the detection and identification of considerable numbers of gamma-ray emitting sources. The large field of view enables INTEGRAL to survey the galactic plane on a regular (~weekly) basis as part of the core programme. The first source catalogue, based on the 1st year of core programme data has been completed and published (Bird et al., 2004). It contained 123 gamma-ray sources (24 HMXB, 54 LMXB, 28 ``unknown'', plus 17 others) - sufficient numbers for a reasonable statistical analysis of their global properties. The detection of previously unknown gamma-ray emitting sources generally exhibiting high intrinsic absorption, is intriguing. The substantial fraction of unclassified gamma-ray sources suggests they must constitute a significant family of objects. In this paper we review the global characteristics of the known galactic sources as well as the unclassified objects. We present Log(N)-Log(S) distributions, angular distributions, and for systems with reliable distance estimates the spatial distributions within the Galaxy and luminosity functions. For the unknown sources, this statistical analysis has shown that they are most likely to be HMXBs containing a highly magnetised neutron star. The lack of X-ray counterparts for these sources indicates a high degree of intrinsic obscuration.

We present the transient nature of INTEGRAL source IGR J17285-2922, identified from a single period of activity during an IBIS/ISGRI Galactic Centre Deep Exposure in September 2003. The source has a maximum detection significance of 14 sigma in the 20-100 keV energy range and exhibits a flux of 6.5 mCrab before it moves out of the ISGRI field of view. The source is visible to at least 150 keV and its spectrum can be fit with a power law slope of Gamma=2.1+/-0.17; a more physical model could not be fit due to poor statistics. Detected characteristics are consistent with the source being a Galactic low-mass X-ray binary harbouring a black hole or neutron star.

The true nature of the progenitor to GRBs remains elusive; one characteristic that would constrain our understanding of the GRB mechanism considerably is gamma-ray polarimetry measurements of the initial burst flux. We present a method that interprets the prompt GRB flux as it Compton scatters off the Earth's atmosphere, based on detailed modelling of both the Earth's atmosphere and the orbiting detectors. The BATSE mission aboard the \textit{CGRO} monitored the whole sky in the 20 keV - 1 MeV energy band continuously from April 1991 until June 2000. We present the BATSE Albedo Polarimetry System (BAPS), and show that GRB 930131 and GRB 960924 provide evidence of polarisation in their prompt flux that is consistent with degrees of polarisation of $\Pi>35$% and $\Pi>50$% respectively. While the evidence of polarisation is strong, the method is unable to strongly constrain the degree of polarisation beyond a systematics based estimation. Hence the implications on GRB theory are unclear, and further measurements essential.

The Burst and Transient Source Experiment (BATSE) aboard CGRO monitored the whole sky through 8 NaI(Tl) crystal scintillators sensitive in the 20 keV - 2 MeV energy band continuously from April 1991 until June 2000. Results are presented on the long term variabilty observed in the brightest Active Galactic Nuclei (AGN) present in the BATSE data archive. This was achieved through the application of the Earth Occultation technique to data flat-fielded using the BATSE Mass Model. Removal of the temporal background variations from the data should allow a more sensitive extraction of source parameters. Analysis of the general trends of the 9-year light curves are also presented.

The BATSE mission aboard CGRO observed the whole sky for 9 years in the 20 keV - 2 MeV energy band. Flat-fielding of the temporal variations in the background present in the data set has been accomplished through a GEANT3 Monte-Carlo simulation - the BATSE Mass Model (BAMM). The Earth Occultation technique (EOT) is used together with a maximum-likelihood imaging approach to construct all-sky maps with ~mCrab sensitivity. Additionally, a non-linear CLEAN algorithm is applied to the all-sky maps to remove large artefacts that are systematic of the EOT. The latest survey results produced through the application of this technique to a subset of the BATSE data set are presented.

The BATSE mission aboard CGRO monitored the whole sky in the 20 keV - 1 MeV energy band continuously from April 1991 until June 2000. Although BATSE had very poor intrinsic angular resolution, the data can be used to survey the entire soft gamma-ray sky with < 1 degree angular resolution and ~ 1 mCrab sensitivity by using the Earth occultation method. This method determines flux by measuring the step in the count rate profile in each BATSE detector as a source rises above or sets below the Earth's limb. A maximum likelihood imaging technique can then be used to build up all-sky maps from the images of the Earth's limb produced by occulting sources. However, since the Earth seen from BATSE has a radius of ~ 70 degrees, the limb images that intersect at the positions of bright point sources have a significant effect over the area of the all-sky map. A method for performing image cleaning on likelihood data has also been developed and is used to effectively remove artefacts from the all-sky maps. This paper describes the 'LIMBO' imaging technique and presents preliminary all-sky maps of 25 - 160 keV emission, the first to be made since the HEAO1-A4 mission of 1978-79.

Orbiting x-ray and gamma ray instruments are subject to large background count rates due to local particle fluxes in the space environment. The ability of an instrument to make calibrated measurements of the flux from a source of interest is highly dependent on accurately determining the background level. We present here a method of calculating the energy dependent background flux for any point in the complete data set recorded by the Burst and Transient Source Experiment (BATSE) in its nine year mission. The BATSE Mass Model (BAMM) uses a Monte Carlo mass modelling approach to produce a data base of the gamma ray background which is then filtered to simulate the background count rate with a 2.048 second time resolution. This method is able to reduce the variations in the background flux by a factor of 8 - 10, effectively `flat-fielding' the detector response. With flat-fielded BATSE data it should be possible to use the Earth occultation technique to produce a hard x-ray all sky survey to the 1-2 mCrab sensitivity limit. BAMM is also capable of estimating the contribution to the spectra measured from gamma ray sources due to the reprocessing of source photons in inactive material surrounding a gamma ray detector. Possible applications of this aspect of the model in the area of Gamma Ray Burst spectral analysis are discussed.

The BATSE mission aboard CGRO can be used to observe hard X-ray sources by using the Earth occultation method. This method relies on measuring a step in the count rate profile in each BATSE detector as a source rises above or sets below the Earth's limb. A major problem in determining the step sizes (and hence the flux) is in extracting the steps from the varying background. A technique for flat-fielding the response of gamma ray detectors has been developed at Southampton. The technique uses a dynamic Monte-Carlo model to simulate the dominant components of the gamma-ray background encountered by the experiment at any point in its orbit. A maximum likelihood imaging method is also being developed that will be used to make a sky survey with all 9 years of the BATSE Continuous data set in the 20 - 500 keV range. An all sky map of 25 - 35 keV emission has been made using 60 days of data and has a 3 sigma flux sensitivity of 14 mCrab.