No Arabic abstract
An array of ten broadband stations was installed on the Popocatepetl volcano (Mexico) for five months between October 2002 and February 2003. 26 regional and teleseismic earthquakes were selected and filtered in the frequency time domain to extract the fundamental mode of the Rayleigh wave. The average dispersion curve was obtained in two steps. Firstly, phase velocities were measured in the period range [2-50] s from the phase difference between pairs of stations, using Wiener filtering. Secondly, the average dispersion curve was calculated by combining observations from all events in order to reduce diffraction effects. The inversion of the mean phase velocity yielded a crustal model for the volcano which is consistent with previous models of the Mexican Volcanic Belt. The overall crustal structure beneath Popocatepetl is therefore not different from the surrounding area, and the velocities in the lower crust are confirmed to be relatively low. Lateral variations of the structure were also investigated by dividing the network into four parts and by applying the same procedure to each sub-array. No well-defined anomalies appeared for the two sub-arrays for which it was possible to measure a dispersion curve. However, dispersion curves associated with individual events reveal important diffraction for 6 s to 12 s periods which could correspond to strong lateral variations at 5 to 10 km depth.
We explore interactions of elastic waves propagating in plates (with soil parameters) structured with concrete pillars buried in the soil. Pillars are 2 m in diameter, 30 m in depth and the plate is 50 m in thickness. We study the frequency range 5 to 10 Hz, for which Rayleigh wave wavelengths are smaller than the plate thickness. This frequency range is compatible with frequency ranges of particular interest in earthquake engineering. It is demonstrated in this paper that two seismic cloaks configurations allow for an unprecedented flow of elastodynamic energy associated with Rayleigh surface waves. The first cloak design is inspired by some approximation of ideal cloaks parameters within the framework of thin plate theory. The second, more accomplished but more involved, cloak design is deduced from a geometric transform in the full Navier equations that preserves the symmetry of the elasticity tensor but leads to Willis equations, well approximated by a homogenization procedure, as corroborated by numerical simulations. The two cloakss designs are strikingly different, and the superior efficiency of the second type of cloak emphasizes the necessity for rigor in transposition of existing cloakss designs in thin plates to the geophysics setting. Importantly, we focus our attention on geometric transforms applied to thick plates, which is an intermediate case between thin plates and semi-infinite media, not studied previously. Cloaking efficiency (reduction of the disturbance of the wave wavefront and its amplitude behind an obstacle) and protection (reduction of the wave amplitude within the center of the cloak) are studied for ideal and approximated cloaks parameters. These results represent a preliminary step towards designs of seismic cloaks for surface Rayleigh waves propagating in sedimentary soils structured with concrete pillars.
The study of volcanic inner density distributions using cosmic muons is an innovative method, which is still in a stage of development. This technique can be used to determine the average density along the muon track, as well as the density distribution within a given volume, by measuring the attenuation of the cosmic muon flux going through it. The aim is to study the volcano domes and magmatic conduit systems within a given time-interval. Our first application will be the Popocatepetl, a large active andesitic stratovolcano built in the Trans-Mexican volcanic arc. Its recent activity includes emplacement of a lava dome, with explosions and frequent scoria and ash emissions. This study is part of a longer-term project of volcanic hazard monitoring that includes other Mexican volcanoes, like the Colima. Muon detector design depends on the volume-of-interest dimensions, as well as on the image-taking frequency required to detect dynamic density variations. Our muon-tracker proposal includes 3 planes, each having 16 independent position-sensitive modules consisting on rectangular aluminum tubes ($10x20x320cm^{3}$) filed with a liquid scintillator. The light collection inside each module is carried out using a wave-length-shifting (WLS) fiber matrix, running along the aluminum-tube length, which is bundled together at the tube extremes. The luminous signal readout is carried out using one SiPM optically coupled to the WLS bundle at each modules end. The main detector characteristics, such as time resolution, surface uniformity, and signal amplitude reconstruction using the time-over-threshold technique, will be presented.
We discuss a methodology to identify observation points for muongraphy of active Colombian Volcanoes and it is found that only Cerro Machin can be studied.
This work presented a block triple-relaxation-time (B-TriRT) lattice Boltzmann model for simulating melting in a rectangular cavity heated from below at high Rayleigh (Ra) number (Ra = 108). The test of benchmark problem shows that present B-TriRT can dramatically reduce the numerical diffusion across the phase interface. In addition, the influences of the location of the heated region are investigated. The results indicate that the location of heated region plays an essential role in melting rate and the full melting occur earliest when the heated region is located in the middle region.
Gravimetric methods are expected to play a decisive role in geophysical modeling of the regional crustal structure applied to geoneutrino studies. GIGJ (GOCE Inversion for Geoneutrinos at JUNO) is a 3D numerical model constituted by ~46 x 10$^{3}$ voxels of 50 x 50 x 0.1 km, built by inverting gravimetric data over the 6{deg} x 4{deg} area centered at the Jiangmen Underground Neutrino Observatory (JUNO) experiment, currently under construction in the Guangdong Province (China). The a-priori modeling is based on the adoption of deep seismic sounding profiles, receiver functions, teleseismic P-wave velocity models and Moho depth maps, according to their own accuracy and spatial resolution. The inversion method allowed for integrating GOCE data with the a-priori information and regularization conditions through a Bayesian approach and a stochastic optimization. GIGJ fits the homogeneously distributed GOCE gravity data, characterized by high accuracy, with a ~1 mGal standard deviation of the residuals, compatible with the observation accuracy. Conversely to existing global models, GIGJ provides a site-specific subdivision of the crustal layers masses which uncertainties include estimation errors, associated to the gravimetric solution, and systematic uncertainties, related to the adoption of a fixed sedimentary layer. A consequence of this local rearrangement of the crustal layer thicknesses is a ~21% reduction and a ~24% increase of the middle and lower crust expected geoneutrino signal, respectively. Finally, the geophysical uncertainties of geoneutrino signals at JUNO produced by unitary uranium and thorium abundances distributed in the upper, middle and lower crust are reduced by 77%, 55% and 78%, respectively. The numerical model is available at http://www.fe.infn.it/u/radioactivity/GIGJ