arxiv: v1 [physics.geo-ph] 9 Oct 2012

Transcription

1 1 The Association of the Moon and the Sun with Large Earthquakes Lyndie Chiou 1, 1 ResearchPipeline.com, Union City, CA, USA Corresponding lyndie@researchpipeline.com arxiv: v1 [physics.geo-ph] 9 Oct 2012 Abstract The role of the moon in triggering earthquakes has been studied since the early 1900s. Theory states that as land tides swept by the moon cross fault lines, stress in the Earth s plates intensifies, increasing the likelihood of small earthquakes. This paper studied the association of the moon and sun with larger magnitude earthquakes (magnitude 5 and greater) using a worldwide dataset from the USGS. Initially, the positions of the moon and sun were considered separately. The moon showed a reduction of 1.74% (95% confidence) in earthquakes when it was 10 hours behind a longitude on earth and a 1.62% increase when it was 6 hours behind. The sun revealed even weaker associations (<1%). Binning the data in 6 hours quadrants (matching natural tide cycles) reduced the associations further. However, combinations of moon-sun positions displayed significant associations. Cycling the moon and sun in all possible quadrant permutations showed a decrease in earthquakes when they were paired together on the East and West horizons of an earthquake longitude (4.57% and 2.31% reductions). When the moon and sun were on opposite sides of a longitude, there was often a small (about 1%) increase in earthquakes. Reducing the bin size from 6 hours to 1 hour produced noisy results. By examining the outliers in the data, a pattern emerged that was independent of earthquake longitude. The results showed a significant decrease (3.33% less than expected) in earthquakes when the sun was located near the moon. There was an increase (2.23%) when the moon and sun were on opposite sides of the Earth. The association with earthquakes independent of terrestrial longitude suggests that the combined moon-sun tidal forces act deep below the Earth s crust where circumferential forces are weaker. Introduction There are many proposed triggers for earthquakes. Examples include underground temperature differences, solar activity and water movements both below and above ground (see, for example, [1 6]). Tides are especially intriguing as a potential trigger because they are induced by the gravitational pull of the moon and the sun which follow calculable orbits. Tides occur both in water and, to a lesser extent, in land. Land tides are a measurable effect wherein the Earth s crust rises and falls as a result of the moon s and sun s gravitational pull. The effect can be up to 20 cm near the Earth s poles [7]. Is it possible that some earthquakes could be triggered by land tides? It seems a reasonable hypothesis that the intersection of land tides with the Earth s fault lines could serve as a trigger for an earthquake. The reverse has proven true: moonquakes can be triggered by land tides caused by the Earth s gravitational influence on the moon [8]. However, the notion that the moon can similarly cause earthquakes is controversial as the Earth has 81 times more mass than the moon [9] and any influence is assumed to be very weak at best with many concluding the effect is non-existent (for example, see references [10, 11]). Despite the controversy, papers continue to be published on both sides of the debate (illustrated by [12 18]), many with conclusions based on minimally-sized datasets corresponding to specific faults. One noticeable exception was reference [19] which reported an analysis using 442,412 earthquakes from magnitude 2.5 to 9. The authors concluded that land tides induced by the moon do trigger earthquakes,

2 2 albeit primarily shallow earthquakes of low magnitude (a 0.5% to 1.0% increase over expected). The moon s greatest influence occurred when it was overhead, corresponding to a rise in the Earth s crust. A short-coming of the analysis was that low magnitude data, particularly below magnitude 4, was not uniformly recorded in the data and potentially lead to a bias in the conclusions. This paper uses a different dataset than the previous study and not only looks at the role of the moon and sun in triggering earthquakes separately, but extends the analysis to examine the role of the sun and moon in conjunction. The analysis is restricted to larger earthquakes, magnitude 5 and greater, since these magnitudes were accurately and comprehensively recorded across the globe. Also, larger magnitude earthquakes can exert a high negative economic impact, making the results more interesting outside of the geological community. Methods The position of the moon and sun as a proxy for land tides A suitable technique to deduce the relationship between land tides and earthquakes would require a dataset that allowed users to extract the land tide value for a given latitude and longitude on a particular date and time. Unfortunately, no such dataset exists. The ability to calculate land tides relies on a detailed map of subsurface properties such as ocean floor depths and underground plate tectonics [20]. For this reason, the longitudinal positions of the moon and sun with respect to earthquakes were used as a proxy for the position of land tides. This approach is not the most ideal since, as mentioned, the actual height of land tides varies not just according to the position of the moon and sun, but also due to localized geographical properties. The assumed correlation of the moon and sun versus land tide heights is an important source of error in the upcoming statistical analysis. This paper measures the occurrence of global earthquakes between 1973 to 2011 relative to the longitudinal positions of the sun and moon and compares the result to a simulation of randomly generated earthquakes. The data was downloaded from the United States Geological Survey (USGS) website [21]. Earthquakes of magnitude 5 and greater were extracted from a date range of January 1, 1973 to July 29, 2011 (the range available when the project was started). As mentioned earlier, lesser magnitude earthquakes, particularly in the magnitudes below 4, were not globally recorded as accurately as those of higher magnitudes. Therefore the analysis was limited to magnitude 5 and above which encompasses larger earthquakes. The resulting dataset consisted of 66,724 earthquakes located around the world. Figure 1 shows a visualization of each earthquake on a world map generated with the data. The following sections describe the method used for establishing the relative difference in position between the moon and sun and each earthquake as well as the histogram method used for the analysis. One Hour Binning To evaluate the possibility of the moon influencing earthquakes, it was necessary to calculate the relative difference in position between the moon and each earthquake. This was done by subtracting the right ascension (RA) coordinates (calculated using Mathematica s ephemerides tables which relied on [22]) for the moon from the longitudinal coordinate of each earthquake, converted to units of time. d rel = Longitude EQ (hrs) RA moon (hrs) (1) The result yielded the relative longitudinal difference between the Earth position of each earthquake and the corresponding moon position (latitude offsets were considered less important since the moon moves primarily circumferentially around the Earth).

3 3 Figure 1. Earthquakes from January 1, 1973 to July 29, 2011, magnitude 5 and greater. The color of the dots denotes magnitude, from blue (magnitude 5) to red (highest magnitude, 9.1). The longitudinal offsets were then binned according to relative time differences in hours. The left hand side of figure 2 shows the entire dataset of calculated longitudinal offsets in a traditional histogram format, binned by hour. The right hand side of figure 2 is an alternate visualization. The 24 wedges represent the 24 hours of possible relative longitudinal distances around the globe from a given earthquake. The height of the wedge reflects the number of earthquakes occurring during that offset hour. The visualization shows the geometrical relationship between the earthquake, the moon s RA coordinate and the number of earthquakes for each distance. Figure 2. Histograms of earthquake-moon longitudinal offsets binned by hour. The left hand side shows the traditional histogram format of the entire moon-earthquake data. The right hand side shows the results referenced to an example earthquake longitude. Note that the results are relative to earthquake longitudes in both histograms, instead of referenced to a specific longitude. The expected mean number of earthquakes/hour was 2780+/ An examination of the histograms shows slight variations in the number of earthquakes versus the

4 4 position of the moon. But does this imply a slight increase in the chance of an earthquake during those moon-earthquake offsets? In a noiseless dataset, the total number of earthquakes would be divided evenly among the 24 hours of possible longitudinal offsets. This would imply 2780 earthquakes per hour of relative offset with a standard deviation of 0 (during the dataset time period, 38.6 years). However, since the dataset is limited to 66,724 data points, it is necessary to run a simulation that generates the equivalent number of random earthquakes and compare the results to the real-world data to determine if the fluctuations bear any statistical significance. This simulation was performed by randomly generating longitudinal offsets occurring throughout the day across the same year range as the real-world data. The distribution along fault lines was implicit as the calculation only used relative distances. In that way we don t need to know the positioning of the fault lines. As mentioned, this was a source of error since local geography was not taken into account. In all, 70 sets of 66,724 earthquakes were generated to achieve a stable group of datasets. The stability was measured by calculating the standard deviation of the hourly bins from the histograms generated from n random datasets. Equations (2) and (3) show how the averaged standard deviation was computed. As more and more datasets were included, the averaged standard deviation converged. stdev bins = stdev( n bin hour ) (2) each hour 24 1 stdev avg = stdev bins (3) n Figure 3 shows the averaged standard deviation for the averaged hourly bins as more and more random datasets are included. 1 Figure 3. Stabilization of the standard deviation of the data. As the number of datasets of randomly generated earthquakes increases, the standard deviation of the height of the averaged bins stabilizes. Each dataset has a full 66,724 points. The resulting simulation predicted a mean of 2780 earthquakes with a standard deviation of ± 51.2 earthquakes. Finally, figure 4 shows the error estimate overlaid on top of the actual measured offset data. The gray overlapping ring shows ±1.96 (standard deviation) which implies a 95% confidence level (the threshold used throughout the paper). A first examination shows that indeed, some moon-earthquake offsets show

5 5 a slight increase (and decrease) compared to what would be expected from totally random earthquakes. A more detailed discussion follows in the Results section of this paper. Figure 4. Relative histograms between earthquake longitude and moon s (left) and sun s (right) longitude (mapped back to Earth). The left shows the same histogram data as figure 2 with an overlay of the expected variation (95% confidence) shown as a gray ring. The actual increase or decrease is listed as a percentage outside each wedge. The inner number shows the relative displacement in hours between the earthquake s and the moon s longitudes. The expected mean was 2,780 earthquakes +/ Six Hour Binning In addition to hourly binning, these methods were applied to six-hour longitudinal bins. Oceanic tides cycle through two high and two low tides lasting six hours over a period of a day. The high tides correspond to the moon sitting nearly overhead and below the observational longitude, while the low tides correspond to the moon near the horizons which are six hours away, hence the choice of six hour bins. While there is about a 1 hour lag between the position of the moon and oceanic tides, there is no such lag for land tides [23]. This allows the data to be re-binned into the following six hours increments which capture the entire phase of a tide: 21 hours to 3 hours (which is equivalent to -3 hours to 3 hours), 3 hours to 9 hours, 9 hours to 15 hours and 15 hours to 21 hours (see figure 5). Moon-Sun Position Combinations Rather than limiting the analysis to just the moon and sun individually, the relative positions of both the moon and sun in conjunction proved to be important. This is the final analytical method used on the data and provided the most interesting set of results.

6 6 Figure 5. The Earth divided into 4 quadrants reflecting the 4 phases of the tide referenced to an earthquake s longitude (shown with a red dot). The longitude indicated is merely an example. Results are relative to all earthquake data longitudes rather than with respect to an absolute longitude. Results Moon and Sun Considered Separately Hourly Bins A priori, it might be guessed that the moon and sun have a low influence on larger magnitude earthquakes. This proved true. The statistical significance of hourly binning of earthquakes with respect to the moon s position produced nothing extraordinary. The same was true for the sun s relative position. Tables 5 and 6 at the end of the paper show the number of earthquakes and Z-values for the moon and sun across the full 24 hours of relative longitudinal distance (the Z-value is defined as how many standard deviations the result lies away from the mean). This data was also presented in figure 4 which shows the moon s results on the left and the sun s results on the right. The results show that the moon and sun have a tiny effect on larger terrestrial earthquakes. Over the course of 38.6 years, the moon showed an increase in the likeliness of earthquakes at 4, 6 and 8 hour longitudinal offsets. Specifically, the moon-earthquake displacements with the largest associations (at the 95% confidence level): 4 hours = 0.741% greater earthquakes 6 hours = 1.62% greater earthquakes 8 hours = 0.501% greater earthquakes Dividing by the 38.6 year duration, this amounted to a combined total of 2.2 extra earthquakes per year, magnitude 5 and greater during those longitudinal offsets. Similar to the moon, the sun showed an effect on terrestrial earthquakes, although the association was even weaker. Earthquakes were more likely than expected when the sun was 17 hours behind a geographical longitude (Z=2.38, 0.746% increase, 95% confidence). On the other hand, certain longitudinal offsets for the sun and moon reduced the chance of an earthquake. In total, the moon was associated with a reduction of 1.8 earthquakes per year during lower risk offsets (9th and 10th hours). The sun was associated with 1.3 fewer earthquakes per year during the longitudinal offsets of 5, 7 and 11 hours.

7 7 Six Hour Bins The associations were non-existent for six-hour bins, defined according to the tidal cycles (figure 5). Tables 1 and 2 show the associations of the sun and moon with the relative positions of earthquakes. In retrospect, the averaging out of signal in the six hour bins should have been obvious due to the haphazard nature of the increases and decreases in the hourly bins. Aside from the minimal influence, one issue in particular stood out in this analysis: there was no clear cycle present. Reference [19] showed a daily cycle similar to the water tide, although the pattern was for shallow, low-magnitude earthquakes which have different origins than high magnitude earthquakes. The next section examines combined sun-moon positions. Surprisingly, this approach resulted in a significant changes in earthquakes for certain moon-sun pairings that were cyclical in nature. Table 1. The Z-value associated with moon-earthquake relative distances with 6-hour binning Moon (quad) EQs Z-value % extra EQs I II III IV The mean and standard deviation of the equivalent size dataset of randomly generated earthquakes was 16,681 ± 125. Table 2. The Z-value associated with sun-earthquake relative distances with 6-hour binning Sun (quad) EQs Z-value % extra EQs I II III IV The mean and standard deviation of the equivalent size dataset of randomly generated earthquakes was 16,681 ± 125. The sun s results produced low statistical fluctuations. Combined Moon-Sun Positions Six Hour Binning While looking at the moon and sun separately revealed little, looking at moon-sun combinations did result in a significant influence. Unlike the expected mean number of earthquakes for the moon and sun considered separately, the combined moon-sun ephemerides caused the mean earthquakes to vary by each bin. The effect was small (<1%) since the bin widths were relatively large, but this was still taken into account when calculating the mean expected values for each of the six-hour bins. Whenever the moon and sun paired together in the horizon, there was a significant decrease in earthquakes. Table 3 lead to the following conclusions for the (moon, sun) quadrants with a 95% confidence level:

8 8 (II,II) = 4.57% fewer earthquakes (IV,IV) = 2.31% fewer earthquakes One physical interpretation is that when the moon and sun are paired together on the horizons they act to horizontally pull the Earth s crust together, helping to seal fault lines. This force behaves differently than when the moon and sun are overhead and pull both sides of a fault upward rather than sideways. Conversely, the statistical increases in earthquakes were shallow and occurred when the moon and sun were not paired in the same quadrant: (II,IV) = 1.08% increase in earthquakes (III,I) = 1.21% increase in earthquakes (III,IV) = 0.345% increase in earthquakes It seems intuitive that when the moon and sun are pulling on opposite sides of a fault they separate the two halves just enough to trigger an earthquake. What about the other permutations? Table 3 shows the measured earthquakes versus expected earthquakes for all possible moon and sun quadrant pairings. Table 3. The Z-value associated with all possible moon and sun quadrant permutations Moon Sun EQs expected EQs Z-value % extra EQs I I I II I III I IV II I II II II III II IV III I III II III III III IV IV I 428I IV II IV III IV IV Based on randomly generated datasets of simulated earthquakes, the standard deviation was ±62.7. The number of expected earthquakes varied by bin and was therefore listed in a separate column in the table. One Hour Binning The final step was to calculate the hourly combined position results. The mean expected earthquakes varied more with hourly binning, amounting to a 4.4% effect as opposed to <1% with six-hour bins. Taking this into account, a large table of 576 bins resulted, of which most combinations produced null or negligible increases. However, 32 positions (about 5.6% of the possible permutations) were statistically

9 9 significant. An initial plot of the moon s and sun s combined positions versus the increase in expected earthquakes showed the apparently random data (left, figure 6). The right hand side of figure 6 showed the results if only the statistically significant results were plotted. Table 7 at the end of the paper lists the values for the non-negligible results (using Z > 2.5 so the chart could fit on one page). A moment should be taken to examine the true significance of the results in the table and figures. Some of the results show very sharp increases and decreases over expected, right next to an hourly bin combination that had no significant difference. For instance, (moon,sun)=(3,4) hrs offset resulted in 13% fewer earthquakes compared to the mean expected number. However, (moon,sun)=(4,4) hrs offset was within the expected range of variation. Clearly, even though the noise instrinsic to the size of the dataset was taken into account, another source of noise must still be present. As mentioned at the start of the paper, the method of looking only at the relative distance between an earthquake and the celestial positions of the moon and sun ignored localized geological properties that affect the heights of land tides. It is likely that this influence added an extra significant source of noise, on top of the noise intrinsic to the dataset. A new approach was thus proposed: treat the statistically significant points listed in Table 7 as outliers on top of a noisy dataset. Examine the outliers for a periodic trend. Periodicity in the result adds validity since the system is cyclical in nature. Then, based on the pattern identified, select the broadest possible bins in order to average out the influence of the extra noise while still preserving the earthquake signal. Figure 6. Earthquakes for various moon/sun offset positions. A plot of the moon/sun positions vs. the number of earthquakes at those positions. The redder the dot, the larger the number of earthquakes, the bluer the dot, the fewer. The plot on the right shows the same data with only the 95% confidence results. The positions which lay within expected variation were removed. Outlier Analysis Ignoring the specific value of the outliers and plotting the fewer-than-expected earthquake bins as blue circles and the higher-than-expected earthquake bins as red squares, a cyclical pattern emerged. Figure 7 shows that whenever the sun was +2.5 or -3.5 hours away from the moon, irrespective of longitude, there was a trend toward a reduction of earthquakes. The axes in figure 7 were labeled such that: x = Longitude EQ (hrs) RA moon (hrs) (4) y = Longitude EQ (hrs) RA sun (hrs) (5)

10 10 Thresholds were selected to identify the region of reduced earthquakes as three hours wide around the central, dotted line, y = x 1/2. This bin definition captured the most blue points and a minimum of red points. Labeled regions in the figure identify the bin definitions. Because the days cycle every 24 hours, some regions are extensions of others, reflected in the labeling. Region 1 accounted for 67% of the fewer-than-expected earthquake moon-sun positions, and only 11% of the higher-than-expected earthquake combinations fell in between these thresholds. Region 2 had only 13% of the lowered combinations and 34% of the higher combinations. Region 3 had 13% lowered combinations and 37% higher combinations. Region 4 had 7% lowered combinations and 17% higher combinations. Figure 7. Significant changes in earthquakes for various moon/sun offset longitudes. On the left, a plot shows the moon/sun relative longitudes. The blue circles show where there were fewer than expected earthquakes and the red squares show where there were higher than expected earthquakes. Because the day cycles every 24 hours, certain regions are extensions of each other and have been labeled to show this relationship. The diagram on the right shows a histogram of the lower than expected earthquakes (blue) and higher than expected earthquakes (red) in the labeled regions. Region 1 (when the sun was +2.5 or -3.5 hours of the moon, independent of longitude) showed reduced earthquakes. Regions 3 and 4 showed increases in high earthquake moon-sun pairings. Do the outliers predict real earthquake trends? The question is answered by summing the earthquakes data in the four regions and again comparing to the 70 randomly generated datasets. These 4 regions encompassed about 16,700 points each as opposed to the 116 points for the hourly bins which reduces random noise fluctuations by a factor of six. The mean expected earthquakes for each bin were recalculated for the four regions. Table 8 shows the results. Indeed, the abundance of outlier lower-than-expected earthquakes correctly identified a region of generally reduced earthquakes. A reduction of 3.33% earthquakes (95% confidence) implied about 14.3 fewer large magnitude earthquakes per year when the sun was +2.5 or -3.5 hours away from the moon, independent of longitude. Conversely, Region 3 (the sun was 8.5 to 14.5 hours behind the moon) showed significantly higher than expected earthquakes, 2.23% at the 95% confidence level. This amounted to 9.6 extra earthquakes per year. Starting in Region 1 and stepping the 6-hour bin through the moon-sun earthquake space an hour at a time revealed a cyclical pattern, seen in figure 8. The minimum earthquake total occured for Region 1

11 11 Table 4. The Z-value associated with the four regions depicted in Figure 7 Region EQs expected EQs Z-value % extra EQs Based on randomly generated datasets of simulated earthquakes, the standard deviation of earthquakes is +/ and the maximum aligned with Region 3. The sinusoidal line shows the estimated mean background for each of the steps, a number which varied due to the ephemerides of the combined moon-sun positions. Figure 8. Cyclical pattern from summing earthquakes by stepping through offsets to y=x-1/2. Stepping through hour-by-hour earthquakes binned +/-3 hours around the line y=x-1/2. For example, Region 1 corresponds to the 0 hours offset from y=x-1/2 (+/-3 hours) and Region 3 corresponds to the 11 hours offset (+/-3 hours). The right hand side shows the raw summed earthquakes plotted with estimated error bars. The sinusoidal line shows the estimated mean background for each of the positions, a number which varied due to the ephemerides of the combined moon-sun positions. Plotting the difference between the mean expected earthquakes and the measured earthquakes produced figure 9. The data for the plot was recorded in Table 8. Fitting an equation to the data gave: EQ extra = 15.3Cos[ 2P i 24 offset y=x 1/2] (6) This equation gives a simple relationship between the measured number of earthquakes over (or less than) the mean expected number of earthquakes which depends strictly on the relative distance between the moon and sun. It is also notable that these conclusions agree with the more narrow conclusions of 6-hour combination binning relative to a specific longitude (Table 3). Discussion The earthquake data, despite coming from a large, accurate record of events, was still noisy. The assumption in the paper was that the method of looking at relative offsets did not account for local properties

12 12 Figure 9. Cyclical pattern from summing earthquakes by stepping through offsets to y=x-1/2. Stepping through hour-by-hour earthquakes binned +/-3 hours around the line y=x-1/2. For example, the 0 hours offset corresponds to Region 1 and the 12 hour offset corresponds to Region 3. that influenced land tide heights and, in turn, influenced the likeliness of an earthquake. The only way that local geography could possibly influence the results would be if specific faults or fault types were over-represented in the outlier moon-sun combinations. Then the local geography would exert an influence. If this was not the case, but the fault lines and types were randomly scattered in the moon-sun earthquake space then local geography could not be a noise source and something else must be influencing the results. A future step to research would be to examine the specific fault and fault type distributions in the data in more detail to verify this statement. Despite this, by examining outliers, a clear association of the positions of the moon and sun and earthquakes was identified. Surprisingly, the clearest trend was for a reduction of earthquakes when the moon and sun were close together in the sky, independent of the longitude of an earthquake. When the moon and sun were on opposite sides of the Earth (Region 3) there was a marked increase in large magnitude earthquakes. The fact that the location of the earthquake was irrelevant may be related to the depth of the earthquakes. Deeper fault lines, closer to the Earth s core, may be less influenced by the circumferential effects of the moon and sun. It is interesting to note that although the association was based on a statistical analysis of a large dataset, the pattern holds for three of the four highest-magnitude earthquakes in the study, including: December 2004 Northern Sumatra earthquake (magnitude 9.1, relative sun RA=0.227, relative moon RA=12.5) March 2005 Indonesia earthquake (magnitude 8.6, relative sun RA=9.17, relative moon RA=18.8) and February 2010 Chile earthquake (magnitude 8.8, relative sun RA=4.09, relative moon RA=15.0) As always, correlation does not necessarily imply causation. One alternate source for the results of the study could be an estimate for a correction to the moon/sun ephemerides. The appearance of an increase or decrease in earthquakes could imply an error in the calculated moon/sun positions at the time of the earthquakes. This paper used Mathematica to generate the RA coordinates which which employed the current standard, VSOP87 (Variations Sculaires des Orbites Plantaires). The precision of this standard is stated as 1 per 4000 years around the year 2000 (see reference [22]), making the possibility of the bias unlikely.

13 13 Acknowledgments The author expresses thanks to the USGS for posting their earthquake data online and Mathematica for making this project possible in a finite amount of time. References 1. Reiter M (2004) Subsurface temperatures and crustal strength changes within the seismogenic layer at arroyo del coyote in the socorro seismic area, central rio grande rift, new mexico. The Geological Society of America 117: Grollimund B, Zoback M (2001) Did deglaciation trigger intraplate seismicity in the new madrid seismic zone? Geology 29: Tsuneji Rikitakea KH (2004) Encyclopedia of Physical Science and Technology pp. http: // 4. Tsuneji Rikitakea KH (2007) Treatise on Geophysics, volume pp. sciencedirect.com/science/article/pii/b S Odintsov KGBKDA K Boyarchuk (2006) Long-period trends in global seismic and geomagnetic activity and their relation to solar activity. Physics and Chemistry of the Earth 31: Taka aki Taira1 FNRMN Paul G Silver1 (2009) Remote triggering of fault-strength changes on the san andreas fault at parkfield. Nature : (1909) Science measures the land tides. New York Times. 8. Roach J (2005) Can the moon cause earthquakes? National Geographic. 9. NASA. Moon fact sheet Knopoff L (1970) Correlation of earthquakes with lunar orbital motions. Earth, Moon, and Planets 2: John E Vidale MJSJDHO Duncan Carr Agnew (1998) Absence of earthquake correlation with earth tides: An indication of high preseismic fault stress rate. Journal of Geophysical Research 103: 24,567-24, B Bodri SI (1988) On the correlation between earth tides and microseismic activity. Physics of the Earth and Planetary Interiors 55: Allen MW (1936) The lunar triggering effect on earthquakes in southern california. Bulletin of the Seismological Society of America April: Simpson JF (1967) Earth tides as a triggering mechanism for earthquakes. Earth and Planetary Science Letters 2: D A Lockner NMB (1999) Premonitory slip and tidal triggering of earthquakes. Journal of Geophysics 104: Kasahara J (2002) Tides, earthquakes and volcanoes. Science 297: N M Beeler DAL (2003) Why earthquakes correlate weakly with the solid earth tides: effects of periodic stress on the rate and probability of an earthquake occurrence. Journal of Geophysics Research 108: 2391.

14 E S Cochran ST J E Vidale (2004) Earth tides can trigger shallow thrust fault earthquakes. Science 306: Laurent Metivier CPCSRMDGP Olivier de Viron (2009) Evidence of earthquake triggering by the solid earth tides. Earth and Planetary Science Letters 278: NASA. Ocean motion and surface currents. tides-observing.htm. 21. USGS. United states geological survey global earthquake search form. gov/earthquakes/eqarchives/epic/epic_global.php. 22. P Bretagnon GF (1988) Planetary theories in rectangular and spherical variables. vsop87 solutions. Astronomy and Astrophysics : (1951) Leveling effect of land-tide on earth and moon. Popular Astronomy 59: 278. Tables

15 15 Table 5. The Z-value associated with moon-earthquake relative distances Moon (hrs) EQs Z-value % extra EQs A positive Z-value reflects more earthquakes than expected while a negative Z-value implies less earthquakes than expected. A simulation was used to calculate the mean number of earthquakes and the mean standard deviation of the equivalent size dataset of randomly generated earthquakes. This produced the values 2780 ± As can be seen in the table, the results had mostly no or low statistical significance.

16 16 Table 6. The Z-value associated with sun-earthquake relative distances Sun (hrs) EQs Z-value % extra EQs A simulation showed the mean number of earthquakes during the time-period to be 2780 ± The sun produced even lower statistically significant results than the moon.

17 17 Table 7. The non-zero values for earthquakes exceeding (or less than) the expected variation in earthquakes at the 95% confidence level. The moon-sun-earthquake relative distances. Moon (hrs) Sun (hrs) EQs Z-value % extra EQs The mean of the equivalent size dataset of randomly generated earthquakes varied by position according to the moon and sun ephemerides (calculated via Mathematica). The means ranged from 111 to 120 earthquakes during the dataset time span (38.6 years). The standard deviation was ± The large number of statistically significant results is attributed to noise introduced by the method of analysis (relying on the position of the sun and moon to gauge the effect of land-tides on earthquake fault lines).

18 18 Table 8. The Z-value associated with the y=x-1/2 region stepped through all 24 hours, as depicted in Figure 8 Region EQs expected EQs Z-value % extra EQs extra EQ/yr Based on randomly generated datasets of simulated earthquakes, the standard deviation of earthquakes is +/- 125.