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.     

Fast inversion of rupture process of the 14 April 2010 Yushu,Qinghai, earthquake

Four results of the rupture process of 14 April 2010 Yushu, Qinghai, earthquake, obtained by inverting the broadband seismic data of Global Seismographic Network (GSN) based on the inversion method of earthquake rupture process, were compared and discussed. It is found that the Yushu earthquake has several basic characteristics as follows: ① There exist two principal sub-events which correspond to two slip-concentrated patches being located near the hypocenter and to the southeast of the epicenter. The rupture of the slip-concentrated patch to the southeast of the epicenter broke though the ground surface; ② The peak slip and peak slip-rate are about 2.1 m and 1.1 m/s, respectively, indicating that the Yushu earthquake is an event with large slip-rate on the fault plane; ③ Overall the Yushu earthquake is a unilateral rupture event with the rupture mainly propagating southeastward. The strong focusing of the seismic energy in the southeast of the epicenter due to the "seismic Doppler effect" reasonably accounts for the tremendous damage in the Yushu city.     

Influences of coupled wind–water processes on suspended sediment grain size: an example from tributaries of the Yellow River / Influences de processus couplés éoliens-hydriques sur la granulométrie des matières en suspension: exemple d'affluents du Fleuve Jaune

Abstract

Some unique coupled wind–water erosion processes exist in the desert-loess transitional zone in the middle Yellow River basin. Based on data from 40 stations on 29 rivers, a study was made on the influence of such processes on suspended sediment grain-size characteristics of the tributaries of the Yellow River. Results show that the percentage of >0.05-mm grain size decreases with the increased annual mean precipitation, but increases with the increase in the annual mean number of sand-dust storm days. The percentage of <0.01-mm grain size increases with the increase in the annual mean precipitation, but decreases with the increase in the annual number of sand-dust storm days. Based on annual mean data from 40 stations, multiple regression equations were established between the percentages of >0.05-mm grain size (r _>0.05) and <0.01-mm grain size (r _<0.01), annual mean precipitation (P _m) and annual mean number of sand-dust storm days (D _ss). On this basis, the relative contributions of the variations in D _ss and P _m to the variation in r _>0.05 and r _<0.01 were estimated. The results indicate that the variation in sand-dust storm frequency exerts greater influences on the variation in grain-size characteristics of suspended load than does the variation in annual mean precipitation. With the increase in the coupled wind–water processes index, expressed by P _m/D _ss, the percentage of >0.05-mm grain size in suspended sediment decreases and the percentage of <0.01-mm grain size increases. With the variation in P _m/D _ss, different combinations of r _>0.05 with r _<0.01 appear, which have some influence on the formation of hyperconcentrated flows. There exist some optimal ratios of coarse to fine fractions in suspended sediment that make sediment concentrations of hyperconcentrated flows the highest. The optimal r _>0.05/r _<0.01 value is related to some range of the index P _m/D _ss. When the P _m/D _ss index falls in this range, the optimum combination of relative coarse with fine sediments in the suspended load appears, and thus results in the peak values of sediment concentration.     

不同上覆水氟浓度对湖泊沉积物氟释放与微生物群落的影响

通过模拟沙湖沉积物-水系统,以沙湖原水氟离子（F^-）浓度为1倍浓度（0.69 mg/L）,设置0.5倍浓度、1倍浓度、1倍浓度灭菌、2倍浓度和4倍浓度共5个实验组,探究不同上覆水F^-浓度背景下沙湖沉积物中F^-的迁移.结果表明,沉积物F^-的释放量随上覆水F^-浓度的增加而呈下降的趋势,其中2倍和4倍浓度组的沉积物由释放F^-转变为吸附F^-.碱性水体有利于沉积物F^-的释放,即随着pH增大,F^-释放量也会增加.微生物的Beta多样性层次聚类分析表明,F^-会抑制Paenisporosarcina与Thiobacillus的相对丰度,但对Fusibacter的生命活动具有促进作用.根据各浓度组间优势菌属相对丰度与环境因子的冗余分析可得,Acinetobacter、Thiobacillus相对丰度与pH呈负相关;Fusibacter相对丰度与F^-浓度呈显著正相关,而Thiobacillus的相对丰度则与F^-表现出弱负相关.通过对各浓度组中重要离子的分析发现,0.5倍组和4倍组中F^-的迁移受Ca²⁺、HCO₃^-和SO₄^2-浓度影响较大;相关性分析表明0.5倍浓度组的Ca²⁺、HCO₃^-浓度与F^-浓度呈极显著正相关,而4倍组中F^-浓度和Ca²⁺、HCO₃^-浓度呈负相关,SO₄^2-浓度和F^-浓度的相关性则是0.5倍组呈负相关,4倍组呈极显著正相关.本文在不同上覆水F^-浓度背景下全面探究了沙湖沉积物中F^-释放和迁移机理,为沙湖及其他含氟地表水的氟污染防治提供理论支持.     

The spatial and temporal changes inV p /V s before strong earthquakes in Kamchatka

A method was developed to obtain from a signal station the spatial and temporal distribution ofV _p /V _s ratios before earthquakes of magnitude>6. It was shown thatV _p /V _s values strongly depend upon the relative positions of the stations, the future large earthquake and the foci of the smaller earthquakes used forV _p /V _s determination. The appearance of a zone of anomalousV _p /V _s values with linear dimensions of the order of 100 km was noted at least 4 years before a deep earthquake of magnitude 7. Similar size anomalous zones were detected one year before some magnitude 6 earthquakes. V _p /V _s values decreased by a small but distinct amount during this time. Additionally, local short term minima inV _p /V _s ratios were observed some months before the major event. The epicenters of the large earthquakes were located within the 100 km size zone where the gradients of theV _p /V _s field were largest.     

Properties of electric turbulence in the polar cap ionosphere

Small-scale (scales of ∼0.5–256 km) electric fields in the polar cap ionosphere are studied on the basis of measurements of the Dynamics Explorer 2 (DE-2) low-altitude satellite with a polar orbit. Nineteen DE-2 passes through the high-latitude ionosphere from the morning side to the evening side are considered when the IMF z component was southward. A rather extensive polar cap, which could be identified using the ɛ-t spectrograms of precipitating particles with auroral energies, was formed during the analyzed events. It is shown that the logarithmic diagrams (LDs), constructed using the discrete wavelet transform of electric fields in the polar cap, are power law (μ ∼ s ^α). Here, μ is the variance of the detail coefficients of the signal discrete wavelet transform, s is the wavelet scale, and index α characterizes the LD slope. The probability density functions P(δE, s) of the electric field fluctuations δE observed on different scales s are non-Gaussian and have intensified wings. When the probability density functions are renormalized, that is constructed of δE/s ^γ, where γ is the scaling exponent, they lie near a single curve, which indicates that the studied fields are statistically self-similar. In spite of the fact that the amplitude of electric fluctuations in the polar cap is much smaller than in the auroral zone, the quantitative characteristics of field scaling in the two regions are similar. Two possible causes of the observed turbulent structure of the electric field in the polar cap are considered: (1) the structure is transferred from the solar wind, which is known to have turbulent properties, and (2) the structure is generated by convection velocity shears in the region of open magnetic field lines. The detected dependence of the characteristic distribution of turbulent electric fields over the polar cap region on IMF B _y and the correlation of the rms amplitudes of δE fluctuations with IMF B _z and the solar wind transfer function (B _y ² + B _z ²)^1/2sin(θ/2), where θ is the angle between the geomagnetic field and IMF reconnecting on the dayside magnetopause when IMF B _z < 0, together with the absence of dependence on the IMF variability are arguments for the second mechanism.     

Earthquake catalogues and modelling strategies. A new testing procedure for the comparison between competing models

In this paper we consider the procedures that, on the basis of an earthquake catalogue, yield the magnitude distribution function F _M (F _M generators). In particular, our attention is focused on the F _M generators that are currently used in the frame of the probabilistic seismic hazard analysis at a site. From an engineering point of view, the behaviour of F _M in the range of strong earthquakes is of crucial importance. On the other hand, in general, the statistical validation of F _M in that range is not feasible because of an insufficient number of strong earthquakes in available catalogues.Our investigation is limited to the problem of the comparison between competing F _M generators. We show that a simple empirical F _M generator (called the CFP generator) can be more reliable than generators based on mathematical models. The numerical experiments that support this statement regards two test-sites and two mathematical models of magnitude distribution, namely the truncated-exponential and the characteristic-type models.The main conclusion is the suggestion of a test (called the CFP-test) that can be summarized as follows. Let X be a site for which, with current criteria, a specific F _M ^X generator based on a mathematical model has been selected. At this point, following the comparison approach described in the paper, the F _M ^X generator is compared with the empirical CFP generator. It can happen (actually it happens in many of our numerical experiments) that this test indicates the CFP generator as more reliable than the selected F _M ^X generator.     

Two notes on the modern theory of inverse problems in geophysics

On the basis of comparison of the approaches to the solution of inverse problems in information theory and geophysics, it is shown that results, obtained in information theory, are suitable to supplement the theory of geophysical inverse problems. The conditions of the existence and uniqueness of the solutions of inverse problems in their practical discrete statement are specified. The terms of ɛ-entropy H _ɛ and informational capacity C _ɛ, characterizing “volumes” of unknown and observed data, are introduced. It is shown, that the instability of the solution of the inverse problem decreases with increase in H _ɛ, (increase in the “complexity” of studied section), if the relation H _ɛ ≤ C _ɛ is maintained.     

Pelagic production and respiration in the Gulf of Papua during May 2004

The metabolic balance between production and respiration in plankton communities of the Gulf of Papua was investigated in May 2004. Water samples taken at 19 stations were allocated to groups on the basis of physico-chemical characteristics. Oxygen consumption and production in flasks incubated in the dark and in the light was determined by micro-Winkler titration. Dark bottle respiration in samples influenced by the estuarine plume averaged 3.09±1.92 (SD) mmol O₂ m⁻³ d⁻¹ and production within surface light bottles averaged 7.63±3.36 (SD)  mmol O₂ m⁻³ d⁻¹. Corresponding values in stations more typical of the central Gulf of Papua were 1.68±1.30 (SD) mmol O₂ m⁻³ d⁻¹ and 1.08±2.25 (SD) mmol O₂ m⁻³ d⁻¹. Despite a shallow (<10 m) euphotic zone within the plume stations, phytoplankton production in the surface layers was sufficiently high to subsidise total water column respiration. Integrating production and respiration over the water column resulted in a calculation of net community production (NCP) of 626±504 (SD) mg C m⁻² d⁻¹, and community respiration (CR) of 712±492 mg C m⁻² d⁻¹ at the plume stations, with an average P:R ratio of 1.97. In the offshore group NCP was 157±450 (SD) mg C m⁻² d⁻¹ and CR was 1620±1576 mg C m⁻² d⁻¹. The average P:R ratio was 1.27. Three of the 7 stations allocated to the offshore group were net heterotrophic. In contrast to earlier studies in the area indicating that the Gulf of Papua waters is heterotrophic [Robertson, A.I., Dixon, P., Alongi, D.M., 1998. The influence of fluvial discharge on pelagic production in the Gulf of Papua, Northern Coral Sea. Estuarine, Coastal and Shelf Science 46, 319–331], our data indicate that in May 2004 the Gulf was in positive metabolic balance, but by only ∼120 mg C m⁻² d⁻¹. We conclude that waters of the Gulf of Papua under riverine influence are net autotrophic, but that within the central Gulf there is a fine metabolic balance alternating between autotrophy and heterotrophy.     

Variations in the amplitude of the chandler wobble

It is shown that within the framework of the Kolmogorov model the “energy” of the pole E(t) = x ₁² + x ₂² can be interpreted as a Markovian process. The exact analytical expression has been obtained for the density of the conditional probability of the quantity E(t) and the problem of the first passage time of the process E(t) has been analyzed. It was shown that the available data on the swing of the function E(t) are not at variance with the Kolmogorov model and a short-period drop of the amplitude of the Chandler wobble in the early 20th century fits this model at Q = 50–200 too; values of Q > 350 are less reasonable.     

匹配滤波方法研究2016年意大利中部MW6.0地震的余震迁移

2016年8月24日意大利中部发生M_W6.0地震,2个月之后,震中附近相继发生M_W5.5、M_W5.9、M_W6.5地震。研究这几次较大地震间的相互触发作用及机制十分必要,然而在大地震之后传统地震目录通常缺失很多余震事件,缺失的余震事件包含着早期余震时空分布和迁移规律的信息,为完善余震目录本文利用匹配滤波方法对M_W6.0地震后80天内的连续数据进行余震检测,得到了数十倍于模板数量的新检测事件,检测事件与模板事件组成的新余震目录完备震级为1.0,提高了地震目录的完备性。依据新余震目录进行余震时空分布研究,结果显示M_W6.0、M_W5.9、M_W6.5地震的早期余震迁移规律不同。M_W6.0地震的早期余震沿着断裂走向同时朝两侧迁移;而M_W5.9、M_W6.5地震的早期余震向南、北迁移表示出不对称的特征。通过拟合余震迁移前端发现,M_W6.0、M_W5.5地震的早期余震朝着随后较大地震的方向迁移,且较符合lgt的特征,表明余震迁移可能与慢滑动有关。     

Statistical features of aftershock distribution size for moderate and large earthquakes in Chinese mainland

Based on data of earthquake sequences with MS≥5.0 in Chinese mainland from 1970 to 2004,for different se-quence types and different rupture modes of the main shock,the relationship between aftershock distribution size Rand the magnitude of the main shock M0 has been studied statistically.Considering the rupture mode of the mainshock,we give the quantitative statistical relationships between R and M0 under 95%confidence level for differentsequence types.Qualitatively,lgR,the logarithm of the aftershock distribution size,is positively correlative to theM0,but the data distribution is dispersed.Viewing from different sequence types,the correlation between R and M0is very weak for isolated earthquake type(IET)sequence,R distributes in the range from 5 to 60 km;For main-shock-aftershock type(MAT),lgR is positively correlative to M0;For multiple main shock type(MMT),the core-lation between lgR and M0 is not very obvious when M0≤6.2 and R distributes in the range from 5 to 70 km,whileit shows a linear correlation when M0≥6.3.The statistical results also show that the occupational ratios of differentsequence types for strike-slip and oblique slip are almost the same.But for dip-slip(mostly are thrust mechanisms),the ratio of MAT is higher than that of IET and MMT.Comparing with previous results,it indicates that,when M0is large enough,R is mainly determined by M0 and there is almost no relationship with the rupture mode of themain shock.     

Detecting remotely triggered temporal changes around the Parkfield section of the San Andreas fault

Detecting temporal changes in fault zone properties at seismogenic depth have been a long-sought goal in the seismological community for many decades.Recent studies based on waveform analysis of repeating earthquakes have found clear temporal changes in the shallow crust and around active fault zones associated with the occurrences of large nearby and teleseismic earthquakes.However,repeating earthquakes only occur in certain locations and their occurrence times cannot be controlled,which may result in inadequate sampling of the interested regions or time periods.Recent developments in passive imaging via auto-and cross-correlation of ambient seismic wavefields (e.g.,seismic noise,earthquake coda waves) provide an ideal source for continuous monitoring of temporal changes around active fault zones.Here we conduct a systematic search of temporal changes along the Parkfield section of the San Andreas fault by cross-correlating relatively high-frequency (0.4-1.3 Hz) ambient noise signals recorded by 10 borehole stations in the High Resolution Seismic Network.After using stretch/compressed method to measure the delay time and the decorrelation-index between the daily noise cross-correlation functions (NCCFs),we find clear temporal changes in the median seismic velocity and decorrelation-index associated with the 2004 M6.0 Parkfield earthquake.We also apply the same procedure to the seismic data around five regional/teleseismic events that have triggered non-volcanic tremor in the same region,but failed to find any clear temporal changes in the daily NCCFs.The fact that our current technique can detect temporal changes from the nearby but not regional and teleseismic events,suggests that temporal changes associated with distance sources are very subtle or localized so that they could not be detected within the resolution of the current technique (～0.2%).     

Variability of parameters of the <Emphasis Type="Italic">F</Emphasis>2-layer maximum in the quiet midlatitude ionosphere under low solar activity: 1. Statistical properties

Results of statistical analysis of the properties of variability of F2-layer maximum parameters (critical frequency foF2 and the height hmF2) in quiet midlatitude ionosphere under low solar activity in the daytime (1000–1500 LT) and nighttime (2200–0300 LT) hours are presented on the basis of Irkutsk station data for 2007–2008. It is found that the distribution density of δfoF2 could be presented as consisting of two distinctly different normal laws of this distribution, one of which corresponds to weak (|δfoF2| < 10%) fluctuations in foF2 and the other corresponds to strong (30% > |δfoF2| > 10%) fluctuations. Weak fluctuations in foF2 to a substantial degree are related to ionospheric variability at times less of than 1–3 h and determine the δfoF2 variability in the daytime hours. Strong fluctuations in foF2 are mainly related to day-to-day variability of the ionosphere at a fixed local time, the variability increasing by approximately a factor of 3 during the transition from day to night and determining the δfoF2 variability in the nighttime hours. The distribution density of ΔhmF2 is close to the normal distribution law. An interpretation of the different character of the distribution densities of δfoF2 and ΔhmF2 is given.     

Application of <Superscript>137</Superscript>Cs fingerprinting technique to interpreting sediment production records from reservoir deposits in a small catchment of the Hilly Loess Plateau,China

According to variations of 137Cs and clay contents, 44 flood couplets were identified in a profile of reservoir deposit with a vertical length of 28.12 m in the Yuntaishan Gully. Couplet 27 at the middle of the profile had the highest average 137Cs content of 12.65 Bq kg-1, which indicated the 1963s' deposits, then 137Cs content decreased both downward and upward in the profile. The second top and bottom couplets had average 137Cs contents of 2.15 Bq kg-1 and 0.92 Bq kg-1, respectively. By integrated analysis of reservoir construction and management history, variations of 137Cs contents over the profile, sediment yields of flood couplets and rainfall data during the period of 1958-1970, individual storms related to the flood couplets were identified. 44 floods with a total sediment yield of 2.36×104 m3 occurred and flood events in a year varied between 1 and 10 times during the period of 1960-1970. 7-10 flood events occurred during the wet period of 1961-1964 with very wet autumn, while only 1-2 events during the dry period of 1965-1969. Average annual specific sediment yield was 1.29×104 t km-2 a-1 for the Yuntaishan Gully during the period of 1960-1970, which was slightly higher than 1.11 ×104 t km-2 a-1 for the Upper Yanhe River Basin above the Ganguyi Hydrological Station and slightly lower than 1.40 ×104 t km-2 a-1 for the nearby Zhifang Gully during the same period. Annual specific sediment yields for the Yuntaishan Gully were correlated to the wet season's rainfalls well.     

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.     

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.     

Effective radius of heating of the lower ionosphere by intense shortwave radiation

The results of a theoretical analysis of the radial distribution of electron temperature T _e in the area of heating of the lower ionosphere by intense shortwave radiation are presented. It was established that effective radius r _eff of heating at a certain height may differ significantly from the characteristic radius of illumination of the ionosphere (a) by radiation at this height. At the boundary of the heating area (r = r _eff ), the characteristic radial scale of T _e changes is less than the corresponding scale of changes in the squared amplitude of the radiation electric field, and it is almost independent of the amplitude value; i.e., the formation of a relatively strong T _e gradient at such a boundary is a common feature of heating of the lower ionosphere by intense shortwave radiation.     

Magnetization of immobilized particle dispersions with two distinct particle sizes

The magnetic properties of specimens containing dispersed magnetic particles, at low concentrations, containing two distinct size ranges, have been measured. The results have been used to evaluate the effectiveness of several parameters which have been used to discriminate the type of particles in a dispersion. The location of the point representing the properties of a specimen on a graph of $\text{M}{}_{\mathrm{<mtext>RS</mtext>}}\text{M}{}_{\mathrm{<mtext>S</mtext>}}$ and $\text{H}{}_{\mathrm{<mtext>R</mtext>}}\text{H}{}_{\mathrm{<mtext>c</mtext>}}$ indicates a possible particle dispersion which would produce a model with equivalent magnetic properties. The behaviour of dispersions of titomagnetic is also considered and these produce ambiguities in the interpretation of the graph. The possibility of modelling a dispersion which has equivalent properties to specimens from a single lava is considered.     

Sulla turbolenza delle correnti marine nello stretto di messina

Riassunto Gli Autori espongono i risultati di un tentativo di valutazione del grado di turbolenza verticale delle correnti marine davanti alle imboccature Nord e Sud dello Stretto di Messina. I valori massimi ottenuti per il coefficiente di diffusione verticale della salinitàK _z sono dell'ordine di 10³–10⁴ cm²/sec.

---

Summary The Authors present the results of some calculations carried out to obtain values of the vertical turbulence (by means of the coefficient of eddy-diffusion of salinityK _z ) of the currents in the Strait of Messina. The maximum values ofK _z are of the order of 10³–10⁴ cm²/sec.

---

---

Comunicazione presentata il 24 Aprile 1957 alla Quinta Assemblea Generale della «Società Italiana di Geofisica e Meteorologia» (Genova: 23–25 Aprile 1957).