Load-Unload Response Ratio and Accelerating Moment/Energy Release Critical Region Scaling and Earthquake Prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

Moment Tensor Analysis of the Acoustic Emission Source in the Rock Damage Process

To further investigate the mechanism of acoustic emission (AE) in the rock fracture experiment, moment tensor analysis was carried out. The AE sources characterized by crack sizes, orientations and fracture modes, are represented by a time-dependent momen

Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

Spatio-Temporal Scanning and Statistical Test of the Accelerating Moment Release (AMR) Model Using Australian Earthquake Data

The Accelerating Moment Release (AMR) preceding earthquakes with magnitude above 5 in Australia that occurred during the last 20 years was analyzed to test the Critical Point Hypothesis. Twelve earthquakes in the catalog were chosen based on a criterion for the number of nearby events. Results show that seven sequences with numerous events recorded leading up to the main earthquake exhibited accelerating moment release. Two occurred near in time and space to other earthquakes preceded by AM R. The remaining three sequences had very few events in the catalog so the lack of AMR detected in the analysis may be related to catalog incompleteness. Spatio-temporal scanning of AMR parameters shows that 80% of the areas in which AMR occurred experienced large events. In areas of similar background seismicity with no large events, 10 out of 12 cases exhibit no AMR, and two others are false alarms where AMR was observed but no large event followed. The relationship between AMR and Load-Unload Response Ratio (LURR) was studied. Both methods predict similar critical region sizes, however, the critical point time using AMR is slightly earlier than the time of the critical point LURR anomaly.

Self-organized generation of transverse waves in diverging cylindrical detonations

Self-organized generation of transverse waves associated with the transverse wave instabilities at a diverging cylindrical detonation front was numerically studied by solving two-dimensional Euler equations implemented with an improved two-step chemical kinetic model. After solution validation, four mechanisms of the transverse wave generation were identified from numerical simulations, and referred to as the concave front focusing, the kinked front evolution, the wrinkled front evolution and the transverse wave merging, respectively. The propagation of the cylindrical detonation is maintained by the growth of the transverse waves that match the rate of increase in surface area of the detonation front to asymptotically approach a constant average number of transverse waves per unit length along the circumference of the detonation front. This cell bifurcation phenomenon of cellular detonations is discussed in detail to gain better understanding on detonation physics.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Transverse evolution of a plasma channel in air induced by a femtosecond laser

We investigate the evolution of filamentation in air by using a longitudinal diffraction method and a plasma fluorescence imaging technique. The diameter of a single filament in which the intensity is clamped increases as the energy of the pump light pulse increases, until multiple filaments appear. (c) 2006 Optical Society of America.

Generation of attosecond pulses in a system with permanent dipole moment

The generation of attosecond pulses in a two-level system with permanent dipole moment is investigated. It is shown due to the presence of permanent dipole moments, that the plateau of the high-order harmonic generation spectrum can be extended to X-ray range. Moreover, attosecond pulses with higher intensity can be synthesized by using both even and odd harmonics because of their quantum interference. (c) 2006 Elsevier B.V. All rights reserved.

Existence of the transverse relaxation time in optically excited bulk semiconductors

Two basic types of depolarization mechanisms, carrier-carrier (CC) and carrier-phonon (CP) scattering, are investigated in optically excited bulk semiconductors (3D), in which the existence of the transverse relaxation time is proven based on the vector property of the interband transition matrix elements. The dephasing rates for both CC and CP scattering are determined to be equal to one half of the total scattering-rate-integrals weighted by the factors (1 - cos chi), where chi are the scattering angles. Analytical expressions of the polarization dephasing due to CC scattering are established by using an uncertainty broadening approach, and analytical ones due to both the polar optical-phonon and non-polar deformation potential scattering (including inter-valley scattering) are also presented by using the sharp spectral functions in the dephasing rate calculations. These formulas, which reveal the trivial role of the Coulomb screening effect in the depolarization processes, are used to explain the experimental results at hand and provide a clear physical picture that is difficult to extract from numerical treatments.

The origin of the transverse relaxation time in optically excited semiconductor quantum wells

The origin of the transverse relaxation time in optically excited semiconductor quantum wells is investigated based on the vector property of the interband transition matrix elements. The dephasing rate due to carrier-carrier (CC) scattering is found to be equal to half of the common momentum relaxation rate. The analytical expression of the polarization dephasing due to CC scattering in two-dimension is established and the dependence of the dephasing rate Gamma(cc) on the carrier density N is determined to be Gamma(cc) = constant (.) N-1/2, which is used to explain the experimental results and provides a promising physical picture. (C) 2004 Elsevier B.V. All rights reserved.

Tunable transverse superresolution with phase-only pupil filters

A new set of pure phase filters for realizing transverse superresolution is presented in this paper. The filters, whose significant features are their ability to tune and their simplicity, consist of one half-wave plate between two quarter-wave plates; the half-wave plate is made of two zones that can rotate with respect to each other. By rotating any zone of the half-wave plate, the central lobe width of the irradiance point spread function (PSF) in the transverse direction can be tunably reduced. At the same time, the axial intensity distribution is analysed in detail.

Transverse superresolution with the radial continuous transmittance filter

The radial continuous transmittance filter is presented to realize transverse superresolution. It consists of two parallel polarizers and a radial birefringent element sandwiched between of them. By adjusting the angle between optical axis of the radial birefringent element and the polarization direction of the polarizers, transverse superresolution can be realized. But transverse superresolution is obtained at the cost of the axial resolution and the increase of the side-lobes in strength. So we then mend such filter, with it not only enhance the transverse resolution but also suppress the influence of the side-lobes and the reduction of the axial resolution. At the same time, the Strehl ratio increases. The advantage of such a filter used in superresolution technique is that it is easy to fabricate because its fabrication does not deal with the variation of the phase. (c) 2005 Elsevier GmbH. All rights reserved.