Source Parameters for 1999 North Anatolian Fault Zone Aftershocks

We develop a data set of aftershock recordings of the 1999, M = 7.4 Izmit and M = 7.2 Duzce (Turkey) earthquakes to study their source parameters. We combined seismograms from 44 stations maintained by several sources (organizations) to obtain a unified data set of events (2.1 ≤ M_w ≤ 5.5). We calculate source parameters of these small earthquakes by two methods that use different techniques to address the difficulty in obtaining source spectra for small earthquakes subject to interference from site response. One method (program NetMoment (NM), Hutchings, 2004) uses spectra of direct S waves in a simultaneous inversion of local high-frequency network data to estimate seismic moment, source corner frequency (f_c), site attenuation (k) and whole-path Q. This approach takes advantage of the source commonality in all recordings for a particular earthquake by fitting a common Brune source spectrum to the data with a and individual k. The second approach (Mayeda et al., 2003) uses the coda method (CM) to obtain “nonmodel-based” source spectra and moment estimates from selected broadband recording sites. We found that both methods do well for events that allow the comparison with seismic moment estimates derived from waveform modeling. Also, source spectra obtained from the two methods are very closely matched for most of the events they have in common. We use an F test to examine the trade-off between k and f_c picks identified by the direct S-wave method. About half of the events could be constrained to have less than a 50% average uncertainty in f_c and k. We used these source spectra solutions to calculate energy and apparent stress and compare these to estimates from the selected “good quality” source spectra from CM. Both studies have values mutually consistent and show a similar increase in apparent stress with increasing moment. This result has added merit due to the independent approaches to calculate apparent stress. We conclude that both methods are at least partially validated by our study, and they both have usefulness for different circumstances of recording local small earthquakes. CM would work well in studies for which there is a broad magnitude range of events and NM works well for local events recorded by band-limited recorders.     

Automatic detection,location and quantification of earthquakes: Application to tsunami warning

We describe a fully automated seismic event detection and location system, providing for real-time estimates of the epicentral parameters of both local and distant earthquakes. The system uses 12 telemetered short-period stations, with a regional aperture of 350 km, as well as two 3-component broad-band stations. Detection and location of teleseismic events is achieved independently and concurrently on the short-period and long-period channels. The long-period data is then used to obtain an estimate of the seismic momentM ₀ of the earthquake through the mantle magnitudeM _m, as introduced byOkal andTalandier (1989). In turn, this estimate ofM ₀ is used to infer the expected tsunami amplitude at Papeete, within 15 minutes of the recording of Rayleigh waves. The performance of the method is discussed in terms of the accuracy of the epicentral parameters and seismic moment obtained in real time, as compared to the values later published by the reporting agencies. Our estimates are usually within 3 degrees of the reported epicenter, and the standard deviation on the seismic moment only 0.19 unit of magnitude for a population of 154 teleseismic events.     

Die Berechnung der Globalstrahlung aus Trübungswert und Bewölkung

Zusammenfassung Es werden einfache halbempirische Formeln zur Berechnung der Globalstrahlung angegeben. Die Abhängigkeit der Globalstrahlung von den meteorologischen Parametern Trübungsgrad (Linke'scher TrübungsfaktorT _g) sowie Bewölkungsart und-bedeckung wurde empirisch aus Messungen in Frankfurt (Main) und Potsdam ermittelt. Die nach diesen Formeln berechnete Globalstrahlungsleistung in ly min^–1 und die 24-stündige Globalstrahlungsenergiesumme in ly wird mit Messergebnissen anderer Autoren in den Tropen, bzw. in Polargebieten, verglichen, wobei sich eine sehr gute Übereinstimmung zeigt.

---

Summary Simple halfempirical formulas for the calculation of the global radiation are presented. The dependence of global radiation on the meteorological parameters: turbidity (turbidity factorT _g afterLinke) and type and amount of cloud was find empirically from the measurements in Frankfurt (Main) and Potsdam. The flux of global radiation (in ly min^–1) and the daily sum of global radiation energy, calculated by these formulas, was compared with measured values of other authors in the tropics and in polar regions, respectively, and shows satisfactory agreement.

     

Estimation of the elastic earth parameters using the SLR LAGEOS 1 and LAGEOS 2 data

We present results for the global elastic parameters (k ₂, k ₃) and (h ₂, l ₂) derived from the analysis of Satellite Laser Ranging (SLR) data. The SLR data for two satellites, LAGEOS 1 and LAGEOS 2, collected over a period of two years, from 1 January 2005 to 1 January 2007, from 18 globally distributed ground stations are analyzed using different approaches. The analysis is performed for the two satellites jointly, and the four elastic parameters are estimated. We carry out a sequential analysis and study the stability of the estimates as a function of the length of the data set used. The adjusted final values for k ₂ (0.3011±0.0001) and k ₃ (0.0896±0.0037) and those for h ₂ (0.6146±0.0006) and l ₂ (0.0883±0.0003) for LAGEOS 1 and LAGEOS 2 tracking data are compared to the estimates derived in other studies.     

Source parameters of some presumed Semipalatinsk underground nuclear explosions

Summary Short-period vertical-componentP-wave spectra of seven presumed Semipalatinsk underground nuclear explosions, recorded by the Swedish seismic station network, are investigated. The events considered have closely spaced foci and cover the magnitude range fromm _b=5.5 tom _b=6.6. Spectra of six of these explosions show pronounced minima, varying from about 1.5 to 1.8 cps, which could be explained as principle minima due toP-pP interference. Supposing a nearsurfaceP-wave velocity at the test area of 4 km/sec, the shot depths are estimated to vary roughly from 750 to 1350 m. In order to obtain an estimate of the yield, the observed spectra are compared withHaskell's theoretical source spectra. For four events, relative yield estimates fit well the predicted values for explosions fired in a granitic medium. The behaviour of the remaining three explosions is discussed in detail.     

Estimation of the elastic earth parameters using slr data for the low satellites STARLETTE and STELLA

In this paper we present estimated values for the global elastic parameters (k ₂, k ₃) and (h ₂, l ₂) derived from the analysis of Satellite Laser Ranging (SLR) data. We analyse SLR data for two low satellites, STARLETTE and STELLA, collected over a period of two years, from 1 January 2005 to 1 January 2007, from 18 globally distributed ground stations. We carry out a sequential analysis for the two satellites jointly, and study the stability of the estimates as a function of the length of the data set used. The adjusted final values of (k ₂, k ₃) and (h ₂, l ₂) for STARLETTE and STELLA are compared to, and are largely found to support, the estimates we previously published based on data for two high satellites LAGEOS 1 and LAGEOS 2. A major discrepancy between the two solutions was only found for the Shida number l ₂.     

A consideration of Hayford's best fitting ellipsoid data using the differential change equations of Vening Meinesz

Summary The relations of the defining parameters of a geodetic system to the deflections of the vertical in this system given byVening Meinesz have been used in the computation of a best fitting ellipsoid. The data was that used byHayford in the computation of the resulting International Ellipsoid. Assuming the undulation (N ₀) at Meades Ranch to be zero, the results are: a=6,378,194 m and l/f=299.9. This a is much nearer to modern values than that of the Internaitonal Ellipsoid. If a flattening of 1/298.24 is fixed andN ₀=+3m,a=6,378,164 m which differs by one meter from a recent determination given byKaula.     

Quantitative estimates of interplate coupling inferred from outer rise earthquakes

Interplate coupling plays an important role in the seismogenesis of great interplate earthquakes at subduction zones. The spatial and temporal variations of such coupling control the patterns of subduction zone seismicity. We calculate stresses in the outer rise based on a model of oceanic plate bending and coupling at the interplate contact, to quantitatively estimate the degree of interplate coupling for the Tonga, New Hebrides, Kurile, Kamchatka, and Marianas subduction zones. Depths and focal mechanisms of outer rise earthquakes are used to constrain the stress models. We perform waveform modeling of body waves from the GDSN network to obtain reliable focal depth estimates for 24 outer rise earthquakes. A propagator matrix technique is used to calculate outer rise stresses in a bending 2-D elastic plate floating on a weak mantle. The modeling of normal and tangential loads simulates the total vertical and shear forces acting on the subducting plate. We estimate the interplate coupling by searching for an optimal tangential load at the plate interface that causes the corresponding stress regime within the plate to best fit the earthquake mechanisms in depth and location.We find the estimated mean tangential load over 125–200 km width ranging between 166 and 671 bars for Tonga, the New Hebrides, the Kuriles, and Kamchatka. This magnitude of the coupling stress is generally compatible with the predicted shear stress at the plate contact from thermal-mechanical plate models byMolnar andEngland (1990), andVan den Buekel andWortel (1988). The estimated tectonic coupling,F _tc , is on the order of 10¹²–10¹³ N/m for all the subduction zones.F _tc for Tonga and New Hebrides is about twice as high as in the Kurile and Kamchatka arcs. The corresponding earthquake coupling forceF _ec appears to be 1–10% of the tectonic coupling from our estimates. There seems to be no definitive correlation of the degree of seismic coupling with the estimated tectonic coupling. We find that outer rise earthquakes in the Marianas can be modeled using zero tangential load.     

Instantaneous photochemical rates in the global stratosphere

Starting with the average actual distribution of ozone (Dütsch [15]) and temperature in the stratosphere, we have calculated the solar intensity as a function of wavelength and the instantaneous rates (molecules cm^–3 sec^–1) for each Chapman reaction and for each of several reactions of the oxides of nitrogen. The calculation is similar to that ofBrewer andWilson [5]. These reaction rates were calculated independently in each volume element in spherical polar coordinates defined by R=1 km from zero to 50, =5° latitude, and ø=15° longitude (thus including day and night conditions). Calculations were made for two times: summer-winter (January 15) and spring-fall (March 22). As input data we take observed solar intensities (Ackerman [1]) and observed, critically evaluated. constants for elementary chemical and photochemical reactions; no adjustable parameters are employed. (These are not photochemical equilibrium calculations.) According to the Chapman model, the instantaneous, integrated, world-wide rate of formation of ozone from sunlight is about five times faster than the rate of ozone destruction, and locally (lower tropical stratosphere) the rate of ozone formation exceeds the rate of destruction by a factors as great as 1000. The global rates of increase of ozone are more than 50 times faster thanBrewer andWilson's [5] estimate of the average annual transfer rate of ozone to the troposphere. The rate constants of the Chapman reactions are believed to be well-enough known that it is highly improbable that these discrepancies are, due to erroneous rate constants. It is concluded that something else besides neutral oxygen species is very important in stratospheric ozone photochemistry. The inclusion of a uniform concentration of the oxides of nitrogen (NO_x as, NO and NO₂) averaging 6.6×10^–9 mole fraction gives a balance between global ozone formation and destruction rates. The inclusion of a uniform mole fraction of NO_x at 28×10^–9 also gives a global balance. These calculations support the hypethesis (Crutzen [10],Johnston [24]) that the oxides of nitrogen are the most important factor in the global, natural ozone balance. Several authors have recently evaluated the natural source strength of NO_x in the stratosphere; the projected fleets of supersonic transports would constitute an artificial source of NO_x about equal to the natural value, thus promising more or less to double an active natural stratospheric ingredient.     

On coda waves of earthquakes

Summary Coda waves viz. the tail portion of an earthquake record have been observed and analysed byCarder, Macelwane and others. They showed that the periods of such waves increase with the increase of epicentral distances.Carder observed that these waves have very little transverse component so that these may be considered as of the type of Rayleigh waves. RecentlyOmote showed that the Coda waves contain three periodsT ₁,T ₂,T ₃ of whichT ₁ increases with epi-central distances as observed by previous observers. ButT ₂,T ₃ remain constant for all earthquakes from different epicentral distances.Omote tried to explain this phenomenon by considering that the surface of the earth consists of several layers andT ₂,T ₃ are free oscillation periods of the surface layers.T ₁ period has been explained bySezawa and also byJeffreys which has been shown byGutenberg. The author has attempted to explain the periodsT ₂,T ₃ by considering passage of cracks at the focal region. The Rayleigh wave character of Coda waves and low velocity of such waves have been explained.     

Hydrocarbon source systems and formation of gas fields in Sichuan Basin

The formation of large and middle gas fields in Sichuan Basin is investigated based on source wntrolling theory and hydrocarbon source systems. It is indicated that Є₁, S_l, P₁, P₂ and T₃ are the main source beds and Є₁/Z₂ d, C₂ h/S₁, P₁/P₂, P₂ ch/P₂, T_1,2/P,T3 x /T3 x are important hydrocarbon source systems in the basin. All these source systems are the prospective formations and exploration spaces of large and middle gas fields. It is also emphasized that hydrocarbon generation intensity is the most important geochemical factor to estimate large and middle gas fields. Project supported by the “85–102” Chinese National Key Science and Technology Project.     

Globalstrahlung und Sonnenaktivität

Summary Using the highest global irradiance values (Sun plus sky) at noon derived from the records gathered during the years 1958 to 1963 at Genova, Monte Capellino and Sauze d'Oulx, it is demonstrated that the thermal energy which the Earth receives from the Sun is dependent from the Sun activity. Precisely, an increase of the thermal energy withR _Z occurs up toR _Z=160–170, whereas for greater values ofR _Z the thermal energy received decreases with the increase ofR _Z.

---

---

Herrn Prof. Dr.Hans Ertel zu seinem 60. Geburtstag gewidmet.     

Earthquake damage detection in the Imperial County Services Building III: Analysis of wave travel times via impulse response functions

The majority of structural health monitoring methods are based on detecting changes in the modal properties, which are global characteristics of the structure, and are not sensitive to local damage. Wave travel times between selected sections of a structure, on the other hand, are local characteristics, and are potentially more sensitive to local damage. In this paper, a structural health monitoring method based on changes in wave travel times is explored using strong motion data from the Imperial Valley Earthquake of 1979 recorded in the former Imperial County Services (ICS) Building, severely damaged by this earthquake. Wave travel times are measured from impulse response functions computed from the recorded horizontal seismic response in three time windows—before, during, and after the largest amplitude response, as determined from previous studies of this building, based on analysis of novelties in the recorded response. The results suggest initial spatial distribution of stiffness consistent with the design characteristics, and reduction of stiffness following the major damage consistent with the spatial distribution of the observed damage. The travel times were also used to estimate the fundamental fixed-base frequency of the structure f₁ (assuming the building deformed as a shear beam), and its changes during this earthquake. These estimates are consistent with previous estimates of the soil–structure system frequency, f_sys, during the earthquakes (f₁<f_sys as expected from soil–structure interaction studies), and with other estimates of frequency (f₁ from ETABS models, and f_sys from ambient vibration tests, and “instantaneous” f₁ from high-frequency pulse propagation).     

Comparison of techniques that use the single scattering model to compute the quality factorQ from coda waves

Using simulated coda waves, the resolution of the single-scattering model to extract codaQ (Q _c ) and its power law frequency dependence was tested. The back-scattering model ofAki andChouet (1975) and the single isotropic-scattering model ofSato (1977) were examined. The results indicate that: (1) The inputQ _c models are reasonably well approximated by the two methods; (2) almost equalQ _c values are recovered when the techniques sample the same coda windows; (3) lowQ _c models are well estimated in the frequency domain from the early and late part of the coda; and (4) models with highQ _c values are more accurately extracted from late code measurements.     

ML of the 1989 Loma Prieta earthquake

The use of regional attenuation in computing the local magnitude, M_L, from strong motion data gathered at distances less than 100 km may lead to systematic underestimates approaching 0·5 magnitude units (Trifunac & Herak, Soil Dynamics and Earthquake Engineering, 1992, 18, 229-41). The use of the attenuation law Att(Δ), for example, with synthetic estimates of Wood-Anderson seismometer response, during the Loma Preita earthquake, leads to estimates of M_L which agree with the surface wave and moment magnitudes, and which are essentially distance-independent.     

Depression in the low-latitude horizontal intensity due to quiet-time ring current

Summary The external field due to plasma within the magnetosphere has been computed as a function ofA _p, which is a measure of solar wind velocity, for very quiet to slightly disturbed conditions using mean daily horizontal intensity from 1932 to 1968 at Alibag. The intensity, corrected for secular change and reduced to a common epoch, showed initially a small increase withA _p followed by a steady depression with further increase in the index. ForA _p7.5, which is representative of conditions over the 33-hour interval during which data relating to low-energy protons were acquired and used byHoffman andBracken [4]²) to compute current distributions, the decrease, computed here from surface data, is 6 . This is in goodagreement with the southward directed field of the quiet-time proton belt 9±5 obtained byHoffman andBracken.     

Soil-structure interaction effects on an instrumented building

The paper highlights the use of fem and bi-directional lumped-mass-storey-stiffness numerical models for the study of the soil–structure interaction (ssi) effects on an instrumented building. Data on the structural response have been obtained through the project for seismic instrumentation of a 16-storey r/c cast-in-place dwelling building (Chisinau, Republic of Moldova) during a series of earthquakes (Gutenberg–Richter M _GR = 5.0−6.7). The effect of soil–structure interaction is clearly observed comparing the responses recorded on foundation and free-field. ssi becomes more pronounced for higher level of ground shaking amplifying the natural period of the structure and slightly suppressing high frequences on the foundation in comparison with the free-field motion.     

Regional moments,energy levels,and a new discriminant

Teleseismic observations of explosions tend to be richer in short-period energy than are earthquakes, thus the effectiveness of them _b M _s discriminant. At regional distances the same basic separation occurs for smaller events in terms ofM _L M ₀ (Woods et al., 1993) andm _b M ₀ (Patton andWalter, 1993). While these studies demonstrate the basic differences in excitation, they suffer in practical application because of the detailed information required in the retrieval ofM ₀ . In this paper, we introduce a new method of discrimination, based on the energy strength (M _E ) from broadband regional records that appears to be effective and efficient. In this method all events are processed as earthquakes, and explosions are distinguished by their stronger energy levels relative to their long-period amplitudes. Results from 29 events recorded by TERRAscope, sampling 15 explosions from NTS and 14 earthquakes from the southwestern United States, are represented, indicating complete separation (45 data points).M _L =3.6 is the smallest event examined to date but the method can probably be extended to even smaller levels in calibrated regions.     

Inversion of source process and related studies of the Taiwan Strait earthquake using genetic algorithm

Applying genetic algorithm to inversion of seismic moment tensor solution and using the data of P waveform from digital network and initial motion directions of P waves of Taiwan network stations, we studied the moment tensor solutions and focal parameters of the earthquake of M=7.3 on 16 September of 1994 in Taiwan Strait and other four quakes of M _L≥5.8 in the near region (21°–26°N, 115°–120°E). Among the five earthquakes, the quake of M=7.3 on September 16, 1994 in Taiwan Strait is the strongest one in the southeastern coast area since Nan’ao earthquake of M=7.3 in 1918. The results show that moment tensor solution of M=7.3 earthquake is mainly double-couple component, and is normal fault whose fault plane is near NW. The strike of the fault plane resembles that of the distributive bands of earthquakes before the main event and fracture pattern shown by aftershocks. The tension stress axis of focal mechanism is about horizontal, near in NE strike, the compressive stress axis is approximately vertical, near in NWW strike. It seems that this quake is controlled by the force of Philippine plate’s pressing Eurasian plate in NW direction. But from the viewpoint of P axis of near vertical and T axis of near horizontal, it is a normal fault of strong tensibility. There are relatively big difference between focal mechanism solution of this quake and those of the four other strong quakes. The complexity of source mechanism solution of these quakes represents the complexity of the process of the strait earthquake sequences. Contribution No. 98A01001, Institute of Geophysics, State Seismological Bureau, China. The subject is supported and helped by Academician Yun-Tai CHEN, Profs. Qing-Yao HONG, Zhen-Xing YAO, Tian-Yu ZHENG, Yao-Lin SHI, Ji-An XU, Bo-Shou HUANG and colleague Mei-Jian AN, Xue-Reng DING, Rui-Feng LIU. De-Chong ZHANG and Ming Li provided the digital data warm-heartedly. Lin-Ying WANG offered us the catalogue of earthquakes in southeastern coastal area in China. Xi-Li WANG and Tong-Xia BAI provided us the issued annual reports data. The authors would like to express their gratitude to all of these people. This paper is sponsored by the National Natural Science Foundation of China and Scientific and Technological Commission of Shantou, Guangdong Province.     

Substitute conductors for electromagnetic response estimates

Various concepts exist to define substitute conductors for empirical response estimates at singular frequencies: Chapman's shell-core model, the Cagniard-Tikhonov apparent resistivity, the Niblett-Bostick and Molochnov transformation, thep ^* —z ^* transformation. They are all interrelated and assign comparable resistivities to the substitute conductor at a given frequency. Applications to synthetic response data of plane and spherical conductors show under which conditions these substitutions come closest to the model and which influence of source dimensions and Earth's sphericity can be expected.p ^* —z ^* transformed global response data forS andDst variations demonstrate how substitute conductors may serve as useful guides in inverse procedures.