GPS for real-time earthquake source determination and tsunami warning systems

We identify the key design aspects of a GPS-based system (and in the future, GNSS-based systems) that could contribute to real-time earthquake source determination and tsunami warning systems. Our approach is based on models of both transient and permanent displacement of GPS stations caused by large earthquakes, while considering the effect of GPS errors on inverted earthquake source parameters. Our main conclusions are that (1) the spatial pattern, magnitude, and timing of permanent displacement of GPS stations can be inverted for the earthquake source and so predict the 3D displacement field of the ocean bottom, thus providing the initial conditions for tsunami models, and (2) there are no inherently limiting factors arising from real-time orbit and positioning errors, provided sufficient near-field GPS stations are deployed. This signal could be readily exploited by GPS networks currently in place, and will be facilitated by the IGS Real-Time Project as it comes to fruition.     

User expectations for multibeam echo sounders backscatter strength data-looking back into the future

With the ability of multibeam echo sounders (MBES) to measure backscatter strength (BS) as a function of true angle of insonification across the seafloor, came a new recognition of the potential of backscatter measurements to remotely characterize the properties of the seafloor. Advances in transducer design, digital electronics, signal processing capabilities, navigation, and graphic display devices, have improved the resolution and particularly the dynamic range available to sonar and processing software manufacturers. Alongside these improvements the expectations of what the data can deliver has also grown. In this paper, we identify these user-expectations and explore how MBES backscatter is utilized by different communities involved in marine seabed research at present, and the aspirations that these communities have for the data in the future. The results presented here are based on a user survey conducted by the GeoHab (Marine Geological and Biological Habitat Mapping) association. This paper summarises the different processing procedures employed to extract useful information from MBES backscatter data and the various intentions for which the user community collect the data. We show how a range of backscatter output products are generated from the different processing procedures, and how these results are taken up by different scientific disciplines, and also identify common constraints in handling MBES BS data. Finally, we outline our expectations for the future of this unique and important data source for seafloor mapping and characterisation.     

Seismic parameters controlling far-field tsunami amplitudes: A review

We present a review of the influence of various parameters of the sources of major oceanic earthquakes on the amplitude of tsunamis at transoceanic distances. We base our computations on the normal mode formalism, applied to realistic Earth models, but interpret our principal results in the simpler framework of Haskell theory in the case of a water layer over a Poisson half-space. Our results show that source depth and focal geometry play only a limited role in controlling the amplitude of the tsunami; their combined influence reaches at most 1 order of magnitude down to a depth of 150 km into the hard rock. More important are the effects of directivity due to rupture propagation along the fault, which for large earthquakes can result in a ten-fold decrease in tsunami amplitude by destructive interference, and the possibility of enhanced tsunami excitation in material with weaker elastic properties, such as sedimentary layers. Modelling of the so-called tsunami earthquakes suggests that an event for which 10% of the moment release takes place in sediments generates a tsunami 10 times larger than its seismic moment would suggest. We also investigate the properties of non-double couple sources and find that their relative excitation of tsunamis and Rayleigh waves is in general comparable to that of regular seismic sources. In particular, landslides involving weak sediments could result in very large tsunamis. Finally, we emphasize that the final amplitude at a receiving shore can be strongly affected by focusing and defocusing effects, due to variations in bathymetry along the path of the tsunami.     

Far-field simulation of the 1946 Aleutian tsunami

We present hydrodynamic far-field simulations of the Aleutian tsunami of 1946 April 1, using both a dislocation source representing a slow earthquake and a dipolar one modelling a large landslide. The earthquake source is derived from the recent seismological study by López and Okal, while the landslide source was previously used to explain the exceptional run-up at Scotch Cap in the near field. The simulations are compared to a field data set previously compiled from testimonies of elderly witnesses at 27 far-field locations principally in the Austral and Marquesas Islands, with additional sites at Pitcairn, Easter and Juan Fernández. We find that the data set is modelled satisfactorily by the dislocation source, while the landslide fails to match the measured amplitudes, and to give a proper rendition of the physical interaction of the wavefield with the shore, in particular at Nuku Hiva, Marquesas. The emerging picture is that the event involved both a very slow earthquake, responsible for the far-field tsunami, and a major landslide explaining the near-field run-up, but with a negligible contribution in the far field.     

Global telescope networks for astronomy education

Ever wanted to give school students or the public a closer look at the stars, without worrying about the vagaries of the UK weather? Paul Roche and Rachel Dodds from the Faulkes Telescope Project explain how you can do just that – without even leaving the classroom!     

One-station estimates of seismic moments from the mantle magnitudeM m: The case of the regional field (1.5°≤Δ≤15°)

We extend to the regional field of distances the procedure of one-station estimation of seismic moments using the mantle magnitudeM _m, as introduced earlier in the case of teleseismic events. A theoretical analysis of the validity of the asymptotic expansion of normal modes in terms of surface waves, which was used in the development ofM _m, upholds the validity of the algorithm for distances as short as 1.5°. This is confirmed by the analysis of a dataset of 149 GEOSCOPE records obtained at distances ranging from 1.5 to 15°, from earthquakes with moments between 10²⁴ and 2.5×10²⁷ dyn-cm. The performance ofM _m as measured in terms of average residual with respect to published values ofM ₀, and standard deviation of the residuals, is not degraded in this distance range, with respect to the teleseismic case. This indicates that the mantle magnitudeM _mcan be reliably used at regional distances, notably for tsunami warning applications.     

Use of the mantle magnitudeM m for the reassessment of the moment of historical earthquakes

The mantle magnitudeM _m is used on a dataset of more than 180 wavetrains from 44 large shallow historical earthquakes to reassess their moments, which in many cases had been previously estimated only on the basis of the earthquake's rupture area. We provide 27 new or revised values ofM _o, based on the spectral amplitudes of surface waves recorded at a number of stations, principally Uppsala and Pasadena. Among them, and most significantly, we document a large low-frequency component to the source of the 1923 Kanto earthquake: the low-frequency seismic moment is 2.9×10²⁸ dyn-cm, in accord with geodetic observations. On the other hand, we revise downwards the seismic moment of the 1906 Ecuador event, which did not exceed 6×10²⁸ dyn-cm.Finally, the study of the 1960 Chilean and 1964 Alaskan earthquakes whose exceptionally large moments are properly retrieved throughM _m measurements, serves proof that this approach performs flawlessly even for the very greatest earthquakes, and is therefore successful in its goal to avoid the saturation effects plaguing any magnitude scale measured at a fixed period.     

M TSU : Recovering Seismic Moments from Tsunameter Records

We define a new magnitude scale, M_TSU, allowing the quantification of the seismic moment M₀ of an earthquake based on recordings of its tsunami in the far field by ocean-bottom pressure sensors (``tsunameters') deployed in ocean basins, far from continental or island shores which are known to affect profoundly and in a nonlinear fashion the amplitude of the tsunami wave. The formula for M_TSU, M_TSU = log₁₀ M₀ − 20 = log₁₀ X (ω) + C_D^TSU + C_S^TSU + C₀, where X (ω) is the spectral amplitude of the tsunami, C_D^TSU a distance correction and C_S^TSU a source correction, is directly adapted from the mantle magnitude M_m introduced for seismic surface waves by Okal and Talandier. Like M_m, its corrections are fully justified theoretically based on the representation of a tsunami wave as a branch of the Earth's normal modes. Even the locking constant C₀, which may depend on the nature of the recording (surface amplitude of the tsunami or overpressure at the ocean floor) and its units, is predicted theoretically. M_TSU combines the power of a theoretically developed algorithm, with the robustness of a magnitude measurement that does not take into account such parameters as focal geometry and exact depth, which may not be available under operational conditions in the framework of tsunami warning. We verify the performance of the concept on simulations of the great 1946 Aleutian tsunami at two virtual gauges, and then apply the algorithm to 24 records of 7 tsunamis at DART tsunameters during the years 1994–2003. We find that M_TSU generally recovers the seismic moment M₀ within 0.2 logarithmic units, even under unfavorable conditions such as excessive focal depth and refraction of the tsunami wave around continental masses. Finally, we apply the algorithm to the JASON satellite trace obtained over the Bay of Bengal during the 2004 Sumatra tsunami, after transforming the trace into a time series through a simple ad hoc procedure. Results are surprisingly good, with most estimates of the moment being over 10²⁹ dyn-cm, and thus identifying the source as an exceptionally large earthquake.