Geology ReportsSearch

Geology topics

W.H.K. Lee

Publications and source records attributed to W.H.K. Lee.

At least 19 recordsLinked to original sources

Seismic guided waves trapped in the fault zone of the Landers, California, earthquake of 1992

A mobile seismic array of seven stations was deployed at 11 sites along the fault trace of the M 7.4 Landers earthquake of June 28, 1992, with a maximum offset of 1 km from the trace. We found a distinct wave train with a relatively long period following the S waves that shows up only when both the stations and the events are close to the fault trace. This wave train is interpreted as a seismic guided wave trapped in a low‐velocity fault zone. To study the distribution of amplitude of the guided waves with distance from the fault trace and also their attenuation with travel distance along the fault zone, we eliminated source and recording site effects by the coda normalization method. The normalized amplitudes of guided waves show a spectral peak at 3–4 Hz, which decays sharply with distance from the fault trace. Spectral amplitudes at high frequencies (8–15 Hz) show an opposite trend, increasing with distance from the fault trace. The normalized amplitudes of guided waves at 3–4 Hz also show a systematic decrease with hypocentral distance along the fault zone, from which we infer an apparent Q of 50. In order to confirm the existence of the guided waves, a dense array of 31 stations was deployed at one of the 11 sites. The resultant records revealed unequivocal evidence for the existence of guided waves associated with the fault zone. By modeling the waveforms as S waves trapped in a low‐velocity waveguide sandwiched between two homogeneous half‐spaces with velocity V s = 3.0 km/s, we infer a waveguide width of about 180 m, a shear velocity of 2.0–2.2 km/s, and a Q of ∼50. Hypocenters of aftershocks with clear guided waves show a systematic distribution both laterally and with depth delineating the extent of the low‐velocity fault zone in three dimensions. We find that the zone extends to a depth of at least 10 km. This zone apparently continues to the south across the Pinto Mountain fault because guided waves are observed at stations north of the Pinto Mountain fault for earthquakes with epicenters south of it. On the other hand, the zone appears to be discontinuous at the fault bend located about 20 km north of the mainshock epicenter; guided waves were observed for stations and epicenters which are located on the same sides of the fault bend but not for those on the opposite sides.

California

Determination of three‐dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes: 1. A homogeneous initial model

Geiger's method of locating local earthquakes has been extended to include the effect of P velocity variation along the ray paths in three dimensions. The crustal structure was modeled by rectangular blocks, and a parameter was assigned to each block describing the perturbation of P wave slowness in the block. On the basis of an initial model, a set of linear equations for the observed first P arrival times was formulated in terms of the source and medium parameters. The source parameters for all the earthquakes in the data set and medium parameters for all the blocks penetrated by the seismic rays were then determined simultaneously by the damped least squares method. A computer program has been written for the case of a homogeneous initial medium model with constant P velocity. It has been tested with artificial data and applied to a set of local earthquakes recorded by a dense seismic array in Bear Valley, California. The resultant velocity distribution in the top 5 km shows a narrow low‐velocity zone of about 5 km/s in the San Andreas fault zone sandwiched between high‐velocity regions of about 6 km/s.

California

Relocation of local earthquakes by seismic ray tracing

Seismic ray tracing is used to relocate a set of local earthquakes recorded by a dense seismic array in Bear Valley, California. The crustal velocity structure is two-dimensional and incorporates most of the known and inferred velocity differences in the region. Relocated hypocenters fall within the San Andreas fault zone, and P residuals, without elevation or station corrections, now appear to be more internally consistent in the rms sense. Ray takeoff angles and azimuths are distinctly different from those of a uniformly layered earth model. Focal mechanism solutions using only traced data are well determined and have few inconsistent data. Zones of convergence and divergence of rays suggest large variations in observed amplitudes across the local network.

California

Rotational motions for teleseismic surface waves

[1] We report the findings for the first teleseismic six degree‐of‐freedom (6‐DOF) measurements including three components of rotational motions recorded by a sensitive rotation‐rate sensor (model R‐1, made by eentec) and three components of translational motions recorded by a traditional seismometer (STS‐2) at the NACB station in Taiwan. The consistent observations in waveforms of rotational motions and translational motions in sections of Rayleigh and Love waves are presented in reference to the analytical solution for these waves in a half space of Poisson solid. We show that additional information (e.g., Rayleigh wave phase velocity, shear wave velocity of the surface layer) might be exploited from six degree‐of‐freedom recordings of teleseismic events at only one station. We also find significant errors in the translational records of these teleseismic surface waves due to the sensitivity of inertial translation sensors (seismometers) to rotational motions. The result suggests that the effects of such errors need to be counted in surface wave inversions commonly used to derive earthquake source parameters and Earth structure.

Geophysical Research Letters

International heat flow commission celebrates 40 years

The outflow of heat from the Earth's interior is, in terms of energy, the most impressive terrestrial phenomenon. Its present rate of about 1021 joules per year is order‐of‐magnitudes greater than the energy dissipation of earthquakes or heat loss from volcanic eruptions. The study of the Earth's internal heat plays an important role in understanding the Earth's origin, internal constitution, and plate tectonics. The thermal structure of the Earth, generally referred to as the geothermics, or also as the (terrestrial) heat flow, has been studied for a long time. However, modern geothermics, one of the fundamental geophysical disciplines, is relatively young. The International Heat Flow Commission (IHFC), under which the academic geothermal research on the international scale is organized, was created only in 1963. The IHFC, operating under the International Association of Seismology and Physics of the Earth's Interior (IASPEI), covers a broad scope of geophysical studies, and links the activities of other associations of the International Union of Geodesy and Geophysics (IUGG).

Eos, Earth and Space Science News

Historical seismograms for unravelling a mysterious earthquake: The 1907 Sumatra Earthquake

History of instrumental seismology is short. Seismograms are available only for a little more than 100 years; high-quality seismograms are available only for the last 50 years and the seismological database is very limited in time. To extend the database, seismograms of old events are of vital importance. Many unusual earthquakes are known to have occurred, but their seismological characteristics are poorly known. The 1907 Sumatra earthquake is one of them (1907 January 4, M = 7.6). Gutenberg and Richter located this event in the outer-rise area of the Sunda arc. This earthquake is known to be anomalous because of its extensive tsunami, which is disproportionate of its magnitude. The tsunami affected the coastal areas over 950 km along the Sumatran coast. We investigated this earthquake using the historical seismograms we could collect from several seismological observatories. We examined the P -wave arrival times listed in the Strassburg Bulletin (1912) and other station bulletins. The scatter of the Observed−Computed traveltime residuals ranges from –30 to 30 s, too large to locate the event accurately. The uncertainty of the epicentre estimated from an S-P grid-search relocation study is at least 1° (∼110 km). We interpreted the Omori seismograms from Osaka, Mizusawa and Tokyo, and the Wiechert seismograms from Göttingen and Uppsala by comparing them with the seismograms simulated from modern broad-band seismograms of the 2002, 2008 and two 2010 Sumatra earthquakes which occurred near the 1907 earthquake. From the amplitude of Rayleigh waves recorded on the Omori seismograms we conclude that the magnitude of the 1907 earthquake at about 30 to 40 s is about 7.8 (i.e. 7.5 to 8.0). The SH waveforms recorded on the Göttingen and Uppsala seismograms suggest that the 1907 earthquake is a thrust earthquake at a shallow depth around 30 km. The most likely scenario is that the 1907 earthquake initiated on the subduction interface, and slowly ruptured up-dip into the shallow sediments and caused the extensive tsunami. Although their quantity and quality are limited, historical seismograms provide key quantitative information about old events that cannot be obtained otherwise. This underscores the importance of preserving historical seismograms.

Geophysical Journal International

Program and abstracts of the Second Tsunami Source Workshop; July 19-20, 2010

In response to a request by the National Oceanic and Atmospheric Administration (NOAA) for computing tsunami propagations in the western Pacific, Eric Geist asked Willie Lee for assistance in providing parameters of earthquakes which may be future tsunami sources. The U.S. Geological Survey (USGS) Tsunami Source Working Group (TSWG) was initiated in August 2005. An ad hoc group of diverse expertise was formed, with Steve Kirby as the leader. The founding members are: Rick Blakely, Eric Geist, Steve Kirby, Willie Lee, George Plafker, Dave Scholl, Roland von Huene, and Ray Wells. Half of the founding members are USGS emeritus scientists. A report was quickly completed because of NOAA's urgent need to precalculate tsunami propagation paths for early warning purposes. It was clear to the group that much more work needed to be done to improve our knowledge about tsunami sources worldwide. The group therefore started an informal research program on tsunami sources and meets irregularly to share ideas, data, and results. Because our group activities are open to anyone, we have more participants now, including, for example, Harley Benz and George Choy (USGS, Golden, Colo.), Holly Ryan and Stephanie Ross (USGS, Menlo Park, Calif.), Hiroo Kanamori (Caltech), Emile Okal (Northwestern University), and Gerard Fryer and Barry Hirshorn (Pacific Tsunami Warning Center, Hawaii). To celebrate the fifth anniversary of the TSWG, a workshop is being held in the Auditorium of Building 3, USGS, Menlo Park, on July 19-20, 2010 (Willie Lee and Steve Kirby, Conveners). All talks (except one) will be video broadcast. The first tsunami source workshop was held in April 2006 with about 100 participants from many institutions. This second workshop (on a much smaller scale) will be devoted primarily to recent work by the USGS members. In addition, Hiroo Kanamori (Caltech) will present his recent work on the 1960 and 2010 Chile earthquakes, Barry Hirshorn and Stuart Weinstein (Pacific Tsunami Warning Center) will present their work on tsunami warning, and Rick Wilson (California Geological Survey) will display three posters on tsunami studies by him and his colleagues.

Open-File Report

Recent advances in rotational seismology

Rotational seismology is an emerging field of study concerned with all aspects of rotational motions induced by earthquakes, explosions, and ambient vibrations. Two recent monographs (Teisseyre et al. 2006; Teisseyre et al. 2008) and a Bulletin of the Seismological Society of America special issue on Rotational Seismology and Engineering Applications (Lee, Celebi et al. 2009) are useful starting points. Rotational seismology is of interest to a wide range of disciplines, including various branches of seismology, earthquake engineering, and geodesy, as well as to physicists using Earth-based observatories for detecting gravitational waves generated by astronomical sources.

Seismological Research Letters

Observing rotational and translational ground motions at the HGSD station in Taiwan from 2007 to 2008

Because of a lack of suitable instruments, rotational ground motions have not been observed until the last decade. Rotational measurements in the near field of earthquakes in Japan ( Takeo, 1998 ) indicate that rotational ground motions are many times larger than expected from the classical elasticity theory. After failing to obtain useful rotational ground motions (using similar rotational sensors as Takeo did), we deployed a far more sensitive rotational velocity sensor (R-1) at the HGSD station in eastern Taiwan. From 7 December 2004 to 12 November 2006, several hundreds of earthquakes were recorded during our Phase 1 operation. This was mostly a learning exercise to solve field operation problems; Phase 1 operations ended when our two R-1 sensors ceased to operate. A K2+R1 instrument was deployed in the spring of 2007 to start our Phase 2 operation. From 8 May 2007 to 17 February 2008, we observed 52 local earthquakes with good rotational velocity signals (with signal-to-noise ratio >∼5), together with excellent translational acceleration signals (with signal-to-noise ratio >∼10). Unfortunately, field operation was interrupted due to flooding of the HGSD station site in mid-February 2008; we just resumed normal operation in June 2008.

Bulletin of the Seismological Society of America

Recording rotational and translational ground motions of two TAIGER explosions in northeastern Taiwan on 4 March 2008

Two explosions were set off on 4 March 2008 at the N3 explosion site in northeastern Taiwan. The code name for the first shot with 3000 kg explosives is N3P and that for the second shot with 750 kg explosives is N3. To record these two explosions, 8 triaxial rotational sensors, 13 triaxial accelerometers, and 12 six-channel, 24 bit dataloggers with Global Positioning System receivers were deployed to continuously record several hours before and after the explosions. These instruments were installed at about 250 m (1 station), 500 m (11 stations), and 600 m (1 station) from the explosions. The 11 stations form a center array with station spacing of about 5 m. Except for one rotational sensor, onscale records were obtained. Although the N3P shot used four times larger amounts of explosives than those used for the N3 shot, the peak ground translational acceleration and rotational velocity at the 13 station sites from the N3P shot are only about 1.5 times larger than those for the N3 shot. We also observed large variations (by tens of percent) of translational accelerations and rotational velocities at the center array with station spacing of about 5 m. The largest peak rotational velocity was observed for the x component: 2.74 and 1.75 mrad/sec at a distance of 254 m from the N3P and N3 shots, respectively. The main purpose of this article is to document our recordings of rotational and translation motions from two explosions in Taiwan and to release the data online for open access. The translational acceleration data from this experiment have been analyzed by Langston et al. (2009) , and we plan to submit an article with analysis of the rotational velocity data in the future.

Bulletin of the Seismological Society of America

Array deployment to observe rotational and translational ground motions along the Meishan Fault, Taiwan: A progress report

Measurements in the near field of earthquakes in Japan ( Takeo, 1998 ) and in Taiwan ( Huang et al. , 2006 ; Liu et al. , 2009 ) indicate that rotational ground motions are many times larger than expected from the classical elasticity theory. The Central Weather Bureau (CWB) deployed four instrumentation sets on the campus of the National Chung-Cheng University, Chia-yi, Taiwan, in order to study in detail earthquakes in the near field for both seismology and earthquake engineering. Both rotational and translational ground motions are being monitored along the active Meishan fault, where a major earthquake occurred in 1906, more than a century ago. The deployed instruments are: (1) a 32-element seismic array in the free field, (2) a 32-element seismic array in a building, (3) a six-channel unit with a broadband seismometer and an accelerometer, and (4) a six-channel unit with an accelerometer and an external rotational velocity sensor.

Bulletin of the Seismological Society of America

Seismic-wave strain, rotation, and gradiometry for the 4 March 2008 TAIGER explosions

Acceleration spatial gradients, horizontal strains, and horizontal rotation were computed using strong-motion array data from the 4 March 2008 TAIGER explosions in northeastern Taiwan and used in conjunction with the original three component acceleration data to perform a gradiometric analysis of the strong ground motion wave train. The analysis yields a complex, frequency-dependent view of the nature of seismic-wave propagation over short propagation distances that imply significant lateral velocity changes in structure. Areal strain and rotation about the vertical axis have equal amplitudes and suggest significant wave scattering within the confines of the river valley where the experiment was performed and/or significant departure from an axisymmetric explosion source. Gradiometry shows that the P wave arrives at the array 35° off-azimuth clockwise from the straight-line path and appears to have been refracted from the northern side of the valley. Large, slowly propagating secondary surface waves initially arrive 45° counterclockwise from the straight-line path but later arrivals are seen to propagate in all directions, including back toward the explosion source. A frequency-dependent radiation pattern for the triple-borehole explosion in comparison to the single-borehole explosion explains the differences in the maximum amplitudes between the sources seen in the acceleration data. The use of seismic strain and rotation with standard particle motion wave fields at a single location allows for a direct view of seismic-wave propagation that illuminates the true nature of the seismogram.

Bulletin of the Seismological Society of America

A glossary for rotational seismology

Some common terms for rotational seismology are compiled from contributions by many authors, including R. DeSalvo, J. R. Evans, E. F. Grekova, C. R. Hutt, H. Igel, C. A. Langston, B. Lantz, E. Majewski, R. Nigbor, J. Pujol, P. Spudich, R. Teisseyre, M. D. Trifunac, R. J. Twiss, and Z. Zembaty. In addition, some glossary terms about earthquakes are excerpted from Aki and Lee (2003), Lee and Wu (2009), and others for the benefit of readers who are not seismologists. It is a first attempt to compile a glossary for rotational seismology (including some terms for continuum mechanics), and it is by no means complete. Several iterations will be required to reach a satisfactory glossary, and this compilation is just a start.

Bulletin of the Seismological Society of America

Suggested readings in continuum mechanics and earthquake seismology

Rotational seismology is a newly emerging field of interest to scientists from widely differing disciplines. We compile here some suggested readings for those wishing to become familiar with disciplines relevant to rotational seismology and its engineering applications. These readings are not exhaustive but contain a number of basic references with occasional annotations. We concentrate on two areas: continuum mechanics and earthquake seismology. Continuum mechanics is the branch of mechanics that deals with the deformation and motion of materials, including various elasticity theories. Although classical elasticity works well for studying earthquakes in the far field, there are alternate elasticity theories that may be more appropriate for studying earthquakes in the near field. The readings in continuum mechanics are intended for seismologists and include (1) fundamental textbooks, (2) intermediate to advanced books, (3) heteromodular theories, and (4) Cosserat theory. A brief introduction to continuum mechanics for seismologists is included as an appendix. In addition, because we wish to attract workers in other disciplines to rotational seismology, we include a short list of suggested readings in earthquake seismology.

Bulletin of the Seismological Society of America

Introduction to the special issue on rotational seismology and engineering applications

Because of a lack of suitable instruments, rotational ground motions have not been observed until the last decade. Rotational measurements in the near field of earthquakes in Japan ( Takeo, 1998 ) indicate that rotational ground motions are many times larger than expected from the classical elasticity theory. After failing to obtain useful rotational ground motions (using similar rotational sensors as Takeo did), we deployed a far more sensitive rotational velocity sensor (R-1) at the HGSD station in eastern Taiwan. From 7 December 2004 to 12 November 2006, several hundreds of earthquakes were recorded during our Phase 1 operation. This was mostly a learning exercise to solve field operation problems; Phase 1 operations ended when our two R-1 sensors ceased to operate. A K2+R1 instrument was deployed in the spring of 2007 to start our Phase 2 operation. From 8 May 2007 to 17 February 2008, we observed 52 local earthquakes with good rotational velocity signals (with signal-to-noise ratio >∼5), together with excellent translational acceleration signals (with signal-to-noise ratio >∼10). Unfortunately, field operation was interrupted due to flooding of the HGSD station site in mid-February 2008; we just resumed normal operation in June 2008.

Bulletin of the Seismological Society of America

State of the art and progress in the earthquake early warning system in Taiwan

We report here the recent progress and future development of real-time seismic monitoring in Taiwan, especially on earthquake early warning system (EWS) developed at the Central Weather Bureau (CWB), using the telemetered signals from strong-motion instruments. Utilizing a quick magnitude determination based on the first 10 sec of signals from a virtual and sub-network which is configured automatically, we were able to reduce the earthquake rapid reporting time to about 30 sec or less. This represents a significant step towards a more realistic earthquake early warning capability. This early warning system has been in operation at CWB since 2002. Comprehensive earthquake reports have been issued mostly in less than 30 sec, with an average of about 22 sec from the origin time. At 3 km/sec for a typical crustal shear-wave velocity, the present operation is not useful if an earthquake occurs less than about 70 km from a city, but the lead time will increase to more than 10 sec for cities at distances greater than 100 km from the earthquake source. In the latter case, a lead time of several seconds will allow pre-programmed emergency response to take place prior to the arrival of strong ground shaking. In order to offer earlier warning within an epicentral distance of less than 70 km, a τ c and Pd method was experimented. We used the Taiwan Strong-Motion Instrumentation Program (TSMIP) accelerograms recorded within epicentral distances of less than 30 km from Mw > 5.0 crustal earthquakes. This method uses the first 3-sec length records from the onset of P-waves and could issue an on-site early warning within 10 sec after the event origin time, reducing the radius of the blind zone from 70 to 25 km from the epicenter.

Book chapter

The IASPEI Seismological Software Library

Since computers became widely available in the early 1960s, seismologists have been using them for data acquisition, processing, and analysis, as well as theoretical computation and modeling. For example, the book by Doornbos (1988) contains a collection of seismological algorithms with the corresponding computer programs available on tape or disk from the World Data Center A for Solid Earth Geophysics. The introduction of personal computers in the early 1980s further revolutionized the use of computers for scientific research. Instead of expensive mainframe computers that required a large staff to operate, inexpensive personal computers allowed creative applications to be implemented by individuals with a shoestring budget.

International Geophysics