Non‐uniqueness with refraction inversion – the Mt Bulga shear zone

The tau‐p inversion algorithm is widely employed to generate starting models with many computer programs that implement refraction tomography. However, this algorithm can frequently fail to detect even major lateral variations in seismic velocities, such as a 50 m wide shear zone, which is the subject of this study. By contrast, the shear zone is successfully defined with the inversion algorithms of the generalized reciprocal method. The shear zone is confirmed with a 2D analysis of the head wave amplitudes, a spectral analysis of the refraction convolution section and with numerous closely spaced orthogonal seismic profiles recorded for a later 3D refraction investigation. Further improvements in resolution, which facilitate the recognition of additional zones with moderate reductions in seismic velocity, are achieved with a novel application of the Hilbert transform to the refractor velocity analysis algorithm. However, the improved resolution also requires the use of a lower average vertical seismic velocity, which accommodates a velocity reversal in the weathering. The lower seismic velocity is derived with the generalized reciprocal method, whereas most refraction tomography programs assume vertical velocity gradients as the default. Although all of the tomograms are consistent with the traveltime data, the resolution of each tomogram is comparable only with that of the starting model. Therefore, it is essential to employ inversion algorithms that can generate detailed starting models, where detailed lateral resolution is the objective. Non‐uniqueness can often be readily resolved with head wave amplitudes, attribute processing of the refraction convolution section and additional seismic traverses, prior to the acquisition of any borehole data. It is concluded that, unless specific measures are taken to address non‐uniqueness, the production of a single refraction tomogram that fits the traveltime data to sufficient accuracy does not necessarily demonstrate that the result is either correct, or even the most probable.     

QLg tomography in Gujarat,Western India

We propose a novel Lg attenuation tomography model (Q_Lg tomography) for the state of Gujarat, Western India, using earthquake data recorded by the Gujarat Seismic Network, operated by the Institute of Seismological Research in Gandhinagar. The waveform dataset consist of 400 3-component recordings, produced by 60 earthquakes with magnitude (M_L) spanning from 3.6 to 5.1, recorded at 60 seismic stations having epicentral distances spanning between 200 and 500 km. Spectral amplitude decays for Lg wave displacement were obtained by generalized inversion at 17 frequencies spanning between 0.9 and 9 Hz. Lg wave propagation efficiency was measured by Lg/Pn spectral ratio categorizing as efficient ratio ≥6 for 86%, intermediate ratio of 3–6 for 10% and inefficient ratio <3 for 4% paths of total 400 ray paths. The earthquake size and quality of waveform recorded at dense network found sufficient to resolve lateral variation of Q_Lg in Gujarat.Average power-law attenuation relationship obtained for Gujarat as Q_Lg(f) = 234f^0.64, which corresponds to high attenuation in comparison to peninsular India shield region and other several regions around the world. Q_Lg tomography resolves the highly attenuating crust of extremely fractured Saurashtra region and tectonically active Kachchh region. The Gujarat average attenuation is also lying in between them. The low attenuation in Cambay and Narmada rift basins and extremely low attenuation in patch of Surendranagar area is identified. This study is the first attempt and can be utilized as pivotal criteria for scenario hazard assessment, as maximum hazard has been reported in highly attenuating tectonically active Kachchh region and in low attenuating Cambay, Narmada and Surendranagar regions. The site and source terms are also obtained along with the Q_Lg inversion. The estimated site responses are comparable with observed local geological condition and agree with the previously reported site amplifications at the same sites. The source terms are comparable with local magnitude estimated from Network. The M_w (Lg) is nearly equivalent to M_L (GSN) and the slight differences are noted for larger magnitude events.     

井间高分辨率纵横波层析成像研究井间油藏

采用视速度-偏振井间波场分离法分离得到纵波和横波波场,使用跨孔地震走时层析成像级联方法计算井间纵波和横波速度结构以及泊松比空间分布.综合纵横波速度、泊松比和测井参数,分析两井间岩性和储层特性.数值模型试验表明该方法垂直分辨率达到主频波长的1/4,并能准确分辨倾斜断层和垂直断层.处理相同纵波资料结果表明本文方法分辨率远远高于以往的井间走时层析成像方法,比井间声波层析成像分辨率更高.使用该方法处理垦71区井间资料,得到同台同源纵波和横波速度结构和泊松比分布.泊松比的低/高变化与测井参数指示的砂/泥岩层基本吻合.根据泊松比分布,参考岩石地震物理学测试数据,可以区分井间介质中的泥岩层、砂岩层,以及砂岩层中的饱水和饱油区,确定储层的连通性,圈定井间的剩余油藏.     

Structural features of the shallow plumbing system of Vulcano Island Italy

In this study, we integrate information gathered from surface geology and tectonics with the results of a shallow (0–2 km b.s.l.) seismic tomography of Vulcano Island (Italy), obtained from the analysis of local earthquakes. The observed low Vp regions correspond to caldera filling products, mainly consisting of pyroclastics, tuffs, lava flows and hyaloclastites. High-velocity anomalies represent intrusive bodies. The striking correspondence between the stratigraphy from deep wells and the calculated velocity structure allows us to reconstruct the geometry and distribution of a main intrusion and to recognize some intra-caldera depressions. The shape and location of the high and low Vp anomalies are consistent with NW–SE and N–S strikes. Eruptive centres younger than 42 kyr, as well as the structural depressions of Vulcano and of the neighbour Lipari Island, align along a N–S direction. The combined interpretation of the available structural data and of the results from the tomography suggests that magmatic reservoirs of Vulcano at shallow depth (>0.5 km) align along a NW–SE strike but their shape is controlled by N–S striking normal faults and/or cracks that accommodate the right-lateral movements of the NW–SE strike-slip fault system.Editorial responsibility: T. Druitt     

Tomographic image of crustal structures across the Chelungpu fault: Is the seismogenic layer structure- or depth-dependent?

In 1999, a large earthquake (M_w = 7.6) occurred along the Chelungpu fault in the fold-and-thrust belt of western Taiwan. To shed more light on the subsurface structures and the seismogenic layers, three-dimensional velocity structures were inverted by using the travel times of both P- and S-waves from 2391 aftershocks recorded by the Central Weather Bureau during the 15 months that followed. From tomography, a typical image of the large-scale thrusting structures in the upper crust across the Chelungpu fault was obtained. In general, high velocities beneath the Western Foothills and Central Ranges are separated from low velocities beneath the Coastal Plain by an east-dipping boundary that is roughly consistent with the Chelungpu fault on the surface. The contrast in velocity on either side of the Chelungpu fault is indicative of about a 7- to 9-km vertical offset in the upper crust. The relocated hypocenter for the Chi-Chi earthquake shifts by 2.2 km toward the northwest, and its focal depth decreases by 0.7 km. A plot of focal depths versus rock velocities where the aftershocks occurred shows earthquakes are more inclined to occur in rock with a velocity of around 5.6 km/s. This strongly suggests the seismogenic layer in the fold-and-thrust belt of Taiwan is more structure-dependent than depth-dependent.     

Laboratory validation of lattice Boltzmann method for modeling pore-scale flow in granular materials

Characteristics of fluid flow through various engineering structures, such as granular filters and asphalt pavements, influence their design life. Numerical simulation of fluid flow is useful for evaluating the hydraulic characteristics of these materials. Among various techniques, the lattice Boltzmann (LB) method is widely accepted due to the ease of implementing boundary conditions and the numerical stability in a wide variety of flow conditions. It has proven to be extremely efficient in the simulation of fluid flow through the complex geometries of granular materials. In this study, two-dimensional and three-dimensional LB models were developed to represent pore-scale monophasic Newtonian incompressible fluid flow in granular materials. Three-dimensional geometries of compacted aggregates and asphalt specimens were generated from X-ray Computed Tomography technique and used as input for the LB model. The accuracy of the models was verified by comparing the results with analytical solutions of simple geometries and hydraulic conductivity measurements on the compacted aggregates and hot mix asphalt specimens. The results of LB simulations were in excellent agreement with those obtained from analytical calculations and laboratory measurements.     

An integrated multi-scale 3D seismic model of the Archaean Yilgarn Craton, Australia

The collection of a range of different seismic data types has greatly improved our understanding of the crustal architecture of Australia's Archaean Yilgarn Craton over the last few years. These seismic data include broadband seismic studies, seismic receiver functions, wide-angle recordings and mine-scale to deep seismic reflection transects. Each data set provides information on the three-dimensional (3D) tectonic model of the Yilgarn Craton from the craton scale through to the mine scale. This paper demonstrates that the integration and rationalisation of these different seismic data sets into a multi-scale 3D geological/seismic model, that can be visualised at once in a single software package, and incorporating all available data sets, significantly enhances this understanding. This enhanced understanding occurred because the integrated 3D model allowed easy and accurate comparison of one result against another, and facilitated the integrated questioning and interrogation across scales and seismic method. As a result, there are feedback questions regarding understanding of the individual seismic data sets themselves, as well as the Yilgarn Craton as a whole.The methodology used, including all the data sets in the model range, had to allow for the wide range of data sets, frequencies and seismic modes. At the craton scale, P-wave, S-wave and surface wave variations constrained the 3D lithospheric velocity model, revealing noticeable large-scale velocity variations within and across the craton. An interesting feature of the data, easily identified in 3D, is the presence of a fast S-wave velocity anomaly (> 4.8 km s^− 1) within the upper mantle. This velocity anomaly dips east and has a series of step-down offsets that coincide approximately with province and terrane boundaries of the Yilgarn Craton.One-dimensional receiver function profiles show variations in their crustal velocity across the craton. These crustal velocity variations are consistent with the larger-scale geological subdivision of the craton, and provide characteristic profiles for provinces and terranes. The receiver function results and the deep seismic reflection data both agree on the depth to the Moho, and both indicate an increase in Moho depth to the east. The 2D seismic refraction results in the south-west of the craton provide crustal thickness information, an indication of middle and lower crustal compositions, and information regarding the broad-scale architectural framework.At the province- and terrane-scale, the deep seismic reflection data and the mine-scale seismic data provide geometric constraints on crustal architecture, in particular the orientation of the region's fault systems as well as variations in the thickness of the granite–greenstone succession. Integration of the results from wide-angle seismic refraction data coincident with the deep seismic reflection data provided additional constraints on likely upper crustal lithologies.The integrated 3D seismic model implies the dominant geodynamic process involved the development of an orogenic belt that developed with a series of contractional (folding and thrusting) events, separated by equally important extensional events. The seismic reflection data in particular suggests that extensional movement on many shear zones was more common than previously thought.The seismic reflection data suggest that the dominant mineral systems involved deeply sourced fluid flowing up crustal-penetrating shear zones. These deeply sourced fluids were further focussed into sites located above fault-breached domal regions in the upper crust.     

Caractérisation de la structure interne et de l'état hydrique de glissements argilo-marneux par tomographie géophysique : l'exemple du glissement-coulée de Super-Sauze (Alpes du Sud,France)

Geophysical methods such as seismic surveying and electrical resistivity imaging appear to be well-adapted to investigate and understand landslide mechanisms. They allow direct and non-intrusive measurement of acoustic velocity and electrical resistivity, two fundamental parameters used to define the physical properties of the reworked moving materials. Both methods were applied at the Super-Sauze site in the French South Alps, where a typical example of an intra-material earthflow-mudslide can be observed. Measurements were taken simultaneously along two orthogonal profiles: one 325 m in length is perpendicular to the axis of the earthflow; the other 235 m in length is located along the axis of the earthflow. The results show a correlation between the seismic and electrical resistivity data, confirming that the simultaneous use of both methods gives complementary information about landslide mechanisms. The seismic data provide information on fracture density variations, whereas the electrical resistivity data provide information on water content variations. To cite this article: G. Grandjean et al., C. R. Geoscience 338 (2006).     

GPS water vapour tomography: preliminary results from the ESCOMPTE field experiment

Water vapour plays a major role in atmospheric processes but remains difficult to quantify due to its high variability in time and space and the sparse set of available measurements. The GPS has proved its capacity to measure the integrated water vapour at zenith with the same accuracy as other methods. Recent studies show that it is possible to quantify the integrated water vapour in the line of sight of the GPS satellite. These observations can be used to study the 3D heterogeneity of the troposphere using tomographic techniques. We develop three-dimensional tomographic software to model the three-dimensional distribution of the tropospheric water vapour from GPS data. First, the tomographic software is validated by simulations based on the realistic ESCOMPTE GPS network configuration. Without a priori information, the absolute value of water vapour is less resolved as opposed to relative horizontal variations. During the ESCOMPTE field experiment, a dense network of 17 dual frequency GPS receivers was operated for 2 weeks within a 20×20-km area around Marseille (southern France). The network extends from sea level to the top of the Etoile chain (700 m high). Optimal results have been obtained with time windows of 30-min intervals and input data evaluation every 15 min. The optimal grid for the ESCOMTE geometrical configuration has a horizontal step size of 0.05°×0.05° and 500 m vertical step size. Second, we have compared the results of real data inversions with independent observations. Three inversions have been compared to three successive radiosonde launches and shown to be consistent. A good resolution compared to the a priori information is obtained up to heights of 3000 m. A humidity spike at 4000-m altitude remains unresolved. The reason is probably that the signal is spread homogeneously over the whole network and that such a feature is not resolvable by tomographic techniques. The results of our pure GPS inversion show a correlation with meteorological phenomena. Our measurements could be related to the land–sea breeze. Undoubtedly, tomography has some interesting potential for the water vapour cycle studies at small temporal and spatial scales.