Low-velocity Channel Research Articles

The characteristics of S-wave velocity and mineralization of the “Double Domes” structure in the eastern Tethys Himalaya

The Tethys Himalayan belt in the southern Qinghai–Tibet Plateau has been intruded by a large amount of leucogranite due to collisional orogeny. In addition, strong tectonic movements since the Cenozoic era have led to the formation of a series of dome structures accompanied by various types of mineralization. The Cuonadong Dome and Yalaxiangbo Dome, forming the “double dome” structure in the eastern part of the Tethys Himalayan Belt, are affected by different geological processes, resulting in differences in their deep structures and affecting the formation of polymetallic minerals. At present, geophysical research on the fine structure of the crust of the “double dome” structure is limited, making it difficult to fully understand the formation of different deep structures in the Tethys Himalayan dome belt. This hinders the progress of research on the genesis of dome structures and large-scale mineralization mechanisms in continental collisional environments. In this study, the ambient noise tomographic method was used to obtain the S-wave velocity structure of the upper crust of the Cuonadong Dome and the Yalaxiangbo Dome, and the following results were found: 1. there is a significant difference in the velocity structures of the two domes, with the core velocity structure of the Yalaxiangbo Dome showing an overall high velocity, extending downward for more than 9 km, while the core of the Cuonadong Dome exhibits low-velocity characteristics, with some high-velocity bodies occurring locally, which may be related to the later intrusion of leucogranite and extensional activity of the Cuona Rift. 2. There are significant differences in the S-wave velocities between the lead–zinc deposits and rare metal deposits in the study area; the lead–zinc deposits occur in basins and graben margins and have large variations in the S-wave velocity, which may be related to the involvement of basin brine in mineralization; below the tungsten–tin–beryllium deposits, the S-wave velocity exhibits high-velocity protrusions, with mineralization occurring at the front ends of the protrusions, which may be caused by crystallization differentiation of leucogranite. 3. The study area has abundant geothermal resources and obvious geothermal structural features. The low-velocity basin in the upper part of the upper crust is a heat storage layer, the low-velocity channel in the middle is a heat-conducting layer, and the lower part is a low-velocity heat source area that continuously supplies heat, forming a special geothermal structural model for the Cuonadong Dome and Yalaxiangbo Dome.

Seismic Evidence for Fluid‐Driven Pore Pressure Increase and Its Links With Induced Seismicity in the Xinfengjiang Reservoir, South China

AbstractXinfengjiang Water Reservoir (XWR), located in the northeast of Guangdong Province, South China, has hosted a large number of earthquakes since its impoundment. In order to understand potential triggering mechanisms of the reservoir‐induced seismicity that could help assess seismic hazard to the Greater Bay Area, the most populated region in China, in this study, we construct a high‐resolution shear‐wave velocity model in shallow crust at depths from surface to ∼10 km based on a high‐density short‐period seismic array recently deployed surrounding the XWR using ambient noise tomography. We image a northwest‐southeast (NW‐SE)‐trending high‐velocity belt at depths of 2–6 km and a relatively low‐velocity channel (LVC) underlying the high‐velocity belt at a depth below 6 km under the XWR. The LVC seems to indicate the existence of a NW‐SE‐trending channel for water infiltration below 6 km under the XWR. Combined with spatiotemporal variation of focal‐depth and frequency of seismicity under the XWR over the past decade, we infer that a NW‐SE‐trending preexisting fault has reactivated due to pore pressure increase via water infiltration from the dam to a recently most active region, Xichang, northwest of XWR. This scenario could possibly explain the occurrence of increasing seismicity in northwest of XWR since 2012. This study highlights the seismic evidence for pore pressure increase due to water infiltration that triggers seismicity and further offers a reference for other forms of fluid‐driven seismicity worldwide.

The shallow structure of Mars at the InSight landing site from inversion of ambient vibrations

Orbital and surface observations can shed light on the internal structure of Mars. NASA’s InSight mission allows mapping the shallow subsurface of Elysium Planitia using seismic data. In this work, we apply a classical seismological technique of inverting Rayleigh wave ellipticity curves extracted from ambient seismic vibrations to resolve, for the first time on Mars, the shallow subsurface to around 200 m depth. While our seismic velocity model is largely consistent with the expected layered subsurface consisting of a thin regolith layer above stacks of lava flows, we find a seismic low-velocity zone at about 30 to 75 m depth that we interpret as a sedimentary layer sandwiched somewhere within the underlying Hesperian and Amazonian aged basalt layers. A prominent amplitude peak observed in the seismic data at 2.4 Hz is interpreted as an Airy phase related to surface wave energy trapped in this local low-velocity channel.

Low-frequency ambient distributed acoustic sensing (DAS): case study from Perth, Australia

SUMMARYAmbient wavefield data acquired on existing (so-called ‘dark fibre’) optical fibre networks using distributed acoustic sensing (DAS) interrogators allow users to conduct a wide range of subsurface imaging and inversion experiments. In particular, recorded low-frequency (<2 Hz) surface-wave information holds the promise of providing constraints on the shear-wave velocity (VS) to depths exceeding 0.5 km. However, surface-wave analysis can be made challenging by a number of acquisition factors that affect the amplitudes of measured DAS waveforms. To illustrate these sensitivity challenges, we present a low-frequency ambient wavefield investigation using a DAS data set acquired on a crooked-line optical fibre array deployed in suburban Perth, Western Australia. We record storm-induced microseism energy generated at the nearby Indian Ocean shelf break and/or coastline in a low-frequency band (0.04−1.80 Hz) and generate high-quality virtual shot gathers (VSGs) through cross-correlation and cross-coherence interferometric analyses. The resulting VSG volumes clearly exhibit surface wave energy, though with significant along-line amplitude variations that are due to the combined effects of ambient source directivity, crooked-line acquisition geometry and the applied gauge length, fibre coupling, among other factors. We transform the observed VSGs into dispersion images using two different methods: phase shift and high-resolution linear Radon transform. These dispersion images are then used to estimate 1-D near-surface VS models using multichannel analysis of surface waves (MASW), which involves picking and inverting the estimated Rayleigh-wave dispersion curves using the particle-swarm optimization global optimization algorithm. The MASW inversion results, combined with nearby deep borehole information and 2-D elastic finite-difference modeling, show that low-frequency ambient DAS data constrain the VS model, including a low-velocity channel, to at least 0.5 km depth. Thus, this case study illustrates the potential of using DAS technology as a tool for undertaking large-scale surface wave analysis in urban geophysical and geotechnical investigations to depths exceeding 0.5 km.

Lithospheric structures of and tectonic implications for the central–east Tibetan plateau inferred from joint tomography of receiver functions and surface waves

SUMMARYWe present an updated joint tomographic method to simultaneously invert receiver function waveforms and surface wave dispersions for a 3-D S-wave velocity (Vs) model. By applying this method to observations from ∼900 seismic stations and with a priori Moho constraints from previous studies, we construct a 3-D lithospheric S-wave velocity model and crustal-thickness map for the central–east Tibetan plateau. Data misfit/fitting shows that the inverted model can fit the receiver functions and surface wave dispersions reasonably well, and checkerboard tests show the model can retrieve major structural information. The results highlight several features. Within the plateau crustal thickness is >60 km and outwith the plateau it is ∼40 km. Obvious Moho offsets and lateral variations of crustal velocities exist beneath the eastern (Longmen Shan Fault), northern (central–east Kunlun Fault) and northeastern (east Kunlun Fault) boundaries of the plateau, but with decreasing intensity. Segmented high upper-mantle velocities have varied occurrences and depth extents from south/southwest to north/northeast in the plateau. A Z-shaped upper-mantle low-velocity channel, which was taken as Tibetan lithospheric mantle, reflecting deformable material lies along the northern and eastern periphery of the Tibetan plateau, seemingly separating two large high-velocity mantle areas that, respectively, correspond to the Indian and Asian lithospheres. Other small high-velocity mantle segments overlain by the Z-shaped channel are possibly remnants of cold microplates/slabs associated with subductions/collisions prior to the Indian–Eurasian collision during the accretion of the Tibetan region. By integrating the Vs structures with known tectonic information, we derive that the Indian slab generally underlies the plateau south of the Bangong–Nujiang suture in central Tibet and the Jinsha River suture in eastern Tibet and west of the Lanchangjiang suture in southeastern Tibet. The eastern, northern, northeastern and southeastern boundaries of the Tibetan plateau have undergone deformation with decreasing intensity. The weakly resisting northeast and southeast margins, bounded by a wider softer channel of uppermost mantle material, are two potential regions for plateau expansion in the future.

A new crustal shear-velocity model in Southwest China from joint seismological inversion and its implications for regional crustal dynamics

SUMMARY Southwest (SW) China is located in a transition site from the active Tibetan Plateau to the stable Yangtze craton, which has complicated tectonic deformation and severe seismic hazards. We combine data from ambient noise, teleseismic body and surface waves, and petroleum wells to better constrain the crustal shear-velocity structure in SW China. We jointly invert the Rayleigh wave dispersion (5–40 s period), Rayleigh wave ZH ratio (20–60 s period), and P-wave receiver function for 114 permanent stations with a stepwise linearized joint inversion method. Compared to previous tomography results, we observe higher shear velocity in the sedimentary rocks within the Sichuan Basin, which is consistent with sonic logging measurements. Our model reveals widespread low-velocity zones in the mid-lower crust, and their boundaries correlate well with major fault systems. Between two main mid-crustal low-velocity channels, a prominent high-velocity region surrounded by earthquakes is observed in the inner zone of the Emeishan large igneous province (ELIP) and around the Anninghe-Zemuhe fault zone. These observations are comparable to regional tomography results using very dense arrays. Based on the results, we suggest that mid-lower crustal ductile flow and upper-crustal rigid fault movement play equally important roles in controlling the regional deformation styles and earthquake distribution in SW China. Our results also resolve thick crust-mantle transition zones beneath the eastern Tibetan Plateau and the inner zone of the ELIP due to ‘top-down’ and ‘bottom-up’ crust-mantle interactions, respectively. Our new model can serve as a reference crustal model of future high resolution model construction in SW China.

Win, win, win: Low cost baffle fish pass provides improved passage efficiency, reduced passage time and broadened passage flows over a low-head weir

The number of low-head barriers to fish migration far outweighs the number of large magnitude barriers and thus the cumulative negative impact on fish communities could also be far greater. Removal of man-made obstructions to fish migration is the most beneficial mitigation measure for reconnecting fragmented rivers but is not always possible and thus fish passes must be constructed. Given the large number of low-head barriers, cheap but effective fish passes must be identified. This study measured passage of brown trout (Salmo trutta L.) at a low-head gauging weir on Eshton Beck, England, before and after low cost baffle (LCB) fish pass construction using passive integrated transponder (PIT) telemetry. The LCB fish pass significantly improved overall passage efficiency from a maximum of 64% to 91%. There was a significant decrease in delay at the obstruction after the LCB fish pass was constructed and fish passed on a greater range of flows (0.08–5.39 m3 s−1) in comparison to before (0.56–1.92 m3 s−1). Fish ascended the fish pass through the low velocity channel (gaps in the baffles) as well as over the baffles, though a higher proportion were detected ascending over baffles at higher flows. It was therefore concluded that similar low-head structures should incorporate this style of fish pass to improve longitudinal connectivity for brown trout and other species with similar passage capabilities.

Water-rich sublithospheric melt channel in the equatorial Atlantic Ocean

The lithosphere–asthenosphere boundary is the most extensive boundary on Earth, separating the mobile plate above from the convecting mantle below, but its nature remains a matter of debate. Using an ultra-deep seismic reflection technique, here we show a systematic seismic image of two deep reflectors that we interpret as the upper and lower limits of the lithosphere–asthenosphere boundary beneath a 40–70-million-year-old oceanic lithosphere in the Atlantic Ocean. These two reflections correspond to 1,260 °C and 1,355 °C isotherms and bound a low-velocity channel, suggesting that the lithosphere–asthenosphere boundary is thermally controlled. We observe a clear age dependency of this sublithospheric channel: its depth increases with age from 72 km where it is 40-Myr-old to 88 km where it is 70-Myr-old, whereas its thickness decreases with age from 18 km to 12 km. We suggest that partial melting, facilitated by water, is the main mechanism responsible for the low-velocity channel. The required water concentration for melting increases with age; nevertheless, its corresponding total mass remains relatively constant, suggesting that most of the volatiles in the oceanic sublithospheric channel originate from a horizontal flux near the ridge axis.

Three-dimensional S-velocity structure of the crust in the southeast margin of the Tibetan plateau and geodynamic implications

The lower crustal flow model is one of several competing models considered to interpret the growth and expansion of southeastern Tibet. However, the dynamic processes involved in the evolution and deformation of the Tibetan plateau remain poorly understood due to the absence of reliable geophysical observations. In this study, 159 earthquakes of magnitude Ms≥6.2 were selected, which were recorded by 50 permanent broadband stations deployed in southwest China, and 1873 pairs of P-wave receiver functions with high signal-to-noise ratio were isolated. On the basis of a linearized inversion algorithm, a two-step inversion procedure was implemented that not only reduced the dependence of the inversion results on the initial models, but also provided a statistical estimation of the solution. Thus, we obtained an accurate 3D image of the S-wave velocity structure of the crust and uppermost mantle in the southeast margin of the Tibetan plateau. Our results reveal that an extensive intracrustal low-velocity zone spreads beneath southwest China, and further suggest that a lower crustal flow coming from eastern Tibet is blocked by the Jinshajiang-Red River fault to the west and the Xiaojiang fault to the east and extends largely by the Sichuan-Yunnan diamond-shaped block. This crustal flow along the southeastern margin of Tibet does not appear to be restricted to two narrow low-velocity channels as was advanced, but rather reflects the southeastward extrusion of crustal material as a way of tectonic escape.

The profile control and displacement mechanism of continuous and discontinuous phase flooding agent

ABSTRACTIn order to further explore the profile control and displacement mechanism of continuous and discontinuous phase flooding agent, the concentration distribution mathematical model of microsphere dispersion system in different branch channels is established, and its particle phase separation is simulated by using microfluidic technology. On this basis, in order to research the migration and plugging characteristics of microspheres, the visualization experiments on micro oil displacement mechanism of polymer solution and microspheres are carried out. And the experiment on the injection, migration, and plugging performance of microspheres in the multi-point core is performed. Results indicate that microspheres are in the axis of the channel due to the effect of fluid shear stress, and preferentially enter the large channel with low resistance and high flow velocity, which results in no particle or few particles in the small aperture and low flow velocity channel. The microspheres have better migration and retention capacity in the core and their migration shows the characteristic of “fluctuating pressure change”. Compared with polymer solution, the alleviation of “entry profile inversion” and the better migration and plugging performance of microsphere dispersion system can realize deep fluid diversion and expand sweep volume.

Low velocity crustal flow and crust–mantle coupling mechanism in Yunnan, SE Tibet, revealed by 3D S-wave velocity and azimuthal anisotropy

We used teleseismic data recorded by a permanent seismic network in Yunnan, SE Tibet, and measured the interstation Rayleigh wave phase velocity between 10 and 60s. A two-step inversion scheme was used to invert for the 3D S-wave velocity and azimuthal anisotropy structure of 10–110km. The results show that there are two low velocity channels between depths of 20–30km in Yunnan and that the fast axes are sub-parallel to the strikes of the low velocity channels, which supports the crustal flow model. The azimuthal anisotropy pattern is quite complicated and reveals a complex crust–mantle coupling mechanism in Yunnan. The N–S trending Lüzhijiang Fault separates the Dianzhong Block into two parts. In the western Dianzhong Block, the fast axis of the S-wave changes with depth, which indicates that the crust and the lithospheric mantle are decoupled. In the eastern Dianzhong Block and the western Yangtze Craton, the crust and the lithospheric mantle may be decoupled because of crustal flow, despite a coherent S-wave fast axis at depths of 10–110km. In addition, the difference between the S-wave fast axis in the lithosphere and the SKS splitting measurement suggests that the lithosphere and the upper mantle are decoupled there. In the Baoshan Block, the stratified anisotropic pattern suggests that the crust and the upper mantle are decoupled.

Two crustal low-velocity channels beneath SE Tibet revealed by joint inversion of Rayleigh wave dispersion and receiver functions

Competing geodynamic models, such as rigid-block extrusion, continuous deformation, and the mid-lower crustal flow, have been proposed to describe the growth and expansion of eastern Tibet. However, the dynamic processes responsible for plateau evolution and deformation remain poorly understood partly due to resolution limitations of previous models of lithospheric structure. On the basis of joint inversion of Rayleigh wave dispersion and receiver functions using data from a newly deployed seismic array, we have obtained a high-resolution 3D image that reveals the distribution of low-velocity zones (LVZs) with unprecedented clarity. The prominent feature of our model is two low-velocity channels that bound major strike-slip faults in SE Tibet and wrap around the Eastern Himalaya Syntaxis, consistent with the clockwise movement of crustal material in this region. Most large earthquakes in this region occurred in the boundaries of the LVZs. We propose that ductile flow within these channels, in addition to shear motion along strike-slip faults, played a significant role in accommodating intensive lithospheric deformation during the eastward expansion of Tibet in the Cenozoic.

Identification of T-Waves in the Alboran Sea

Analyses of seismograms from ~1,100 north-Moroccan earthquakes recorded at stations of the Red Sismica de Andalucia (Southern Spain) reveal the systematic presence of late phases embedded in the earthquake codas. These phases have distinctive frequency contents, similar to the P and S spectra and quite different to the frequency contents of the earthquake codas. They are best detected at near-shore stations. Their amplitudes decay significantly with distance to the shoreline. The delays with respect to the P-wave onsets of the preceding earthquakes are consistently around 85 s. Late phases are only detected for earthquakes located in a small region of about 100 × 60 km centered at 35.4°N, 4.0°W near the northern coast of Morocco. Several hypotheses could, in principle, explain the presence of these late phases in the seismograms, for example, the occurrence of low-energy aftershocks, efficient wave reflections, or Rayleigh waves generated along the source-station paths. However, we conclude that the most-likely origin of these phases corresponds to the incidence of T-waves (generated by conversion from elastic to acoustic energy in the north-Moroccan coast) in the southern coast of the Iberian Peninsula. T-waves are thought to be generated by energy trapping in low-velocity channels along long oceanic paths; in this case, we demonstrate that they can be produced in much shorter paths as well. Although T-waves have been already documented in other areas of the Mediterranean Sea, this is the first time that they have been identified in the Alboran Sea.

Seismic visibility of a deep subduction channel – insights from numerical simulation of high-frequency seismic waves emitted from intermediate depth earthquakes

Abstract. Return flow in a deep subduction channel (DSC) has been proposed to explain rapid exhumation of high pressure–low temperature metamorphic rocks, entirely based on the fossil rock record. Supported by thermo-mechanical models, the DSC is envisioned as a thin layer on top of the subducted plate reaching down to minimum depths of about 150 km. We perform numerical simulations of high-frequency seismic wave propagation (1–5 Hz) to explore potential seismological evidence for the in situ existence of a DSC. Motivated by field observations, for modeling purposes we assume a simple block-in-matrix (BIM) structure with eclogitic blocks floating in a serpentinite matrix. Homogenization calculations for BIM structures demonstrate that effective seismic velocities in such composites are lower than in the surrounding oceanic crust and mantle, with nearly constant values along the entire length of the DSC. Synthetic seismograms for receivers at the surface computed for intermediate depth earthquakes in the subducted oceanic crust for models with and without DSC turn out to be markedly influenced by its presence or absence. While for both models P and S waveforms are dominated by delayed high-amplitude guided waves, models with DSC exhibit a very different pattern of seismic arrivals compared to models without DSC. The main reason for the difference is the greater length and width of the low-velocity channel when a DSC is present. Seismic velocity heterogeneity within the DSC or oceanic crust is of minor importance. The characteristic patterns allow for definition of typical signatures by which models with and without DSC may be discriminated. The signatures stably recur in slightly modified form for earthquakes at different depths inside subducted oceanic crust. Available seismological data from intermediate depth earthquakes recorded in the forearc of the Hellenic subduction zone exhibit similar multi-arrival waveforms as observed in the synthetic seismograms for models with DSC. According to our results, observation of intermediate depth earthquakes along a profile across the forearc may allow to test the hypothesis of a DSC and to identify situations where such processes could be active today.

Nature of the Vrancea seismic zone (Eastern Carpathians) – New constraints from dispersion of first-arriving P-waves

The Vrancea region of the southeastern Carpathians is one of the most active seismic zones in Europe and it is well-known for its strong intermediate depth earthquakes. Seismic tomography had revealed a high-velocity body beneath Vrancea and the Moesian platform that extends to a depth of at least 350 km and can be interpreted as descending lithosphere. The strong earthquakes occur within the northeastern part of this high-velocity body, in a very limited seismogenic volume at intermediate depth (70–180 km). Several geodynamic models have been proposed for this area. They can be split into two main categories, in terms of the nature of the high-velocity anomaly, which may (a) be associated with descending relic oceanic lithosphere beneath the bending zone of the SE-Carpathians, either attached or already detached from the continental crust; or (b) it may represent continental lithosphere that has been delaminated, after continental collision and orogenic thickening. Based on currently available information, it appears difficult to distinguish between these two types of models. In this paper we attempt to shed more light on the nature of the seismic anomaly, as well as that of the origin of the intermediate depth seismicity in the Vrancea zone, by investigating the waveform character of P-waves excited by local earthquakes beneath this area, and in particular the dependence of group arrival times on frequency. We present observations of such a dispersion from stations situated at the bending zone of the SE-Carpathians. On the other hand, signals from the same earthquakes, but observed at reference stations outside of the anomalous zone do not show that frequency dependence. A natural explanation for these observations is that it is caused by the presence of a low-velocity channel at the top of the seismic anomaly, which is too thin to be resolved by classical seismic tomographic techniques. Similar observations of dispersed first-arriving P-waves have been made above subduction zones around the world, in which low-velocity layers with a thickness of several kilometers are known to exist. This suggests that a tabular slab of subducted oceanic crust is present within the seismic anomaly under the Vrancea region, and that the anomaly consists of subducted oceanic lithosphere rather than continental lithosphere, at least at depths shallower than the seismically active zone.

Lithospheric expression of geological units in central and eastern North America from full waveform tomography

The EarthScope TA deployment has provided dense array coverage throughout the continental US and with it, the opportunity for high resolution 3D seismic velocity imaging of both lithosphere and asthenosphere in the continent. Building upon our previous long-period waveform tomographic modeling in North America, we present a higher resolution 3D isotropic and radially anisotropic shear wave velocity model of the North American lithospheric mantle, constructed tomographically using the spectral element method for wavefield computations and waveform data down to 40 s period. The new model exhibits pronounced spatial correlation between lateral variations in seismic velocity and anisotropy and major tectonic units as defined from surface geology. In the center of the continent, the North American craton exhibits uniformly thick lithosphere down to 200–250 km, while major tectonic sutures of Proterozoic age visible in the surface geology extend down to 100–150 km as relatively narrow zones of distinct radial anisotropy, with Vsv>Vsh. Notably, the upper mantle low velocity zone is present everywhere under the craton between 200 and 300 km depth. East of the continental rift margin, the lithosphere is broken up into a series of large, somewhat thinner (150 km) high velocity blocks, which extend laterally 200–300 km offshore into the Atlantic Ocean. Between the craton and these deep-rooted blocks, we find a prominent narrow band of low velocities that roughly follows the southern and eastern Laurentia rift margin and extends into New England. We suggest that the lithosphere along this band of low velocities may be thinned due to the combined effects of repeated rifting processes and northward extension of the hotspot related Bermuda low-velocity channel across the New England region. We propose that the deep rooted high velocity blocks east of the Laurentia margin represent the Proterozoic Gondwanian terranes of pan-African affinity, which were captured during the Rodinia formation but left behind after the opening of the Atlantic Ocean. Our results suggest that recurring episodes of tectonic events that are well exposed at the surface also leave persistent scars in the continental lithosphere mantle, marked by isotropic and radially anisotropic velocity anomalies that reach as deep as 100–150 km.

Waveguide effects within the Philippine Sea slab beneath southwest Japan inferred from guided SP converted waves

SUMMARY High-sensitivity, dense seismic observations and numerical simulations revealed a distinct guided SP converted wave of oceanic mantle earthquakes that occurred within the subducting Philippine Sea (PHS) plate and a waveguide effect of the PHS Plate beneath southwest Japan. Based on seismograms, we identify the oceanic mantle earthquakes of the double seismic zone and find a new type of later phase for these earthquakes. The primary features of these later phases are as follows: (1) these phases are observed at the northeastern inland stations of the convex slab axis, (2) the vertical and radial components are dominant and (3) the apparent velocity is approximately 8.0 km s –1 . To examine the origin of these later phases, we simulate a2 -DP-SV seismic wavefield considering the PHS slab geometry. Based on the numerical models, we conclude that the later phases are guided SP converted waves trapped within the oceanic crust, which acts as low-velocity layer. A series of numerical simulations indicates that two factors are necessary to reproduce the observed SP guided waves: (i) a low-velocity channel must exist between stations and sources and (ii) the oceanic crust must only be in contact with the lower crust near the stations.

High resolution Rayleigh wave phase velocity tomography in northern North China

SUMMARY This study presents the Rayleigh wave phase velocity tomographic results in northern North China. The data are from 190 broad-band and 10 very broad-band stations of the North China Seismic Array and 50 permanent stations nearby. All available teleseismic vertical component time-series are used to extract the phase velocity dispersion curves of the fundamental mode Rayleigh wave by interstation method. Tomographic maps are obtained at periods of 10, 15, 25 and 60 s with a grid spacing of 0.25°× 0.25°. The short-period phase velocity maps show good correlation with the geological and tectonic features. To be specific, lower velocities correspond to North China Basin and depression area whereas higher velocities are associated with Taihangshan and Yanshan uplifts. At 25 s, there are obvious low-velocity anomalies in Jizhong depression and Beijing–Tianjin–Tangshan region, indicating that ascendant low velocity channel may be formed beneath these areas and induce the velocity difference in the upper crust. The phase velocity map at 60 s reflects the upper-mantle information in the study area. High-velocity anomalies are observed at Yanshan blocks north to Zhangjiakou–Bohai seismic belt, suggesting that the materials are stable beneath these areas or the asthenosphere is at deeper location. Low-velocity anomalies are mainly south to the seismic belt, implying the asthenosphere is shallower and the materials are transformed by the open stretching rift trending NNE, resulting in many NNE-directed fault belts. These structural differences at depth may be controlled by the fault activity and strong tectonic movements.

Study on the Critical Heat Flux Mechanism Model under the Condition of Low Velocity and Subcooled Boiling in the Narrow Channel

Based on the bubble dynamics, the generation and separation of bubble plays a leading role in the heat transfer of boiling condition. Before the bubble separation, micro liquid layer below bubbles is evaporated to dryness and heat transfer deterioration is the ultimate cause of boiling crisis, proposed a critical heat flux mechanism in the conditions of low velocity and narrow channel. A critical heat flux mechanism model is established. By calculation, in the low mass flow range the mechanism model has a high precision to predict CHF in a larger pressure and export dry degree ranges. Using the model the relationship of parameters is studied.

Finding a realistic velocity distribution based on iterating forward modelling and tomographic inversion

Tomography is like a photograph that was taken by a camera with blurred and defective lenses that deform the shapes and colours of objects. Reporting quantitative parameters derived from tomographic inversion is not always adequate because tomographic results are often strongly biased. To quantify the results of tomographic inversion, we propose a forward modelling and tomographic inversion (FM&TI) approach that aims to find a more realistic solution than conventional tomographic inversion. The FM&TI scheme is based on the assumption that if two tomograms derived from the inversion of observed and synthetic data are identical, the synthetic structure may appear to be closer to the real unknown structure in the ground than the inversion result. However, the manual design of the synthetic velocity distribution is usually time-consuming and ambiguous. In this study, we propose an approach that automatically searches for a probabilistic model. In this approach, a synthetic model is iteratively updated while taking into account the bias of the model in previous stages of the FM&TI performance. Here, we present an example of synthetic modelling and real data processing for an active source refraction data set corresponding to a marine profile across the subduction zone in Chile at about 32°S latitude. A key feature of the model is a low-velocity channel above the subducted oceanic crust, which was defined in the synthetic model and expected in the real case. The conventional first arrival traveltime tomography was barely able to resolve this channel. However, after several iterations of the FM&TI modelling, we succeeded in reconstructing this channel clearly. In the paper, we briefly discuss the nature of this low-velocity subduction channel, and we compare the results with other studies.