Deterministic signals in precipitation in Arkansas,Louisiana, and Mississippi

Evidence for a signal near 19 years is found in 51 out of 65 rain gauge records in the tri-state region of Arkansas, Louisiana, and Mississippi, While another term with period 10 to 11 years is found in 39 instances. These are identified as the 18.6-year luni-solarM _n and solar cycleS _c signals in climate variability. In terms of both amplitude and phase the luni-solar wave trains are highly nonstationary; in particular, the waves commonly change phase by 180°. Examples of theM _n andS _c terms in corn yield for Arkansas and Missouri are presented, and the implications for economic science briefly discussed.     

Modeling of the runup heights of the Hokkaido-Nansei-Oki tsunami of 12 July 1993

The Hokkaido-Nansei-Oki earthquake (M _w 7.7) of July 12, 1993, is one of the largest tsunamigenic events in the Sea of Japan. The tsunami magnitudeM _t is determined to be 8.1 from the maximum amplitudes of the tsunami recorded on tide gauges. This value is larger thanM _w by 0.4 units. It is suggested that the tsunami potential of the Nansei-Oki earthquake is large forM _w . A number of tsunami runup data are accumulated for a total range of about 1000 km along the coast, and the data are averaged to obtain the local mean heightsH _n for 23 segments in intervals of about 40 km each. The geographic variation ofH _n is approximately explained in terms of the empirical relationship proposed byAbe (1989, 1993). The height prediction from the available earthquake magnitudes ranges from 5.0–8.4 m, which brackets the observed maximum ofH _n , 7.7 m, at Okushiri Island.     

Case 25 application of the concentration parameter of seismoactive faults to Southern California

In this paper we evaluate the present state of the seismic regime in Southern California using the concentration parameter of seismogenic faults (K _sf ,Sobolev andZavyalov, 1981). The purpose of this work is to identify potential sites for large earthquakes during the next five or ten years. The data for this study derived from the California Institute of Technology's catalog of southern California earthquakes, and spanned the period between 1932 to June 1982. We examined events as small asM _L 1.8 but used a magnitude cutoff atM _L =3.3 for a detailed analysis. The size of the target earthquakes (M _M ) was chosen as 5.3 and 5.8.The algorithm for calculatingK _sf used here was improved over the algorithm described bySobolev andZavyalov (1981) in that it considered the seismic history of each elementary seismoactive volume. The dimensions of the elementary seismoactive volumes were 50 km×50 km and 20 km deep. We found that the mean value ofK _sf within 6 months prior to the target events was 6.1±2.0 for target events withM _L 5.3 and 5.41.8 for targets withM _L 5.8. Seventy-three percent of the targets withM _L 5.8 occurred in areas whereK _sf was less than 6.1. The variance of the time between the appearance of areas with lowK _sf values and the following main shocks was quite large (from a few months to ten years) so this parameter cannot be used here for accurate predictions of occurrence time.Regions where the value ofK _sf was below 6.1 at the end of our data set (June, 1982) are proposed as the sites of target earthquakes during the next five to ten years. The most dangerous area is the area east of San Bernardino whereK _sf values are presently between 2.9 and 3.7 and where there has been no earthquake withM _L 5.3 since 1948.     

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.     

Partial breaking of a mature seismic gap: The 1987 earthquakes in New Britain

To better understand the mechanics of subduction and the process of breaking a mature seismic gap, we study seismic activity along the western New Britain subduction segment (147°E–151°E, 4°S–8°S) through earthquakes withm _b 5.0 in the outer-rise, the upper area of subducting slab and at intermediate depths to 250 km, from January 1964 to December 1990. The segment last broke fully in large earthquakes of December, 28, 1945 (M _s =7.9) and May 6, 1947 (M _s =7.7.), and its higher seismic potential has been recognized byMcCann et al., (1979). Recently the segment broke partially in two smaller events of February, 8, 1987 (M _s =7.4) and October 16, 1987 (M _s =7.4), leaving still unbroken areas.We observe from focal mechanisms that the outer-rise along the whole segment was under pronounced compression from the late 60's to at least October 1987 (with exception of the tensional earthquake of December 11, 1985), signifying the mature stage of the earthquake cycle. Simultaneously the slab at intermediate depths below 40 km was under tension before the earthquake of October 16, 1987. That event, with a smooth rupture lasting 32 sec, rupture velocity of 2.0 km/sec, extent of approximately 70 km and moment of 1.2×10²⁷ dyne-cm, did not change significantly the compressive state of stress in the outer-rise of that segment. The earthquake did not fill the gap completely and this segment is still capable of rupturing either in an earthquake which would fill the gap between the 1987 and 1971 events, or in a larger magnitude event (M _s =7.7–7.9), comparable to earthquakes observed in that segment in 1906, 1945 and 1947.     

Relative seismic moment tensor determination for Vrancea intermediate depth earthquakes

The method of relative seismic moment tensor determination proposed byStrelitz (1980) is extended a) from an interactive time domain analysis to an automated frequency domain procedure, and b) from an analysis of subevents of complex deep-focus earthquakes to the study of individual source mechanism of small events recorded at few stations.The method was applied to the recovery of seismic moment tensor components of 95 intermediate depth earthquakes withM _L=2.6–4.9 from the Vrancea region, Romania. The main feature of the obtained fault plane solutions is the horizontality ofP axes and the nonhorizontal orienaation ofT axes (inverse faulting). Those events with high fracture energy per unit area of the fault can be grouped unambiguously into three depth intervals: 102–106 km, 124–135 km and 141–152 km. Moreover, their fault plane solutions are similar to ones of all strong and most moderate events from this region and the last two damaging earthquakes (November 10, 1940 withM _W=7.8 and March 4 1977 withM _W=7.5) occurred within the third and first depth interval, respectively. This suggests a possible correlation at these depths between fresh fracture of rocks and the occurrence of strong earthquakes.     

A semi-empirical method for predicting hydrological drought magnitudes in the Canadian prairies

Abstract

A hydrological drought magnitude (M _T ) expressed in standardized terms is predicted on annual, monthly and weekly time scales for a sampling period of T years in streamflow data from the Canadian prairies. The drought episodes are considered to follow the Poisson law of probability and, when coupled with the gamma probability distribution function (pdf) of drought magnitude (M) in the extreme number theorem, culminate in a relationship capable of evaluating the expected value, E(M _T ). The parameters of the underlying pdf of M are determined based on the assumption that the drought intensity follows a truncated normal pdf. The E(M _T ) can be evaluated using only standard deviation (σ), lag-1 autocorrelation (ρ) of the standardized hydrological index (SHI) sequence, and a weighting parameter Φ (ranging from 0 to 1) to account for the extreme drought duration (L _T ), as well as the mean drought duration (L_m ), in a characteristic drought length (L_c ). The SHI is treated as standard normal variate, equivalent to the commonly-used standardized precipitation index. A closed-form relationship can be used for the estimation of first-order conditional probabilities, which can also be estimated from historical streamflow records. For all rivers, at the annual time scale, the value of Φ was found equal to 0.5, but it tends to vary (in the range 0 to 1) from river to river at monthly and weekly time scales. However, for a particular river, the Φ value was nearly constant at monthly and weekly time scales. The proposed method estimates E(M _T ) satisfactorily comparable to the observed counterpart. At the annual time scale, the assumption of a normal pdf for drought magnitude tends to yield results in close proximity to that of a gamma pdf. The M _T , when transformed into deficit-volume, can form a basis for designing water storage facilities and for planning water management strategies during drought periods.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Sharma, T.C. and Panu, U.S., 2013. A semi-empirical method for predicting hydrological drought magnitudes in the Canadian prairies. Hydrological Sciences Journal, 58 (3), 549–569.     

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.     

Stability of premonitory seismicity pattern and intermediate-term earthquake prediction in Central Italy

The algorithm CN makes use of normalized functions. Therefore the original algorithm, developed for the California-Nevada region, can be directly applied, without adjustment of the parameters, to the determination of the Time of Increased Probability (TIP) of strong earthquakes for Central Italy. The prediction is applied to the events with magnitudeMM ₀=5.6, which in Central Italy have a return period of about six years. The routinely available digital earthquake bulletins of the Istituto Nazionale di Geofisica (ING), Rome, permits continuous monitoring. Here we extend to November 1994 the first study made by Keilis-Boroket al. (1990b). On the basis of the combined analysis of seismicity and seismotectonic, we formulate a new regionalization, which reduces the total alarm time and the failures to predict, and narrows the spatial uncertainty of the prediction with respect to the results ofKeilis-Borok et al. (1990b).The premonitory pattern is stable when the key parameters of the CN algorithm and the duration of the learning period are changed, and when different earthquake catalogues are used.The anlysis of the period 1904–1940, for whichM ₀=6, allows us to identify self-similar properties between the two periods, in spite of the considerably higher seismicity level of the earlier time interval compared with the recent one.     

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.     

Probabilistic assessment of earthquake recurrence in the January 26, 2001 earthquake region of Gujrat, India

Kutch region of Gujrat is one of the most seismic prone regions of India. Recently, it has been rocked by a large earthquake (M _w = 7.7) on January 26, 2001. The probabilities of occurrence of large earthquake (M≥6.0 and M≥5.0) in a specified interval of time for different elapsed times have been estimated on the basis of observed time-intervals between the large earthquakes (M≥6.0 and M≥5.0) using three probabilistic models, namely, Weibull, Gamma and Lognormal. The earthquakes of magnitude ≥5.0 covering about 180 years have been used for this analysis. However, the method of maximum likelihood estimation (MLE) has been applied for computation of earthquake hazard parameters. The mean interval of occurrence of earthquakes and standard deviation are estimated as 20.18 and 8.40 years for M≥5.0 and 36.32 and 12.49 years, for M≥6.0, respectively, for this region. For the earthquakes M≥5.0, the estimated cumulative probability reaches 0.8 after about 27 years for Lognormal and Gamma models and about 28 years for Weibull model while it reaches 0.9 after about 32 years for all the models. However, for the earthquakes M≥6.0, the estimated cumulative probability reaches 0.8 after about 47 years for all the models while it reaches 0.9 after about 53, 54 and 55 years for Weibull, Gamma and Lognormal model, respectively. The conditional probability also reaches about 0.8 to 0.9 for the time period of 28 to 40 years and 50 to 60 years for M≥5.0 and M≥6.0, respectively, for all the models. The probability of occurrence of an earthquake is very high between 28 to 42 years for the magnitudes ≥5.0 and between 47 to 55 years for the magnitudes ≥6.0, respectively, past from the last earthquake (2001).     

Probabilistic Assessment of Earthquake Recurrence in Northeast India and Adjoining Regions

Northeast India and adjoining regions (20°–32° N and 87°–100° E) are highly vulnerable to earthquake hazard in the Indian sub-continent, which fall under seismic zones V, IV and III in the seismic zoning map of India with magnitudes M exceeding 8, 7 and 6, respectively. It has experienced two devastating earthquakes, namely, the Shillong Plateau earthquake of June 12, 1897 (M _w 8.1) and the Assam earthquake of August 15, 1950 (M _w 8.5) that caused huge loss of lives and property in the Indian sub-continent. In the present study, the probabilities of the occurrences of earthquakes with magnitude M ≥ 7.0 during a specified interval of time has been estimated on the basis of three probabilistic models, namely, Weibull, Gamma and Lognormal, with the help of the earthquake catalogue spanning the period 1846 to 1995. The method of maximum likelihood has been used to estimate the earthquake hazard parameters. The logarithmic probability of likelihood function (ln L) is estimated and used to compare the suitability of models and it was found that the Gamma model fits best with the actual data. The sample mean interval of occurrence of such earthquakes is estimated as 7.82 years in the northeast India region and the expected mean values for Weibull, Gamma and Lognormal distributions are estimated as 7.837, 7.820 and 8.269 years, respectively. The estimated cumulative probability for an earthquake M ≥ 7.0 reaches 0.8 after about 15–16 (2010–2011) years and 0.9 after about 18–20 (2013–2015) years from the occurrence of the last earthquake (1995) in the region. The estimated conditional probability also reaches 0.8 to 0.9 after about 13–17 (2008–2012) years in the considered region for an earthquake M ≥ 7.0 when the elapsed time is zero years. However, the conditional probability reaches 0.8 to 0.9 after about 9–13 (2018–2022) years for earthquake M ≥ 7.0 when the elapsed time is 14 years (i.e. 2009).     

Crustal structure of the Peninsular India

Summary Earthquake parameters for the forty aftershocks of the main Koyna earthquake of 10 December, 1967, have been determined. Depths of the foci of the earthquakes have been found to vary between 2 to 17 km. The velocities for the phasesP _g ,P ^*,P _n have been observed to be 5.78±0.00, 6.58±0.04, 8.19±0.02 km/sec, and forS _g ,S ^*,S _n to be 3.42±0.00, 3.92±0.01 and 4.62±0.01 km/sec respectively. A two-layered crustal model has been interpreted for the Peninsular shield with the average thickness of the granitic layer as 20 km and that for the basaltic layer as 18.7 km. A plot of the epicenters suggests a NNE to SSW orientation of the fault.     

The time of occurrence and the magnitude of the largest aftershock over India

Summary The time of occurrence and the magnitude of the largest aftershock in relation to the main shock have been studied for India and its neighbourhood based on the USCGS data during the years 1963–1971. It is found that the largest aftershock occurs within 2 hours after the main shock in about 50% of the cases and frequency of occurrencen(t) of the largest aftershocks decreases hyper-bolically with the intervalt after the main event and could be represented by a law of the formn(t)=At ^–h whereA andh are constants. The probability of occurrence of the largest aftershock within 2 hours of the main shock is found to be higher over island are regions of the world. The difference (M ₀–M ₁) of the magnitude of the largest aftershockM ₁ to that of the main shockM ₀ as a measure of aftershock activity does not show any marked regional variation over India and its neighbourhood, as was reported by Mogifor Japan. Examination of the values ofM ₁/M ₀ and the constantb in Gutenberg-Richter's frequency magnitude relationship reveals a range of variation in both; high values ofM ₁/M ₀ have been found to be associated with high values ofb in many tectonic earthquakes and thus not, restricted to reservoir associated seismic activity.     

The <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 3.1–4.7 earthquakes in the southern Baltic Sea and adjacent areas in 2000, 2001 and 2004

The area south and east of the Baltic Sea has very minor seismic activity. However, occasional events occur as illustrated by four events in recent years, which are analysed in this study: near Wittenburg, Germany, on May 19, 2000, M _w = 3.1, near Rostock, Germany, on July 21, 2001, M _w = 3.4 and in the Kaliningrad area, Russia, two events on September 21, 2004 with M _w = 4.6 and 4.7. Locations, magnitudes (M _L and M _w) and focal mechanisms were determined for the two events in Germany. Synthetic modeling resulted in a well-confined focal depth for the Kaliningrad events. The inversion of macroseismic observations provided simultaneous solutions of the location, focal depth and epicentral intensity. The maximum horizontal compressive stress orientations obtained from focal mechanism solutions, approximately N–S for the two German events and NNW–SSE for the Kaliningrad events, show a good agreement with the regionally oriented crustal stress field.     

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.     

Rayleigh wave magnitude scaleM s

Summary TheM _s Rayleigh wave magnitude formula is revised for purposes of eliminating the variable effects of near distances and propagation pats on the values computed from standard long period seismograms. The improved formulation employs a revised distance correction function and period dependent path correction that normalisesM _s to large teleseismic distance, 20 second values. The method for computing the path corrections is described. The magnitude scale presented here givesM _s values which are within ±0.1 magnitude units of the Gutenberg and Prague magnitude formulae.     

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.     

2017年3月渤海地震序列微震检测与发震构造分析

地震序列发震构造研究是区域地震活动性和地震危险性分析的重要基础.2017年3月渤海海域发生地震序列活动,该序列发生在郯城-庐江断裂带与张家口-渤海地震带的交汇部位,区域构造较为复杂.然而在渤海海域,连续运行的固定地震监测仪器难以布设,导致地震监测能力相对较弱.本文首先采用模板匹配方法对序列遗漏地震进行检测,再使用波形互...