Reader TTS

ABSTRACT On September 1st, 2018 a devastating explosion occurred on the facility of an oil refinery near Ingolstadt, Germany. We analyzed data of 400 permanent and temporary seismic stations

and find strong seismo-acoustic signals on more than 80 seismic stations. The infrasound signal is detectable on seismic stations within 10–350 km from the source, with 40 km spatial

resolution. We confirm the explosion site both by the seismic and seismo-acoustic arrivals. Apart from seismic P- and S-waves, we identified three separate acoustic phases with celerities of

332, 292, and 250 m/s, respectively, each of which has a particular spatial pattern of positive detections at the ground. Seismo-acoustic amplitudes are strongly affected by the type of

seismic installation but still allow insight into regional infrasound attenuation. Our observations likely represent tropospheric, stratospheric, and thermospheric phases. We performed 3D

acoustic ray tracing to validate our findings. Tropospheric and thermospheric arrivals are to some extent reproduced by the atmospheric model. However, ray tracing does not predict the

observed acoustic stratospheric ducts. Our findings suggest that small-scale variations had considerable impact on the propagation of infrasound generated by the explosion. SIMILAR CONTENT

BEING VIEWED BY OTHERS MYSTERIOUS SEISMOACOUSTIC SIGNALS OF EASTERN HELWAN QUARRY BLASTS 2022 Article Open access 27 May 2025 RAYLEIGH WAVE ATTENUATION AND PHASE VELOCITY MAPS OF THE GREATER

ALPINE REGION FROM AMBIENT NOISE Article Open access 25 November 2024 USING DENSE SEISMO-ACOUSTIC NETWORK TO PROVIDE TIMELY WARNING OF THE 2019 PAROXYSMAL STROMBOLI ERUPTIONS Article Open

access 14 July 2021 INTRODUCTION Dense seismic networks such as the transportable USArray or the European AlpArray1 enable scientist to study the propagation of seismic waves in

unprecedented spatial detail. While they are designed to record and analyze seismic wavefields with great precision2,3,4, recent studies highlight the ability of seismic networks to also

detect and track acoustic signals over distances up to hundreds of kilometers5,6. As dedicated infrasound installations are rare anywhere on earth, this renders seismic arrays and dense

regional seismic networks particularly useful to study the infrasound signals generated by strong explosive events, both of natural7,8,9,10 and man-made origin6,11,12,13. Data from an

accidental rocket explosion recorded by the transportable USArray highlighted the potential wealth of information that could be gained from spatially dense infrasound records14. Here we

report on a devastating explosion that occurred on September 1st, 2018 near Ingolstadt, Germany. The explosion has created both seismic and acoustic waves that were detected by seismic

stations up to 400 km distance. The dense spatial station coverage of the AlpArray seismic network allows to study the propagation of both the seismic and the acoustic wave field in great

detail. Here, we focus in particular on the seismo-acoustic observations, that reveal a unique number of various seismo-acoustic phases detected on regional scale and potentially allow for

improved insight into local and regional infrasound propagation and into seismo-acoustic coupling. OBSERVATIONS SEISMO-ACOUSTIC DATA According to news reports the explosion occurred on

September 1st, 2018 between 03:00 and 03:15 UTC at the site of an oil refinery operated by BayernOil near the village of Vohburg in Bavaria, Germany (Coordinates: 48.766, 11.595). Local

newspaper report that the explosion was audible in several kilometers distance and caused severe damage to buildings in neighboring villages15. According to reports, hot fuel escaped through

a crack from a highly pressurized (25 bar) container and eventually got ignited. Structural damage near the explosion indicates that the explosion occurred above ground. The explosion site

is located 12 km east of Ingolstadt and 10 km south of the Gräfenberg broadband array16 and is well-covered by temporary broadband stations of the AlpArray seismic network1 with station

spacings of approximately 40 km (see Fig. 1). The AlpArray stations provide excellent spatial coverage for this event, in particular as Bavaria in southern Germany is only sparsely

instrumented with permanent sensors. Additionally, most AlpArray stations provide data in real time, so timely analysis of the explosion data was possible. We acquired data from a total of

400 seismic broadband stations within 400 km radius of the explosion site (Fig. 1). Waveform data and associated metadata were downloaded from the European Integrated Data Archive (EIDA). We

focused exclusively on vertical component data, with 100–200 Hz sampling rate. All data were corrected for the instrument response to obtain ground velocity records. We manually inspected

the data for impulsive signals that are potentially related to the explosion, and travel at seismic velocities (2–8 km/s) as well as acoustic velocities (≈0.33 km/s). Figure 2 shows a record

section of all available data up to 400 km from the explosion. We identified both seismic P and S phases up to a distance of 400 km (the maximum extent of our study area). They are

traveling with an apparent velocity of approximately 5.5–5.9 km/s and 3.5 km/s, respectively. Additionally, three separated phase arrivals are visible at times between 50 and 1300 seconds

after the explosion and up to 350 km distance. They appear most prominently in the 1–2 Hz frequency range (see Fig. S1, electronic supplement, for a comparison of different frequency bands).

These phases travel with velocities of 310–340 m/s (1st phase; visible at 0–200 km distance), 270–300 m/s (2nd phase; visible at 150–325 km distance) and 240–260 m/s (3rd phase; visible at

250–350 km distance) (see also Table 1). There is considerable offset between neighboring phases (≈120 s between phases 1 and 2; ≈160–180 s between phases 2 and 3), and within certain

distance ranges (e.g. 150–175 km and 250–350 km) more than one phase is visible. These velocities strongly suggest they represent seismo-acoustic signals that were generated by the refinery

explosion. The seismo-acoustic explosion signal is detected by a large number (>80) of seismic stations, complicating any manual analysis of single-station data. Thus, we defined

objective criteria to assess if the signal is positively detected at a given station. Based on the observations described above we defined three separate arrival time windows, calculated

using the observed effective velocities: Phase 1, 310–340 m/s; Phase 2, 270–300 m/s; Phase 3, 240–260 m/s. For each time window we determine the maximum absolute value of the vertical peak

ground velocity (PGV) in the 1–2 Hz frequency band, and the corresponding celerity (=total travel time from source to receiver divided by horizontal distance). Background noise is calculated

as three times the standard deviation (3_σ_) of seismic amplitudes (in the 1–2 Hz band) within a 60 s window that ends 20 s before the explosion. If the measured vertical PGV exceeds the

background noise by a factor of three we label the respective station as a positive detection. This way we identify the seismo-acoustic signal on in total 84 different seismic stations

within 350 km from the explosion site. 35 of these are temporary AlpArray installations. None of the stations shows high SnR signals of more than one individual phase. Note that potentially

more stations may show weaker but still manually recognizable signals of the explosion and secondary phases. These are however not considered in this study. Table 1 summarizes the number of

detections for the three different phases, as well as the calculated celerities. Figure S2 displays the distribution of all individual arrival times and calculated celerities as read from

the PGV in each time window. SPATIAL DISTRIBUTION OF DETECTIONS The positive detections of the seismo-acoustic explosion signal show pronounced spatial patterns that are significantly

different among the three separate phases. Figures 3 and S3 show the spatial distributions of amplitudes for all three phases. Phase 1 is predominantly observed in southerly to northwesterly

directions (135–315° azimuth from the source) at distances up to 200 km (see Figs 3 and S3a). A branch of weak signals also potentially stretches out northwards past stations of the

Gräfenberg array up to distances of almost 250 km. Seismo-acoustic phase 2 is recorded at distances between 150–325 km, and in narrow corridors towards North (330–45°) and West (240–300°)

and is absent elsewhere (Figs 3 and S3b). Phase 3 is recorded in a broad azimuth range (350–115°) from North to South-East of the explosion site at distances of 250–350 km (see Figs 3 and

S3c). Note, that within the distance ranges of phases 2 and 3 (beyond 150 km distance) the spatial station coverage is greatly reduced into north-western directions (beyond the extent of the

AlpArray seismic network) and north-eastern directions (AlpArray stations existing, but no real-time data available). LOCATION: SEISMIC VS. ACOUSTIC In order to independently locate and

time the source of the detected signals we identified 9 stations that show the earliest arrivals of both the seismic and the seismo-acoustic phase. For these stations we manually picked the

first onset of the seismic P-arrival and the first onset of the seismo-acoustic arrival (see Figs S4 and S5, supplement and Table 2). Both sets of onset readings were independently used for

a simple grid search routine that determines the epicenter and origin time of the explosion by minimizing the root-mean-squared residuals between observed and theoretical arrival times of

all stations simultaneously. Grid points are spaced by 1 km and time is searched in steps of 0.1 s for seismic waves and 1 s for acoustic waves. We obtained an estimate of the propagation

velocities through linear regression along the picked first onset, assuming that the velocity is constant and isotropic. This resulted in a propagation velocity of 5.9 km/s for the seismic

P-waves and 335 m/s for the acoustic waves (see Fig. S2). Figure 4 compares the location results based on seismic and acoustic phases. Using the seismic phases we locate the source 1 km

(i.e. one grid point) north of the true origin. Using the acoustic phases we locate the source 1.5 km (i.e. one grid point along diagonal) south-west of the true origin. This clearly

confirms the refinery site as the epicenter of both the seismic and the seismo-acoustic signal. Relative to the mean total travel time, the remaining travel time residuals for the best

epicenter solution are similar for the two approaches (1.5 s/108 s ≈ 1.4% acoustic; 0.07 s/7.6 s ≈ 0.9% seismic). The extension of the 10%-error ellipse is small (<1 km) both for the

seismic and acoustic result (see Fig. 4). Based on the acoustic picks, we determine the origin time of the explosion to 03:11:44 UTC, and based on the seismic picks the the best-fitting

origin time is 03:11:46 UTC. Thus we define the origin time to the mean value of 03:11:45 UTC ± 1 s. The explosion is listed in the event bulletin of the International Seismological Centre

(ISC, event ID 61268997717) with an origin time of 03:11:45 UTC, in agreement with our findings. The hypocenter is located at the surface and 2.5 km south of the refinery site. The local

magnitude is M_l_ = 2.0 and the body wave magnitude is m_b_ = 1.9. FIRST ONSET POLARITIES On the same stations as used for the localization we also analyzed the polarity of the first onset

for both the P-wave and the seismo-acoustic wave. For this purpose data was not filtered between 1–2 Hz as for most of the other measurements, but with a wide passband of 0.1–25 Hz to avoid

artificial disturbance of the waveforms (note that the lower corner frequency of 0.1 Hz is far below any frequencies contained in the signal. Without any filtering, the onsets are hardly

visible on several stations due to strong noise. For consistency we applied the filter even to the low-noise stations). All P-wave onsets show a clear positive vertical first motion

(upwards), regardless of the azimuth between station and source (see Fig. S4, supplement). In contrast, the onsets of the seismo-acoustic waves are generally less pronounced and do not show

a regular pattern (see Fig. S5, supplemental), with some stations showing clear downward first motion and others suggesting upward first motion, with no correlation to the azimuth.

SEISMO-ACOUSTIC AMPLITUDES AND ATTENUATION The dense spatial coverage around the explosion site allows some potential insight into the attenuation of the seismo-acoustic signal with

distance. For this analysis we only considered stations with a signal-to-noise ratio of larger than three (in the 1–2 Hz frequency band). Figure 5 displays the maximum vertical amplitude of

the seismo-acoustic signal with distance. Our data show large variations in seismo-acoustic amplitudes, but on average within 200 km from the source the observed seismo-acoustic amplitudes

attenuate by approximately 1–2 orders of magnitude. Notably, we observe that seismo-acoustic amplitudes are generally higher on the temporary AlpArray installations compared to the permanent

seismic stations. Within 25–50 km from the source, the PGV induced by the seismo-acoustic signal at AlpArray stations is (on average) more than 70 times larger than average PGV measured on

permanent stations (e.g. the Gräfenberg array). Within 50–150 km from the source the seismo-acoustic mean PGV measured on temporary AlpArray stations is still 10 times higher than on

permanent stations at the same distance from the source. Factors that impact the seismo-acoustic amplitudes are reviewed in the discussion section. ACOUSTIC RAYPATH MODELING The previous

observations clearly indicate that we are observing a seismo-acoustic phenomenon, and that the refinery explosion created strong infrasound signals. Thus, similar to the approach in6 we

modeled the propagation of infrasound rays from the source up to 400 km distance, in order to compare modeled surface bounce points with our observations. We used the GeoAc ray tracing

suite18,19 to perform 3D ray tracing in a proper atmospheric model. We extracted the required atmospheric parameters (wind speed, wind direction, temperature, etc.) from an hourly updated

forecast model provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). We retrieved the 3D ECMWF model for 2018/09/01 03:00 UTC on a grid ranging from 45° to 52°

longitude and from 7° to 17° latitude with 0.125° resolution. For distances up to 400 km from the source thermospheric arrivals are expected, so we extended the ECMWF model to higher

altitudes using the MSISE-0020 and HWM1421 empirical models. To guarantee a smooth transition between ECMWF and HWM14/MSISE-00 models we mixed the models between 67 and 77 km altitude. Below

67 km altitude we used the pure ECMWF model and above 77 km altitude the pure model derived from MSISE00 and HWM14. A linear transition between the models is implemented in the 67 km–77 km

range. The resulting merged input model for GeoAc (calculated using the barometric formula) measures height above sea level (with z = 0 being sea level). However, as infrasound only

propagates in air above ground, we subtracted the local topography at each grid point from the height column such that for each grid point the height measures start at ground level (z = 0

being ground level instead of sea level). Figure 6 shows the surface bounce points resulting from the 3D acoustic ray tracing, colored by the expected celerity, and compares the modeling

results to our measurements. The acoustic ray tracing predicts a high-celerity tropospheric duct resulting in several bounce points within 0–150 km distance from the source, in the 130–300°

azimuth range. Notably, the model does not predict any stratospheric ducts and consequently no surface bounce points with medium celerities. Thermospheric reflections with low celerity are

expected at all azimuths, yet at different distances depending on azimuth. All calculated thermospheric reflections show celerities increasing with distance. Please refer to the discussion

section for a detailed comparison of measurements and modeling. The ray tracing results predict precise travel times for several seismic stations that co-locate with calculated tropospheric

bounce points (see Table 2). We deduct this theoretical travel time from the observed first onset pick to independently determine the origin time of the explosion. Measurements from six

stations within 20 and 86 km distance result in a mean origin time of 03:11:41.9 UTC ± 2.8 s (see Table 2 for details), slightly earlier than estimated by the seismic and acoustic grid

search, but still within the estimated uncertainties. Note, however, that we could not use the same set of stations as for the grid search, because the modeling did not reproduce any bounce

points north of the explosion. All stations with modeled travel times available are exclusively south/west of the explosion, potentially also introducing spatial bias. DISCUSSION

SEISMO-ACOUSTIC DETECTIONS The seismic and seismo-acoustic signals we have observed are clearly associated with the refinery explosion near Ingolstadt, Germany on September 1st, 2018. The

origin times determined from seismic and acoustic grid search as well as the ray tracing agree within their respective uncertainties. Additionally, all observed acoustic phases show a stable

and steady move-out from the source with distance, evidence by well-aligned acoustic arrival times and little scatter in measured celerities (see Fig. S2). Thus, we are confident that all

signal-to-noise ratio based detections for phases 1 and 2 and also most of the detections for phase 3 are in fact due to the refinery explosion and not due to local noise sources.

Additionally, we searched the publicly accessible seismic event bulletin of the Swiss Seismological Service (SED), the earthquake service for Bavaria (LMU), the monthly E-Mail event bulletin

for Austria (ZAMG) and the ISC bulletin for seismic events that occurred close in time to the Ingolstadt explosion. Except the explosion itself (registered in the ISC event bulletin) there

are no seismic events registered within 400 km from the explosion site and within hours after the explosion. Consequently we can rule out that e.g. the late arrivals of phase 3 are due to

other seismic events. Closer inspection of the seismic record sections also shows no indications for potential secondary infrasound sources (e.g. from topography22). Each of the observed

seismo-acoustic phases displays a different celerity, indicating that they each follow an individual propagation path. The values of 332 m/s (phase 1), 292 m/s (phase 2) and 250 m/s (phase

3) suggest that the three phases correspond to three separate infrasound ducts at different heights within the atmosphere23. The measured celerities match the expected values for

tropospheric (phase 1), stratospheric (phase 2) and thermospheric (phase 3) waveguides. Thus we conclude that we have observed all three types of infrasound pathways on the seismo-acoustic

recordings of the explosion. Below we will discuss the comparison of our measurements with acoustic ray tracing. The dense spatial sampling of the infrasound field by the AlpArray seismic

network reveals large zones in which no or only very weak seismo-acoustic arrivals were measured. This observation should, however, not be confused with classical shadow zones known for

infrasound propagation23, as it’s still influenced by the seismo-acoustic coupling. Notably, the infrasound array IS26 of the CTBTO International Monitoring System, which is located near the

border tripoint of Germany, Austria, and the Czech Republic, (near station GEC2 in Figs 1 and 3) falls into an area where no or only very weak seismo-acoustic signals were detected in this

study. However, the acoustic signal was in fact detected by the very sensitive infrasound array IS2624. This demonstrates that acoustic energy may be present even though it does not map into

the seismic recordings, because of the required seismo-acoustic coupling which may only be effective above a certain acoustic amplitude threshold25. LOCATION AND ORIGIN TIME We obtained the

origin time of the explosion by three independent approaches: Seismic grid search, acoustic grid search and modeled travel times, all based on the manual onset picks. The results from

seismic and acoustic grid search agree within 2 seconds difference and determine the origin time on average to 03:11:45 UTC ±1 s. Due to the uncertainty of the underlying velocity estimates

we do not speculate if the different origin times may have a physical explanation. The determination of the origin time depends on the seismic and acoustic velocities used during the grid

search. Any false assumption may map into both the location and origin time. We tested seismic velocities between 5–6 km/s and acoustic velocities between 330–350 m/s and compared the

resulting output of the grid search. The best fitting location is only weakly affected by choosing a different velocity, and varies around 1 km from ground truth. This indicates that the

local velocity structure (in both ground and atmosphere) is rather homogeneous and that assuming a constant value is reasonable. Consequently, changing the propagation velocities only

affects the origin time, with slower propagation resulting in earlier origin time. However, ground truth video records are available that rule out any origin time earlier than 03:11:44, UTC

(K. Koch, personal communication). Using the modeling results, that is subtracting modeled acoustic travel times from the first onset of the acoustic signal, results in origin times

considerably earlier than estimated by the grid search approach (see Table 2). However, since bounce points are only modeled West and South of the explosion site, the results strongly rely

on correct measures of the acoustic velocities and trajectories. As video material contradicts any origin time earlier than 03:11:44, we conclude that the ray tracing likely overestimated

the true travel times, and that origin times based on ray tracing may not be reliable in this case. SIGNAL POLARITIES AND AMPLITUDES For the nine seismic stations closest to the explosion

site the first motion polarity of the the P wave shows a clear upwards motion (see Fig. S4), regardless of the azimuth. This is consistent with an isotropic source model as expected for

explosive sources. The first motion polarity of the seismo-acoustic signals is less conclusive, yet it indicates a downward motion on several stations (see Fig. S5). This polarity would be

in agreement with acoustic pressure pulses arriving from above ground, pushing down the earth surface and thereby locally inducing downward first motion. Still, local seismo-acoustic

coupling conditions in the vicinity of a station might potentially obscure clear first motion polarities, which could be why we observe a less conclusive pattern for acoustic waves compared

to the seismic waves. In principle, the high spatial density of our measurements could allow to determine the attenuation of the infrasound waves on regional scale, especially in potential

tropospheric ducts. However, seismic sensors measure ground velocity and the exact relation of ground velocity to air pressure changes is unknown and might be highly site-specific. We

observe large scatter in the maximum amplitudes of the seismo-acoustic signal (Fig. 5). For any given distance range the measured seismo-acoustic amplitudes generally scatter within 1–2

orders of magnitude. Differences in seismo-acoustic amplitudes may be due to different attenuation properties along the acoustic raypath, in particular caused by potentially anisotropic

attenuation due to winds26. Differences in seismo-acoustic amplitudes may also be due to local coupling conditions. The ray tracing results due not indicate strong differences in amplitudes,

and the fluctuations we observe are likely due to specific seismo-acoustic coupling conditions at the individual seismic stations. Generally, air-to-ground coupling is poorly understood,

especially when it comes to relating seismic amplitudes to acoustic amplitudes25,26,27. Potential mechanisms for air-to-ground coupling are direct vertical forcing of the ground due to air

pressure pulses28,29 and so called air-coupled Rayleigh waves30. The latter are generated if Rayleigh waves velocities in near surface layers are equal to the speed of sound or to the sweep

velocity of the acoustic wave field30. The arrival times of the seismo-acoustic signals in this study are consistent with purely acoustic propagation and do not indicate any air-coupled

Rayleigh wave as origin. Thus, we focus on the direct surface forcing by a pressure pulse. We suggest that this mechanism is responsible both for the primary creation of seismic waves at the

explosion site as well as the creation of seismo-acoustic signals at each individual seismic station. However, very few studies exist that attempt to investigate this coupling mechanism in

detail, in particular in the context of near-surface explosions. Generally, soft surface conditions will allow for more efficient seismo-acoustic coupling than stiff surface conditions.

Several studies are available that investigate the efficiency of seismo-acoustic coupling after bolide/fireball explosions at high altitudes27,31. Such studies highlight that the coupling

efficiency is affected by frequency31 and even more so by local site conditions27,32. To our knowledge no study investigates the seismo-acoustic coupling efficiency at the frequencies

observed in this study (1–2 Hz). Additionally, there is still considerable discrepancy among different studies about the efficiency of air-to-ground coupling, ranging from efficiencies as

low as 10−5 30,31 to values of 0.527, again likely due to site effects. Still, Edwards _et al_. provide a detailed analysis of air-to-ground coupling, and point out that in principle a

scaling relation between seismic ground motion (in m/s) and air pressure changes (in Pa) can be established27, if the properties of the shallow subsurface are well known. However, such

information (e.g. Vs30) is unfortunately not available for the stations we used in this study. Other studies suggest that energy transfer from air to ground requires a soft and

unconsolidated surface layer of at least 10 m thickness32. In this study, seismo-acoustic amplitudes measured on temporary seismic stations are notably higher than amplitudes measured on

permanent seismic stations at the same distance (see Fig. 5). Within 25–50 km from the source, the seismo-acoustic amplitude is on average more than 70 times higher on temporary AlpArray

stations than on permanent stations. At larger distances, temporary stations still on average record 10 times larger seismo-acoustic amplitudes than permanent stations. Permanent seismic

stations are commonly installed on hard rock in underground vaults and are carefully shielded from external disturbances, such as pressure variations. In particular the seismic stations of

the Gräfenberg array just 10 km north of the explosion site are installed into solid limestone formations16. All permanent stations within 50 km from the explosion are part of the Gräfenberg

array. Hence, we expect weak seismo-acoustic coupling, as evidenced in Fig. 5. In contrast, most of the temporary AlpArray stations are installed in near-surface structures such as

buildings and within the top-most potentially soft surface layer. In particular, all temporary stations south of the explosion site are inside the soft Molasse unit. Consequently, these

stations are potentially more sensitive to atmospheric perturbations, as they are both more exposed and situated in soft surface structures. Additionally, the AlpArray stations are generally

less carefully shielded against pressure variations33,34,35,36. Both effects will allow for more efficient seismo-acoustic coupling into temporary stations than into permanent stations.

This may explain the large variation of seismo-acoustic amplitudes in Fig. 5. Thus, we only point out that despite the amplitude variations, between 0–200 km from the source (that is within

the range of phase 1 = tropospheric ducting) seismo-acoustic amplitudes on average decay by approximately one order of magnitude per 100 kilometers (Fig. 5). This holds even if some of

highest or lowest amplitudes are solely due to site effects. For distances larger than 200 km from the source, we cannot observe any clear trend. Generally, we do not account for potential

azimuthal anisotropy. The large scatter in seismo-acoustic amplitudes renders any detailed analysis of the acoustic attenuation unreliable, yet we attempt to relate our observations to the

findings of other studies. Le Pichon _et al_. derived an attenuation formula for infrasound from a synthetic dataset, assuming the presence of a stratospheric duct but no tropospheric

duct37. At 180 km distance from the source they predict that without a tropospheric duct amplitudes are attenuated by 4–5 orders of magnitude due to geometrical spreading and dissipation of

direct waves within the troposphere. Our seismo-acoustic observations indicate that within the troposphere attenuation may be significantly less if an efficient infrasound duct is present.

However, the unknown and highly site-dependent seismo-acoustic coupling prohibits an in-depth analysis of the acoustic attenuation in this study. To conclude we highlight that several

studies suggest using the dependence of air-to-ground coupling efficiency on local material properties to turn seismo-acoustic observations into a tool to locally infer such properties of

the near surface layers32,38. This requires simultaneous acquisition of pressure and seismic data, though. In turn, if site properties such as Vs30 were known for all sites in this study, it

could be possible to reconstruct the acoustic wave field from seismic records alone27. INFRASOUND PROPAGATION Above we argued that phases 1, 2 and 3 likely correspond to infrasound ducting

within the troposphere, stratosphere and thermosphere, respectively. We performed acoustic ray tracing to test our hypothesis and to explain the distinct spatial detection patterns that we

observe for the three different phases. Phase 1 is predominantly observed at azimuths between 135–315° at distances up to 200 km (see Figs 3 and S3a). These detections coincide well with

expected surface bounce points and expected celerities predicted by the numerical modeling (see Fig. 6), confirming that we observe infrasound waves guided within the troposphere. The

corresponding turning height for the tropospheric reflections is less than 1 km. However, ray path modeling for several stations West/South of the explosion site apparently overestimates the

acoustic travel times, resulting in unrealistically early origin times for the explosion (see Table 2). This indicates, that the model may underestimate celerities within the tropospheric

duct or produce wrong trajectories. It should be noted that very shallow turning heights as observed here should generally interpreted with caution, being on the limits of what ray tracing

can resolve. Towards south-west ray tracing also predicts tropospheric waves only up to 80 km from the source, while we still observe them at 160 km distance. Additionally, the modeling does

not reproduce the narrow corridor of high-celerity detections north of the explosion site. Weather balloon data is available from a launch site 80 km north of the explosion site, however

only for 2018-09-01 00:00 UTC (3 hours before the explosion). Using wind and temperature data measured by the balloon for the ray tracing we are able to reproduce some tropospheric

reflections towards north, with very shallow turning heights of <1 km (see Fig. S6). Despite the three hours offset in time, this indicates that sharp but small-scale variations in the

atmosphere that are not captured by the smooth ECMWF model might be responsible for the tropospheric detections north of the explosion site. We observe phase 2 in the azimuth range 240–300°

(west) and 330–45° (north) (see Figs 3 and S3b) and at distances of 150–325 km from the explosion. Note however that within this distance range there is an azimuthal gap with little to no

stations between 300–330° (north–west). Thus, we cannot state for sure if the two patches of stratospheric detections are separated or potentially connected all the way from west to north.

In any case, we observe clear evidence for stratosphere-guided waves that are not at all reproduced by ray tracing (Fig. 6). The potential presence of a stratospheric duct is usually

evaluated by the Veff-ratio, which is defined as the ratio between the effective sound speed at stratospheric altitude and the sound speed at ground level37. Stratospheric reflections are

expected when the effective sound speed at height is faster than at ground, that is for Veff-ratio > 1. Figure 7 shows a map of the effective sound speed ratio Veff-ratio 15 minutes

before the explosion (03:00 UTC), based on the ECMWF atmospheric model. Not knowing the exact stratospheric boundary, we chose the maximum effective sound speed within 35–60 km altitude as a

measure for the stratosphere. Figure 7 displays the Veff-ratio for our study area as if each respective point on the map was the hypothetical turning point of a stratospheric ray

propagating radially outwards from the explosion site. Note that in order to have a detection at the surface the corresponding turning point would be approximately halfway in between the

detection point and the explosion site. The map clearly indicates that for the study area the ECMWF model only predicts Veff-ratios smaller than one (except at the edge of the study area

towards south-west. However, the corresponding potential surface bounce points would lie beyond our study area). This is because the effective sound speed at ground at the explosion site is

always slightly larger than the effective sound speed at stratospheric heights in any direction and at any distance. Consequently, ray tracing should not predict any stratospheric

reflections. Note, that this considerations about the spatial distribution of the Veff-ratio are based on a simplified 1D propagation of the rays in each direction. The actual modeling was

done with fully three-dimensional propagation, which may slightly deviate from the simplified 1D path. Our stratospheric detections north of the explosion follow a similar corridor as the

tropospheric detections. As the atmospheric conditions are generally close to allowing stratospheric reflections with Veff-ratio close to 1, we speculate that small-scale atmospheric

perturbations are responsible for both the tropospheric and stratospheric reflections towards north. Such perturbations might be particularly effective if located close to the infrasound

source39,40. Similarly, the stratospheric detections towards west could be due to small-scale perturbations that are not included in the atmospheric model or due to gravity waves. Phase 3 is

predominantly detected in an azimuth range from 350–115° (Figs 3 and S3c), but it shows up only weakly on Fig. 2 and should be interpreted with caution. Based on the measured celerity and

distances we suggest that phase 3 corresponds to thermospheric arrivals. The acoustic modeling, however, predicts thermospheric arrivals at all azimuths and at slightly different distances

than observed (Fig. 6). For northwards and eastwards directions we detect phase 3 about 50 km closer to the source than expected. Few detections towards north-east even fall 100 km closer

than predicted by modeling. Note, that the modeled celerity of all thermosperic arrivals increases with distance. This feature is for example reproduced by a group of measurements towards

the East (in Austria) that show a similar increase with distance for the measured celerity and generally a good match between modeled and measured celerity, despite the 50 km offset in

bounce points. This supports the interpretation of phase 3 as thermospheric. Generally, however comparing measurements to model results suggests that the thermospheric conditions on that day

were not properly reproduced by the atmospheric model. Given that the effective sound velocity in the thermosphere always exceeds the effective sound velocity at ground level, thermospheric

phases are always predicted by modeling. However, they are only rarely observed even on microbarographs due to the relatively low density of the atmosphere at high altitudes that strongly

attenuates especially higher frequencies41. Gibbons _et al_. observed weak seismo-acoustic thermospheric arrivals on seismic arrays in Scandinavia12. Beside low celerity values they confirm

that thermospheric arrivals consist of lower frequencies26. Assink _et al_. report similar effects after the 2017 North Korea nuclear test13. We cannot clearly identify such effects in our

measurements, but we note that the potential thermospheric arrivals are no longer present for frequencies higher than 3 Hz, while tropo- and stratospheric detections are still well-visible

(Fig. S1, supplemental). Compared to a similar study by6 our observations are less in agreement with predictions by acoustic ray tracing, even though we calculated the infrasound propagation

in the fully 3D atmospheric model. The discrepancy between observation and model calls for explanation. Small scale disturbances at low altitudes may not be captured by the smooth ECMWF

models, but only in local measurements from e.g. weather balloons. Such measurements, however, are only taken every 12 hours and may not properly represent the state of the atmosphere for

the time of interest. Regarding the discrepancy among modeled and observed stratospheric reflections we point out that the explosion occurred close in time to autumnal equinox. Studies

highlight that at this time of the year there is little stratospheric wind, enhancing turbulence and small-scale wave activity in the middle atmosphere (stratosphere-mesosphere)13. In fact,

the ECMWF model we applied does generally show weak stratospheric winds. This may allow for disturbances that are not accounted for in the atmospheric models, but only in local measurements

e.g. from Lidar42. Such measurements, however, are only taken every 12 hours and may not properly represent the state of the atmosphere for the time of interest. During calm stratospheric

conditions small-scale structures and wind shear layers can occur39, causing unexpected reflections and infrasound arrivals, such as our tropospheric and stratospheric detections north of

the explosion site. Assink _et al_. point out that high frequency signals (>1 Hz) such as the ones we observe are more likely to be affected by small-scale structures13. In turn this

demonstrates the usefulness of dense spatial sampling of the infrasound wave field. Theoretically, the discrepancy between observations and models could be utilized to improve the underlying

atmospheric models. Dense seismo-acoustic measurements may then eventually allow to infer atmospheric properties on several hundreds of kilometers scale43. CONCLUSIONS AND OUTLOOK We

observed strong seismic and seismo-acoustic signals generated by an explosion inside an oil refinery near Ingolstadt, Germany. Seismic phases are detectable up to 400 km from the source

(extension of our study area), and seismo-acoustic phases are detectable up to 350 km from the source. The dense seismic station coverage of 40 km allowed to study the propagation of seismic

waves and infrasound waves with high spatial resolution. We detected tropospheric, stratospheric, and potentially thermospheric reflections. Our work demonstrates the usefulness of dense

seismic networks to study both seismic and infrasonic wave phenomena in great detail. We highlight that such networks may also be valuable for atmospheric studies, in particular as dedicated

infrasound sensors are usually rare. We demonstrate that on regional scale there might be considerable discrepancy between the modeled and the true state of the atmosphere. Especially

temporary seismic installations are very sensitive to impulsive atmospheric perturbations, as e.g. generated by explosions. However, highly site-specific seismo-acoustic coupling complicates

the interpretation of seismo-acoustic amplitudes. Ultimately, with a more precise understanding of seismo-acoustic coupling, it might be possible to infer the state of the atmosphere from

seismic records. Even more so, if all seismic stations from the temporary AlpArray network (or comparable future networks) were equipped with additional pressure sensors - as was done in

parts for the USArray - explosive events like the one from this work would provide an unprecedented wealth of information, potentially enabling us to invert for atmospheric properties over

hundreds of kilometers14. DATA AND RESOURCES Seismic waveform data used in this study is available for download at the European Integrated Data Archive (EIDA) at

http://www.orfeus-eu.org/data/eida/index.html (last accessed August 2019). This study is based in part on data from the AlpArray seismic network44 which at the time of publication was not

publicly available. Please visit http://www.alparray.ethz.ch (last accessed August 2019) for more details on data access and for more information about the AlpArray seismic network. The ISC

earthquake catalog search engine is provided at: http://www.isc.ac.uk/iscbulletin/search/catalogue/ (last accessed August 2019). The Swiss seismic event bulletin is available at:

http://seismo.ethz.ch/en/research-and-teaching/products-software/earthquake-catalogues/ (last accessed August 2019). The Bavarian seismic event bulletin is available at:

http://erde.geophysik.uni-muenchen.de/ (last accessed August 2019). The Austrian event bulletin is available via E-Mail request ([email protected]). ECMWF forecast models are available to

registered users through the ECMWF Meteorological Archival and Retrieval System (MARS) or a python API: https://www.ecmwf.int/en/forecasts/datasets/archive-datasets (last accessed August

2019). The HWM14 model is accessible via a python interface: https://github.com/rilma/pyHWM14 (last accessed August 2019). The MSISE-00 model is accessible via a python interface:

https://github.com/scivision/msise00 (last accessed August 2019). The GeoAc ray tracing code is kindly provided on GitHub: https://github.com/LANL-Seismoacoustics/GeoAc (last accessed August

2019). All waveform data processing was done using the ObsPy toolbox45 versions 1.1.046 and 1.1.147. Codes, maps and graphs were created using44 Python 2.7.12 and Matplotlib 1.5.1. For this

study we used seismic data from several permanent seismic networks and we appreciate the continuous operation of these seismic networks by the responsible institutions (see references): BE

net48, BW net49, C4 net, CH net50, FR net51, GE net52, GR net53, GU net54, HU net55, IU net56, IV net57, MN net58, NI net59, OE net60, OX net61, RF net62, SI net, SK net63, 533 SL net64, ST

net65, SX net66 and TH net67. REFERENCES * Hetenyi, G. _et al_. The alparray seismic network: A large-scale european experiment to image the alpine orogen. _Surv. Geophys._ 39, 1009–1033,

https://doi.org/10.1007/s10712-018-9472-4 (2018). Article  ADS  PubMed  PubMed Central  Google Scholar  * Bowden, D. C., Tsai, V. C. & Lin, F.-C. Amplification and attenuation across

usarray using ambient noise wavefront tracking. _J. Geophys. Res. Solid Earth_ 122, 10086–10101, https://doi.org/10.1002/2017JB014804 (2017). Article  ADS  Google Scholar  * Buehler, J. S.,

Mancinelli, N. J. & Shaerer, P. M. S-to-rayleigh wave scattering from the continental margin observed at usarray. _Geophys. Res. Lett._ 45, 4719–4724,

https://doi.org/10.1029/2017GL076812 (2018). Article  ADS  Google Scholar  * Ekstrom, G. & Nettles, M. Observations of seismometer calibration and orientation at usarray stations,

2006–2015. _Bull. Seismol. Soc. Am._ 108, 2008–2021, https://doi.org/10.1785/0120170380 (2018). Article  Google Scholar  * Walker, K. T., Shelby, R., Hedlin, M. A., de Groot-Hedlin, C. &

Vernon, F. Western u.s. infrasonic catalog: Illuminating infrasonic hot spots with the usarray. _J. Geophys. Res._ 116, B12305, https://doi.org/10.1029/2011JB008579 (2011). Article  ADS 

Google Scholar  * Schneider, F. M. _et al_. Seismo-acoustic signals of the baumgarten (austria) gas explosion detected by the alparray seismic network. _Earth Planet. Sci. Lett._ 502,

104–114, https://doi.org/10.1016/j.epsl.2018.08.034 (2018). Article  ADS  CAS  Google Scholar  * de Groot-Hedlin, C. D. & Hedlin, M. A. H. Infrasound detection of the chelyabinsk meteor

at the usarray. _Earth Planet. Sci. Lett._ 402, 337–345, https://doi.org/10.1016/j.epsl.2014.01.031 (2014). Article  ADS  CAS  Google Scholar  * de Groot-Hedlin, C. D. & Hedlin, M. A. H.

A new automated approach to detecting and locating seismic events using data from a large network. _Bull. Seismol. Soc. Am._ 108, 2032–2045, https://doi.org/10.1785/0120180072 (2018).

Article  Google Scholar  * Hedlin, M. A. H., Ritsema, J., de Groot-Hedlin, C. D. & Hetland, E. A. A multidisciplinary study of the 17 january 2018 bolide terminal burst over southeast

michigan. _Seismol. Res. Lett._ 89, 2183–2192, https://doi.org/10.1785/0220180157 (2018). Article  Google Scholar  * Matoza, R. S. _et al_. Local, regional, and remote seismo-acoustic

observations of the april 2015 vei 4 eruption of calbuco volcano, chile. _J. Geophys. Res. Solid Earth_ 123, 3814–3827, https://doi.org/10.1002/2017JB015182 (2018). Article  ADS  Google

Scholar  * Ottemöller, L. & Evers, L. G. Seismo-acoustic analysis of the buncefield oil depot explosion in the uk, 2005 december 11. _Geophys. J. Int._ 182, 1123–1134,

https://doi.org/10.1111/j.1365-246X.2007.03701.x (2008). Article  ADS  Google Scholar  * Gibbons, S. J. _et al_. The european arctic: A laboratory for seismoacoustic studies. _Seismol. Res.

Lett._ 86, 917–928, https://doi.org/10.1785/0220140230 (2015). Article  Google Scholar  * Assink, J., Averbuch, G., Shani-Kadmiel, S., Smets, P. & Evers, L. A seismo-acoustic analysis of

the 2017 north Korean nuclear test. _Seismol. Res. Lett._ 89, 2025–2033, https://doi.org/10.1785/0220180137 (2018). Article  Google Scholar  * Vergoz, J., le Pichon, A. & Millet, C. The

antares explosion observed by the usarray: An unprecedented collection of infrasound phases recorded from the same event. In le Pichon, A., Blanc, E. & Hauchecorne, A. (eds) _Infrasound

Monitoring for Atmospheric Studies_, 349–386, https://doi.org/10.1007/978-3-319-75140-59 (Springer, Cham, 2019). * Donaukurier.

https://www.donaukurier.de/lokales/ingolstadt/Brand-Raffinerie-Irsching-DKmobil-Bayernoil-Katastrophe-Ursache-der-Explosion-steht-noch-nicht-fest;art599,3993333, last accessed August 2019

(2019). * Buttkus, E. B. (ed.) _Ten Years of the Gräfenberg Array_, http://www.schweizerbart.de//publications/detail/isbn/9783510960514/Geologisches_Jahrbuch_Reihe_E_Heft (Schweizerbart

Science Publishers, Stuttgart, Germany, 1986). * International Seismological Centre. On-line bulletin, http://www.isc.ac.uk, internatl.seismol. cent., thatcham, united kingdom,

https://doi.org/10.31905/D808B830 (2016). * Blom, P. & Waxler, R. Impulse propagation in the nocturnal boundary layer: Analysis of the geometric component. _J. Acoust. Soc. Am._ 131,

3680–3690, https://doi.org/10.1121/1.3699174 (2012). Article  ADS  PubMed  Google Scholar  * Blom, P. Interaction of the Cyclonic Winds with the Infrasonic Signal Generated by a Large

Maritime Storm. Dissertation, University of Mississippi, ProQuest/UMI (UMI No. 3567512) (2013). * Picone, J. M., Hedin, A. E., Drob, D. P. & Aikin, A. C. Nrlmsise-00 empirical model of

the atmosphere: Statistical comparisons and scientific issues. _J. Geophys. Res. Space Phys._ 107, 1468, https://doi.org/10.1029/2002JA009430 (2002). Article  ADS  CAS  Google Scholar  *

Drob, D. P. _et al_. An update to the horizontal wind model (hwm): The quiet time thermosphere. _Earth Space Sci._ 2, 301–319, https://doi.org/10.1002/2014EA000089 (2015). Article  ADS 

Google Scholar  * Shani-Kadmiel, S., Assink, J. D., Smets, P. S. M. & Evers, L. G. Seismoacoustic coupled signals from earthquakes in central italy: Epicentral and secondary sources of

infrasound. _Geophys. Res. Lett._ 45, 427–435, https://doi.org/10.1002/2017GL076125 (2017). Article  ADS  Google Scholar  * Negraru, P. T., Golden, P. & Herrin, E. Infrasound propagation

in the “zone of silence”. _Seismol. Res. Lett._ 81, 614–624, https://doi.org/10.1785/gssrl.81.4.614 (2010). Article  Google Scholar  * Koch, K. & Pilger, C. Investigation of

seismo-acoustic signals from the baumgarten and ingolstadt explosions in December 2017 and september 2018. _Geophys. Res. Abstr. EGU Gen. Assembly 2019_ 21, EGU2019–9144–1 (2019b). Google

Scholar  * Arrowsmith, S. J., Johnson, J. B., Drob, D. P. & Hedlin, M. A. H. The seismoacoustic wavefield: A new paradigm in studying geophysical phenomena. _Rev. Geophys._ 48, RG4003,

https://doi.org/10.1029/2010RG000335 (2010). Article  ADS  Google Scholar  * Whitaker, R. W. & Mutschlecner, J. P. A comparison of infrasound signals refracted from stratospheric and

thermospheric altitudes. _J. Geophys. Res. Atmospheres_ 113, D08117, https://doi.org/10.1029/2007JD008852 (2008). Article  ADS  Google Scholar  * Edwards, W. N., Eaton, D. W., McCausland, P.

J., ReVelle, D. O. & Brown, P. G. Calibrating infrasonic to seismic coupling using the stardust sample return capsule shockwave: Implications for seismic observations of meteors. _J.

Geophys. Res._ 112, B10306, https://doi.org/10.1029/2006JB004621 (2007). Article  ADS  Google Scholar  * Ben-Menahem, A. & Singh, S. J. _Seismic Waves and Sources_. (Springer, New York,

NY, 1981). * de Angelis, S., Fee, D., Haney, M. & Schneider, D. Detecting hidden volcanic explosions from mt. cleveland volcano, Alaska with infrasound and ground-coupled airwaves.

_Geophys. Res. Lett._ 39, L21312, https://doi.org/10.1029/2012GL053635 (2012). Article  ADS  Google Scholar  * Tauzin, B., Debayle, E., Quantin, C. & Coltice, N. Seismoacoustic coupling

induced by the breakup of the 15 february 2013 chelyabinsk meteor. _Geophys. Res. Lett._ 40, 3522–3526, https://doi.org/10.1002/grl.50683 (2013). Article  ADS  Google Scholar  * Brown, P.

G., ReVelle, D. O., Tagliaferri, E. & Hildebrand, A. R. An entry model for the tagish lake fireball using seismic, satellite and infrasound records. _Meteorit. & Planet. Sci._ 37,

661–675, https://doi.org/10.1111/j.1945-5100.2002.tb00846.x (2002). Article  ADS  CAS  Google Scholar  * Langston, C. A. Seismic ground motions from a bolide shock wave. _J. Geophys. Res._

109, B12309, https://doi.org/10.1029/2004JB003167 (2004). Article  ADS  Google Scholar  * Fuchs, F. _et al_. Site selection for a countrywide temporary network in austria: noise analysis and

preliminary performance. _Adv. In Geosci._ 41, 25–33, https://doi.org/10.5194/adgeo-41-25-2015 (2015). Article  ADS  Google Scholar  * Fuchs, F., Kolínský, P., Gröschl, G. & Bokelmann,

G. & the AlpArrayWorking Group. Alparray in austria and slovakia: technical realization, site description and noise characterization. _Adv. Geosci._ 43, 1–13,

https://doi.org/10.5194/adgeo-43-1-2016 (2016). Article  Google Scholar  * Molinari, I. _et al_. Swiss-alparray temporary broadband seismic stations deployment and noise characterization.

_Adv. Geosci._ 43, 15–29, https://doi.org/10.5194/adgeo-43-15-2016 (2016). Article  Google Scholar  * Vecsey, L. _et al_. Data quality control and tools in passive seismic experiments

exemplified on the czech broadband seismic pool mobnet in the alparray collaborative project. _Geosci. Instrumentation, Methods Data Syst._ 6, 505–521, https://doi.org/10.5194/gi-6-505-2017

(2017). Article  ADS  Google Scholar  * le Pichon, A., Ceranna, L. & Vergoz, J. Incorporating numerical modeling into estimates of the detection capability of the ims infrasound network.

_J. Geophys. Res._ 117, D05121, https://doi.org/10.1029/2011JD016670 (2012). Article  ADS  Google Scholar  * Tanimoto, T. & Wang, J. Low-frequency seismic noise characteristics from the

analysis of co-located seismic and pressure data. _J. Geophys. Res. Solid Earth_ 123, 5853–5885, https://doi.org/10.1029/2018JB015519 (2018). Article  ADS  Google Scholar  * Chunchuzov, I.

_et al_. Study of the wind velocity-layered structure in the stratosphere, mesosphere, and lower thermosphere by using infrasound probing of the atmosphere. _J. Geophys. Res. Atmospheres_

120, 8828–8840, https://doi.org/10.1002/2015JD023276 (2015). Article  ADS  Google Scholar  * Jones, R. M. & Bedard, A. J. Infrasonic ray tracing applied to small-scale atmospheric

structures: Thermal plumes and updrafts/downdrafts. _The J. Acoust. Soc. Am._ 137, 625–632, https://doi.org/10.1121/1.4906175 (2015). Article  ADS  PubMed  Google Scholar  * Sutherland, L.

C. & Bass, H. E. Atmospheric absorption in the atmosphere up to 160 km. _The J. Acoust. Soc. Am._ 115, 1012–1032, https://doi.org/10.1121/1.1631937 (2004). Article  ADS  CAS  Google

Scholar  * le Pichon, A. _et al_. Comparison of co-located independent ground-based middle atmospheric wind and temperature measurements with numerical weather prediction models. _J.

Geophys. Res. Atmospheres_ 120, 8318–8331, https://doi.org/10.1002/2015JD023273 (2015). Article  ADS  Google Scholar  * Gibbons, S., Kvaerna, T. & Nasholm, P. Characterization of the

infrasonic wavefield from repeating seismo-acoustic events. In le Pichon, A., Blanc, E. & Hauchecorne, A. (eds) _Infrasound Monitoring for Atmospheric Studies_,

https://doi.org/10.1007/978-3-319-75140-510 (Springer, Cham, 2019). * AlpArray Seismic Network. Alparray seismic network (aasn) temporary component, alparray working group, datacite link:

http://data.datacite.org/10.12686/alparray/z3_2015 -project webpage: http://www.alparray.ethz.ch/home, https://doi.org/10.12686/alparray/z32015 (2015). * Krischer, L. _et al_. Obspy: a

bridge for seismology into the scientific python ecosystem. _Comput. Sci. & Discov._ 8, 014003, https://doi.org/10.1088/1749-4699/8/1/014003 (2015). Article  ADS  Google Scholar  * The

ObsPy Development Team. (2017 October 27) ObsPy 1.1.0, Zenodo, https://doi.org/10.5281/zenodo.165135 (2017). * The ObsPy Development Team. (2019 February 17) ObsPy 1.1.1, Zenodo,

https://doi.org/10.5281/zenodo.1040770 (2019). * Royal Observatory of Belgium. Belgian seismic network. international federation of digital seismograph networks. other/seismic network,

https://doi.org/10.7914/SN/BE (1985). * Department of Earth and Environmental Sciences, Geophysical Observatory, University of Munchen. Bayernnetz, international federation of digital

seismograph networks, other/seismic network, https://doi.org/10.7914/SN/BW (2001). * Swiss Seismological Service (SED) at ETH Zurich. National seismic networks of switzerland, eth zurich,

https://doi.org/10.12686/sed/networks/ch (1983). * RESIF. Resif-rlbp french broad-band network, resif-rap strong motion network and other seismic stations in metropolitan france, resif -

reseau sismologique & geodesique francais, https://doi.org/10.15778/RESIF.FR (1995). * GEOFON Data Centre. Geofon seismic network, deutsches geoforschungszentrum gfz,

https://doi.org/10.14470/TR560404 (1993). * Federal Institute For Geosciences And Natural Resources (BGR). German regional seismic network (grsn), https://doi.org/10.25928/mbx6-hr74 (1976).

* University of Genova. Regional seismic network of north western italy, international federation of digital seismograph networks, https://doi.org/10.7914/SN/GU (1967). * Kovesligethy Rado

Seismological Observatory (Geodetic and Geophysical Institute, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences (MTA CSFK GGI KRSZO)). Hungarian national

seismological network, deutsches geoforschungszentrum gfz, https://doi.org/10.14470/UH028726I (1988). * Albuquerque Seismological Laboratory (ASL)/USGS. Global seismograph network (gsn -

iris/usgs), international federation of digital seismograph networks, https://doi.org/10.7914/SN/IU (1988). * INGV Seismological Data Centre. Rete sismica nazionale (rsn), istituto nazionale

di geofisica e vulcanologia (ingv), Italy, https://doi.org/10.13127/SD/X0FXnH7QfY (1997). * MedNet project partner institutions. Mediterranean very broadband seismographic network (mednet),

istituto nazionale di geofisica e vulcanologia (ingv), italy, https://doi.org/10.13127/SD/fBBBtDtd6q (1988). * OGS (Istituto Nazionale di Oceanografia e di Geofisica Sperimentale) and

University of Trieste. North-east italy broadband network, international federation of digital seismograph networks, https://doi.org/10.7914/SN/NI (2002). * ZAMG - Zentralanstalt für

Meterologie und Geodynamik. Austrian seismic network. international federation of digital seismograph networks. other/seismic network, https://doi.org/10.14470/TR560404 (1987). * OGS

(Istituto Nazionale di Oceanografia e di Geofisica Sperimentale). North-east italy seismic network, international federation of digital seismograph networks, https://doi.org/10.7914/SN/OX

(2016). * University of Trieste. Friuli venezia giulia accelerometric network, international federation of digital seismograph networks, https://doi.org/10.7914/SN/RF (1993). * ESI SAS

(Earth Science Institute Of The Slovak Academy Of Sciences). National network of seismic stations of slovakia, deutsches geoforschungszentrum gfz, https://doi.org/10.14470/FX099882 (2004). *

Slovenian Environment Agency. Trentino seismic network, international federation of digital seismograph networks, https://doi.org/10.7914/SN/SL (2001). * Geological Survey-Provincia

Autonoma di Trento. Trentino seismic network, international federation of digital seismograph networks, https://doi.org/10.7914/SN/ST (1981). * Leipzig University. Sxnet saxon seismic

network. international federation of digital seismograph networks, https://doi.org/10.7914/SN/SX (1988). * Institut fuer Geowissenschaften, Friedrich-Schiller-Universitaet Jena. Thüringer

seismisches netz (tsn). international federation of digital seismograph networks. other/seismic network, https://doi.org/10.7914/SN/TH (2009). Download references ACKNOWLEDGEMENTS This work

was funded by the Austrian Science Fund FWF through projects number P26391 and P30707. We thank an anonymous reviewer and Karl Koch for his very careful and insightful reviews which greatly

improved the manuscript. We are grateful to ORFEUS and EIDA for providing the tools to access the seismic waveform data. We greatly appreciate the fruitful discussions about our findings at

the CTBTO Infrasound Technology Workshop 2018 in Vienna, Austria. We acknowledge the use of data from the AlpArray network - please visit the project homepage

http://www.alparray.ethz.ch/en/seismic_network/backbone/data-policy-and-citation/ (last accessed August 2019) for a full list of people contributing to the AlpArray seismic network. AlpArray

seismic network team list: The AlpArray seismic network was made possible by the AlpArray Working Group: György Hetényi, Rafael Abreu, Ivo Allegretti, Maria-Theresia Apoloner, Coralie

Aubert, Simon Besançon, Maxime Bés de Berc, Götz Bokelmann, Didier Brunel, Marco Capello, Martina Čarman, Adriano Cavaliere, Jérôme Chéze, Claudio Chiarabba, John Clinton, Glenn Cougoulat,

Wayne C. Crawford, Luigia Cristiano, Tibor Czifra, Ezio D’alema, Stefania Danesi, Romuald Daniel, Anke Dannowski, Iva Dasović, Anne Deschamps, Jean-Xavier Dessa, Cécile Doubre, Sven Egdorf,

ETHZ-SED Electronics Lab, Tomislav Fiket, Kasper Fischer, Wolfgang Friederich, Florian Fuchs, Sigward Funke, Domenico Giardini, Aladino Govoni, Zoltán Gráczer, Gidera Gröschl, Stefan

Heimers, Ben Heit, Davorka Herak, Marijan Herak, Johann Huber, Dejan Jarić, Petr Jedlička, Yan Jia, Hélène Jund, Edi Kissling, Stefan Klingen, Bernhard Klotz, Petr Kolínský, Heidrun Kopp,

Michael Korn, Josef Kotek, Lothar Kühne, Krešo Kuk, Dietrich Lange, Jürgen Loos, Sara Lovati, Deny Malengros, Lucia Margheriti, Christophe Maron, Xavier Martin, Marco Massa, Francesco

Mazzarini, Thomas Meier, Laurent Métral, Irene Molinari, Milena Moretti, Anna Nardi, Jurij Pahor, Anne Paul, Catherine Péquegnat, Daniel Petersen, Damiano Pesaresi, Davide Piccinini, Claudia

Piromallo, Thomas Plenefisch, Jaroslava Plomerová, Silvia Pondrelli, Snježan Prevolnik, Roman Racine, Marc Régnier, Miriam Reiss, Joachim Ritter, Georg Rümpker, Simone Salimbeni, Marco

Santulin, Werner Scherer, Sven Schippkus, Detlef Schulte-Kortnack, Vesna Šipka, Stefano Solarino, Daniele Spallarossa, Kathrin Spieker, Josip Stipčević, Angelo Strollo, Bálint Süle,

Gyöngyvér Szanyi, Eszter Szűcs, Christine Thomas, Martin Thorwart, Frederik Tilmann, Stefan Ueding, Massimiliano Vallocchia, Luděk Vecsey, René Voigt, Joachim Wassermann, Zoltán Wéber,

Christian Weidle, Viktor Wesztergom, Gauthier Weyland, Stefan Wiemer, Felix Wolf, David Wolyniec, Thomas Zieke, Mladen Živčić, Helena Žlebčíková. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS

* University of Vienna, Department of Meteorology and Geophysics, Vienna, 1090, Austria Florian Fuchs, Petr Kolínský & Götz Bokelmann * Helmholtz Centre Potsdam, GFZ German Research

Centre for Geosciences, Potsdam, 14473, Germany Felix M. Schneider * University of Innsbruck, Department of Atmospheric and Cryospheric Sciences, Innsbruck, 6020, Austria Stefano Serafin

Authors * Florian Fuchs View author publications You can also search for this author inPubMed Google Scholar * Felix M. Schneider View author publications You can also search for this author

inPubMed Google Scholar * Petr Kolínský View author publications You can also search for this author inPubMed Google Scholar * Stefano Serafin View author publications You can also search

for this author inPubMed Google Scholar * Götz Bokelmann View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS F.F. performed the seismic and

seismo-acoustic waveform analysis and mainly wrote the manuscript. F.M.S. compiled the atmospheric model and computed the acoustic ray tracing and contributed to the manuscript. P.K.

analyzed seismic waveform data and determined the origin time. S.S. provided the ECMWF model and atmospheric interpretations. G.B. supervised the study, initiated research ideas and

contributed to the manuscript. All authors reviewed the manuscript. CORRESPONDING AUTHOR Correspondence to Florian Fuchs. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no

competing interests. ADDITIONAL INFORMATION PUBLISHER’S NOTE: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

SUPPLEMENTARY INFORMATION ELECTRONIC SUPPLEMENT RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use,

sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative

Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated

otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds

the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and

permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Fuchs, F., Schneider, F.M., Kolínský, P. _et al._ Rich observations of local and regional infrasound phases made by the AlpArray seismic

network after refinery explosion. _Sci Rep_ 9, 13027 (2019). https://doi.org/10.1038/s41598-019-49494-2 Download citation * Received: 15 January 2019 * Accepted: 27 August 2019 * Published:

10 September 2019 * DOI: https://doi.org/10.1038/s41598-019-49494-2 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a

shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative