Low-Frequency Atmospheric Gravity Waves from Vertical Tectonic Deformation During Two Recent Chilean Megathrust Events : the 2010 Maule ( Mw 8 . 8 ) , and 2014 Iquique ( Mw 8 . 2 ) Earthquakes

RESEARCH ARTICLE Low-Frequency Atmospheric Gravity Waves from Vertical Tectonic Deformation During Two Recent Chilean Megathrust Events: the 2010 Maule (Mw8.8), and 2014 Iquique (Mw8.2) Earthquakes Takeshi Mikumo, Takuo Shibutani, Makiko Iwakuni and Nobuo Arai Kyoto University, Uji, Kyoto 611-0001, Japan, Disaster Prevention Research Institute, Kyoto University, Uji, Kyoto 611-0001, Japan Japan Weather Association, Tokyo 170-6045, Japan Disaster Mitigation Research Center, Nagoya University, Nagoya 464-8601, Japan


INTRODUCTION
The 2010 earthquake (Mw = 8.8) took place off the west coast of the Maule region in southwestern Chile (-35.909°N,-72.733°E, due to Universidad de Chile) on February 27, 2010, 06:34 UTC, as shown in Fig. (1a).This region is located just north to the source region of the giant 1960 Valdivia, Chile earthquake (Mw = 9.5).Far-field seismic observations [1 -5] provided the general features of the fault rupture process, and various geodetic measurements including a radar technique clarified the spatial extent of the rupture zone [6 -10].These various measurements including GPS, InSAR, and near-field land-level changes clearly suggested that the rupture zone over the Nazca-South America plate boundary extends for about 500 km in the NNE-SSW direction nearly parallel to the trench axis and sea-shore line, and for about 200 km in the E-W direction.Synthetic inversion of all these data has suggested a maximum coseismic slip of 16 m on the dipping fault plane in its northern segment and about 10 m over two patches in its southern segment, and this slip pattern yields maximum ground uplift of 4.0 m and subsidence of 0.85 m at the two sides of the rupture zone [10].The second 2014 Iquique earthquake (Mw = 8.2) occurred also in the west coast of northern Chile (-19.642°N,-70.817°E at a depth of 21.9 km, due to USGS), as also shown in Fig. (1a).It took on April 1, 2014, 23:46 UTC, about 1600 km north of the 2010 Maule earthquake.This region has been identified as the northern Chile seismic gap [11], [12], since a megathrust earthquake of 1877 with a magnitude of 8.7 -8.9 [13].Aftershocks of this mainshock event extend nearly parallel to the trench axis for about 150-200 km.Teleseismic P-waves and W-phase inversions from global seismic recordings suggested a maximum slip of ~ 2 -6.7 m on a thrust fault plane dipping with 14°, indicating a dip of the subducting plate in this region [14], while the inversion of Tsunami data suggested a slip of 5 -7 m around a depth of 23 km in a more concentrated area of 60 km by 40 km [15].This mechanism should have caused some amount of uplift and subsidence near the bottom of the oceanic trench and a part of the inland region, although their amounts have not been estimated from geodetic surveys, unlike the case of the Maule earthquake.
With these circumstances, low-frequency atmospheric pressure waves were identified during 6.5 hours and 4.7 hours after the first and second great earthquakes, respectively, by sensitive microbarographs at several global International Monitoring System (IMS) stations located in South America and its surroundings, as shown in 1b.It is to be mentioned here that similar observations of low-frequency atmospheric waves have been recorded after three megathrust events, since the 1964 Alaskan earthquake (Mw = 9.2) [16,17], the 2004 Sumatra-Andaman earthquake (Mw = 9.2) [18,19], and the 2011 great Off-Tohoku earthquake (Mw = 9.0) [20,21].In the present article, we try to find the evidence of observations for low-frequency gravity waves propagating from each of the source regions of the above two great earthquakes through the atmosphere to some of these microbarograph stations.We also calculate synthetic waveforms, incorporating various seismic source characteristics described above and based on a standard atmospheric sound velocity structure up to a height of 220 km.Our main purpose here is to demonstrate that the observed atmospheric gravity waves have actually been generated from some part of the source regions of the two great earthquakes, and also try to elucidate, if possible, the overall source characteristics for exciting the low-frequency atmospheric waves from these megathrust events.

OBSERVATIONS OF PROPAGATING ATMOSPHERIC WAVES
The IMS stations shown in Fig. (1b) that recorded low-frequency atmospheric pressure waves from the two earthquakes are summarized in Tables (1 and 2), respectively.These stations are located in South America and its surroundings in the distance range from the epicenter between 1,813 km and 7,684 km in the case of the 2010 Maule earthquake Fig. (2a), and between 1,569 km and 6,481 km in the case of the 2014 Iquique earthquake Fig. (2b).Four of these stations are located within the South America continent, another four are on islands in the Atlantic Ocean northeast or east of it, and the other three are located on islands in the southwestern Pacific Ocean west of the continent.However, observed records at some of these stations are perturbed by quite large atmospheric noise and probably by high oceanic waves.Accordingly, in the present analysis, we selected each four stations with relatively less perturbed and coherent waves from the two events.27,2010.The p-p amplitudes are on the filtered records for periods between 14.0 min and 4.5 min.Arrival times are the approximate arrival of waves on the above filtered records, and a question mark means uncertain data for its time, and not adopted here.Record interval between the initial and end times is 7,200 sec (120 min) on all records.In the upper two traces of (Figs. 3, 4), are shown each two components of the original records of pressure perturbation over a time interval of 7,200 sec (2 hours), which have been observed at IS08 station during the Maule earthquake, and at IS02 station during the Iquique event, respectively.The long time length has been taken to look at how are atmospheric noise and possible arrival of target signals during the two hours.Because these records include high-frequency noise and long-term inclinations, as shown in each of these upper traces, we apply a second-order Butterworth bandpass-filter to each two channels of the above original records, with a unit amplitude over a frequency range between 1.19 and 3.70 mHz (or between 14.0 and 4.5 min), which could cover the frequency range of gravity waves propagating in the lower atmosphere [22].The above frequency range does not cover the propagating acoustic modes, because the modes are heavily contaminated by various noise in the present case.These filtered traces are aligned in the lower two traces of each of the above figures.The next step is to compare above filtered channels at each station with the synthetic waveforms given in Section 4, in order to discuss the source characteristics of the above two megathrust earthquakes, after estimating approximate propagation phase velocity of the gravity waves radiated from these events.

WAVEFORM MODELING
We perform waveform modeling for the above filtered records of pressure perturbation, to confirm that the observed waves have been actually excited somewhere within the extensive source regions of these great earthquakes, and propagated through the lower to middle atmosphere to the observation stations.To do this, we incorporate such a realistic thermal structure as the US Standard Atmosphere, formerly called the ARDC atmospheric thermal model.The structure of the atmosphere is composed mainly of four layers.The temperature decreases first up an altitude of 15 km in the troposphere, increases up to 45 -50 km in the stratosphere, then decreases again at 80 -90 km in the mesosphere, and finally the temperature increases rapidly in the thermosphere up to 220 km.Although somewhat modified models, such as CIRA [23], MISISE [24], and NRL-C2Z [25] have been presented to date, we refer to the US standard structure as a more global scale structure to discuss lower frequency waves.For the ARDC standard atmosphere, the phase and group velocities of several acoustic and gravity modes have been calculated [26,27], where the air density is assumed to decrease exponentially with altitude.In addition, Harkrider [22]  atmosphere to a point source and a receiver both located on the ground surface, by approximating the above thermal structure by 39 horizontally stratified, isothermal layers terminated either with the isothermal half-space or with the free surface at an altitude of 220 km.In the present article, we assume the former case in view of our previous test calculations [18].
The propagation of acoustic-gravity waves in the lower to middle atmosphere can be described by the following atmospheric transfer function, A(ω) exp[− iϕ A (ω)], where, A Aj (ω) is the dynamic response of the atmosphere to a surface point source and receiver, k is the wave number, r is the distance from the source to the receiver, and C j (ω) is the phase velocity, with j indicating the mode number.These functions have been calculated for several acoustic and gravity modes [18 -21], and a few examples of the Green's functions from a point source to a far-field station have been shown in the case of the 2004 Sumatra-Andaman earthquake [18].Beside the atmospheric transfer function, the observed waveforms should include the effects of the source time function, the dimension of the source region including the velocity of its expansion, and the instrumental response as indicated below.
If we assume a linear system coupled between the solid earth, the sea or land surface, and the atmosphere, then the pressure perturbation that would be observed at any station in the far-field, compared with the source dimension and the wavelength, can be written in the time domain as in the following equations [18,21], where S(ω), D(ω), A(ω), and B(ω) are the source time spectrum, the source finiteness spectrum, the atmospheric transfer function defined as in equation (1), and the barograph response, respectively.ϕ (ω)s are their phase spectra, respectively.Here, c = 2(2/π) 1/2 is simply a numerical constant derived from Harkrider's formulations [22], and (r/Rsinθ) 1/2 is the approximate curvature correction for the effect of energy spreading over a spherical earth instead of a flat surface, with the distance r to the recording station from the source, the Earth's radius R, and the spherical colatitude of the station θ.S(ω) for pressure perturbation includes the sound velocity c and air density ρ near the ground surface and the rise time of vertical displacement τ.D(ω) includes the source dimensions L i and W i of the assumed rectangular source, the expanding velocity ν, and β the azimuth from the source to the station with respect to the direction of source expansion.These functions, including the parameters not explicitly mentioned here, have been calculated for various cases, and are shown in [18,21].

COMPARISON BETWEEN THE OBSERVED AND SYNTHETIC WAVEFORMS
Our next step is to incorporate the various source parameters of these great earthquakes into the calculation of synthetic waveforms, and then to compare them with the corresponding observed waveforms of pressure perturbation from the two different earthquakes, respectively.This is to elucidate the generation mechanism of low-frequency atmospheric pressure waves and also the possible features of coseismic vertical crustal deformation associated with these megathrust earthquakes.In this sense, we compare the two megathrust earthquakes that took place in the same subduction zone over the Nazca-South America plate boundary.

2010 Maule Earthquake
For this earthquake, the location, the entire source dimension, and total coseismic slip obtained from three different types of geodetic measurements, are shown in Moreno et al. [10].The overall source dimension extending for about 500 km in the NNE-SSW direction nearly parallel to the trench axis may be separated into the northern and southern segments with respect to the location of the epicenter, and the width of about 200 km in the E-W direction may be separated by two zones of ground uplift and subsidence, respectively, as shown in Fig. (2a).We denote these four zones as U1, U2, D1, and D2, respectively, which are specified with their approximate lengths L1 = 220 km and L2 = 280 km and their widths W1 = 150 km and W2 = 50 km.Average tectonic ground uplifts in the two zones enclosed by L1 and W1 and by L2 and W1 are successively varied within possible values between 0.5 m and 3.5 m, as indicated in Table (3).Average subsidences within two zones enclosed by L1 and W2 and by L2 and W2 are fixed as -1.0 m and -0.5 m.Also assumed is the average velocity of source spreading from the epicenter as 2.5 km/s.We confirmed by a number of trial calculations that these fixed values did not have large effects on the synthetic waveforms, as shown in our previous studies [18,21].Some of the above estimates are taken by referring to the geodetic measurements given in Table S2 of Moreno et al. [10].The rest of the important parameter to be estimated is the rise time of vertical tectonic displacement in the source region, which could directly affect the calculated waveforms.These parameters are varied successively between 1.5 and 3.5 min in the zones of uplift, while fixing it as 2.0 min in the zone of subsidence as a first approximation.The synthetic waveforms for 5 possible cases with different rise times have been calculated, and the results are given for station IS41 in an additional PDF file [see Graph 1].The five cases in Syn2 to Syn7 give generally similar waveforms, expecting that a good fit between the filtered and synthetic waveforms would be obtained at this station.3 in most of these cases, whereas there remain some discrepancies between the filtered waveforms and the synthetics S2 (from Case Syn7 in Table 3 in these figures.This suggests that these latter cases S2 do not provide satisfactory explanations to the observations.From the first case, however, we would be able to infer that the average uplift in the northern and southern segments may be approximately 3.0 m and 1.5 m, respectively, and the corresponding rise time in the two zones may be around 3.0 min (See Table 3).There are some differences in their absolute amplitudes, however, for which we believe they are attributed to prevailing atmospheric conditions around these stations.We would say, however, that there are more or less generally satisfactory agreements between the filtered and synthetic waveforms S1 composed mainly of gravity modes, at least, for a limited time interval around the first part of about 22 min.The waveforms S1 in (Fig. 5b) provides the best case.

2014 Iquique Earthquake
For the 2014 Iquique earthquake, seismic source parameters have been estimated from various observations up to this time, from aftershock distribution, teleseismic waves and Tsunami observations, by several investigators mentioned below.The overall source dimension appears to extend for about 150 km [2] to 200 km (USGS) in the N-S direction inferred from aftershock distribution, which appears to be almost parallel to the trench axis and the shore line in this region.Its width in the E-W direction may be about 100 km (USGS-NEIC).Preliminary inversions of teleseismic [2] and Tsunami [14] observations, and Tsunami and GPS observations [15], yield approximate spatial slip distribution on its dipping fault, but coseismic vertical displacements on land have not been estimated from geodetic surveys.We assume the maximum length of the uplifted zone as 200 km along the NNW-SSE direction (L1 + L2) with its maximum width (W1) to be 100 km beneath the oceanic region, and its maximum width of a possible subsided zone as 50 km (W2), as shown in Fig. (2b).On the other hand, inversion of Tsunami data suggested a more concentrated source area of 60 km by 40 km near the hypocenter, in the total dimension of 200 km × 160 km, with a maximum slip of 5 to 7 m at a depth around 23 km [14].We also assume the rupture starting point at two different positions; one is located at the northwestern periphery of the source region (U1), and the other is assumed to lie at the northwestern central position of the entire source region.For these two possible cases, we refer them to unilateral (U) and bilateral (B) rupture extension, respectively, with a constant rupture velocity of 2.5 km/sec.We calculate synthetic waveforms for the two assumed cases.For the assumed bilateral case (B), we denote the divided four zones, U1, U2, D1, and D2, similar to the case for the 2010 Maule earthquake, with L1 = L2 = 100 km, W1 = 100 km and W2 = 50 km.For these subdivided zones, we assign average uplifts in the two zones enclosed by L1and W1 and by L2 and W1, as 3 m and 2 m, respectively, and average subsidence in the zones enclosed by L1 and W2 and by L2 and W2 are assumed as -1.0 and -0.5 m, respectively, similar to the case of the Maule earthquake.The maximum positive displacement of the initial water elevation has been estimated as 1.54 m from the inversion analysis of Tsunami data [14].The rise times of vertical displacement in these zones are assumed here as 3.0, 2.0, 2.0, and 2.0 min, respectively.For the assumed unilateral case (U), on the other hand, L1 = 200 km and L2= 0 km, and the other parameters are fixed as in the same as in the case of (B).All these assumed parameters are given in Table 4. 2.0 2.0 2.0 2.0 Syn40 -Syn44 are assumed for unilateral cases, while Syn51 -Syn53 are assumed for bilateral cases, respectively.D1U and D2U indicate the average uplifted displacements at zones U1 and U2, and D1D and D2D indicate the average subsided displacements at zones D1 and D2, respectively.T1U, T2U, T1D and T2D stand for the assigned rise time τ s in these cases.Mark ** (Syn40) indicates a most probable case, and * (Syn41) follows this.The bilateral cases are assumed for comparison, and not finally adopted.

Fig. (6a).
Filtered and synthetic waveforms at station IS41 for the case of the 2014 Iquique earthquake.The lower two traces (SUL and SBL) indicate the corresponding synthetic waveforms from Syn40 (unilateral case) and Syn51 (bilateral case), given in (Table 4).The time shown there indicates the absolute time UTC, (00:25-01:25).The third (SUL) and fourth traces (SBL) in these figures indicate the corresponding synthetic waveforms for two different cases for assumed unilateral (U) and bilateral (B) rupture extension, respectively.Detailed comparisons between the filtered and synthetic waveforms are described in the figure captions.It appears that the filtered waveforms are more or less generally well simulated, by the synthetic waveforms (SUL) for the assumed unilateral case (U) from Syn40, at least, for a limited time interval of the first part of 22 min or so.The waveforms SUL in Fig. (6b) shows the best case.On the other hand, this is not the case (SBL) for the bilateral rupture (B) from Syn51.Comparing the results for the two cases (U) and (B), it may be concluded that the unilateral case (U) in which the rupture spreads southwards from the northwestern periphery of the source region, appears to be more likely for the 2014 earthquake, unlike the case of the 2010 Maule event.

RESULTS AND DISCUSSION
As described in the previous section, we calculated synthetic waveforms to be observed at each four IMS stations, during two Chilean megathrust events, the 2010 Maule and 2014 Iquique earthquakes, respectively, incorporating their several probable source parameters, together with the standard atmospheric temperature structure up to a height of 220 km.The comparison between these synthetics and the corresponding low-frequency waveforms observed for gravity waves appears to provide satisfactory agreement, at least, for their initial time interval of about 22 min or so.There still remain, however, some discrepancies between their absolute amplitudes.There leaves a possibility that some of the discrepancy might be due to our inadequate estimate of the source parameters, including average vertical ground displacements and its time constants in adopted or assumed source models.Another possible reason might be attributed to prevailing atmospheric conditions during the time of these earthquakes which could be somewhat deviated from the adopted standard atmospheric model.For these possible reasons, we did not make further attempts to fit exactly the calculated amplitudes to their corresponding observations.We would like to mention here that these are the fourth and fifth megathrust earthquakes, from which lowfrequency atmospheric waves including gravity modes have been emitted, since the first, unusual observation after the 1964 Alaskan earthquake (Mw = 9.2) [16,17], the Sumatra-Andaman earthquake (Mw=9.2) [18,19], and the 2011 Off-Tohoku earthquake (Mw=9.0)[21].
It is to be noteworthy that the present 2010 Maule, and 2014 Iquique earthquakes in Chile, even with slightly smaller magnitudes of Mw = 8.8 and 8.2 as compared with the three previous megathrust earthquakes, appear to have also generated low-frequency gravity waves with somewhat smaller amplitudes that could propagate up to long distances.This is the reason why the two Chilean earthquakes in 2010 and 2014 also yielded quite large-scale uplift and subsidence of the sea bottom over a wide area in the Nazca-South America subduction zone.
Previous theoretical studies on Tsunami generation [28] and numerical simulations [29] clearly show that if the wavelength of the sea-bottom deformation is much longer than the water depth and if the deformation takes place in a few minutes or so, then the swelling and depression of the sea-surface behaves almost exactly like the sea-bottom deformation.These conditions are met exactly in the source regions of the two great earthquakes off the western coast of South America mentioned here.The rapid change of the sea-surface could excite atmospheric disturbances just above it.For a part of land area, if extended to some part of the source region, the ground surface would provide more direct upward pressure increase into the air.In any case, the pressure perturbation of the air directly above it has been estimated as less than a few to several Pa or even less, if the density of the air near the sea or ground surface, ρ =1.193 x 10 -3 gr/cm 3 , the sound velocity c =320 m/s near the ground surface, and also their dissipation effects during propagation are taken into account.The disturbances over a wide source region favor the condition that low-frequency atmospheric waves could be efficiently excited and propagated through the lower to middle atmosphere to long distances.The frequency response of the sensors used at most of microbarograph stations, particularly at IMS stations, however, drops rapidly towards lower frequencies.Nevertheless, it was barely successful to observe atmospheric gravity waves this time.It is hoped in the future that the response could extend to cover lower frequencies to be able to observe lowfrequency acoustic and gravity waves to reveal the lower to upper atmospheric temperature structures as well as to discuss the source characteristics of large earthquakes from atmospheric pressure waves.For the above five megathrust earthquakes, there is a variety in their source dimension and the amount of vertical tectonic deformation depending on their dip angle of the fault plane.Nevertheless, the rise time of their vertical displacement, which is directly related to the amplitude and waveforms of acoustic-gravity waves, ranged between 1.5 and 3.5 min.This reflects partly the total rupture time of the entire fault plane, but might suggest the possible maximum time for tectonic deformation that could take place in subduction zones.
It is to be also mentioned here that after this Maule, Chile earthquake, Tsunami arrival has been reported at a number of seashore sites particularly along Chile, Peru, and Ecuador in South America [30].It has also been found that the Tsunami waves arrived at two IMS stations, Easter Islands and Marquesas Islands, both are located in southwestern Pacific region and included in Table (1).The Tsunami arrival was with a height of 0.35 -0.18 m at 3 hours and a half later than the gravity waves, indicating the Tsunami speed in these cases was between 220 and 239 m/sec.This means that the acoustic arrival might be used as early warning for Tsunami arrival [20,21], although the Tsunami amplitudes this time were not high enough to issue urgent warning at far-field stations.After the 2014 Iquique earthquake, on the other hand, Tsunami waves arrived at several sites within 4 hours along the western coast of Chile up to a height of 0.7 -2.0 m [31].

CONCLUSION
Low-frequency atmospheric gravity waves have been recorded after two recent Chilean megathrust events, during 6.5 hours after the 2010 Maule (Mw8.8), and during 4.7 hours after the 2014 Iquique (Mw8.2) earthquakes, respectively, with apparent phase velocity of about 319 m/s and 337 m/s at several microbarograph stations in South America and its surrounding regions.To compare with these observations, we constructed synthetic waveforms to be recorded at these stations, incorporating seismic source characteristics of these earthquakes and a sound velocity structure up to a height of 220 km.The comparison suggests that the first part of these waves, at least, within 22 min have been actually excited by vertical tectonic deformation at the source regions of the two megathrust events, respectively.From this comparison, the time constant of the vertical deformation appears to be in the range between 2 and 3 min.

Fig. (
Fig. (1a).Index map.South America continent, and the location of two source regions of the 2010 Maule earthquake (M) and the 2014 Iquique earthquake (Q).
Fig. (1b).Location of IMS stations in and around South America.

Fig. (
Fig. (2a).Source region of the 2010 Maule earthquake.The fault rupture started from the epicenter located in the mid-position.Notation MFZ indicates the location of the Mendana Fracture Zone off northern Chile.Other notations in the figure refer to the text.

(Table 1 )
Fig. (2b).Source region of the 2014 Iquique earthquake.The fault rupture is assumed to start from different positions; For the unilateral case (U), the rupture is assumed to start from the northern periphery (solid asterisk), while for the bilateral case (B), the rupture is assumed to start from the mid-position (open asterisk).Other notations in the figure refer to the text.

Fig. ( 4 ).
Fig. (4).An example of the recorded and filtered waveforms of pressure perturbation during the time (02:00-04:00) at an IMS station IS02 after the 2014 Iquique earthquake.The upper two traces (CH1 and CH2) provides the original records of pressure perturbation for the time interval of 7, 200 sec (the initial time indicated at the lowest left end), before and after the arrival of target atmospheric signals.The lower two traces (F1 & F2) are their filtered waveforms for the period range between 14 min and 4.5 min (or between 1.19 and 3.70 mHz) (See Chapter 3).
also calculated the dynamic response of the

Graph 1 .Table 3 .
Fig. (5a).Filtered and synthetic waveforms at station IS41 for the case of the 2010 Maule earthquake.Each of the upper two traces (F1 and F2) reproduces the filtered waveforms for the time interval of 3,600 sec (one hour).The lower two traces (S1 and S2) indicate the corresponding synthetic waveforms.The two synthetics S1 and S2 come from Syn7 and Syn2 in Table (3), respectively.The time shown below indicates the absolute time UTC, (a) IS41 (07:30-08:30).

Fig. (
Fig. (5b).Filtered and synthetic waveforms at station IS08 for the case of the 2010 Maule earthquake.Each of the upper two traces (F1 and F2) reproduces the filtered waveforms for the time interval of 3,600 sec (one hour).The lower two traces (S1 and S2) indicate the corresponding synthetic waveforms.The two synthetics S1 and S2 come from Syn7 and Syn2 in Table (3), respectively.The time shown below indicates the absolute time UTC, (b) IS08 (08:03:20-09:03:20).

Fig. (
Fig. (5d).Filtered and synthetic waveforms at all 4 stations for the case of the 2010 Maule earthquake.Each of the upper two traces (F1 and F2) reproduces the filtered waveforms for the time interval of 3,600 sec (one hour).The lower two traces (S1 and S2) indicate the corresponding synthetic waveforms.The two synthetics S1 and S2 come from Syn7 and Syn2 in Table(3), respectively.The time shown in each bracket below indicates the absolute time UTC, (d) IS21 ( 13:00:00-14:00:00).

Table 4 .
Source parameters incorporated in synthetic waveforms for the 2014 Iquique earthquake.

Table 2 . Location of IMS stations and summary of observations for the 2014 Iquique earthquake.
April 1, 2014.The p-p amplitudes are on the filtered records for periods between 14 min and 4.5 min.Arrival times are for approximate arrival on the filtered records, and a question mark means uncertain data, and not adopted here.Record interval between the initial and end times is 7200 sec (120 min) for all records.