Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.     

A traveltime reciprocity discrepancy in the Podvin & Lecomte time3d finite difference algorithm

We have found that the extensively used finite difference scheme time3d produces time fields which are 'asymmetric' in the sense that traveltimes computed to the right of the source are faster than traveltimes computed to the left. All finite difference schemes will, as they are approximations to the wave equation, to some extent fail to obey reciprocity perfectly. We show, however, that the errors in time3d may be significant—and unnecessarily large. An asymmetry in the point source initialization has been identified, and after correction time3d produces time fields with an improved reciprocity.     

3D behaviour of contamination in landfill sites using 2D resistivity/IP imaging: case studies in Portugal

Aiming at defining a valid spatial contamination model, resistivity and induced polarization (IP) measurements were used to investigate contamination plumes in the vicinity of two municipal landfills (Ovar and ílhavo). Previous geophysical surveys and underground water samples confirmed the contamination. However 2D resistivity/IP surveys enabled in obtaining a more accurate spatial model. The Ovar survey consisted of two profiles with nine Wenner soundings each; the ílhavo survey was carried out along two individual lines using a Wenner standard pseudo-section. In both situations, negative IP values were found associated with positive IP values, which can be explained mainly by 2D or 3D geometric effects caused by the presence of the conductive plumes. The data were modelled using a 2D inversion program (RES2DINV) and the resulting resistivity and chargeability distributions were displayed as pseudo-sections. The resistivity and chargeability pseudo-sections define the contamination plumes and the sedimentary structure. These case studies illustrate the advantages of 2D resistivity/IP surveys for the mapping of shape and dimension of contamination associated with landfills.     

利用三维地震在戈壁沙漠区解决复杂地质构造问题

在分析戈壁沙漠地区地球物理特征的基础上，针对性的提出了解决煤田地质构造问题的方法，包括三维地震野外施工方法的选取、施工因素的选择、资料处理手段；并以获得的实际地质成果证实三维地震在解决复杂地质构造问题的可行性。     

Seismic Hazards Assessment for Radioactive Waste Disposal Sites in Regions of Low Seismic Activity

A comparative tectonic quiescence and lack of earthquakes make the stable centres of continents attractive for siting long-term radioactive waste storage facilities. The low rates of deformation in such regions, however, make it difficult to characterize their long-term seismotectonic behaviour, leading to uncertain estimates for the very low probability hazard estimates required by society. In an attempt to overcome the deficiency of both contemporary seismicity and paleoseismic data in central Canada, we have used earthquake histories from regions with similar seismotectonic characteristics from around the world. Substituting space for time, we estimate a long-term rate per 10⁶ km² of 0.004 magnitude ≥6 earthquakes per annum, of which 33–100% might rupture to the surface.