The air flow and heat transfer in gravel embankment in permafrost areas

A comparative numerical investigation of transient temperature profile and pore-air velocities in horizontal rock block embankments are conducted using the “gravels model”, in which the embankment is composed of stones and air, and the “porous media model” respectively. As the velocities from the “gravels model” directly reflect the true flow of air and winter-time convection, in this paper it can be concluded that computational results from the “gravels model” are superior to the “porous media model”. In addition, the “gravels model” has the advantages of reflecting the effect of the dimensions and collocation of gravels upon the temperature fields. Therefore, the computation of the gravels embankment is mainly based on the gravels model. Simulation results show that in summer, a clockwise circulation of the pore-air extends throughout most of the embankment. However its motion is very weak that results in relatively straight horizontal isotherm lines. And heat transfer is mainly maintained through conduction. But in winter, the pore-air velocities are higher and multiple vortexes are formed in the embankment. Natural convection then becomes the dominant influence on the isotherm shapes within the embankment. The isotherms are complex and alternative upward and downward flowing plumes exist. The winter-time convection can further reduce the temperature of the foundation soil beneath the gravel embankment. In addition, the effects of the gravel dimensions within the embankment have been analyzed and compared in the gravels model. It shows that in winter, large stones, e.g. 200 mm, lead to stronger vortexes than those of small stones, say 60 mm. Consequently, the zone of low-temperature beneath the large-stone embankment extends deeper into the ground.     

Numerical analysis for cooling effect of open boundary ripped-rock embankment on Qinghai-Tibetan railway

At present, the Qinghai-Tibetan railway is being built, and it will pass across more than 550-km perma-frost regions. Therefore, the key to the stability of therailway embankment lies in solving the permafrost problem. Because global warming and existence of railway tend to degrade the permafrost in these re-gions[1], more difficulties and problems are induced in the construction and maintenance of railway. In the area where the mean annual air temperature is higher than a certain value, the …     

Effect of the oceanic lithosphere velocity on free convection in the asthenosphere beneath mid-ocean ridges

The paper presents results obtained in experiments on a horizontal layer heated from below in its central part and cooled from above; the layer models the oceanic asthenosphere. Flow velocity and temperature profiles are measured and the flow structure under boundary layer conditions is determined (at Rayleigh numbers Ra > 5 × 10⁵). The flow in the core of a plane horizontal layer heated laterally and cooled from above develops under conditions of a constant temperature gradient averaged over the layer thickness. The flow core is modeled by a horizontal layer with a moving upper boundary and with adiabatic bounding surfaces under conditions of a constant horizontal gradient of temperature. Exact solutions of free convection equations are found for this model in the Boussinesq approximation. Model results are compared with experimental data. Temperature and flow velocity ranges are determined for the boundary layer regime. Based on the experimental flow velocity profiles, an expression is found for the flow velocity profile in a horizontal layer with a mobile upper boundary heated laterally and cooled from above. Free convection velocity profiles are obtained for the asthenosphere beneath a mid-ocean ridge (MOR) with a mobile lithosphere. An expression is obtained for the tangential stress at the top of the asthenosphere beneath an MOR and the total friction force produced by the asthenospheric flow at the asthenosphere-lithosphere boundary is determined.     

Study on wall-type gravel drains as liquefaction countermeasure for underground structures

One of the most dramatic causes of damage to engineering structures during earthquakes has been the development of soil liquefaction beneath and around the structures. In order to dissipate the excess pore water pressures near structures, gravel drains are usually employed. In this study, the use of recycled concrete crushed stones as gravel drain materials is addressed. In order to investigate the performance of wall-type gravel drains, two series of shaking table tests were performed. The test results showed that gravel drains, when appropriate grain size distribution is considered, effectively dissipate the excess pore water pressure underneath the structure, and consequently reduce the magnitude of uplift. To supplement the laboratory tests, finite element analyses were also performed. For specified structure, ground and earthquake conditions, there is a critical width of gravel drain at which no uplift of structure will occur. The results of the model tests and the finite element analyses were then employed in developing design charts for determining the critical width of gravel drain to prevent buoyant rise of structure when the surrounding soil mass liquefies.     

Long wavelength thermal convection between non-conducting boundaries

Non-linear Rayleigh-Bénard convection in a fluid layer is considered as a model of convection in the Earth's upper mantle. Previous studies have shown that when the temperature is held fixed at one of the boundaries of the layer, convection takes place in cells of width of the order of the layer depth or less. We investigate the effects of a different thermal boundary condition, in which the flux of heat is held fixed on both layer boundaries; then if this flux is just greater than that required for the onset of convection, motion takes place on horizontal scales much greater than the layer depth. An analytical treatment of the equations, based on an expansion in the depth-to-width ratio of the cells, shows that cells of a definite horizontal scale are the fastest growing according to linearised theory, but that these cells are unstable to ones of larger wavelength than themselves. Thus the dominant wavelength lengthens with time. The results hold whether the heat flux is generated internally of comes from beneath the layer. These results produce flow patterns similar to those found when the heat flux is much greater than the critical value. The results have important consequences for the understanding of mantle convection.     

Testing density-dependent groundwater models: two-dimensional steady state unstable convection in infinite,finite and inclined porous layers

This study proposes the use of several problems of unstable steady state convection with variable fluid density in a porous layer of infinite horizontal extent as two-dimensional (2-D) test cases for density-dependent groundwater flow and solute transport simulators. Unlike existing density-dependent model benchmarks, these problems have well-defined stability criteria that are determined analytically. These analytical stability indicators can be compared with numerical model results to test the ability of a code to accurately simulate buoyancy driven flow and diffusion. The basic analytical solution is for a horizontally infinite fluid-filled porous layer in which fluid density decreases with depth. The proposed test problems include unstable convection in an infinite horizontal box, in a finite horizontal box, and in an infinite inclined box. A dimensionless Rayleigh number incorporating properties of the fluid and the porous media determines the stability of the layer in each case. Testing the ability of numerical codes to match both the critical Rayleigh number at which convection occurs and the wavelength of convection cells is an addition to the benchmark problems currently in use. The proposed test problems are modelled in 2-D using the SUTRA [SUTRA––A model for saturated–unsaturated variable-density ground-water flow with solute or energy transport. US Geological Survey Water-Resources Investigations Report, 02-4231, 2002. 250 p] density-dependent groundwater flow and solute transport code. For the case of an infinite horizontal box, SUTRA results show a distinct change from stable to unstable behaviour around the theoretical critical Rayleigh number of 4π² and the simulated wavelength of unstable convection agrees with that predicted by the analytical solution. The effects of finite layer aspect ratio and inclination on stability indicators are also tested and numerical results are in excellent agreement with theoretical stability criteria and with numerical results previously reported in traditional fluid mechanics literature.     

The long wavelength component of mantle convection

Theoretical models of mantle convection predict that the major temperature fluctuations within the mantle are confined to narrow horizontal boundary layers and vertical plumes. These fluctuations result in heterogeneities in seismic body wave velocities which could, in principle, be detected by seismic tomographic techniques. However, recent tomographic images of the spatial variations of temperature in the mantle are considerably “out of focus” in that only the longest wavelength components can be resolved. To assess this partial recovery of the total tomographic image, theoretical temperature fields have been generated with a numerical model of high Rayleigh number mantle convection and then Fourier analysed in two dimensions. Upon re-synthesizing the model temperature fields, the Fourier series expansions were truncated at various levels of resolution. The truncated expansions, containing only the long wavelength components of the model temperature fields, are compared to both the complete field and the tomographic images of the mantle. At the current level of resolution it seems unlikely that seismic tomography could distinguish between layered and whole-mantle convection. Estimates, based on current tomographic data, of long wavelength fluctuations of temperature and surface topography are predicted, in the case of whole-mantle convection, to represent approximately 10% of the total temperature variation, and approximately 50% of the total topographic relief. Thus topography at the core-mantle boundary may be more accurately inferred from seismic tomography than may the characteristic lateral temperature fluctuation which drives the convective circulation.     

冬季亚洲中强震与夏季我国东部主雨带的准同纬性

根据1954～2003年50年亚洲中强震与我国地温资料,分析了冬季震中纬度—震涡中心纬度—冬季我国东部强地热涡中心纬度—夏季我国东部主雨带纬度4者之间的相互关系,发现冬季中强震与夏季主雨带之间具有较好的准同纬性. 这是每年汛期降水预报的一个指标. 文中也讨论了应用此预报指标时的一些特例.     

Two-dimensional mantle flow beneath a rigid accreting lithosphere

This study considers two-dimensional mantle flow beneath a rigid lithosphere. The lithosphere which forms the upper boundary of a convecting region moves with a prescribed uniform horizontal velocity, and thickens with distance from the accreting plate boundary as it cools. Beneath the lithosphere, the mantle deforms viscously by diffusion creep and is heated radiogenically from within. Solutions for thermal convection beneath the lithosphere are obtained by finite-difference methods. Two important conclusions have resulted from this study: (1) convective patterns of large aspect ratio are stable beneath a rigid moving lithosphere; (2) even for a lithosphere velocity as small as 3 cm/yr. and a Rayleigh number as large as 10⁶, mantle circulation with large aspect ratio is driven dominantly by the motion of the lithosphere rather than by temperature gradients within the flow. Gravity, topography and heat flow are determined and implications for convection in the upper mantle are discussed.     

地幔对流与岩石层板块的相互耦合及影响──（Ⅱ）地幔混合对流理论及其应用

观测到的板块运动包含两种能量分布几乎相等的运动形态：极型场和环型场．纯粹的由热驱动的地幔自由对流不能预期和解释环型运动的产生．本文提出地幔混合对流理论，既考虑了热驱动的自由对流，也考虑了由板块自身激发的强迫对流．根据板块处于动力学平衡状态的观测事实，建立了相应的模型．数值结果表明，根据混合对流模型所预期的板块速度场，既能产生极型场，也能产生环型场，而且在空间分布特征及功率谱分布上与观测资料符合相当好．地幔物质的上升流动基本和洋脊对应，而下降流动和俯冲带对应．     

Ventifacts distribution in Qatar

This study attempts to investigate the distribution of ventifacts in Qatar. It is believed that ventifacts are confined to the areas within about 5 km of the Miocene or Mio-Pliocene Hofuf formations and the spreads of continental gravels derived from them. Three hypotheses were formulated: (1) Ventifacts in Qatar are confined to areas within about 5 km of the Hofuf formations and the spreads of continental gravels derived from them. The distribution of ventifacts within these areas varies according to the nature of the ground surface; (2) The most active ventifaction areas are where the continental gravels merge with the Eocene limestone because of the increase in saltation particle speed in these areas where bedrock or bare limestone is exposed; (3) The unit area ratio of ventifact to non-ventifact pebbles varies inversely with the total amount of pebbles. To test these hypotheses, nine land class categories were identified in the three major Hofuf formations. Line transects were carried out from randomly selected stations near the middle of the Hofuf formations. Along each transect systematic sampling was carried out at 200 m intervals. The data were processed using a WANG MVP 2200 computer with software developed for the project. It was found that ventifacts tend to concentrate on the outer edges of the continental gravels in areas of limestone outcrop and limestone pavement. Higher areas have big gravel counts and a low ratio of ventifacts while the low-lying plains have small gravel counts and a higher ratio of ventifacts. In certain areas ‘ventifact fields’ were found where the density of ventifacts was as high as 30 per m². Many of the ventifacts in these fields were buried beneath the surface suggesting that the ventifaction predates the present site conditions. Other high ventifact density areas were discovered where the ventifacts have collected in shallow depressions or hollows on the limestone plateaux. Water action has washed these ventifacts, a high proportion of which are dreikanters, into the hollows, where they have been partially buried in fine alluvial silts. These ‘ventifact graveyards’ are generally only a few metres wide but contain large numbers of fine specimens.     

地幔对流与岩石层板块的相互耦合及影响──（Ⅱ）地幔混合对流理论及其应用

观测到的板块运动包含两种能量分布几乎相等的运动形态：极型场和环型场．纯粹的由热驱动的地幔自由对流不能预期和解释环型运动的产生．本文提出地幔混合对流理论，既考虑了热驱动的自由对流，也考虑了由板块自身激发的强迫对流．根据板块处于动力学平衡状态的观测事实，建立了相应的模型．数值结果表明，根据混合对流模型所预期的板块速度场，既能产生极型场，也能产生环型场，而且在空间分布特征及功率谱分布上与观测资料符合相当好．地幔物质的上升流动基本和洋脊对应，而下降流动和俯冲带对应．     

Particle path characteristics at the large gravel‐bed river Danube: results from a tracer study and numerical modelling

A tracer study performed on a 3 km long reach of the Danube River in Austria is presented. Forty artificial stones of three different sizes (intermediate b‐axis: 25 mm, 40 mm, 70 mm) were produced and a coded radio acoustic transmitter was implanted. The measurement system had to be improved to be applicable to large rivers with water depths up to 12 m. The positions of the stones were observed approximately once a week, depending on hydrology, over a period of at least one year by radio‐tracking from a boat, including a 15 year flood event. Transport paths and velocities, as well as the incipient motion of bedload transport, could be monitored for the first time on a large gravel‐bed river. The particle paths were found to be mostly bankline‐parallel, even though the stones passed a 30° river bend. The median of the transverse particle displacement was found to be 4% of the longitudinal displacement. Calculations considering both transverse slope and transverse flow velocities showed transverse transport to be 6·6% of the longitudinal transport indicating that marginal lateral transport is mainly influenced by morphology. A three‐dimensional (3D) numerical model using a stochastic particle tracing approach was validated with the data, indicating that the observed positions are well reproduced by the model. Within the observation period, 74% of all stones passed the reach. With more than 1000 detections, particle transport could be characterized by a mean travel velocity of about 10 m per day (variable for the different grain sizes); single tracer stones were transported up to 1000 m during a single flood event. Size‐selective behaviour could be shown and the incipient motion of the large 70 mm gravel was detected at lower discharges than predicted by commonly used uniform bedload transport formulae. Copyright © 2012 John Wiley & Sons, Ltd.     

Particle transport in gravel-bed rivers: Revisiting passive tracer data

Data from tracer experiments were compiled and analysed in order to explore the role of geomorphological, hydrological and sedimentological constraints on fluvial gravel transport in gravel-bed rivers. A large data set from 217 transport episodes of tagged stones were compiled from 33 scientific papers. Our analyses showed that while magnitude of peak discharge is a major control on gravel transport and mobility, tracer travel distances show some scale dependence on the morphological configuration of the channel. Our results also highlight differences in the way tracers are displaced between step–pool and riffle and pool channels. The riffle–pool sequence seems to be a more efficient trap for travelling gravels than the step–pool pair. In addition, in step–pool channels there are clear differences in tracer transport between observations of first displacements after tracer seeding (unconstrained-stone conditions), and second and subsequent observations of tracer displacements (constrained-stone conditions). The comparison between tracer experiments under constrained conditions and those under unconstrained conditions also highlights the importance of bed state and structures in gravel mobility. The results of this study confirm that sediment transport in gravel-bed rivers is a complex process, whereby sedimentological and geomorphological controls are superimposed on the hydraulic forcing. © 2018 John Wiley & Sons, Ltd.     

Crust and upper mantle structure of the Ailao Shan-Red River fault zone and adjacent regions

Using arrival data of the body waves recorded by seismic stations, we reconstructed the velocity structure of the crust and upper mantle beneath the southeastern edge of the Tibetan Plateau and the northwestern continental margin of the South China Sea through a travel time tomography technique. The result revealed the apparent tectonic variation along the Ailao Shan-Red River fault zone and its adjacent regions. High velocities are observed in the upper and middle crust beneath the Ailao Shan-Red River fault zone and they reflect the character of the fast uplifting and cooling of the metamorphic belt after the ductile shearing of the fault zone, while low velocities in the lower crust and near the Moho imply a relatively active crust-mantle boundary beneath the fault zone. On the west of the fault zone, the large-scale low velocities in the uppermost mantle beneath western Yunnan prove the influence of the mantle heat flow on volcano, hot spring and magma activities, however, the upper mantle on the eas     

The velocity structure of the upper mantle in the transition zone from the East European Platform to Western Europe from seismic noise data

The surface wave tomography from ambient seismic noise recorded at stations in Western Europe (WE) and on the East European Platform (EEP) revealed the structure of the crust and upper mantle in the transitional zone from the Precambrian platform to the younger geological units in Western Europe. The Tornquist-Teisseyre Line separating these structures is clearly traced as a transition zone from the high velocities beneath EEP to the low velocities beneath WE in the crust and upper mantle, which extends to a depth of 150?C170 km. Below 200 km the relationship between the velocities beneath EEP and WE becomes the opposite. A similar relationship between the velocities in the upper mantle down to a depth of 300 km is observed on the southern boundary, where EEP borders on the northern segment of the Alpine-Himalayan seismic belt.     

Exhalation of radon and thoron: the question of the effect of thermal gradients in soil

Experimental measurements show no evidence for diurnal variation in exhalation of radon and thoron from soil due to convection induced by thermal gradients in the top few decimeters of soil as suggested by several authors. Estimates based on convective calculations indicate that even in the unlikely event of either vertical temperature gradients large enough to cause vertical instability, or horizontal gradients sufficient to cause significant convection, any effect would be too small to be detected. These same calculations suggest that it is difficult to conceive of cases involving typical thermal gradients in unfractured porous media such as soil where thermally induced convection would play an important role in transport of radon or thoron.     

The effect of horizontal resolution on gravity waves simulated by the GFDL “SKYHI” general circulation model

To examine the effects of horizontal resolution on internal gravity waves simulated by the 40-level GFDL SKYHI general circulation model, a comparison is made between the 3° and 1° resolution models during late December. The stratospheric and mesospheric zonal flows in the winter and summer extratropical regions of the 1° model are much weaker and more realistic than the corresponding zonal flows of the 3° model. The weaker flows are consistent with the stronger Eliassen-Palm flux divergence (EPFD).The increase in the magnitude of the EPFD in the winter and summer extratropical mesospheres is due mostly to the increase in the gravity wave vertical momentum flux convergence (VMFC). In the summer extratropical mesosphere, the increase in the resolvable horizontal wavenumbers accounts for most of the increase in the gravity wave VMFC. In the winter extratropical mesosphere, the increase of VMFC associated with large-scale eastward moving components also accounts for part of the increase in the gravity wave VMFC.The gravity waves in the summer and winter mesosphere of the 1° model are associated with a broader frequency-spectral distribution, resulting in a more sporadic time-distribution of their VMFC. This broadening is due not only to the increase in resolvable horizontal wavenumbers but also occurs in the large-scale components owing to wave-wave interactions. It was found that the phase velocity and frequency of resolvable small-scale gravity waves are severely underestimated by finite difference approximations.     

A simple model of convection in the terrestrial planets

Abstract

In this paper we study analytically the simplest fluid mechanical model which can mimic the convective behavior which is thought to occur in the solid mantles of the terrestrial planets. The convecting materials are polycrystalline rocks, whose creep behavior depends very strongly on temperature and probably also on pressure. As a simple model of this situation, we consider the flow of a Newtonian viscous fluid, whose viscosity depends strongly on temperature (only), and in fact has an infinite viscosity below a certain temperature, and a constant viscosity above this temperature. This model would also be directly relevant to the convection of a melt beneath its own solid phase (e.g. water below ice, though in that case there are other physical complications).

As a consequence of this assumption, there is a vigorous convection zone overlain by a stagnant lid, as also observed in analogous laboratory experiments (Nataf and Richter, 1982). The analysis is then very similar to that of Roberts (1979), but the extension to variable viscosity introduces important differences, most notably that the boundary between the lid and the convecting zone is unknown, and not horizontal. The resulting buoyancy induced stresses near this boundary are much larger than the stresses produced by buoyancy in the side-wall plumes, and mean that the dynamics of this region, and hence also the heat flux, are independent of the rest of the cell. We give a first order approximation for the Nusselt number-Rayleigh number relationship.     

利用Rayleigh波相速度和群速度联合反演青藏高原东北缘S波速度结构

利用青藏高原东北缘地区固定和流动地震台网2007年8月到2012年1月期间记录的远震波形,运用小波变换频时分析方法分别测定了1216和653条周期从15到140 s的台站间基阶Rayleigh相速度和群速度频散曲线.通过对上述频散进行反演,重构了青藏高原东北缘分辨率高达0.5°×0.5°的2-D相速度和群速度分布图.然后通过对所提取到的每个格网点Rayleigh波相速度和群速度频散进行联合反演,得到了研究区下方一维S波速度结构.最后通过线性插值,得到了青藏高原东北缘下方地壳上地幔三维S波结构.结果表明,印度板块向北俯冲已经达到班公-怒江缝合带附近;在柴达木盆地北部祁连山下面我们发现了亚洲板块,且其没有表现出明显的向南俯冲的迹象;在两大板块中间,我们观测到延伸到250 km深度的低速异常,该低速异常可能是地幔物质底辟上涌现象造成的.