530 likes | 632 Views
Statistical Analysis of Electrostatic Turbulences over Seismic Regions. T. Onishi and J.J. Berthelier Centre d'Etude des Environnements Terrestre et Planétaires (CETP) Saint-Maur, France. EuroPlaNet Strategic Workshop on Earthquakes: Ground-based and Space Observations
E N D
Statistical Analysis of Electrostatic Turbulences over Seismic Regions T. Onishi and J.J. Berthelier Centre d'Etude des Environnements Terrestre et Planétaires (CETP) Saint-Maur, France EuroPlaNet Strategic Workshop on Earthquakes: Ground-based and Space Observations Graz, Austria, 1-2 June, 2007
OUTLINE • Purpose of this study • Frequency Classification of power spectra of ELF/VLF emissions • Interferences from other instruments (ISL) • Selection of Earthquake data • Preliminary results. • Conclusion and Future Work
TYPICAL CHARACTERISTICS OF WAVES DETECTED BY ICE Log(μV2/Hz) Ordinary ELF hiss Electrostatic Turbulance ELF hiss below Cross-over freq. ΩH+
Purpose of study ICE level-1 VLF spectra data with characteristic frequencies (ΩH+, etc…) calculated from IAP data. Purpose: characterize the shape of the frequency spectra to determine emissions with different origin, propagation condition etc… Track the characteristics of these emissions to search for changes linked with seismic activity. ΩH+
Purpose of study ICE level-1 VLF spectra data with characteristic frequencies (ΩH+, etc…) calculated from IAP data. Purpose: characterize the shape of the frequency spectra to determine emissions with different origin, propagation condition etc… Track the characteristics of these emissions to search for changes linked with seismic activity. It is easy to do it on just one spectral plot. But, there are ten of thousands of them. ΩH+
Automatic identification of characteristic frequencies---- Procedure ---- Savitzky-Golay smoothing ”pwr_smooth” “Minimum” filter in time Digitalization of “pwr_smooth” ”pwr_smooth_bin” Detection of characteristic frequencies on “pwr_smooth” from ”pwr_smooth_bin”
Automatic identification of characteristic frequencies---- “minimum” filter in time domain ----
Automatic identification of characteristic frequencies2. Smoothing in frequency domain : Savitzky-Golay filter
Automatic identification of characteristic frequencies2. Smoothing in frequency domain
Automatic identification of characteristic frequencies2. Smoothing in frequency domain : Digital filter Reduces number of candidate points and makes it easy to pick one
Automatic identification of characteristic frequencies2. Smoothing in frequency domain
So the analytical tool is ready! We can start the statistical study, using actual earthquake data!
So the analytical tool is ready! We can start the statistical study, using actual earthquake data! To begin with, let us see the electrostatic turbulence at low frequency.
So the analytical tool is ready! We can start the statistical study, using actual earthquake data! To begin with, let us see the electrostatic turbulence at low frequency. But…. We have a problem !!!.
Interferences due to the swept Langmuir probe ISL • Parasites are present in the form of modulation Why? Burst mode
Interferences due to the swept Langmuir probe ISL • Parasites are present in the form of modulation Why? • Knowledge of the waveform of a parasite is not enough to separate parasites from the natural emissions. (Nonlinear effects) Only Burst Mode Burst mode
Interferences due to the swept Langmuir probe ISL • Parasites are present in the form of modulation Why? • Knowledge of the waveform of a parasite is not enough to separate parasites from the natural emissions. (Nonlinear effects) Only Burst Mode • Removal of parasite signals is very critical for the analysis of low frequency emissions (i.e. Electrostatic turbulences) Detection and removal of parasites Burst mode
Parasite position inside a packet Why modulation ? VLF power spectra and waveform with Blackmann-Harris window ULF potential variations (S1 and S2) VLF power spectra obtained from VLF waveform VLF waveform and time-average
Parasite position inside a packet Why modulation ? VLF power spectra and waveform with Blackmann-Harris window ULF potential variations (S1 and S2) VLF power spectra obtained from VLF waveform VLF waveform and time-average
Parasite position inside a packet Why modulation ? VLF power spectra and waveform with Blackmann-Harris window ULF potential variations (S1 and S2) VLF power spectra obtained from VLF waveform VLF waveform and time-average
Parasite position inside a packet Why modulation ? VLF power spectra and waveform with Blackmann-Harris window ULF potential variations (S1 and S2) VLF power spectra obtained from VLF waveform VLF waveform and time-average
Parasite position inside a packet Why modulation ? VLF power spectra and waveform with Blackmann-Harris window ULF potential variations (S1 and S2) VLF power spectra obtained from VLF waveform VLF waveform and time-average
Parasite position inside a packet Why modulation ? VLF power spectra and waveform with Blackmann-Harris window ULF potential variations (S1 and S2) VLF power spectra obtained from VLF waveform VLF waveform and time-average
Parasite position inside a packet Why modulation ? VLF power spectra and waveform with Blackmann-Harris window ULF potential variations (S1 and S2) VLF power spectra obtained from VLF waveform VLF waveform and time-average
Why modulation ? • Voltage sweep of Langmuir probe is performed every ~1.0112 second. • Exact duration of one packet is 0.0512 second. • Relative time position voltage drop inside a packet is periodic in every 4 times.
We understand why parasites show a modulation. Now how can we remove them?
How to detect a parasite position Potential peak position is not precise enough to define the corresponding potential drop. Peak of the sweep is detected first. A point where a potential decreases by 1/e are detected on both S1. Parasite position is defined as the mid point of these two.
Removal of parasite signals from VLF waveform data and construction of clean VLF power spectra
Removal of parasite signals from VLF waveform data and construction of clean VLF power spectra
Removal of parasite signals from VLF waveform data and construction of clean VLF power spectra 40 spectra averaged with parasites 40 spectra averaged except those with parasites
After the voltage change in a potential sweep Sweep voltage is reduced from 7.6V to 3.8V from Orbit 2154.0. Parasite effect is also reduced. But……
After the voltage change in a potential sweep But with different contour levels, the parasite effect is evident. Parasite effects may be small. But changes due to EQs may be smaller. Such tiny changes can be masked by parasites. Therefore, parasite removal is important for low frequency analysis!
Now we have a tool and clean data. We can start the statistical study with actual EQ data.
Now we have a tool and clean data. We can start the statistical study with actual EQ data. But… which EQ data can we use?
Now we have a tool and clean data. We can start the statistical study with actual EQ data. But… which EQ data can we use? A bunch of earthquakes often occurs in the sameregion and at about the same time (few days difference)
Now we have a tool and clean data. We can start the statistical study with actual EQ data. But… which EQ data can we use? A bunch of earthquakes often occurs in the same region and at about the same time (few days difference) If we should find an anomaly before one earthquake, How do we know if it is a precursor to the earthquake Or a post-seismic phenomenon of another earthquake?
Earthquake selection • Defining M as the magnitude of the main earthquake, following conditions are checked in EQ selection. • All EQs of magnitude smaller than M-2 are ignored. • No EQs within the dobrovolny distance of the main EQ in the preceding 10 days. • Number of EQs selected is in total 664. • M > 7.0 : 4 • 6.0 < M < 6.9 : 27 • 5.0 < M < 5.9 : 633
First attempt to statistical study of electrostatic turbulenceSelection Criteria • Earthquakes of magnitude 5.4 ≤ M ≤ 5.9 are used. There are 169 earthquakes. • Orbit data are used if … • In Burst mode, • ap-index ≤ 15, • -40 < Latitude < 40, • No MTB activation, • Within 200km from the epicenter In total, 59 orbits for 32 earthquakes remained.
First attempt to statistical study of electrostatic turbulenceFirst Result : Power Spectre at 20Hz Time (day) relative to EQ
First attempt to statistical study of electrostatic turbulenceFirst Result : Power Spectre at 20Hz Example with 4 EQ data sets and 7 orbits
First attempt to statistical study of electrostatic turbulenceFirst Result : Power Spectre at 20Hz Time (day) relative to EQ
First attempt to statistical study of electrostatic turbulenceFirst Result : Frequency for log(power) = -2.0 Time (day) relative to EQ
First attempt to statistical study of electrostatic turbulence First Result : Frequency for log(power) = -2.0 Time (day) relative to EQ
Better conditions: Ap-index < 10. Magnitudes 5.0 < M and Satellite distance < 100km Time (day) relative to EQ