Showing papers in "arXiv: Geophysics in 2010"

A review of wildland fire spread modelling, 1990-present, 1: Physical and quasi-physical models

Abstract: In recent years, advances in computational power and spatial data analysis (GIS, remote sensing, etc) have led to an increase in attempts to model the spread and behaviour of wildland fires across the landscape. This series of review papers endeavours to critically and comprehensively review all types of surface fire spread models developed since 1990. This paper reviews models of a physical or quasi-physical nature. These models are based on the fundamental chemistry and/or physics of combustion and fire spread. Other papers in the series review models of an empirical or quasi-empirical nature, and mathematical analogues and simulation models. Many models are extensions or refinements of models developed before 1990. Where this is the case, these models are also discussed but much less comprehensively.

Scaling of plate-tectonic convection with pseudoplastic rheology

Abstract: The scaling of plate-tectonic convection is investigated by simulating thermal convection with pseudoplastic rheology and strongly temperature-dependent viscosity. The effect of mantle melting is also explored with additional depth-dependent viscosity. Heat-flow scaling can be constructed with only two parameters, the internal Rayleigh number and the lithospheric viscosity contrast, the latter of which is determined entirely by rheological properties. The critical viscosity contrast for the transition between plate-tectonic and stagnant-lid convection is found to be proportional to the square root of the internal Rayleigh number. The relation between mantle temperature and surface heat flux on Earth is discussed on the basis of these scaling laws, and the inverse relationship between them, as previously suggested from the consideration of global energy balance, is confirmed by this fully dynamic approach. In the presence of surface water to reduce the effective friction coefficient, the operation of plate tectonics is suggested to be plausible throughout the Earth history.

The dynamo bifurcation in rotating spherical shells

Abstract: We investigate the nature of the dynamo bifurcation in a configuration applicable to the Earth's liquid outer core, i.e. in a rotating spherical shell with thermally driven motions. We show that the nature of the bifurcation, which can be either supercritical or subcritical or even take the form of isola (or detached lobes) strongly depends on the parameters. This dependence is described in a range of parameters numerically accessible (which unfortunately remains remote from geophysical application), and we show how the magnetic Prandtl number and the Ekman number control these transitions.

Source geometry from exceptionally high resolution long period event observations at Mt Etna during the 2008 eruption

Abstract: During the second half of June, 2008, 50 broadband seismic stations were deployed on Mt Etna volcano in close proximity to the summit, allowing us to observe seismic activity with exceptionally high resolution. 129 long period events (LP) with dominant frequencies ranging between 0.3 and 1.2 Hz, were extracted from this dataset. These events form two families of similar waveforms with different temporal distributions. Event locations are performed by cross-correlating signals for all pairs of stations in a two-step scheme. In the first step, the absolute location of the centre of the clusters was found. In the second step, all events are located using this position. The hypocentres are found at shallow depths (20 to 700 m deep) below the summit craters. The very high location resolution allows us to detect the temporal migration of the events along a dike-like structure and 2 pipe shaped bodies, yielding an unprecedented view of some elements of the shallow plumbing system at Mount Etna. These events do not seem to be a direct indicator of the ongoing lava flow or magma upwelling.

Retrospective Evaluation of the Five-Year and Ten-Year CSEP-Italy Earthquake Forecasts

Abstract: On 1 August 2009, the global Collaboratory for the Study of Earthquake Predictability (CSEP) launched a prospective and comparative earthquake predictability experiment in Italy. The goal of the CSEP-Italy experiment is to test earthquake occurrence hypotheses that have been formalized as probabilistic earthquake forecasts over temporal scales that range from days to years. In the first round of forecast submissions, members of the CSEP-Italy Working Group presented eighteen five-year and ten-year earthquake forecasts to the European CSEP Testing Center at ETH Zurich. We considered the twelve time-independent earthquake forecasts among this set and evaluated them with respect to past seismicity data from two Italian earthquake catalogs. In this article, we present the results of tests that measure the consistency of the forecasts with the past observations. Besides being an evaluation of the submitted time-independent forecasts, this exercise provided insight into a number of important issues in predictability experiments with regard to the specification of the forecasts, the performance of the tests, and the trade-off between the robustness of results and experiment duration. We conclude with suggestions for the future design of earthquake predictability experiments.

Solitary Wave Benchmarks in Magma Dynamics

Abstract: We present a model problem for benchmarking codes that investigate magma migration in the Earth's interior. This system retains the essential features of more sophisticated models, yet has the advantage of possessing solitary wave solutions. The existence of such exact solutions to the nonlinear problem make it an excellent benchmark problem for combinations of solver algorithms. In this work, we explore a novel algorithm for computing high quality approximations of the solitary waves and use them to benchmark a semi-Lagrangian Crank-Nicholson scheme for a finite element discretization of the time dependent problem.

Climate Change and Its Causes, A Discussion About Some Key Issues

Abstract: This article discusses the limits of the Anthropogenic Global Warming Theory advocated by the Intergovernmental Panel on Climate Change. A phenomenological theory of climate change based on the physical properties of the data themselves is proposed. At least 60% of the warming of the Earth observed since 1970 appears to be induced by natural cycles which are present in the solar system. A climatic stabilization or cooling until 2030-2040 is forecast by the phenomenological model.

Segmentation of Fault Networks Determined from Spatial Clustering of Earthquakes

Abstract: We present a new method of data clustering applied to earthquake catalogs, with the goal of reconstructing the seismically active part of fault networks. We first use an original method to separate clustered events from uncorrelated seismicity using the distribution of volumes of tetrahedra defined by closest neighbor events in the original and randomized seismic catalogs. The spatial disorder of the complex geometry of fault networks is then taken into account by defining faults as probabilistic anisotropic kernels, whose structures are motivated by properties of discontinuous tectonic deformation and previous empirical observations of the geometry of faults and of earthquake clusters at many spatial and temporal scales. Combining this a priori knowledge with information theoretical arguments, we propose the Gaussian mixture approach implemented in an Expectation-Maximization (EM) procedure. A cross-validation scheme is then used and allows the determination of the number of kernels that should be used to provide an optimal data clustering of the catalog. This three-steps approach is applied to a high quality relocated catalog of the seismicity following the 1986 Mount Lewis ($M_l=5.7$) event in California and reveals that events cluster along planar patches of about 2 km$^2$, i.e. comparable to the size of the main event. The finite thickness of those clusters (about 290 m) suggests that events do not occur on well-defined euclidean fault core surfaces, but rather that the damage zone surrounding faults may be seismically active at depth. Finally, we propose a connection between our methodology and multi-scale spatial analysis, based on the derivation of spatial fractal dimension of about 1.8 for the set of hypocenters in the Mnt Lewis area, consistent with recent observations on relocated catalogs.

Pressure Solution in Sedimentary Basins: Effect of Temperature Gradient

Abstract: Pressure solution is an important process in sedimentary basins, and its behaviour depends mainly on the sediment rheology and temperature distribution The compaction relation of pressure solution is typically assumed to be a viscous one and is often written as a relationship between effective stress and strain rate A new derivation of viscous compaction relation is formulated based on more realistic boundary conditions at grain contacts A nonlinear diffusion problem with a moving boundary is solved numerically and a simple asymptotic solution is given to compare with numerical simulations Pressure solution is significantly influenced by the temperature gradient Porosity reduction due to pressure solution is enhanced in an environment with a higher thermal gradient, while porosity decreases much slowly in the region where the thermal gradient is small Pressure solution tends to complete more quickly at shallower depths and earlier time in higher temperature environment than that in a low one These features of pressure solution in porous sediments are analysed using a perturbation method to get a solution for the steady state Comparison with real data shows a reasonably very good agreement

Source Mechanism of Long Period events recorded by a high density seismic network during the 2008 eruption on Mt Etna

Abstract: 129 Long Period (LP) events, divided into two families of similar events, were recorded by the 50 stations deployed on Mount Etna in the second half of June 2008. During this period lava was flowing from a lateral fracture after a summit strombolian eruption. In order to understand the mechanisms of these events, we perform moment tensor inversions. Inversions are initially kept unconstrained to estimate the most likely mechanism. Numerical tests show that unconstrained inversion leads to reliable moment tensor solutions because of the close proximity of numerous stations to the source positions. However, single forces cannot be accurately determined as they are very sensitive to uncertainties in the velocity model. Constrained inversions for a crack, a pipe or an explosion then allow us to accurately determine the structural orientations of the source mechanisms. Both numerical tests and LP event inversions emphasise the importance of using stations located as close as possible to the source. Inversions for both families show mechanisms with a strong volumetric component. These events are most likely generated by cracks striking SW-NE for both families and dipping 70\deg SE (Family 1) and 50\deg NW (Family 2). For Family 1 events, the crack geometry is nearly orthogonal to the dike-like structure along which events are located, while for Family 2 the location gave two pipe-like bodies which belong to the same plane as the crack mechanism. The orientations of the cracks are consistent with local tectonics, which shows a SW-NE weakness direction. The LP events appear to be a response to the lava fountain occurring on the 10th of May, 2008 as opposed to the flank lava flow.

Constraining the shape of a gravity anomalous body using reversible jump Markov chain Monte Carlo

Abstract: Typical geophysical inversion problems are ill-posed, non-linear and non-unique. Sometimes the problem is trans-dimensional, where the number of unknown parameters is one of the unknowns, which makes the inverse problem even more challenging. Detecting the shape of a geophysical object underneath the earth surface from gravity anomaly is one of such complex problems, where the number of geometrical parameters is one of the unknowns. To deal with the difficulties of non-uniqueness, ill-conditioning and nonlinearity, a statistical Bayesian model inference approach is adopted. A reversible jump Markov chain Monte Carlo (RJMCMC) algorithm is proposed to overcome the difficulty of trans-dimensionality. Carefully designed within-model and between-model Markov chain moves are implemented to reduce the rate of generating inadmissible geometries, thus achieving good overall efficiency in the Monte Carlo sampler. Numerical experiments on a 2-D problem show that the proposed algorithm is capable of obtaining satisfactory solutions with quantifiable uncertainty to a challenging trans-dimensional geophysical inverse problem. Solutions from RJMCMC appear to be parsimonious for the given prior, in the sense that among the models satisfactorily represent the true model, models with higher posterior probabilities tend to have fewer number of parameters. The proposed numerical algorithm can be readily adapted to other similar trans-dimensional geophysical inverse applications. Keywords: trans-dimensional geophysical inversion, reversible jump Markov chain Monte Carlo; gravity anomaly.

A Unified Approach to Mechanical Compaction, Pressure Solution, Mineral Reactions and the Temperature Distribution in Hydrocarbon Basins

Abstract: In modelling sediment compaction and mineral reactions, the rheological behaviour of sediments is typically considered as poroelastic or purely viscous. In fact, compaction due to pressure solution and mechanical processes in porous media is far more complicated. A generalised model of viscoelastic compaction and the smectite to illite mineral reaction in hydrocarbon basins is presented. A one-step dehydration model of the mineral reaction is assumed. The obtained nonlinear governing equations are solved numerically and different combinations of physical parameters are used to simulate realistic situations in typical sedimentary basins. Comparison of numerical simulations with real data has shown very good agreement with respect to both the porosity profile and the mineral reaction.

Evidence for Anthropogenic Surface Loading as Trigger Mechanism of the 2008 Wenchuan Earthquake

Abstract: Two and a half years prior to China's M7.9 Wenchuan earthquake of May 2008, at least 300 million metric tons of water accumulated with additional seasonal water level changes in the Minjiang River Valley at the eastern margin of the Longmen Shan. This article shows that static surface loading in the Zipingpu water reservoir induced Coulomb failure stresses on the nearby Beichuan thrust fault system at 7 earthquakes. Finally, the daily event rate of the micro-seismicity (M>0.5) correlates well with the static stress perturbations, indicating destabilization.

NUV/Blue spectral observations of sprites in the 320-460 nm region: N2 (2PG) Emissions

Abstract: : A near-ultraviolet (NUV) spectrograph (320-460 nm) was flown on the EXL98 aircraft sprite observation campaign during July 1998. In this wavelength range video rate (60 fields/sec) spectrographic observations found the NUV/blue emissions to be predominantly N2(2PG). The negligible level of N+2 (1NG) present in the spectrum is confirmed by observations of a co-aligned, narrowly filtered 427.8 nm imager and is in agreement with previous ground-based filtered photometer observations. The synthetic spectral fit to the observations indicates a characteristic energy of 1.8 eV, in agreement with our other NUV observations.

Fractal Models of Earthquake Dynamics

Abstract: Our understanding of earthquakes is based on the theory of plate tectonics Earthquake dynamics is the study of the interactions of plates (solid disjoint parts of the lithosphere) which produce seismic activity Over the last about fifty years many models have come up which try to simulate seismic activity by mimicking plate plate interactions The validity of a given model is subject to the compliance of the synthetic seismic activity it produces to the well known empirical laws which describe the statistical features of observed seismic activity Here we present a review of two such models of earthquake dynamics with main focus on a relatively new model namely The Two Fractal Overlap Model

Integrated Geophysical Measurements on a Test Site for Detection of Buried Steel Drums

Abstract: Geophysical methods are increasingly used to detect and locate illegal waste disposal and buried toxic steel drums. This study describes the results of a test carried out in clayey-sandy ground where 12 empty steel drums had previously been buried at 4-5 m below ground level. This test was carried out with three geophysical methods for steel-drum detection: a magnetometric survey, electrical resistivity tomography with different arrays, and a multifrequency frequency-domain electromagnetic induction survey. The data show that as partially expected, the magnetometric and electromagnetic induction surveys detected the actual steel drums buried in the subsurface, while the electrical resistivity tomography mainly detected the changes in some of the physical properties of the terrain connected with the digging operations, rather than the actual presence of the steel drums.

Enhancing the Google imagery using a wavelet filter

Abstract: In some previous papers we proposed the use of free software for a processing of the Google satellite imagery. Here we discuss the use of a wavelet filter for the same purposes. This filter is a tool included in a freely downloadable software (Iris), well-known for the processing of astronomical images. Combining the image obtained after applying the wavelet filter, with an image created with Gimp and AstroFracTool, the visibility of the landforms, as obtained from Google Maps, is strongly increased. Among several possible examples, we proposed a crater, a paleochannel and the Great Bend of the Nile.

Comment on "NOx production in laboratory discharges simulating blue jets and red sprites" by H. Peterson et al

Abstract: Peterson et al. [2009] intend to determine the NOx production of sprites and blue jets from laboratory experiments. Sprite discharges are known to be physically similar to streamer discharges at standard temperature and pressure. However, streamer experiments that simulate geophysical discharges require voltage pulses of one polarity, in contrast to the oscillating voltages of Peterson et al. Furthermore, Peterson et al. make errors in their comparison of current densities, sprite volumes and gas densities. Finally, there is a calculation error of 6 orders of magnitude in an important result.

Ionized Nitrogen Mono-hydride Bands are Identified in the Pre-solar and Carbonado Diamond Spectra

Abstract: None of the well established Nitrogen related IR absorption bands, common in synthetic and terrestrial diamonds, have been identified in the pre-solar diamond spectra. In the carbonado diamond spectra only the single nitrogen impurity (C centre) is identified and the assignments of the rest of the nitrogen-related bands are still debated. It is speculated that the unidentified bands in the Nitrogen absorption region are not induced by Nitrogen but rather by Nitrogen-hydrides because in the interstellar environment Nitrogen reacts with Hydrogen and forms NH+; NH; NH2; NH3. Among these Hydrides the electronic configuration of NH+ is the closest to Carbon. Thus this ionized Nitrogen-mono-hydride is the best candidate to substitute Carbon in the diamond structure. The bands of the substitutional NH+ defect are deduced by red shifting the irradiation induced N+ bands due to the mass of the additional Hydrogen. The six bands of the NH+ defects are identified in both the pre-solar and the carbonado diamond spectra. The new assignments identify all of the nitrogen-related bands in the spectra, indicating that pre-solar and carbonado diamonds contain only single nitrogen impurities.

The Critical State and the Steady-State solution in Granular Solid Hydrodynamics

Abstract: Central South University, Changsha 410083, China(Dated: June 29, 2010)The transition from partially elastic to perfectly plastic behavior is considered employing“granularsolid hydrodynamics.” For two reasons, this unjamming step is rendered coherent and transparent:First, the underlying physics is identified, being the relaxation of the elastic strain to its steady-statesolution. Second, many results are obtained as simple, analytic expressions.

Hot-Cold Spots in Italian Macroseismic Data

Abstract: The site effect is usually associated with local geological conditions, which increase or decrease the level of shaking compared with standard attenuation relations. We made an attempt to see in the macroseismic data of Italy some other effects, namely, hot/cold spots in the terminology of Olsen (2000), which are related to local fault geometry rather than to soil conditions. We give a list of towns and villages liable to amplify (+) or to reduce (-) the level of shaking in comparison with the nearby settlements. Relief and soil conditions cannot always account for the anomalous sites. Further, there are sites where both (+) and (-) effects are observed depending on the earthquake. The opposite effects can be generated by events from the same seismotectonic zone and along the same direction to the site. Anomalous sites may group themselves into clusters of different scales. All isolated anomalous patterns presented in this paper can be used in hazard analysis, in particular, for the modeling and testing of seismic effects.

The nonlinear anomalous lattice elasticity associated with the high-pressure phase transition in spodumene: A high precission static compression study

Abstract: The high-pressure behavior of the lattice elasticity of spodumene, LiAlSi2O6, was studied by static compression in a diamond-anvil cell up to 9.3 GPa. Investigations by means of single-crystal XRD and Raman spectroscopy within the hydrostatic limits of the pressure medium focus on the pressure ranges around similar to 3.2 and similar to 7.7 GPa, which have been reported previously to comprise two independent structural phase transitions. While our measurements confirm the well-established first-order C2/c-P2(1)/c transformation at 3.19 GPa (with 1.2% volume discontinuity and a hysteresis between 0.02 and 0.06 GPa), both unit-cell dimensions and the spectral changes observed in high-pressure Raman spectra give no evidence for structural changes related to a second phase transition. Monoclinic lattice parameters and unit-cell volumes at in total 59 different pressure points have been used to re-calculate the lattice-related properties of spontaneous strain, volume strain, and the bulk moduli as a function of pressure across the transition. A modified Landau free energy expansion in terms of a one component order parameter has been developed and tested against these experimentally determined data. The Landau solution provides a much better reproduction of the observed anomalies than any equation-of-state fit to data sets truncated below and above P (tr), thus giving Landau parameters of K (0) = 138.3(2) GPa, K' = 7.46(5), lambda (V) = 33.6(2) GPa, a = 0.486(3), b = -29.4(6) GPa and c = 551(11) GPa.

Gambling scores in earthquake prediction analysis

Abstract: The number of successes 'n' and the normalized measure of space-time alarm 'tau' are commonly used to characterize the strength of an earthquake prediction method and the significance of prediction results. To evaluate better the forecaster's skill, it has been recently suggested to use a new characteristic, the gambling score R, which incorporates the difficulty of guessing each target event by using different weights for different alarms. We expand the class of R-characteristics and apply these to the analysis of results of the M8 prediction algorithm. We show that the level of significance 'alfa' strongly depends (1) on the choice of weighting alarm parameters, (2) on the partitioning of the entire alarm volume into component parts, and (3) on the accuracy of the spatial rate of target events, m(dg). These tools are at the disposal of the researcher and can affect the significance estimate in either direction. All the R-statistics discussed here corroborate that the prediction of 8.0<=M<8.5 events by the M8 method is nontrivial. However, conclusions based on traditional characteristics (n,tau) are more reliable owing to two circumstances: 'tau' is stable since it is based on relative values of m(.), and the 'n' statistic enables constructing an upper estimate of 'alfa' taking into account the uncertainty of m(.).

Temperature induced pore fluid pressurization in geomaterials

Abstract: Summary The theoretical basis of the thermal response of the fluid-saturated porous materials in undrained condition is presented. It has been demonstrated that the thermal pressurization phenomenon is controlled by the discrepancy between the thermal expansion of the pore fluid and of the solid phase, the stress-dependency of the compressibility and the non-elastic volume changes of the porous material. For evaluation of the undrained thermo-poro-elastic properties of saturated porous materials in conventional triaxial cells, it is important to take into account the effect of the dead volume of the drainage system. A simple correction method is presented to correct the measured pore pressure change and also the measured volumetric strain during an undrained heating test. It is shown that the porosity of the tested material, its drained compressibility and the ratio of the volume of the drainage system to the one of the tested sample, are the key parameters which influence the most the error induced on the measurements by the drainage system. An example of the experimental evaluation of undrained thermoelastic parameters is presented for an undrained heating test performed on a fluid-saturated granular rock.

Crater-like landform in Bayuda desert (a processing of satellite images)

Abstract: The paper is proposing a survey of a region in Sudan, the Bayuda desert, using the satellite images as obtained from Google Maps The images reveal the ring granitic structure of the region enclosed by a bend of river Nile To enhance the features of the landform, images are processed with a method based on fractional calculus, able to increase the rendering of edges without deteriorating the overall quality of images Besides the ringed structure of the region, a huge crater-like structure with a diameter of 10 km is evidenced

Arabia: from craters to stone circles

Abstract: The Arabia Shield has a volcanic nature inside. A region of the Western Saudi Arabia is in fact covered with vast fields of lava known as harraat. These lands are spotted by many stone circles and other quite interesting archaeological remains of the Neolithic period, such as the "desert kites", the hunters used to guide the game across the harrah in some corrals. With Google Maps, we can observe both sceneries, the volcanic nature of the land and a portrait of Arabia during the Neolithic times.

Earthquake prediction analysis: The M8 algorithm

Abstract: The quality of space-time earthquake prediction is usually characterized by a two-dimensional error diagram (n,tau), where 'n' is the rate of failures-to-predict and 'tau' is the normalized measure of space-time alarm. The most interesting space measure for analysis of a prediction strategy is the rate of target events m(dg) in a sub-area 'dg'. In this case the quantity H=1-(n+tau) determines the prediction capability of the strategy. The uncertainty of m(dg) causes difficulties in estimating 'H' and the statistical significance, 'alfa', of prediction results. We investigate this problem theoretically and show how the uncertainty of the measure can be taken into account in two situations, viz., the estimation of 'alfa' and the construction of a confidence zone for (n,tau)-parameters of the random strategies. We use our results to analyse the M8 earthquake prediction algorithm.

Symbolic landforms created by ancient earthworks near Lake Titicaca

Abstract: Interesting landforms created by an ancient network of earthworks are shown, using Google satellite imagery enhanced by an image processing This network covers a large part of the land near the Titicaca Lake Satellite images clearly display the slopes of hills criss-crossed with terrace walls and the surfaces of the plains covered with raised fields, indicating that this was once a highly productive agricultural place for the south central Andes Some of the landforms are rather remarkable, having a clear symbolic function Among them, there are structures which seem to represent birds, where ponds are their eyes

A comparison of new calculations of 10be production in the earths polar atmosphere by cosmic rays with 10be concentration measurements in polar ice cores between 1939-2005 - a troubling lack of concordance paper #1

Abstract: Using new calculations of 10Be production in the Earths atmosphere which are based on direct measurements of the 11-year solar modulation effects on galactic cosmic rays and spacecraft measurements of the cosmic ray energy spectrum, we have calculated the yearly average production of 10Be in the Earths atmosphere by galactic and solar cosmic rays since 1939. During the last six 11-year cycles the average amplitude of these production changes is 36%. These predictions are compared with measurements of 10Be concentration in polar ice cores in both the Northern and Southern hemisphere over the same time period. We find a large scatter between the predicted and measured yearly average data sets and a low cross correlation ~0.30. Also the normalized regression line slope between 10Be production changes and 10Be concentration changes is found to be only 0.4-0.6; much less than the value of 1.0 expected for a simple proportionality between these quantities, as is typically used for historical projections of the relationship between 10Be concentration and solar activity. The distribution of yearly averages in the 10Be concentration level in the data from the Dye-3 ice core in Greenland for the time period 1939-1985, contains a "spike" of high concentration one year averages which is not seen in the production calculations. These and other features suggest that galactic cosmic ray intensity changes which affect the production of 10Be in the Earths atmosphere are not the sole source of the 10Be concentration changes and confirm the importance of other effects, for example local and regional climatic effects, which could be of the same magnitude as the 10Be production changes.

Radiation belt electron precipitation in the upper ionosphere at middle latitudes before strong earthquakes

Abstract: In this article we present examples of a wider study of space-time correlation of electron precipitation event of the Van Allen belts with the position and time of occurrence of strong (M>6.5) earthquakes. The study is based on the analysis of observations of electron bursts (EBs) with energies 70 - 2350 keV at middle geographic latitudes, which were detected by DEMETER satellite (at an altitude of ~700 km). The EBs show a relative peak-to-background increase usually < 100, they have a time duration ~0.5 - 3 min, energy spectrum with peaks moving in higher energies as the satellite moves towards the equator, and highest energy limit <~500 keV. The EBs are observed in the presence of VLF waves. The flux-time profile of the EBs varies in East Asia and Mediterranean Sea at the similar geographic latitudes, due to the differentiation of the magnitude of the earth's magnetic field. The most important result of our study is the characteristic temporal variation of electron precipitation variation which begins with incremental rate several days / a few weeks before major earthquakes, then presents a maximum, and finally become weaker, with a minimun or a pause of the electron precipitation some hours before the onset of the earthquake.