The ionospheric response over the UK to major bombing raids during World War II

. The Earth’s ionosphere is subject to disturbance from above (via solar variability and space-weather effects) and from below (such as tectonic activity, thunderstorms and sudden stratospheric warmings). Identifying the relative contribution of these effects remains challenging, despite recent advances in spacecraft monitoring near-Earth space. Man-made explosions provide a quantiﬁable proxy for natural terrestrial sources, enabling their impact on ionospheric variability to be studied. In this paper, the contribution of ground-based disturbances to ionospheric variability is investigated by considering the response of the ionospheric F2 layer over Slough, UK, to 152 major bombing raids over Europe during World War II, using a superposed epoch analysis. The median response of the F2 layer is a signiﬁcant decrease in peak electron concentration ( ∼ 0 . 3 MHz decrease in fo F2). This response is consistent with wave energy heating the thermo-sphere, enhancing the (temperature-dependent) loss rate of O + ions. The analysis was repeated for a range of thresholds in both time of bombing before the (noon) ionospheric measurement and tonnage of bombs dropped per raid. It was found that signiﬁcant ( ∼ 2–3 σ ) deviations from the mean oc-curred for events occurring between approximately 3 and 7 h ahead of the noon ionospheric measurements and for raids using a minimum of between 100 and 800


Introduction
The source of much of the observed variability within the Earth's ionosphere remains poorly understood.In this study we examine unique ionospheric measurements made above Slough, UK, for the duration of the Second World War (WWII) in order to determine whether any of the observed variability could be attributed to the major bombing campaigns across Europe.
Production of ionisation in the Earth's upper atmosphere is predominantly through photo-ionisation by solar extreme ultra violet (EUV) and X-ray radiation, while loss rates are very sensitive to the temperature and composition of the neutral thermosphere.As a result, the long-term average behaviour of the ionosphere is closely tied to solar activity and is well understood.Transient space-weather phenomena such as coronal mass ejections, high-speed solar wind streams and energetic particle events can temporarily perturb the ambient ionospheric conditions by enhancing ionospheric production through impulsive brightening of the solar atmosphere (so-lar flares), enhancing ionospheric loss rates (through heating of the thermosphere, affecting neutral composition and loss rates), and through direct enhancement of ionisation by precipitation of energetic particles.While the details of such processes are still the subject of ongoing research, once again, the underlying physics is broadly understood.Despite this, there are further sources of ionospheric variability that remain unaccounted for and it has been suggested (e.g.Rishbeth and Müller-Wodarg, 2006) that the source of this variability is from the lower atmosphere.Sources such as earthquakes (e.g.Astafyeva et al., 2013, and references therein), thunderstorms (Davis and Johnston, 2005;Yu et al., 2015) and explosions (e.g.Rishbeth, 1991) have been cited as potential causes of ionospheric variability, with a variety of proposed mechanisms including pressure waves, gravity waves, infrasound and modulation of the global electric circuit.
There have been a number of case studies into the impact of terrestrial explosions on the upper atmosphere (e.g.Rishbeth, 1991), most notably the events surrounding the explosion at the Flixborough chemical plant in 1974 (Jones and Spracklen, 1974;Krasnova et al., 2003), while Pohotelov et al. (1991) considered the impact on the F2 and Es ionospheric layers during a 32-day period of bombing in the Gulf War.In this current paper we make use of historical records to identify large bombing raids over mainland Europe during WWII and, using a superposed epoch, or composite analysis (Chree, 1913), look for any consistent response in ionospheric measurements made at the Radio Research Station at Slough in the UK.Using historical records, reasonable quantitative estimates of the type and tonnage of explosives for each raid can be made, enabling the raids to be subdivided by size.

Early ionospheric measurements
The Radio Research Station (latterly the Radio and Space Research Station and ultimately the Appleton Laboratory) located at Ditton Park near Slough (Gardiner et al., 1982) conducted routine measurements of the Earth's ionosphere from 1933 to 1996.These measurements continue today at the Rutherford Appleton Laboratory near Chilton, UK.This sequence represents the longest continuous set of ionospheric measurements in the world.The technique used exploits the fact that a transmitted radio pulse is returned from an ionised atmospheric layer when the radio frequency of the pulse, f , matches the local plasma frequency.From this, the local electron concentration, N, can be determined via the relation f ≈ 9 √ N. By transmitting a sequence of radio pulses over a range of radio frequencies, it is therefore possible to construct a height profile of ionospheric electron concentration.Such measurements are most often presented in the form of an ionogram, a plot of "virtual height" (estimated from the time of flight of the signal assuming a vacuum) against radio frequency.From such records the virtual height and the peak frequency (and therefore electron concentra-tion) returned from each ionospheric layer can be determined and tabulated.While modern digital soundings automatically identify a comprehensive set of such parameters, this was a time-consuming task for the earliest analogue measurements, and so only the peak frequency of the ionospheric F layer (denoted foF2) was initially routinely scaled.While other parameters were scaled intermittently, foF2 represents the most comprehensive parameter scaled from these data that exists in a digital form.The original photographic prints of these early ionograms are held by the UK Solar System Data Centre, from which additional information can potentially be gleaned (Davis et al., 2013).

Bombing raids during World War II
Looking for a signature in the UK ionospheric records from Allied bombing campaigns over Europe between 1943 and 1945 may not seem like the most obvious of studies, but there are several reasons as to why a signature from such raids may be detectable over others.While the bombing of London by the Luftwaffe between September 1940 and May 1941 (popularly known as the "London Blitz") would have generated explosions at a closer proximity to the ionospheric measurements being made above Slough, this bombing was more or less continuous, making it difficult to separate the impact of wartime raids from those of natural seasonal variability.In addition, it is well documented that the Luftwaffe did not possess any four-engined long-range bombers (e.g.Beker, 1969).The mainstay of the Luftwaffe was the Heinkel 111, a twin-engined bomber capable of carrying 4400 lb (1667 kg) of bombs.Using external racks, the aircraft could carry 7900 lbs (3600 kg) of bombs, but the external racks blocked the internal bomb bay and significantly impaired the aircraft's performance (Regnat, 2004).
In contrast, the Allied airforces' use of four-engined bombers enabled them to carry much larger bombing loads, including individual high-explosive (HE) bombs of increasing capacity.A regular Avro Lancaster load designed for bombing of heavily industrialised cities (Mason, 1989) consisted of 1 × 4000 lb (1667 kg) amatol-filled ("Cookie"), 3 × 1000 lb (455 kg) minol-or tritonal-filled, impact-fused highcapacity (HC) bomb short-finned, short-delay, tail-armed HE bombs, and up to six additional compartments filled with 4 or 30 lb incendiary bombs.An alternative configuration used for the bombing of factories, railyards and dockyards consisted of 14 × 1000 lb (1667 kg) medium case (MC), general purpose (GP) short-finned HE bombs.The other mainstay of the RAF during this period was the Handley Page Halifax bomber which had a maximum bomb load of 13 000 lbs (5897 kg).Typical loads consisted of six 500 lb, six 1000 lb and two 2000 lb HE bombs or six 500 lb and four 2000 lb HE bombs.Unlike the Lancaster, the Halifax was not able to carry the 4000 lb or larger bombs.A third aircraft used by the RAF in combined or individual raids was the De Havilland Mosquito, which was capable of carrying a single 4000 lb bomb.The USAAF B-17 "Flying Fortress" bomber was able to carry a bomb load between 4500 lb (2000 kg) and 8000 lb (3600 kg) depending on the range of the mission.This aircraft too was capable of carrying the larger 4000 lb bombs.The use of additional aircraft, such as the Wellington (4500 lb bomb load) and Short Stirling (8000-14 000 lb bomb load) bombers, was gradually phased out during this period of the war.
Amatol, minol and tritonal represent various mixtures of tri-nitro-toluene (TNT) with aluminium or ammonium nitrate.Originally formulated to extend limited supplies of TNT, these mixtures provided similar or even enhanced explosive energy compared with TNT alone (Maienschein, 2002).Torpex was 50 % more powerful than TNT.One metric ton of TNT has an explosive energy of 4.184 × 10 9 J.
The "Cookie", used by the RAF, was the first "blockbuster" bomb.The minimum height at which an aircraft could safely drop the 4000 lb "Cookie" without being damaged by the resulting shockwave was considered to be 6000 feet (1800 m).Even so, there are anecdotal accounts of aircraft being damaged despite following this instruction (Nelmes and Jenkins, 2002).Later versions, designed to penetrate and destroy reinforced underground bunkers, were even bigger.These so-called "earthquake" bombs included the 12 000 lb (5450 kg) "Tall Boy" and the 22 000 lb (10 900 kg) "Grand Slam".
While the bombing of London continued during 1943-1945, it was mostly via V1 and V2 rockets which, while individually devastating, did not compare with the explosive energy associated with that of a heavy bombing raid.

Method
While there was much military activity throughout Europe and beyond during the Second World War, some of the largest bombing raids over Berlin and other European targets took place between 1943 and 1945.The dates of these heavy raids are listed in Table 1, along with an estimate of the combined weight of high explosives used.Where no reference to the total mass of HE bombs could be found for a particular raid, an estimate was made from the recorded number of aircraft, weighted by the mean ratio of HE to total bomb load from raids for which this information was known (0.667).These dates were selected to coincide with the latter part of the war, during which intensive bombing of London was less prevalent and the raids over mainland Europe were more intense (in terms of tonnage of high explosives dropped); the length of time spent bombing was usually much shorter; and more raids occurred during the day (Middlebrook and Everitt, 1985).The times listed in this table mostly came from Berlin air-raid records (Demps, 2004).These were augmented with timings gleaned from various eye-witness accounts (see Appendix A), in particular from the archives of 550 and 410 RAF squadrons whose records have been made available online.For aircraft records, the times of the first and last recorded bombings were used.For times taken from the times of air-raid alarms, the start and end of these warnings were used.Where no detail about the length of each raid was given, the start and end times are identical.The Berlin air-raid records were recorded in local time, while military records are most likely recorded in GMT.Where known, this time difference has been taken into account, but for some this may introduce an uncertainty of 1 h into the analysis.The start times from these raids were used as "trigger" times in a superposed epoch or Chree analysis (Chree, 1913).This type of analysis is useful in geophysics for detecting a faint signal in a noisy data sequence (e.g.Davis et al., 1997).By calculating the median response of a dataset around these trigger times, any small but repeatable signal is reinforced, while any signal not associated with the trigger events is expected to occur at random and therefore cancel out when averaged.Medians were used in order to minimise the influence of outliers in the analysis.
The ionospheric data used in this study contain a strong seasonal variation introduced by annual and solar cycle variations in ion production and loss.In order that the seasonal distribution did not dominate the signal in any superposed epoch analysis, a 30-day running median was subtracted from these data.Data from the resulting parameter, δfoF2, were then aligned according to the trigger times and combined in the superposed epoch study.Since it is not possible for the ionosphere to be influenced by a given raid prior to its occurrence, the study was further constrained to ensure that each trigger event was aligned with the first subsequent ionospheric data point within a given time window.The length of time between a raid and the subsequent ionospheric measurement, along with the minimum tonnage considered for a "large" raid, are subjective variables in this analysis.In order to address this, the analysis was repeated for a range of time windows, from 0 to ≤ 23 h, and the size of a bombing raid was defined by the minimum estimated total tonnage of HE used, from ≥ 100 to ≥ 1000 t in steps of 100 t.

Results
The results of a typical superposed epoch analysis using the estimated start time of each bombing raid are shown in Fig. 1.This analysis uses thresholds of ≤ 10 h of the ionospheric observation for raids using ≥ 700 t of HE for which there are 14 events.The median ionospheric response is shown as the black line in the top panel.Here the standard error in the median is shown as a grey shaded area around the line, while the dashed and dotted lines represent the 95th and 99th percentiles obtained by repeating the analysis 10 000 times using a similar number of random trigger times drawn from the same epoch for which no major bombing raid has been identified.It can be seen that the ionosphere is significantly weaker (1.9 standard deviations from the mean and around Table 1.Dates, times, locations and estimated tonnage of major bombing raids conducted over Europe between 1943 and 1945.Times for Berlin raids were taken from the duration of air-raid warnings (Demps, 2014).Other times were taken from a variety of eye-witness accounts (see Appendix A).Raids for which the total tonnage of HE bombs was estimated from the type and number of aircraft involved are marked with a * .Noon values of foF2 and 30-day running median foF2 values over Slough corresponding to time = 0 in the superposed epoch analysis are presented in columns 6 and 7. Of the 152 events considered in this study, there are 29 days for which there is currently no noon foF2 value available.the 99th percentile level when compared with the range calculated from random trigger events) at time = 0 (the day of the bombing raids).In order to ensure that geomagnetic disturbances did not contribute to the ionospheric response, all events for which Ap exceeded 48 (Kp > 5) were not consid-ered.This is reflected in the low median values in Ap presented in the lower panel.

Date
In order to investigate the sensitivity of this analysis to the arbitrary thresholds applied to the data, the analysis was repeated for a range of thresholds in both the length of time The grey shaded area represents the standard error in these median values, while the dashed and dotted lines represent the 95th and 99th percentiles of 10 000 repeated analyses using random control days in which no bombing raids were identified.At time = 0 (within 24 h of the trigger events) δfoF2 is depleted (1.9 standard deviations from the mean).Arbitrary threshold values were used, corresponding to > 700 t of high explosives per raids occurring within 10 h before the noon ionospheric measurements.
between the raid and the ionospheric measurement (from 0 to 23 h) and the size of a raid as defined by the minimum amount of explosives used (from 100 to 1000 t).The significance of each of these analyses was estimated by calculating the mean and standard deviation of all response times other than at time = 0 and calculating the number of standard deviations the time = 0 measurement was from this mean.The result of these analyses is shown in Fig. 2. The results are plot-ted on a grid with maximum time between raid and measurement along the y axis and minimum tonnage of the raid along the x axis.The significance of the response at time = 0 for each analysis is shown as a colour contour.As the thresholds used are varied, the number of raids included in each analysis will vary.This is reproduced in a similar grid also presented in Fig. 2, with the colour axis representing the number of events used in each analysis.Two results are clear from this for a range of thresholds.Events preceding the ionospheric measurements were considered for time windows from 0 to 23 h ahead of the noon ionospheric measurement and the minimum tonnage of HE bombs used in each raid was also varied from 100 to 1000 t.The most significant ionospheric response occurred for events occurring within 5 h before noon, in which a minimum of 300 t of HE bombs were dropped.Panel (b) presents the number of events contributing to each analysis.The significance of the response decreases if the threshold is extended beyond 5 h ahead of the ionospheric measurement.This indicates that events occurring at larger time offsets are not contributing to the observed median response.The significance of the result 5 h ahead decreases as the threshold on the minimum tonnage of HE bombs increases.This is likely due to the decreasing number of events contributing to each analysis.
analysis.Firstly a significant (> 2σ ) ionospheric response is obtained for a broad range of trigger thresholds (3 ≤ time window ≤ 7 h & 100 ≤ minimum tonnage ≤ 800 t).Secondly, the most significant result (2.99σ ) is seen for a time window of ≤ 5 h and a minimum tonnage per raid of 300 t.The results of this analysis are presented in Fig. 3, for which there are 20 events.Consideration of the number of events contributing to each study suggests that 14 or more events are required before the random noise is reduced to a level where a significant signal can be detected.The ionospheric responses for the analyses that contain the most events (with thresholds exceeding 15 h) are not significant, indicating that these events do not contribute to the response observed within 5 h of the ionospheric observations.

Discussion
From the analysis undertaken in this study, no minimum threshold in HE tonnage is resolved.In order to investigate this, data for smaller raids need to be included.In addition to the major raids considered here, there are many more smaller-scale raids involving fewer or smaller aircraft.For example, Mosquito aircraft were used in many hundreds of bombing raids throughout this period (Middlebrook and Everitt, 1985).Given the fact that the current list of major bombing raids used in this study is by no means comprehensive, it is likely that information about many hundreds of smaller raids would need to be included in order to tease out the signal of such raids from the background ionospheric variability in which further large raids were occurring.As such we consider this beyond the scope of the current study.
For the ionosphere at the altitude of the F2 region (∼ 200-300 km) above the UK to respond to bombing raids conducted at distances up to 1000 km away, the bombing must have generated pressure waves that were capable of propagating to ionospheric altitudes.A sound wave travelling this distance in the lower atmosphere would arrive within an hour.The speed of sound is temperature-dependent and the temperature decreases with altitude in the troposphere and mesosphere before increasing in the thermosphere.Since the thermosphere represents the most significant fraction of the vertical profile, it is likely that a sound wave propagating vertically as well as horizontally would arrive even sooner.One potential mechanism therefore is of a pressure wave propagating upwards in all directions.At higher altitudes its amplitude increases until it breaks in the upper atmosphere, depositing its energy as heat.A very rough estimate of the anticipated thermospheric temperature rise can be obtained by considering the specific heat capacity of the atmosphere which can be expressed as where Q is the energy input in joules (4.184 × 10 12 for 1000 metric tonnes of TNT), C p is the molar-specific heat capacity of N 2 (∼ 29.1 J mol −1 K −1 ), n is the number of moles of gas m −3 (at ionospheric altitudes, the number density of the atmosphere is ∼ 10 16 m −3 or 1.66×10 −8 moles m −3 ) and T is the change in temperature (K).Assuming the energy is equally distributed throughout a cylinder of atmosphere 1000 km in radius and 300 km in height, this gives a temperature rise of ∼ 9 K.
The dominant ion species in the mid-latitude ionospheric F region is O + , whose loss rate is temperature-dependent (Rees, 1989).However, the dominant mechanism by which O + ions are lost is through reaction with N 2 and O 2 molecules in the reactions The overall loss rate, β, for O + ions can therefore be expressed as where [N 2 ] and [O 2 ] are the concentrations of N 2 and O 2 molecules respectively and k 1 and k 2 are the rate coefficients for the two reactions.These rate coefficients are also temperature-dependent (Rees, 1989).The combined loss rate for O + ions is therefore dependent on both reaction rates and the concentration of thermospheric species.Müller-Wodarg et al. (1998) modelled the ionospheric and thermospheric response to localised thermospheric cooling (≤ 40 K) during a total solar eclipse.They predicted an 8 % increase in foF2 (∼ 0.2 MHz) due to the contraction of the atmosphere and an increase in the [O] / [N 2 ] ratio caused, in part, by a reduction in the concentration of N 2 .It is reasonable to assume that the atmospheric expansion due to energy from localised bombing raids would have an analogous, if opposite, effect on the ionosphere and thermosphere.A rise in the background thermospheric temperature would result in an enhanced loss rate, with the equilibrium between production and loss being established at a lower peak electron concentration.Such equilibrium would be reached within minutes of perturbation, well within the resolution of the ionospheric data.Grandin et al. (2015) studied the impact on foF2 of high-speed streams on Earth.They found that a thermospheric temperature increase of 20-50 K may result in a decrease in foF2 by 0.5-1.0MHz.If the bombing resulted in the generation of shock waves or atmospheric gravity waves, their horizontal propagation speed would need to be of the order of 300 km h −1 , while the vertical velocity component would need to be around 100 km h −1 in order to affect the atmosphere above Slough.There is evidence that turbulence generated in the lower thermosphere by space shuttle launches can propagate 1000 km horizontally within 8 h (Kelley et al., 2009).While this example was specific to the lower thermosphere at altitudes between 100 and 115 km, it nevertheless has a similar time constant to that observed for the ionospheric response to bombing in the current study.Such a mechanism may therefore contribute to the observed effect.
Infrasonic waves generated by explosions are launched preferentially in a vertical direction (e.g.Blanc, 1985).Any impact on the upper atmosphere overhead would then require horizontal transport to move that region over Slough.For the scale sizes involved this would require winds of the order of 300 km h −1 to blow consistently in a westward direction for the impact to be observed within 3 h, as suggested by the data.For this to happen continuously throughout the epoch being studied is unlikely.Whatever the cause of the ionospheric depletion, the impact appears to last less than 24 h, since only the subsequent noon value is affected.
One metric ton of TNT has an explosive energy of 4.184× 10 9 J, which is of the same order of energy as a cloud to ground lightning stroke.While a ground-based explosion and a lightning stroke are somewhat different in location and duration, it is not unfeasible that wave energy generated by lightning could also influence the ionosphere in a similar way.Since the occurrence of lightning has distinctive diurnal and seasonal cycles, it is feasible that this mechanism could contribute to the observed seasonal anomaly in ionospheric F-region electron concentrations (Rishbeth and Müller-Wodarg, 2006).

Conclusions
We have shown that the influence of major bombing raids over Europe during the latter half of WWII caused a significant depletion of the ionosphere above Slough, UK.This depletion is consistent with an enhanced ionospheric loss rate resulting from thermospheric heating via the dissipation of wave energy in the upper atmosphere.While the list of bombing raids used in this analysis is by no means complete, it is nevertheless sufficient to reduce the background noise in a composite analysis to a level where the ionospheric response to bombing is statistically significant.No lower limit to the minimum mass of explosives required to generate such a response was revealed in this study, though raids using ≥ 100 < 800 t of HE were observed to deplete the ionosphere above Slough if they took place between 3 and 7 h before the ionospheric measurement.It is possible that raids occurring outside of these thresholds could still produce an effect, but the current study contains insufficient data to test this.
One metric ton of TNT has an explosive energy of 4.184× 10 9 J, which is of the same order of energy as a cloud to ground lightning stroke.Since the occurrence of lightning has distinctive diurnal and seasonal cycles, it is feasible that a similar mechanism could contribute to the observed seasonal anomaly in ionospheric F-region electron concentrations (Rishbeth and Müller-Wodarg, 2006).
The duration of the ionospheric effect appears to be limited to within 24 h.This is currently restricted by the resolution of the digitised data.While hourly foF2 data were digitised from the original analogue ionograms, these data could be re-examined to increase the temporal resolution and investigate the behaviour of other parameters such as the D, E, F1 and sporadic-E layers.The ionospheric sporadic-E layer has already been shown to be enhanced by terrestrial lightning (Davis and Johnson, 2005;Johnson and Davis, 2006;Davis and Lo, 2008).Examining the response of this layer to terrestrial explosions would also provide further information on mechanism(s) involved in this process.The digitisation of these records is therefore highly desirable for a range of different research topics.
Data availability.The ionospheric data used in this study are publicly available via the UK Solar System Data Centre (https://www.ukssdc.ac.uk/).The majority of the information about the times, dates and munitions used in the bombing raids was obtained from the published references within this paper.Additional individual sources are cited in Appendix A.

Figure 1 .
Figure1.A superposed epoch analysis of the ionospheric response to major bombing raids over Europe.The black line in the top panel shows the median response in δfoF2 (foF2 with a 30-day median subtracted) to 60 bombing raids used as a trigger event in the analysis.The grey shaded area represents the standard error in these median values, while the dashed and dotted lines represent the 95th and 99th percentiles of 10 000 repeated analyses using random control days in which no bombing raids were identified.At time = 0 (within 24 h of the trigger events) δfoF2 is depleted (1.9 standard deviations from the mean).Arbitrary threshold values were used, corresponding to > 700 t of high explosives per raids occurring within 10 h before the noon ionospheric measurements.

Figure 2 .
Figure2.Panel (a) presents the relative significance of the ionospheric response at time zero (estimated standard deviations from the mean) for a range of thresholds.Events preceding the ionospheric measurements were considered for time windows from 0 to 23 h ahead of the noon ionospheric measurement and the minimum tonnage of HE bombs used in each raid was also varied from 100 to 1000 t.The most significant ionospheric response occurred for events occurring within 5 h before noon, in which a minimum of 300 t of HE bombs were dropped.Panel (b) presents the number of events contributing to each analysis.The significance of the response decreases if the threshold is extended beyond 5 h ahead of the ionospheric measurement.This indicates that events occurring at larger time offsets are not contributing to the observed median response.The significance of the result 5 h ahead decreases as the threshold on the minimum tonnage of HE bombs increases.This is likely due to the decreasing number of events contributing to each analysis.

Figure 3 .
Figure3.The same as Fig.1, except for the most significant response using 20 bombing raids within threshold values corresponding to > 300 t of high explosives per raids occurring within 5 h before the noon ionospheric measurements.The ionospheric response at time = 0 is around the 99th percentile of 10 000 control studies using random dates on which no bombing raids are catalogued within the study.