Prediction of Great Japan Earthquake 11 March of 2011 by analysis of microseismic noise. Does Japan approach the next Mega-EQ?
|
|
- Hugh Reed
- 6 years ago
- Views:
Transcription
1 Prediction of Great Japan Earthquake 11 March of 2011 by analysis of microseismic noise. Does Japan approach the next Mega-EQ? Alexey Lyubushin Institute of Physics of the Earth Russian Academy of Sciences Moscow Ver. 05 May 2018
2 When and Were the prediction of Tohoku EQ was published before the earthquake Lyubushin A.A. Multi-fractal Properties of Low-Frequency Microseismic Noise in Japan, Book of abstracts of 7th General Assembly of the Asian Seismological Commission and Japan Seismological Society, 2008 Fall meeting, Tsukuba, Japan, November 2008, p.92. Lyubushin A.A. Synchronization Trends and Rhythms of Multifractal Parameters of the Field of Low-Frequency Microseisms Izvestiya, Physics of the Solid Earth, 2009, Vol. 45, No. 5, pp Lyubushin A.A. The Statistics of the Time Segments of Low-Frequency Microseisms: Trends and Synchronization Izvestiya, Physics of the Solid Earth, 2010, Vol., No. 6, pp Lyubushin A.A. Synchronization of multi-fractal parameters of regional and global low-frequency microseisms European Geosciences Union General Assembly 2010, Vienna, of May, 2010, Geophysical Research Abstracts, Vol. 12, EGU , 2010, Lyubushin A., Multifractal Parameters of Low-Frequency Microseisms // V. de Rubeis et al. (eds.), Synchronization and Triggering: from Fracture to Earthquake Processes, GeoPlanet: Earth and Planetary Sciences 1, DOI / _15, Springer-Verlag Berlin Heidelberg, 2010, 8p., Chapter 15, pp Lyubushin A.A. Synchronization phenomena of low-frequency microseisms. European Seismological Commission, nd General Assembly, September 06-10, 2010, Montpelier, France. Book of abstracts, p.124, session ES6. Lyubushin A.A. Cluster Analysis of Low-Frequency Microseismic Noise Izvestiya, Physics of the Solid Earth, 2011, Vol. 47, No. 6, pp (submitted April 26, 2010) Publications after the Tohoku earthquake and about increasing seismic danger at the Nankai Trough. Lyubushin A.A. (2011) Seismic Catastrophe in Japan on March 11, 2011: Long-Term Prediction on the Basis of Low-Frequency Microseisms Izvestiya, Atmospheric and Oceanic Physics, 2011, Vol., No. 8, pp Lyubushin, A. (2012) Prognostic properties of low-frequency seismic noise. Natural Science, Vol.4,No.8A, doi: 10./ns Lyubushin A.A. (2013) Mapping the Properties of Low-Frequency Microseisms for Seismic Hazard Assessment. Izvestiya, Physics of the Solid Earth, 2013, Vol.49, No.1, pp ISSN , DOI: 10.11/S Lyubushin A.A. (2013) Spots of Seismic Danger Extracted by Properties of Low-Frequency Seismic Noise European Geosciences Union General Assembly 2013, Vienna, of April, 2013, Geophysical Research Abstracts, Vol. 15, EGU , Lyubushin, A. (2013) How soon would the next mega-earthquake occur in Japan? Natural Science, Vol.5, No.8A1, 1-7. doi: 10./ns A Lyubushin A.A. (2014) Dynamic estimate of seismic danger based on multifractal properties of low-frequency seismic noise. Natural Hazards, January 2014, Volume 70, Issue 1, pp DOI /s
3 Asian Seismological Commission, Tsukuba, Japan, Nov 2008 The earliest prediction (at the middle of 2008): the magnitude of oncoming catastrophe will be more than 8.3.
4
5 The latest prediction: dangerous time for waiting catastrophe begins at the middle of 2010
6 Positions of 78 seismic stations of broadband network F-net in Japan. N, deg , M = , M = E, deg 148
7 Main sources of lowfrequency microseisms Atmospheric Pressure Variations Oceanic Waves Low-Frequency Microseismic Noise: Slow Earthquakes, Creep 2 min Periods 1000 min
8 Multi-fractal singularity spectrum Measure of the random signal variability at the time interval [ t δ/2, t + δ/2 ] Multi-fractal singularity spectrum F (α) and its parameters: α - support width and α - generalized Hurst exponent. 1.0 F (α) - fractal dimensionality of the set of time moments t for which h(t) = α µ(t,δ ) δ h(t) δ α = α max α min t δ α min α α max α
9 Minimum normalized entropy of squared orthogonal wavelet coefficients distribution Minimum normalized entropy : N En = p log( p ) / log( N) min k = 1 k N En 1, pk = ck / c j, j= 1 c j k - orthogonal wavelet coefficients, minimum is taken by wavelets from Daubechies family.
10 Data transform Initial LHZ-records, t = 1 sec t = 1 min Seismic Noise, t = 1 min Downsampling Detrending Donoho share Log(Variance) Kurtosis Estimates within adjacent daily time intervals, t = 1 day α Linear Predictability Index Spectral Exponent α α min Smoothness Index Normalized Entropy
11 Graphics of seismic noise of the length 1 day with sampling time step 1 minute after removing tidal trends High α Low normalized entropy En Low danger, high variability of stochastic behavior types, a lot of low-frequency spikes because of mutual movements of non-consolidated small blocks of the Earth's crust, energy is not accumulated. Low α High normalized entropy En High danger, the behavior of the seismic noise is much more uniform, small blocks of the Earth's crust are consolidated, the energy could be accumulated. α = 0.249, En = α = 0.821, En = α = 0.154, En = α = 0.757, En = α = 0.207, En = α = 0.809, En = Time, minutes Time, minutes
12 Maps of multi-fractal singularity spectra support width α. Low α values ( blue and violet regions ) indicate synchronization and danger α α N, deg N, deg E, deg E, deg From the beginning of 1997 up to 25 of Sept 2003: the area of future seismic catastrophe is characterized by relatively low α and it is not split into North and South parts. From 26 of Sept 2003 up to 10 of March 2011: the area of future seismic catastrophe is characterized by relatively low α and the previous area of low α values is split into North and South parts. From March 14, 2011 by April, 2018: the North part of the relatively low α values before was realized as the area of Great Japan Earthquake 11 of March 2011, M=9.0, whereas the South part is still characterized by low α values.
13 α The sequence of mean maps of multi-fractal singularity spectra support width of lowfrequency seismic noise waveforms for 16 years of observations at broad-band seismic network F-net at Japan Small black star hypocenter of earthquake , M=8.3 Large black star hypocenter of Tohoku megaearthquake , M=
14
15 Monthly maps of singularity spectrum support width α, corresponding to 2017
16 Monthly maps of singularity spectrum support width α, corresponding to 2018
17 Maps of low-frequency seismic noise wavelet-based normalized entropy, i.e. normalized entropy of the noise waveforms squared wavelet coefficients for the best orthogonal wavelet which is found for each station within each daily time window from the minimum of entropy En En N, deg N, deg E, deg E, deg From the beginning of 1997 up to 25 of Sept 2003: the area of future seismic catastrophe is characterized by relatively high values of normalized entropy and it is not split into North and South parts. From 26 of Sept 2003 up to 10 of March 2011: the area of future seismic catastrophe is characterized by relatively high values of normalized entropy and the previous area of high entropy values is split into North and South parts. From March 14, 2011 by April, 2018: the North part of the relatively high values of normalized entropy before was realized as the area of Great Japan Earthquake 11 of March 2011, M=9.0, whereas the South part is still characterized by high entropy values.
18 En The sequence of mean maps of minimum normalized entropy of squared orthogonal wavelet coefficients of low-frequency seismic noise waveforms for 16 years of observations at broad-band seismic network F-net at Japan. Small black star hypocenter of earthquake , M=8.3 Large black star hypocenter of Tohoku megaearthquake , M=
19
20 Averaged monthly maps of minimum normalized entropy of squared wavelet coefficients, corresponding to 2017
21 Averaged monthly maps of minimum normalized entropy of squared wavelet coefficients, corresponding to 2018
22 Maps of low-frequency seismic noise waveforms kurtosis κ κ N, deg 45 N, deg E, deg E, deg From the beginning of 1997 up to 25 of Sept From 26 of Sept 2003 up to 10 of March From March 14, 2011 by April, Kurtosis = ( x) ( x) , x - increments after removing tidal trends within each daily window
23 κ The sequence of mean maps of seismic noise waveforms kurtosis for 16 years of observations at broad-band seismic network F-net at Japan. Small black star hypocenter of earthquake , M=8.3 Large black star hypocenter of Tohoku megaearthquake , M=
24
25 Averaged monthly maps of seismic noise waveforms kurtosis, corresponding to 2017
26 Averaged monthly maps of seismic noise waveforms kurtosis, corresponding to 2018
27
28
29 Averaged maps of annual correlation coefficient between multi-fractal singularity spectrum support width α and generalized Hurst exponent α* within moving time window of the length 5 days with mutual shift 7 days N, deg E, deg Before Tohoku EQ the region of future event was extracted by relatively high values of correlation coefficient What is observed for time period March 14, 2011 April, 2018
30 Averaged maps of annual correlation coefficient between kurtosis and generalized Hurst exponent within moving time window of the length 5 days with mutual shift 7 days N, deg N, deg E, deg E, deg Before Tohoku EQ the region of future event was extracted by relatively high values of correlation coefficient. What is observed for time period March 14, 2011 April, 2018
31 Averaged maps of annual correlation coefficient between minimum wavelet-based normalized entropy and generalized Hurst exponent within moving time window of the length 5 days with mutual shift 7 days Before Tohoku EQ the region of future event was extracted by relatively low values of correlation coefficient (high absolute values) What is observed for time period March 14, 2011 April, 2018
32 Cluster analysis of daily median values of 3 multifractal parameters within moving time window ( t ) * ( t) Let s consider moving time window of the length L =5 days and let ξ = ( α, α, α ) be 3D vector within current time window, t = 1,..., L, t is time ( t ) index, numerating vectors. Our purpose is investigating clustering properties of clouds of 3D vectors ξ with each 1-year time window. In particular, we are interesting what is the best number of clusters. Before making cluster procedure a preliminary operation of normalizing and winzorizing was performed within each time window for each scalar component ( t ) ( ) ξ k of vectors ξ t ( ). Here k = 1,2,3 is the index numerating scalar components of the vector ξ t 1 L ( t ) 2 1 L. Let ξk = L ξ, ( ( t ) 2 t= 1 k σ ) k = ( L 1) ξ t 1 k ξ be sample = ( t ) estimates of mean values and variance of scalar components of the 3D vector ξ. Let s perform iterations which consist in coming to values ( t ) ( t ) ( ) ζ k = ( ξ k ξk ) / σ k and clipping values ζ t k exceeding over and under thresholds ± 4σ k. These iterations are stopped when the values ξ k and σ became stable k and equal to the following values: ξ = k 0, σ 1. After this preliminary operation at each current time window we have a cloud consisting of L 3D vectors ( t) ζ. k = Let s split some cloud into given probe number q of clusters using standard k-means cluster procedure. Let Γ r, r = 1,..., q be clusters, zr ζ / n ζ Γ r mean r ( ) vector of the cluster Γ r, n r be a number of vectors ζ t q within cluster Γ r, n = L. K-means procedure minimizes sum with respect to positions of clusters centers z r. Let the best number of clusters r= 1 r S z z z z1,..., zq min q 2 ( 1,..., q) = ζ r min z1,..., zq r= 1 ζ Γ r J ( q) = min S( z,..., z ) * q was solved from maximum of pseudo-f-statistics: 1 q =. We try probe number of clusters within range 2 q. The problem of selecting PFS( q) = σ ( q) σ ( q) max, σ = J ( q) /( L q) q 0 σ ( q) ν z z, ν n / L, z ζ / L q L 2 2 ( t) 1 = r r 0 r = r 0 = r= 1 t= 1
33 Pseudo-F-statistics map as an estimate of current seismic danger (a1)-(a3) plots of daily median values of multifractal singularity spectrum support width α, * generalized Hurst exponent α and minimum Hölder-Lipschitz exponent α min from all 78 stations of broadband seismic network F-net in Japan, green lines present running average within time windows of the length 57 days; (b) plot of the best numbers of clusters for the sequence of clouds consisting of 5 daily 3D * vectors ( α, α, α min ) from moving time window of the length 5 days with mutual shift 3 days. The best number of clusters is defined from the maximum of pseudo-f-statistics. Vertical red line indicates time moment of Tohoku mega-earthquake on March 11, 2011, M = 9.1. Two-dimensional diagram (c) presents dependence of pseudo- F-statistics on the probe number of clusters which is varying from 2 up to within each time window. Plot (d) presents mean value of pseudo-f-statistics averaged by all probe numbers of clusters in dependence on right-hand end of moving time window of the length 5 days. Plot (e) presents the sequence of time moments of strong earthquakes M 7 in the rectangular domain with coordinates 28 ο N Latitude 48 ο N; 128 ο E Longitude 156 ο E which is a rather broad vicinity of Japan islands.
34
35 Time Prediction: Squared correlation between mean values of multi-fractal parameters and of microseisms from all F-net stations estimated within 1 year moving time window. 0.8 Strong earthquake is only a foreshock , M = , M = Mean correlation since Best fitted line to the points of minimum correlations 0.4 Blue arrows were plotted at the April of The 2-nd arrow indicates the middle of 2010 as the lower estimate of catastrophe time moment Right-hand end of moving time window of the length 1 year
36
37 Dangerous spots of high coherence of GPS noise at Nankai Trough Alexey Lyubushin, Institute of Physics of the Earth, Russian Academy of Sciences, Moscow, Network of 13 GPS stations in Japan Map of kernel estimate of probability density functions is presented for nodes of the regular net sized x, which are realized the spatial maximum of the frequency-dependent maximal values of the multiple spectral function of coherence of GPS time series from the nearest 10 workable stations for all 3 components of GPS time series with sampling time step 5 minutes within moving time window of the length 5 days. Time interval Maximum danger at the vicinity of the point N and 1E This result is an independent confirmation of the earlier conclusion made by the analysis of seismic noise Data were downloaded from the site of Nevada Geodetic Laboratory: ftp://gneiss.nbmg.unr.edu/rapids_5min/kenv/ E-005 1E-006 1E-007 1E-008 1E-009 1E-010 1E-011 1E-012 1E Periods, minutes Examples of multiple spectral coherence functions in dependence on periods from 10 nearest workable GPS station within time window of the length 5 days. Multiple coherence function: is based on using canonical coherences: The spectral matrices were calculated in a moving time window of the length 14 samples (5 days) with mutual shift 288 samples (1 day) for each node of regular grid from 10 nearest workable stations using model of vector autoregression of 5 th order. The GPS station is considered workable in the current time window if its registration interval includes the considered time window and the number of missing values does not exceed a predetermined maximum allowable proportion of the total length equal to 0.1. The missed values are filled using information about records from neighbor time interval of the same length as the length of gaps. Before calculating the spectral matrix in each time window the trend is removed by polynomial of 4 th order and ±3σ winsorizing was performed.
38 Conclusions. During preparing huge seismic catastrophe geological medium become more simple and homogenous. This manifests in the following properties of lowfrequency microseisms noise waveforms: 1) Multi-fractal singularity spectrum support width α decreases, i.e. the noise structure become more simple, less multi-fractal. 2) Minimum value of normalized entropy of squared orthogonal wavelet coefficients increases. 3) According to the results of seismic noise processing from Japan broadband network F-net the region of Nankai Trough could be the place of the next mega-earthquake. 4) In Japan, since 2003, the regime of natural fluctuations of seismic hazard has been established with a period of about 2 years. The previous 2-year cycle preceded the last strong earthquake Kumamoto on in the south of Japan. It follows that the next 2018 will be a year of increased seismic danger for Japan.
Spots of Seismic Danger Extracted by Properties of Low-Frequency Seismic Noise
Spots of Seismic Danger Extracted by Properties of Low-Frequency Seismic Noise Alexey Lyubushin Institute of Physics of the Earth, Moscow European Geosciences Union General Assembly, Vienna, 07 12 April
More informationSeismic Catastrophe in Japan on March 11, 2011: Long-Term Prediction on the Basis of Low-Frequency Microseisms
ISSN -8, Izvestiya, Atmospheric and Oceanic Physics,, Vol. 6, No. 8, pp. 9 9. Pleiades Publishing, Ltd.,. Original Russian Text A.A. Lyubushin,, published in Geofizicheskie protsessy i biosfera,, Vol.,
More informationPreliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code.
Preliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code. J. Zahradnik 1), F. Gallovic 1), E. Sokos 2) G-A. Tselentis 2) 1) Charles University
More informationSCIENCE OF TSUNAMI HAZARDS
SCIENCE OF TSUNAMI HAZARDS ISSN 8755-6839 Journal of Tsunami Society International Volume 32 Number 4 2013 A METHOD FOR THE ESTIMATION OF TSUNAMI RISK ALONG RUSSIA s FAR EAST G.V. Shevchenko, D. E. Zolotukhin,
More informationSimulation of Strong Ground Motions for a Shallow Crustal Earthquake in Japan Based on the Pseudo Point-Source Model
6 th International Conference on Earthquake Geotechnical Engineering -4 November 25 Christchurch, New Zealand Simulation of Strong Ground Motions for a Shallow Crustal Earthquake in Japan Based on the
More informationGemona valley (Italy)
Gemona valley (Italy) Situation of the valley Gemona is a city located in the north-eastern part of Italy. It is mainly built on an alluvial fan. The area were Gemona is located is a seismologically active
More informationInternational Centre for Global Earth Models (ICGEM)
International Centre for Global Earth Models (ICGEM) 1 International Centre for Global Earth Models (ICGEM) http://icgem.gfz-potsdam.de/ Franz Barthelmes, Elmas Sinem Ince, Sven Reißland Overview The ICGEM
More informationLong-term properties of annual maximum daily rainfall worldwide
European Geosciences Union General Assembly 2011 Vienna, Austria, 03 08 April 2011 Session: HS7.4/AS4.9/CL3.4 Hydrological change versus climate change Long-term properties of annual maximum daily rainfall
More informationSubject Index. Block maxima, 3 Bootstrap, 45, 67, 149, 176 Box-Cox transformation, 71, 85 Brownian noise, 219
Subject Index Entries in this index are generally sorted with page number as they appear in the text. Page numbers that are marked in bold face indicate that the entry appears in a title or subheading.
More informationJCR (2 ), JGR- (1 ) (4 ) 11, EPSL GRL BSSA
Dun Wang ( ) In collaboration with: Hitoshi Kawakatsu, Jim Mori, Kazuki Koketsu, Takuto Maeda, Hiroshi Tsuroka, Jiancang Zhunag, Lihua Fang, and Qiang Yao School of Geosciences, China University of Geosciences
More informationPROBABILISTIC SEISMIC HAZARD MAPS AT GROUND SURFACE IN JAPAN BASED ON SITE EFFECTS ESTIMATED FROM OBSERVED STRONG-MOTION RECORDS
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 3488 PROBABILISTIC SEISMIC HAZARD MAPS AT GROUND SURFACE IN JAPAN BASED ON SITE EFFECTS ESTIMATED FROM
More informationSOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION
Synopses of Master Papers Bulletin of IISEE, 47, 19-24, 2013 SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION Fabricio Moquete Everth* Supervisor:
More informationSub-Gaussian Model Based LDPC Decoder for SαS Noise Channels
Sub-Gaussian Model Based LDPC Decoder for SαS Noise Channels Iulian Topor Acoustic Research Laboratory, Tropical Marine Science Institute, National University of Singapore, Singapore 119227. iulian@arl.nus.edu.sg
More informationSource modeling of hypothetical Tokai-Tonankai-Nankai, Japan, earthquake and strong ground motion simulation using the empirical Green s functions
Source modeling of hypothetical Tokai-Tonankai-Nankai, Japan, earthquake and strong ground motion simulation using the empirical Green s functions Y. Ishii & K. Dan Ohsaki Research Institute, Inc., Tokyo
More informationREPORT ON THE TOHOKU AREA PASIFIC OFFSHORE EARTHQUAKE
REPORT ON THE TOHOKU AREA PASIFIC OFFSHORE EARTHQUAKE GENERAL PERSPECTIVE The Highest Magnitude Ever Recorded The 2011 off the Pacific Coast of Tohoku Earthquake (hereafter, the 2011 Tohoku- Pacific Earthquake
More informationLong-period Ground Motion Characteristics of the Osaka Sedimentary Basin during the 2011 Great Tohoku Earthquake
Long-period Ground Motion Characteristics of the Osaka Sedimentary Basin during the 2011 Great Tohoku Earthquake K. Sato, K. Asano & T. Iwata Disaster Prevention Research Institute, Kyoto University, Japan
More informationThe Illustrated Wavelet Transform Handbook. Introductory Theory and Applications in Science, Engineering, Medicine and Finance.
The Illustrated Wavelet Transform Handbook Introductory Theory and Applications in Science, Engineering, Medicine and Finance Paul S Addison Napier University, Edinburgh, UK IoP Institute of Physics Publishing
More informationSeparation of a Signal of Interest from a Seasonal Effect in Geophysical Data: I. El Niño/La Niña Phenomenon
International Journal of Geosciences, 2011, 2, **-** Published Online November 2011 (http://www.scirp.org/journal/ijg) Separation of a Signal of Interest from a Seasonal Effect in Geophysical Data: I.
More informationIf you are looking for the book Geodetic Deformation Monitoring: From Geophysical to Engineering Roles: IAG Symposium Jaén, Spain, March 7-19,2005
Geodetic Deformation Monitoring: From Geophysical To Engineering Roles: IAG Symposium Jaén, Spain, March 7-19,2005 (International Association Of Geodesy Symposia) READ ONLINE If you are looking for the
More informationSHale gas. in relations to technological activity at The Geysers geothermal field, California
SHale gas Exploration and Exploitation induced Risks Temporal static stress drop variations in relations to technological activity at The Geysers geothermal field, California Monika Staszek, Beata Orlecka-Sikora,
More informationCharacterization of Induced Seismicity in a Petroleum Reservoir: A Case Study
Characterization of Induced Seismicity in a Petroleum Reservoir: A Case Study Edmond Sze, M. Nafi Toksöz, and Daniel R. Burns Earth Resources Laboratory Dept. of Earth, Atmospheric and Planetary Sciences
More informationThree Fs of earthquakes: forces, faults, and friction. Slow accumulation and rapid release of elastic energy.
Earthquake Machine Stick-slip: Elastic Rebound Theory Jerky motions on faults produce EQs Three Fs of earthquakes: forces, faults, and friction. Slow accumulation and rapid release of elastic energy. Three
More informationCoseismic slip model
Figure 3 - Preliminary highly smoothed model of coseismic slip for the 11 March UCL Institute for Risk & Disaster Reduction Magnitude 9.0 (JMA scale) earthquake Occurred at 02:46:23 pm local time near
More informationAVERAGE AND VARIATION OF FOCAL MECHANISM AROUND TOHOKU SUBDUCTION ZONE
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 414 AVERAGE AND VARIATION OF FOCAL MECHANISM AROUND TOHOKU SUBDUCTION ZONE Shunroku YAMAMOTO 1 Naohito
More informationGEO-VIII November Geohazard Supersites and Natural Laboratories Progress Report. Document 9
GEO-VIII 16-17 November 2011 Geohazard Supersites and Natural Laboratories Progress Report Document 9 This document is submitted to GEO-VIII for information. Geohazard Supersites and Natural Laboratories
More informationby Pentti Varpasuo* Fortum Nuclear Services Ltd Rajatorpantie 8, Vantaa, Fortum
THE INFORMATION REPORT OF EARTHQUAKES IN NORTHWEST OF RUSSIA CLOSE TO THE CITY OF KALININGRAD ON SEPTEMBER 21, 2004 AND THE ASSESSMENT OF THEIR EFFECTS AT OL3 NPP SITE by Pentti Varpasuo* Fortum Nuclear
More informationA pilot project on the comphrensive diagnosis of earthquake precursors on Sakhalin Island: Experiment results from 2007
A pilot project on the comphrensive diagnosis of earthquake precursors on Sakhalin Island: Experiment results from 2007 Alexey Romanov, Youry Urlichich, Sergey Pulinets, Alexander Romanov and Victor Selin
More informationTECTONIC MOVEMENTS MONITORING OF SUDETIC MARGINAL FAULT USING SHORT GPS BASELINES
POSTER PAPER SESSION TECTONIC MOVEMENTS MONITORING OF SUDETIC MARGINAL FAULT USING SHORT GPS ASELINES ernard Kontny Department of Geodesy and Photogrammetry Agricultural University of Wroclaw, Poland Abstract
More information4-D Seismic Tomography for the Complex System of Strong Earthquakes: Formulation of a Problem
4-D Seismic Tomography for the Complex System of Strong Earthquakes: Formulation of a Problem Tatyana A. Smaglichenko and Ingi Th. Bjarnason Abstract Geodynamic processes are acting in the Earth s interior
More informationRadiative Climatology of the North Slope of Alaska and the Adjacent Arctic Ocean
Radiative Climatology of the North Slope of Alaska and the Adjacent Arctic Ocean C. Marty, R. Storvold, and X. Xiong Geophysical Institute University of Alaska Fairbanks, Alaska K. H. Stamnes Stevens Institute
More informationGROUND MOTION SPECTRAL INTENSITY PREDICTION WITH STOCHASTIC GREEN S FUNCTION METHOD FOR HYPOTHETICAL GREAT EARTHQUAKES ALONG THE NANKAI TROUGH, JAPAN
GROUND MOTION SPECTRAL INTENSITY PREDICTION WITH STOCHASTIC GREEN S FUNCTION METHOD FOR HYPOTHETICAL GREAT EARTHQUAKES ALONG THE NANKAI TROUGH, JAPAN Masayuki YOSHIMI 1, Yasuto KUWAHARA 2, Masayuki YAMADA
More informationThe Non-volcanic tremor observation in Northern Cascadia. Hsieh Hsin Sung 3/22
The Non-volcanic tremor observation in Northern Cascadia Hsieh Hsin Sung 3/22 Reference Kao, H., S. J. Shan, H. Dragert, and G. Rogers (2009), Northern Cascadia episodic tremor and slip: A decade of observations
More informationA Modified Earthquake Model Based on Generalized Barabási Albert Scale-Free
Commun. Theor. Phys. (Beijing, China) 46 (2006) pp. 1011 1016 c International Academic Publishers Vol. 46, No. 6, December 15, 2006 A Modified Earthquake Model Based on Generalized Barabási Albert Scale-Free
More informationLong term properties of monthly atmospheric pressure fields
European Geosciences Union General Assembly 2012 Vienna, Austria, 22 27 April 2012 Session HS7.4/AS4.17/CL2.10: Climate, Hydrology and Water Infrastructure Long term properties of monthly atmospheric pressure
More informationSpatial and temporal variability of wind speed and energy over Greece
European Geosciences Union General Assembly 2014 Vienna, Austria, 27 April 02 May 2014 Session HS5.4/CL4.8/ERE1.16 : Hydropower and other renewable sources of energy for a sustainable future: modeling
More informationChapter Introduction
Chapter 4 4.1. Introduction Time series analysis approach for analyzing and understanding real world problems such as climatic and financial data is quite popular in the scientific world (Addison (2002),
More informationData Sharing at the National Research Institute for Earth Science and Disaster Resilience
Data Sharing at the National Research Institute for Earth Science and Disaster Resilience Katsuhiko SHIOMI Network Center for Earthquake, Tsunami and Volcano, NIED 2017.12.05 International workshop on
More informationAN OVERVIEW AND GUIDELINES FOR PROBABILISTIC SEISMIC HAZARD MAPPING
CO 2 TRACCS INTERNATIONAL WORKSHOP Bucharest, 2 September, 2012 AN OVERVIEW AND GUIDELINES FOR PROBABILISTIC SEISMIC HAZARD MAPPING M. Semih YÜCEMEN Department of Civil Engineering and Earthquake Studies
More informationV r : A new index to represent the variation rate of geomagnetic activity
Earthq Sci (2010)23: 343 348 343 Doi: 10.1007/s11589-010-0731-9 V r : A new index to represent the variation rate of geomagnetic activity Dongmei Yang 1, Yufei He 1 Chuanhua Chen 2 and Jiadong Qian 3 1
More informationThe geomagnetic field variations and earthquake activity
The geomagnetic field variations and earthquake activity Yu. Ya. Ruzhin 1), V. A. Novikov 2) 1) Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation (IZMIRAN), Russian Academy
More informationDiscriminating Between Explosions and Earthquakes by the SDFA Method. Raquel Romes Linhares 1
American Review of Mathematics and Statistics December 2017, Vol. 5, No. 2, pp. 1-5 ISSN: 2374-2348 (Print), 2374-2356 (Online) Copyright The Author(s).All Rights Reserved. Published by American Research
More informationFinite data-size scaling of clustering in earthquake networks
Finite data-size scaling of clustering in earthquake networks Sumiyoshi Abe a,b, Denisse Pastén c and Norikazu Suzuki d a Department of Physical Engineering, Mie University, Mie 514-8507, Japan b Institut
More informationCrustal deformation by the Southeast-off Kii Peninsula Earthquake
Crustal deformation by the Southeast-off Kii Peninsula Earthquake 51 Crustal deformation by the Southeast-off Kii Peninsula Earthquake Tetsuro IMAKIIRE, Shinzaburo OZAWA, Hiroshi YARAI, Takuya NISHIMURA
More informationVertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes
Vertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes N.J. Gregor Consultant, Oakland, California, USA N.A. Abrahamson University of California, Berkeley, USA K.O. Addo BC
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationOccurrence of quasi-periodic slow-slip off the east coast of the Boso peninsula, Central Japan
LETTER Earth Planets Space, 9, 11 1, Occurrence of quasi-periodic slow-slip off the east coast of the Boso peninsula, Central Japan Shinzaburo Ozawa, Hisashi Suito, and Mikio Tobita Geographical Survey
More informationRAPID SOURCE PARAMETER DETERMINATION AND EARTHQUAKE SOURCE PROCESS IN INDONESIA REGION
RAPI SOURCE PARAETER ETERINATION AN EARTHQUAKE SOURCE PROCESS IN INONESIA REGION Iman Suardi Supervisor: Yui YAGI EE06012 Akio KATSUATA Osamu KAIGAICHI ABSTRACT We proposed a magnitude determination method
More informationCredibility of climate predictions revisited
European Geosciences Union General Assembly 29 Vienna, Austria, 19 24 April 29 Session CL54/NP4.5 Climate time series analysis: Novel tools and their application Credibility of climate predictions revisited
More informationTitle. Author(s)Fujii, Yoshiaki; Kodama, Jun-ichi; Fukuda, Daisuke. CitationProceedings ISRM Congress 2015, 2015: 513. Issue Date
Title Giant Earthquakes are Occurring at Lunar Phases Spec Author(s)Fujii, Yoshiaki; Kodama, Jun-ichi; Fukuda, Daisuke CitationProceedings ISRM Congress 2015, 2015: 513 Issue Date 2015-05- Doc URL http://hdl.handle.net/2115/5113
More informationJackson County 2013 Weather Data
Jackson County 2013 Weather Data 61 Years of Weather Data Recorded at the UF/IFAS Marianna North Florida Research and Education Center Doug Mayo Jackson County Extension Director 1952-2008 Rainfall Data
More informationDepth-dependent slip regime on the plate interface revealed from slow earthquake activities in the Nankai subduction zone
2010/10/11-14 Earthscope Workshop Depth-dependent slip regime on the plate interface revealed from slow earthquake activities in the Nankai subduction zone Kazushige Obara, ERI, Univ. Tokyo Recurrence
More informationEarthquake Source. Kazuki Koketsu. Special Session: Great East Japan (Tohoku) Earthquake. Earthquake Research Institute, University of Tokyo
2012/9/24 17:20-17:35 WCEE SS24.4 Special Session: Great East Japan (Tohoku) Earthquake Earthquake Source Kazuki Koketsu Earthquake Research Institute, University of Tokyo 1 Names and features of the earthquake
More informationSeismic Activity and Crustal Deformation after the 2011 Off the Pacific Coast of Tohoku Earthquake
J-RAPID Symposium March 6-7, 2013 Seismic Activity and Crustal Deformation after the 2011 Off the Pacific Coast of Tohoku Earthquake Y. Honkura Tokyo Institute of Technology Japan Science and Technology
More informationStatistical Properties of Marsan-Lengliné Estimates of Triggering Functions for Space-time Marked Point Processes
Statistical Properties of Marsan-Lengliné Estimates of Triggering Functions for Space-time Marked Point Processes Eric W. Fox, Ph.D. Department of Statistics UCLA June 15, 2015 Hawkes-type Point Process
More informationCOMPLEX TRACE ANALYSIS OF SEISMIC SIGNAL BY HILBERT TRANSFORM
COMPLEX TRACE ANALYSIS OF SEISMIC SIGNAL BY HILBERT TRANSFORM Sunjay, Exploration Geophysics,BHU, Varanasi-221005,INDIA Sunjay.sunjay@gmail.com ABSTRACT Non-Stationary statistical Geophysical Seismic Signal
More informationGLOBAL CHANGES OF THE SEISMICITY OF THE EARTH
GLOBAL CHANGES OF THE SEISMICITY OF THE EARTH Valentin Ulomov Schmidt Institute of Physics of the Earth, Russian Academy of Sciences, (Russia, Moscow), ulomov@ifz.ru Introduction Beginning in 2010 was
More informationSpatio-temporal variation in slip rate on the plate boundary off Sanriku, northeastern Japan, estimated from small repeating earthquakes
Spatio-temporal variation in slip rate on the plate boundary off Sanriku, northeastern Japan, estimated from small repeating earthquakes T. Matsuzawa, N. Uchida, T. Igarashi *, N. Umino, and A. Hasegawa
More informationA GLOBAL SURGE OF GREAT EARTHQUAKES FROM AND IMPLICATIONS FOR CASCADIA. Thorne Lay, University of California Santa Cruz
A GLOBAL SURGE OF GREAT EARTHQUAKES FROM 2004-2014 AND IMPLICATIONS FOR CASCADIA Thorne Lay, University of California Santa Cruz Last 10 yrs - 18 great earthquakes: rate 1.8/yr; rate over preceding century
More informationAPPENDIX G-7 METEROLOGICAL DATA
APPENDIX G-7 METEROLOGICAL DATA METEOROLOGICAL DATA FOR AIR AND NOISE SAMPLING DAYS AT MMR Monthly Normals and Extremes for Honolulu International Airport Table G7-1 MMR RAWS Station Hourly Data Tables
More informationSTRONG GROUND MOTION PREDICTION FOR HUGE SUBDUCTION EARTHQUAKES USING A CHARACTERIZED SOURCE MODEL AND SEVERAL SIMULATION TECHNIQUES
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 655 STRONG GROUND MOTION PREDICTION FOR HUGE SUBDUCTION EARTHQUAKES USING A CHARACTERIZED SOURCE MODEL
More informationAspects of risk assessment in power-law distributed natural hazards
Natural Hazards and Earth System Sciences (2004) 4: 309 313 SRef-ID: 1684-9981/nhess/2004-4-309 European Geosciences Union 2004 Natural Hazards and Earth System Sciences Aspects of risk assessment in power-law
More informationDevelopment of Ground Motion Time Histories for Seismic Design
Proceedings of the Ninth Pacific Conference on Earthquake Engineering Building an Earthquake-Resilient Society 14-16 April, 2011, Auckland, New Zealand Development of Ground Motion Time Histories for Seismic
More informationChuichi Arakawa Graduate School of Interdisciplinary Information Studies, the University of Tokyo. Chuichi Arakawa
Direct Numerical Simulations of Fundamental Turbulent Flows with the Largest Grid Numbers in the World and its Application of Modeling for Engineering Turbulent Flows Project Representative Chuichi Arakawa
More informationMultifractal Analysis of Seismicity of Kutch Region (Gujarat)
P-382 Multifractal Analysis of Seismicity of Kutch Region (Gujarat) Priyanka Midha *, IIT Roorkee Summary The geographical distribution of past earthquakes is not uniform over the globe. Also, in one region
More informationAyaka Shoda University of Tokyo. M. Ando A, K. Okada, K. Ishidoshiro B, Y. Aso, K. Tsubono Kyoto University A, KEK B
Ayaka Shoda University of Tokyo M. Ando A, K. Okada, K. Ishidoshiro B, Y. Aso, K. Tsubono Kyoto University A, KEK B Table of contents 1. Introduction 2. Torsion-bar Antenna 3. Simultaneous observational
More informationEffect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L02611, doi:10.1029/2007gl032129, 2008 Effect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami S. Koshimura, 1 Y.
More informationSchool of Computer Science and Engineeering,VIT University,Vellore. Emai: Received on: Accepted on:
ISSN: 975-766X CODEN: IJPTFI Available Online through Research Article www.ijptonline.com AN APPROACH TO PREDICT EARHQUAKE DAMAGE USING REGRESSION TECHNIQUE Ankit Panchariya, J.Jabanjalin Hilda, C.Srimathi,
More informationUpper crust structure under CTBTO station «Petropavlovsk-Kamchatsky» by endogenic microseismic activity
Upper crust structure under CTBTO station «Petropavlovsk-Kamchatsky» by endogenic microseismic activity Chebrov V.N., Kugaenko Yu.A., Saltykov V.A. Geophysical Survey of Russia, Kamchatkan Branch, Petropavlovsk-Kamchatsky
More informationWavelet Analysis of CHAMP Flux Gate Magnetometer Data
Wavelet Analysis of CHAMP Flux Gate Magnetometer Data Georgios Balasis, Stefan Maus, Hermann Lühr and Martin Rother GeoForschungsZentrum Potsdam, Section., Potsdam, Germany gbalasis@gfz-potsdam.de Summary.
More informationSeismic Source Mechanism
Seismic Source Mechanism Yuji Yagi (University of Tsukuba) Earthquake Earthquake is a term used to describe both failure process along a fault zone, and the resulting ground shaking and radiated seismic
More informationSTUDYING THE IMPORTANT PARAMETERS IN EARTHQUAKE SIMULATION BASED ON STOCHASTIC FINITE FAULT MODELING
STUDYING THE IMPORTANT PARAMETERS IN EARTHQUAKE SIMULATION BASED ON STOCHASTIC FINITE FAULT MODELING H. Moghaddam 1, N. Fanaie 2* and H. Hamzehloo 1 Professor, Dept. of civil Engineering, Sharif University
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 20, 30 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Seismic Spectra & Earthquake Scaling laws. Seismic Spectra & Earthquake Scaling laws. Aki, Scaling law
More informationArthur Frankel, William Stephenson, David Carver, Jack Odum, Robert Williams, and Susan Rhea U.S. Geological Survey
Probabilistic Seismic Hazard Maps for Seattle: 3D Sedimentary Basin Effects, Nonlinear Site Response, and Uncertainties from Random Velocity Variations Arthur Frankel, William Stephenson, David Carver,
More informationInvestigation of long period amplifications in the Greater Bangkok basin by microtremor observations
Proceedings of the Tenth Pacific Conference on Earthquake Engineering Building an Earthquake-Resilient Pacific 6-8 November 2015, Sydney, Australia Investigation of long period amplifications in the Greater
More informationSeismic Sequences Branching Structures: Long-Range Interactions and Hazard Levels
Open Journal of Earthquake Research, 201,, 189-20 http://www.scirp.org/journal/ojer ISSN Online: 219-931 ISSN Print: 219-923 Seismic Sequences Branching Structures: Long-Range Interactions and Hazard Levels
More informationFORECASTING THE M6.9 KYTHERA EARTHQUAKE OF 8 JANUARY 2006: A STEP FORWARD IN EARTHQUAKE PREDICTION?
European Geosciences Union General Assembly 2006 Vienna, Austria,, 02 07 April 2006 FORECASTING THE M6.9 KYTHERA EARTHQUAKE OF 8 JANUARY 2006: A STEP FORWARD IN EARTHQUAKE PREDICTION? Andreas Tzanis Department
More informationMagnitude 8.3 SEA OF OKHOTSK
A powerful earthquake in Russia's Far East was felt as far away as Moscow, about 7,000 kilometers (4,400 miles) west of the epicenter, but no casualties or damage were reported. The epicenter was in the
More informationIGF PAS activity in WP4
SHale gas Exploration and Exploitation induced Risks This project has received funding from the European Union s Horizon 2020 research and innovation programme under grant agreement No 640896. IGF PAS
More informationRISKY HIGH-RISE BUILDINGS RESONATING WITH THE LONG-PERIOD STRONG GROUND MOTIONS IN THE OSAKA BASIN, JAPAN
RISKY HIGH-RISE BUILDINGS RESONATING WITH THE LONG-PERIOD STRONG GROUND MOTIONS IN THE OSAKA BASIN, JAPAN K. Miyakoshi 1 and M. Horike 2 ABSTRACT : 1 Earthquake Engineering Group, Geo-Research Institute,
More informationData analysis methods in the neuroscience Zoltán Somogyvári
Data analysis methods in the neuroscience Zoltán Somogyvári Wigner Research Centre for Physics of the Hungarian Academy of Sciences Methods applicable to one time series The Fourier transformation The
More informationThe Japanese University Joint Seismic Observations at the Niigaka-Kobe Tectonic Zone
Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 2*,**/ pp. +-- +.1 * The Japanese University Joint Seismic Observations at the Niigaka-Kobe Tectonic Zone The Japanese University Group of the Joint Seismic Observations
More informationAn autocorrelation method to detect low frequency earthquakes within tremor
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L16305, doi:10.1029/2008gl034560, 2008 An autocorrelation method to detect low frequency earthquakes within tremor Justin R. Brown, 1 Gregory C. Beroza, 1 and David
More informationComputational Seismology: A Practical Introduction. Heiner Igel Computational Seismology 1 / 18
Computational Seismology: A Practical Introduction Heiner Igel Computational Seismology 1 / 18 Introduction Heiner Igel Computational Seismology 2 / 18 Goals of the course Introduction Understand methods
More informationModels for models. Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research
Models for models Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research Outline Statistical models and tools Spatial fields (Wavelets) Climate regimes (Regression and clustering)
More informationData S.A. 30s. Xe= Gravity anom. sigd=0.06 corlength=1000s sigp=1. dg microgal. dg microgal. Data S.A. 30s. Xe=
Gravity anom. sigd=.6 corlength=4s sigp=1. Gravity anom. sigd=.6 corlength=6s sigp=1. Xe=.594 Λ=4s Xe=.595 Λ=6s -2. -15. -1. -5.. Tim -2. -15. -1. -5.. Tim Gravity anom. sigd=.6 corlength=8s sigp=1. Gravity
More informationCase study of Japan: Reference Frames in Practice
Case study of Japan: Reference Frames in Practice Basara Miyahara and Koji Matsuo Sponsors: 1 Outline Introduction - Geodetic reference frame of Japan - GEONET : continuous GNSS observation system Geometric
More informationChapter 2 Analysis of Solar Radiation Time Series
Chapter 2 Analysis of Solar Radiation Time Series Abstract In this chapter, various kinds of analysis are performed by using solar radiation time series recorded at 12 stations of the NREL database. Analysis
More informationBROADBAND STRONG MOTION SIMULATION OF THE 2004 NIIGATA- KEN CHUETSU EARTHQUAKE: SOURCE AND SITE EFFECTS
Third International Symposium on the Effects of Surface Geology on Seismic Motion Grenoble, France, 30 August - 1 September 2006 Paper Number: 105 BROADBAND STRONG MOTION SIMULATION OF THE 2004 NIIGATA-
More information2. Athens Water Supply & Sewerage Company (EYDAP), - URL:
Preseismic electric field strange attractor like precursor analysis applied on large (Ms > 5.5R) EQs, which occurred in Greece during December 1 st, 2007 - April 30 th, 2008. Thanassoulas 1, C., Klentos
More informationSeismological Study of Earthquake Swarms in South-Eastern Puerto Rico
Seismological Study of Earthquake Swarms in South-Eastern Puerto Rico Principal Investigator: Lillian Soto-Cordero [1] Co-Investigators: Victor Huérfano 1, Leonardo Cano [2], Robert Watts [3], and Christa
More informationEffect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics
Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological
More informationP217 Automated P- and S-Wave Picking of Micro Earthquakes Re-Corded by a Vertical Array
P217 Automated P- and S-Wave Picking of Micro arthquakes Re-Corded by a Vertical Array T. Fischer* (Geophysical Institute, Czech Academy of Science), A. Bouskova (Geophysical Institute, Czech Academy of
More informationTHE ANTHROPOGENIC FACTOR IN THE VARIATIONS OF MAGNETIC FIELD. O.D. Zotov
THE ANTHROPOGENIC FACTOR IN THE VARIATIONS OF MAGNETIC FIELD O.D. Zotov Geophysical Observatory Borok, IPE, RAS, Borok, Russia, e-mail: ozotov@inbox.ru Abstract. This research addresses to the problem
More informationSeismic hazard modeling for Bulgaria D. Solakov, S. Simeonova
Seismic hazard modeling for Bulgaria D. Solakov, S. Simeonova Bulgarian seismic network and foreign stations used in epicenter location Sismicity in Bulgaria and surroundings (M 4.) Epicentral map for
More informationLong-period ground motion characterization by cross wavelet transform
Long-period ground motion characterization by cross wavelet transform *Tsoggerel Tsamba 1) and Masato Motosaka 2) 1), 2) International Research Institute of Disaster Science, Tohoku University, Sendai980-8579,Japan
More informationTohoku-oki event: Tectonic setting
Tohoku-oki event: Tectonic setting This earthquake was the result of thrust faulting along or near the convergent plate boundary where the Pacific Plate subducts beneath Japan. This map also shows the
More informationComparison of Short-Term and Time-Independent Earthquake Forecast Models for Southern California
Bulletin of the Seismological Society of America, Vol. 96, No. 1, pp. 90 106, February 2006, doi: 10.1785/0120050067 Comparison of Short-Term and Time-Independent Earthquake Forecast Models for Southern
More informationSOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION
SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION Kimiyuki Asano 1 and Tomotaka Iwata 2 1 Assistant Professor, Disaster Prevention
More informationEffects of Hurricanes on Ambient Noise in the Gulf of Mexico
Effects of Hurricanes on Ambient Noise in the Gulf of Mexico Mark A. Snyder Naval Oceanographic Office 1002 Balch Blvd. Stennis Space Center, MS 39522-5001 U.S.A. Abstract - Long-term omni-directional
More informationInterpretation of the Omori law
Interpretation of the Omori law Anatol V. Guglielmi Schmidt Institute of Physics of the Earth, Russian Academy of Sciences, ul. B. Gruzinskaya 10, Moscow, 13995 Russia e-mail: guglielmi@mail.ru This paper
More information