Correspondence: Response of a gravimeter to an instantaneous step in gravity.

Introduction

Montagner et al.[1] presented evidence that a signal observed prior to P-waves on gravimeters operating in Japan during the 2011 Tohoku-Oki earthquake was caused by gravitational field changes propagating at light speed. In this comment, I explore the expected response of gravimeters attached to the Earth’s surface and subjected to an instantaneous acceleration step, , where is a Heaviside step function. I argue that at the time of the step, inertial accelerations from the elastic response of the Earth exactly cancel the gravitational step; in order to observe the gravity step using a gravimeter, one must wait until the Earth deforms and begins to re-equilibrate. This conclusion is an application of the principle of equivalence[2] between gravitational and inertial mass. In the case of the signal reported[1] for Tohoku-Oki, I estimate that inclusion of the Earth’s elastic response should significantly decrease the expected acceleration at the time of the P-wave.

To determine a complete solution for the response of a gravimeter attached to the Earth’s surface, I would need to (1) calculate the spatial distribution of gravitational changes caused by mass redistribution in an earthquake, (2) calculate the motion at the gravimeter due to unbalanced elastic and gravitational forces, and (3) sum the gravity change with the inertial accelerations determined in step 2. Finding the Earth’s elastic response is a complex problem that can be approached by summing the responses of the Earth’s normal modes. When the Earth response problem is decomposed into a sum of modes, the time response of each mode simplifies into a set of equivalent linear single-degree-of-freedom (sdof) oscillator problems driven by a step function forcing function of the appropriate amplitude. Alternatively, I could use a finite-element code to calculate both the change in gravity and also the motion of the Earth caused by those changes in gravitational body forces. Performing those calculations would transform this brief communication into a long treatise. However, I can provide insight into this problem by introducing a whimsical gedanken experiment.

Consider the response of a gravimeter attached to the Earth’s surface, and that the gravitational constant is suddenly transformed to zero, or . You might guess that when this happened, we would all instantly become weightless. However, what would actually happen is that the Earth (previously compressed by gravity) would instantaneously begin to expand with an outward acceleration equal to gravity. That is, initially our weight would be virtually unchanged until the Earth’s outward acceleration decreased as it started to approach its new uncompressed weightless state; only then would we truly feel weightless.

This gedanken experiment illuminates the basic principle of my comment about using gravimeters to detect changes in gravity. Gravimeters are vertical accelerometers that measure the sum of the vertical acceleration of the Earth’s surface plus the vertical gravitational acceleration. When the Earth is in gravitational/elastic equilibrium, then there is no inertial acceleration of the Earth and then the gravimeter measures only the gravitational acceleration.

A more quantitative description of the dynamics of an elastic body that responds to gravitational changes comes from considering the case of a linear sdof oscillator that is attached to some stationary point; for a deformable body, it is the center of mass. The solution of this sdof provides the time history of motions of the Earth described by (Earth’s fundamental radial mode). It also provides a solution of the gedanken problem (at least for the first seconds when the response is linear). Figure 1 shows parameters used for this calculation. I assume that the downward gravitational force on is , and that prior to the gravity change, the Earth is in static equilibrium, and that the Earth’s elastic response is linear for small perturbations from equilibrium.

Fig. 1

Schematic showing the parameters used to describe the acceleration of an elastically supported mass in the presence of a downward gravitational acceleration with an instantaneous step of . is upward motion with respect to elastic/gravitational equilibrium

The equation of the generalized motion of the Earth’s surface iswhere , and . is the characteristic frequency at which the Earth oscillates about its new equilibrium (about 8.1 × 10−4 Hz), and , where for . The motions for small times can be determined by recognizing that the third term in (1) is negligible at short times ( is initially zero). That is, for short times,

Since the system is initially at rest (i.e., equilibrium), , and then . The full solution of (1) with at-rest initial conditions iswhere . Taking time derivatives of (3),

And

If the damping is small , then andand

The total vertical acceleration of the mass (inertial acceleration plus gravity) is then

Notice that as , the inertial response of the Earth, , exactly cancels the gravitational change . This is because inertial accelerations from unbalanced elastic forces are indistinguishable from unbalanced gravitational accelerations. Notice that this cancelling effect is lessened if the damping is increased; there is a second term in (5) that increases as . This term is from viscous forces that are not part of the static equilibrium problem, and this term grows linearly with time. If the Earth was heavily damped, then the gravitational forces would be easier to observe with an accelerometer. However, since of the Earth is large (>500), this damping term is very small.

Assuming that the elastic response of the Earth is approximately given by , I used (6) to estimate that the acceleration signal after 12 s would be about . Even after 120 s, it only grows to .

Although this simple sdof gedanken problem is a gross oversimplification of the Earth’s response, it does allow us to guess that the acceleration response at a point will include the inertial effects of changing the gravimeter’s initial equilibrium position to a new equilibrium position appropriate for the gravity perturbation. The gravimeter will then oscillate about that new position until the motions damp out (from inelastic and radiation damping). Nevertheless, the initial acceleration of the gravimeter will be equal and opposite to the gravity change. Considering the Tohoku earthquake, there are only two plausible length/time scales to determine this oscillation period; the first is the Earth’s dimension (as in ), and the second is the source–observer distance. For the Tohoku earthquake, the travel time of elastic deformations between the source and observer seems to be the natural scale. I speculate that the effective period will be on the order of 500 s, which is four times the S-wave travel time. for this deformation may be about 500. That is, I hypothesize that the net acceleration on a gravimeter is approximatelyWhere , is the hypocentral distance and is the shear wave velocity. If the gravity signal is changing with time as , then the expected gravimeter signal isWhere is the convolution operator. For short times after the earthquake origin, a Taylor expansion of the cosine can be used to approximate (10) as

Therefore, the expected gravimeter signal for the Tohoku earthquake that is shown in Fig. 4a of Montagner et al.[1], should be modified as

This means that the gravimeter response is very sluggish with time and it is poorly suited to rapidly detect earthquakes as is required for seismic early warning. Examination of Fig. 4b from Montagner et al.[1] shows that their predicted gravitational signal significantly overpredicts the accelerations recorded by the gravimeter. Perhaps modification of their predicted acceleration with the convolution of (11) will produce better agreement with the observation.

Data sharing is not applicable to this article as no data sets were generated or analyzed during the current study.