Approximate determination of an Origin—Earthquake-hunting in Japan—Determinations by direction of motion—Direction indicated by destruction of buildings—Direction determined by rotation—Cause of rotation—The use of time observations—Errors in such observations—Origin determined by the method of straight lines—The method of circles, the method of hyperbolas, the method of co-ordinates—Haughton’s method—Difference in time between sound, earth, and water waves—Method of Seebach. One of the most practical problems which can be suggested to the seismologist is the determination of the district or districts in any given country from which earthquake disturbances originate. With a map of a country before us, shaded with tints of different intensity to indicate the relative frequency of seismic disturbance in various districts, we at once see the localities where we might dwell with the least disturbance, and those we should seek if we wish to make observational seismology a study. Before erecting observatories for the systematic investigation of earthquakes in a country, it would be necessary for us, in some way or other, to examine the proposed country to find out the most suitable district. The special problem of determining approximately the origin or origins of a set of earthquakes would be given to us. Having made this preliminary investigation, the next point is, by means of observatories so arranged that Approximate determination of origins.—In 1880 I obtained a tolerably fair idea of the distribution of seismic energy throughout Japan, by compiling the facts obtained from some hundreds of communications received from various parts of the country respecting the number of earthquakes that had been felt. The communications were replies to letters sent to various residents in the country and to a large number of public officers. By taking these records, in conjunction with the records made by instruments, it was ascertained that in Japan alone there were certainly 1,200 shocks felt during the year, that is to say, three or four shocks per day. The greater number of these shocks were felt along the eastern coast, commencing at Tokio, in the south, and going northwards to the end of the main island. These shocks were seldom felt on the west coast. It appeared as if the central range of mountains formed a barrier to their progress. Similarly, ranges of mountains to the south-west of Tokio prevented the shocks from travelling southwards. Proceeding in this way the conclusion was arrived at that the west coast, the southern part of Japan, and the islands of Shikoku and Kiushiu, had their own local earthquakes. Earthquake-hunting.—These preliminary enquiries having shown that the northern part of Japan was a better district for seismological observations than the southern half, the next step was to subject the northern half to a closer analysis. This analysis was commenced by sending The actual results obtained in three months by this method of working are shown in the accompanying map, which embraces the northern half of the main island of Nipon and part of Yezo. The shaded portion of the map indicates the mountainous districts, which are traversed by ranges varying in height from about 2,000 and 7,000 feet. The dotted lines show the boundaries of the more important groups of earthquakes which were recorded. I. is the western boundary of earthquakes, which at places to the eastward are usually felt somewhat severely. Some of these have been felt the most severely at or near Hakodadi, whilst farther south their effects have been weak. Occasionally the greatest effect has been near to Kameishi. Sometimes these earthquakes terminate along the western boundaries of III. or IV., not being able to pass the high range of mountains which separate the plain of Musashi from Kofu. Fig. 29.—Northern Japan. Mountainous districts shaded with oblique lines. II. is the boundary of a shock confined to the plain which surrounds Kofu. These earthquakes are evidently quite local. Many of the disturbances have evidently originated beneath the ocean, having come in upon the land in the direction of the arrows a or b. III. This line indicates the boundary of a group of shocks which are often experienced in Tokio. These may come in the directions d, e, or f. It is probable that some of them originate to the eastward of Yokohama, on or near to the opposite peninsula. IV. V. and VI. The earthquakes bounded by these lines probably originate in the directions c or d. VII. The earthquakes bounded by this line probably come from the direction e. VIII. This line gives us the boundary of earthquakes which may come from the direction b. The above boundaries sometimes do not extend so far to the westward as they are shown. At other times, groups like V. and VI. extend farther to the south-west. These earthquake boundaries, which so clearly show the effects of high mountains in preventing the extension of motion, have been drawn up, not from single earthquakes, but from a large series of earthquakes which have been plotted upon blank maps, and are now bound together to form an atlas. To give an idea of the material upon which I have been working, I may state that between March 1 and March 10, 1882, I received records of no less than thirty-four distinct shocks felt in districts between Hakodate and Tokio, and for each of these it is quite possible to draw a map. In addition to the boundaries of disturbances given in the accompanying map, other boundaries might be drawn for shocks which were more local in their character. The groups which contained the greatest number of shocks are III., IV., V., VI., and VII. By Determination of earthquake origins from the direction of motion.—If we assume that an earthquake is propagated from a centre as a series of waves, in which normal vibrations are conspicuous, and obtain at two localities, not in the same straight line with the origin, Direction determined from destruction of buildings.—When an observer first sees a town that has been partially shattered by an earthquake, all appears to be confusion, and it is difficult to imagine that in such apparent chaos we are able to discover laws. If, however, we take a general view of this destruction and compare together similarly built buildings, it is possible to discover that similar and similarly situated structures have suffered in a similar manner. By carefully analysing the destruction we are enabled to infer the direction in which the destroying forces have acted. It was chiefly by observing the cracks in buildings, and the direction in which bodies were overthrown or projected, that Mallet determined the origin of the Neapolitan earthquake. From the observations given in Chapter VII. it would appear that, with destructive earthquakes, walls which are transverse to the direction of motion are most likely to be overturned, whilst, with small earthquakes, these walls are the least liable to be fractured. From a critical examination of the general nature of the damage done on the buildings of a town, earthquake observers have shown that the direction of a shock may often be approximately determined. The direction in The rotation of bodies.—It has often been observed that almost all large earthquakes have caused objects like tombstones, obelisks, chimneys, &c., to rotate. One of the most natural and at the same time most simple explanations is to suppose that during the shock there had been a twisting, or backward and forward screw-like motion in the ground. Amongst the Italians and the Mexicans earthquakes producing an effect like this are spoken of as ‘vorticosi.’ In the Calabrian earthquake, not only were bodies like obelisks twisted on their bases, but straight rows of trees seem to have been left in interrupted zigzags. These latter phenomena have been explained upon the assumption of the interference of direct waves and reflected waves, the consequence of which being that points in close proximity might be caused to move in opposite directions. Reflections such as these would be most likely to occur near to the junction of strata of different elasticity, and it may be remarked that it is often near such places that much twisting has been observed. Another way in which it is possible for twisting to have taken place would be by the interference of the normal and transverse waves which probably always exist in an earthquake shock, or by the meeting of the parts of the normal wave itself, one having travelled in a direct line from the origin, whilst the other, travelling through variable material, has had its direction changed. Mallet, however, has shown that the rotation may have been in many cases brought about without the supposition of any actual twisting motion of the earth— Although the above explanation is simple, and also in many cases probably true, it hardly appears sufficient to account for all the phenomena which have been observed. Thus, for instance, if the stones in the Yokohama cemetery, at the time of the earthquake of 1880, had been twisted in consequence of the cause suggested by Mallet, we should most certainly have found that some stones had turned in one direction whilst others had been twisted in another. By a careful examination of the rotated stones, I found that every stone—the stones being in parallel lines—had revolved in the same direction, namely in a direction opposite to that of the hands of a watch. As it would seem highly improbable that the centre of greatest friction in all these stones of different sizes and shapes should have been at the same side of their centres of gravity, an effect like this could only be explained by the conjoint action of two successive shocks, the direction of one being transverse to the other. Although fully recognising the sufficiency of two transverse shocks to produce the effects which have been observed in Yokohama, I will offer what appears to me to be the true explanation of this phenomenon: it was first suggested by my colleague, Mr. Gray, and appears to be simpler than any with which I am acquainted. Fig. 30. If any columnar-like object, for example a prism which the basal section is represented by a b c d (see fig. 30), receives a shock at right angles to b c, there will be a tendency for the inertia of the body to cause it to overturn on the edge b c. If the shock were at right angles to d c, the tendency would be to overturn on the edge d c. If the shock were in the direction of the diagonal c a, the tendency would be to overturn on the point c. Let us, however, now suppose the impulse to be in some direction like e g, where g is the centre of gravity of the body. For simplicity we may imagine the overturning effect to be an impulse given through g in an opposite direction—that is, in the direction g e. This force will tend to tip or make the body bear heavily on c, and at the same time to whirl round c as an axis, the direction of turn being in the direction of the hands of a watch. If, however, the direction of impulse had been e' g, then, although the turning would still have been round c, the direction would have been opposite to that of the hands of a watch. To put these statements in another form, imagine g e' to be resolved into two components, one of them along g c and the other at right angles, g f. Here the component of the direction g c tends to make the body tip on c, whilst the other component along g f causes revolution. Similarly g e may be resolved into its two components g c and g f', the latter being the one tending to cause revolution. From this we see that if a body has a rectangular section, so long as it is acted upon by a shock which is parallel to its sides or to its diagonals, there ought not to Further, it will be observed that the tendency of the turning is to bring a stone, like the one we are discussing, broadside on to the shock; therefore, if a stone with a rectangular cross section has turned sufficiently, the direction of a shock will be parallel to one of its faces, but if it has not turned sufficiently it will be more nearly parallel to its faces in their new position than it was to its faces when in their original position. If a stone receives a shock nearly parallel with its diagonal, on account of its instability it may turn either positively or negatively according as the friction on its base or some irregularity of surface bearing most influence. Similarly, if a stone receives a shock parallel to one of its faces, the twisting may be either positive or negative, but the probability is that it would only turn slightly; whereas in the former case, where the shock was nearly parallel to a diagonal, the turning would probably be great. Determination of direction from instruments.—When speaking about earthquakes it was shown, as the result of many observations, that the same earthquake in the space of a few seconds, although it may sometimes have only one direction of motion, very often has many One remarkable fact which has been pointed out in reference to artificial earthquakes produced by exploding charges of gunpowder or dynamite, and also with regard to certain earthquakes, is that the greatest motion of the ground is inwards, towards the point from which the disturbance originated. Should this prove the rule, it gives a means of determining, not only the direction of earthquake, but the side from which it came. Determination of earthquake origins by time observations.—The times at which an earthquake was felt at a number of stations are among the most important observations which can be made for the determination of an earthquake origin. The methods of making time observations, and the difficulties which have to be overcome, have already been described. When determining the direction from which a shock has originated, or determining the origin of the shock by means of time observations, it has been usual to assume that the velocity of propagation of the shock has been uniform from the origin. The errors involved in this assumption appear to as follows:— 1. We know from observations on artificial earthquakes that the velocity of propagation is greater between stations near to the origin of the shock than it is between more remote stations; and also the velocity of propagation varies with the initial force which produced the disturbance. If our points of observation are sufficiently close together as compared with their distance from the 2. We have reasons for believing that the transit velocity of an earthquake is dependent on the nature of the rocks through which it is propagated. Errors which arise from causes of this description will, however, be practically eliminated if our observation points are situated on an area sufficiently large, so that the distribution of the causes tending to alter the velocity of a shock balance each other. It must be remarked, that causes of this description may also produce an alteration in the direction of our shock. Other errors which may sometimes enter into our results, when determining the origin of shocks by means of observations on velocities, are the assumptions that the disturbance has travelled along the surface from the epicentrum and not in a direct line from the centrum. Again, it is assumed that the origin is a point, whereas it may possibly be a cavity or a fissure. Lastly, if we desire extreme accuracy, we must make due allowance for the sphericity of the earth and the differences of elevation of the observing stations. I. The method of straight lines.—Given a number of pairs of points a0, a1, b0, b1, c0, c1, &c., at each of which the shock was felt simultaneously, to determine the origin. Theoretically if we bisect the line which joins a0 and a1 by a line at right angles to a0, a1, and similarly bisect the lines b0, b1, c0, c1, all these bisecting lines a0, a1, b0, b1, c0, c1, &c., ought to intersect in a point, which point will be the epicentrum or the point above the origin. This method will fail, first, if a0, a1, b0, b0, c0, c1 form a continuous straight line, or if they form a series of parallel lines. Hopkins gives a method based on a principle similar to the one which is here employed—namely, given that a shock arrives simultaneously at three points to determine, the centre. In this case, the relative positions of the three points, where the time of arrival was simultaneous, must be accurately known, and these three points must not lie in a straight line, or the method will fail. For practical application the problem must be restricted to the case of three points which do not lie nearly in the same straight line. II. The method of circles.—Given the times t0, t1, t2, &c., at which a shock arrived at a number of places a0, a1, a2, &c., to determine the position from which the shock originated. Suppose a0 to be the place which the shock reached first, and that it reached a1, a2, a3, &c., successively afterwards. Let t1 - t0 = a With a1, a2, a3, &c. as centres, describe circles with radii proportional to the known qualities a, b, c, &c., and also a circle which passes through a0 and touches these circles. The centre of the last circle will be the epicentrum. The radii proportional to a, b, c, &c. may be represented by the quantities ax, bx, cx, &c., where x is the velocity of propagation of the shock. It will be observed that the times at which the shock arrived at three places might alone be sufficient. If, instead of taking the times of arrival of a shock, the arrival of a sea wave be taken, the result would be a closer approximate to the absolute truth. It will be observed that this method is not a direct one, but is one of trial. If, however, an imaginary case In the following table the first column gives the times at which the sea wave arrived at each of these places in Iquique time; in the second column the difference between these times and the time at which it reached Huanillos is given; in the third column the distances through which a sea wave, propagated at the rate of 350 feet per second, could travel during the intervals noted in the second column is given.
The distances marked in the third column are used as radii of the circles drawn round the places to which they respectively refer. The centre of the circle drawn to touch the circles of the first column, and at the same time to pass through Huanillos, is marked c. The position from which the shock originated appears therefore to have occurred very near to a place lying in Long. 7° 15' W. and Lat. 21° 22' S. Fig. 31. The actual operations which were gone through in making the accompanying map were as follows. First, the places with which we had to deal were represented on a map in orthographical projection, the centre of projection III. The method of hyperbolas.—The method which I call that of hyperbolas is only another form of the method of circles. It is, however, useful in special cases, as, for instance, where we have the times of arrival of earthquakes at only two stations. Between Tokio and Yokohama, at which places I frequently obtain tolerably accurate time records, the method has been applied on several occasions with advantage. In the preceding example let us suppose that the only time records which we had were for Huanillos and Mejillones, and that the wave was felt at the IV. The method of co-ordinates.—Given the times at which a shock arrived at five or more places, the position of which we have marked upon a map, or chart, to determine the position on the map of the centre of the shock, its depth, and the velocity of propagation. Commencing with the place which was last reached by the shock, call these places p, p1, p2, p3, and p4, and let the times taken to reach these places from the origin be respectively t, t1, t2, t3, and t4. Through p draw rectangular co-ordinates, and with a scale measure the co-ordinates of p1, p2, p3, and p4, and let these respectively be a1, b1; a2, b2; a3, b3; a4, b4. Then if x, y, and z be the co-ordinates of the origin of the shock, d, d1, d2, d3, and d4, the respective distances of p, p1, p2, p3, and p4 from this origin, and v the velocity of the shock, we have
Because we know the actual times at which the waves arrived at the places p, p1, p2, p3, p4, we know the values t—t1, t—t2, t—t3, t—t4. Call these respectively m, p, q, and r. Suppose t known, then
Subtracting equation No. 1 from each of the equations 2, 3, 4, and 5, we obtain,
Now let v2 = u, and 2v2 t = w. Then
We have here four simple equations containing the four unknown quantities x, y, u, and w. x and y determine the origin of the shock. The absolute velocity v equals v u. From v and w we obtain t. Substituting x, y, v, and t in the first equation we obtain z. We have here assumed that the points of observation have all about the same elevation above sea level. If it is thought necessary to take these elevations into account, a sixth equation may be introduced. If the propagation of the wave is considered as a horizontal one, as would be done when calculating the position of the epicentrum or point above the origin, by means of the times of arrival of a sea wave, the ordinate z of the first five equations would be omitted. Working in this way the resulting four equations, viz. 2a1 x + 2b1 y + um2 - wm2 = a12 + b12 &c. &c. &c. remained unchanged. Applying this method to the same example as that
From this data we find the co-ordinates x and y of this origin to be 85·8 and 56·7; whilst the velocity of propagation = 45 feet per second. Measuring these ordinates upon the map, we obtain a centre lying very near Long. 71° 5' W. and Lat. 21° 22' S., a position which is very near to that which has already been obtained by other methods. If instead of Huanillos we substitute the ordinates and time of arrival of the sea wave for Pabalon de Pica, another point for the origin will be obtained lying farther out at sea. To obtain the best result, the method to be taken will evidently be, first to reject those places at which it seems likely that some mistake has been made with the time observations, and then with the remaining places to form as many equations as possible, and from these to obtain a mean value. This is a long and tedious process, and as the time observations of this particular earthquake are probably one and all more or less inaccurate, it is hardly worth while to follow the investigation farther. In this example, as in the preceding ones, it will be In these determinations it will also have been observed that it is assumed that the disturbance has radiated from a centre, and, therefore, approached the various stations in different directions. If we assume that we have three stations very near to each other as compared with their distances from the origin, so that we can assume that the wave fronts at the various stations were parallel, the determination of the direction in which the wave advanced appears to be much simplified. The determination of the direction in which a wave has passed across three stations was first given by Professor Haughton. Haughton’s method.—Given, the time of an earthquake shock at three places, to determine its horizontal velocity and coseismal line. The solution of this is contained in the formula When a, b, and c are three stations at which a shock is observed at the times t1, t2, and t3; a, b, and c are the distances between a, b, and c, and ? is the angle made by the coseismal lines x a x, y b y, and the line a b, which are assumed to be parallel. This I applied in the case of the Iquique earthquake, but owing to the smallness of the angles between the three stations a, b, and c, the result was unsatisfactory. The problem ought to be restricted, first, to places which are a long distance away from a centre, and, secondly, to places which are not nearly in a straight line. This problem may be solved more readily by geometrical methods. The difference in time of the arrival of two disturbances.—In the various calculations which have been made to determine an origin based on the assumption of a known or of a constant velocity, we have only dealt with a single wave, which may have been a disturbance in the earth or in the water. A factor which has not yet been employed in this investigation is the difference in time between the arrival of two disturbances; one propagated, for instance, through the earth, and the other, for example, through the ocean. The difference in the times of the arrival of two waves of this description is a quantity which is so often recorded that it is well not to pass it by unnoticed. To the waves mentioned we might also add sound waves, which so frequently accompany destructive earthquakes, and, in some localities, as, for instance, in Kameishi, in North Japan, are also commonly associated with earthquakes of but small intensity. It was by observing the difference in time between the shaking and the sound in different localities that Signor Abella was enabled to come to definite conclusions regarding the origin of the disturbances which affected the province of Neuva Viscoya in the Philippines, in 1881; the places where the interval of time was short, or the places where the two phenomena were almost simultaneous, being, in all probability, nearer to the origin than when the intervals were comparatively large. I myself applied the method with considerable success A method similar to the above was first suggested by Hopkins. It depended on the differences of velocity with which normal and transversal waves are propagated.[86] Seebach’s method.—To determine the true velocity of an earthquake, the time of the first shock, and the depth of the centre. Fig. 32. Let the straight line m, m1, m2, m3 represent the surface of the earth shaken by an earthquake. For small earthquakes, to consider the surface of the earth as a plane will not lead to serious errors. If an earthquake originates at c, then to reach the surface at m it traverses a distance h in the time t. To reach the surface at m1 it traverses a distance h + x1 in a time t2. If v equals the velocity of propagation, then t = h/v, t1 = h + x1/v, t2 = h + x2/v, &c. Seebach now says that if we have given the position of m or epicentrum of the shock, and draw through it rectangular axes like m m3 and m t3, and lay down The intersection of the asymptote with the ordinate axis is the time point of the first shock, which, because the scale for time and for space were taken as equal, gives the absolute position of the centrum. This intersection is shown by dotted lines. Knowing the position of the centrum, we can directly read from our diagram how far the disturbance has been propagated in a given time. |