The influence of local topography for wind direction on Mars: two examples of dune fields in crater basins

In this work, we perform an analysis of large dark dunes within Moreux Crater and Herschel Crater on Mars using High Resolution Imaging Science Experiment (HiRISE) and Context Camera (CTX) data sets. These data allow us to conduct a detailed analysis of dune morphology and slip faces, concluding that the studied dune fields are influenced by topographically‐controlled complex wind directions. Our morphological analysis reveals that inside Moreux Crater in particular, the topographic setting dominates the wind flow direction, leading to the development of a sand transport pathway encircling the central peak of the crater. The dune fields in Herschel Crater are also affected by winds controlled by variable topography as suggested by the presence of complex dunes and dune fields. Our analysis indicate that the studied dune systems is not the result of paleo‐wind regimes. Furthermore, we perform thermal inertia measurements using thermal emission spectrometer (TES) data, which indicate that the studied dune fields consist of medium sand 250–500 µm in diameter. Copyright © 2012 John Wiley & Sons, Ltd.     

Diffusive approximation for unsteady mud flows with backwater effect

The adoption of the Diffusive Wave (DW) instead of the Full Dynamic (FD) model in the analysis of mud flood routing within the shallow-water framework may provide a significant reduction of the computational effort, and the knowledge of the conditions in which this approximation may be employed is therefore important. In this paper, the applicability of the DW approximation of a depth-integrated Herschel–Bulkley model is investigated through linear analysis. Assuming as the initial condition a steady hypocritical decelerated flow, induced by downstream backwater, the propagation characteristics of a small perturbation predicted by the DW and FD models are compared. The results show that the spatial variation on the initial profile may preclude the application of DW model with a prescribed accuracy. Whenever the method is applicable, the rising time of the mud flood must satisfy additional constraints, whose dependence on the flow depth, along with the Froude number and the rheological parameters, is deeply analyzed and discussed.     

Boundary conditions effect on linearized mud-flow shallow model

The occurrence of roll-waves in mud-flows is investigated based on the formulation of the marginal stability threshold of a linearized onedimensional viscoplastic (shear-thinning) flow model. Since for this kind of non-Newtonian rheological models this threshold may occur in a hypocritical flow, the downstream boundary condition may have a nonnegligible effect on the spatial growth/decay of the perturbation. The paper presents the solution of the 1D linearized flow of a Herschel and Bulkley fluid in a channel of finite length, in the neighbourhood of a hypocritical base uniform flow. Both linearly stable and unstable conditions are considered. The analytical solution is found applying the Laplace transform method and obtaining the first-order analytical expressions of the upstream and downstream channel response functions in the time domain. The effects of both the yield stress and the rheological law exponent are discussed, recovering as particular cases both power-law and Bingham fluids. The theoretical achievements may be used to extend semi-empirical criteria commonly employed for predicting roll waves occurrence in clear water even to mud-flows.     

Flow behaviour and runout modelling of a complex debris flow in a clay‐shale basin

Identification of debris‐flow hazard areas necessitates the knowledge of the flow thickness and the runout distance. Both have been investigated using a numerical runout model. On the Faucon stream (South French Alps), representative of clay‐shale basins, results of various rheological tests and numerical experiments are presented and discussed. The calibration of the model was undertaken using the results of both geomorphological surveys and sedimentological analyses. Rheological tests using either a parallel‐plate rheometer, a coaxial rheometer, slump tests, and an inclined plane were carried out on several samples. Results have shown that the flow behaviour could be described by an Herschel‐Bulkley constitutive equation. The rheological responses of several natural suspensions collected from surficial deposits (sandstones, moraines, weathered black marls) were also investigated. In order to model the runout of the flow, the model BING was used. The model describes well the influence of each type of sediment on the behaviour (runout distance, deposit thickness) of the flow, although the velocities were significantly overestimated. Different risk scenarios are tested and discussed. Copyright © 2005 John Wiley & Sons, Ltd.     

Comparing single-phase,non-Newtonian approaches with experimental results: Validating flume-scale mud and debris flow in HEC-RAS

Post-fire debris flows and tailing impoundment failures destroy lives and property. These geologic hazards – and other similar processes – fall on a continuum between classic Newtonian flood analyses and geotechnical stability analyses. The US Army Corps of Engineers (USACE) is developing a non-Newtonian library (DebrisLib) that includes a suite of rheological and clastic approaches to hyper-concentrated, mudflow and debris flow dynamics. The Hydrologic Engineering Center (HEC) has implemented these non-Newtonian methods into the widely used, public-domain open-channel hydraulics and morphodynamic software, HEC-RAS (river analysis system). This work presents part of the verification and validation of these non-Newtonian approaches, applying several rheological equations to published laboratory results high-concentration flume experiments. This study tested the linear Bingham model as well as the turbulent and Bagnold quadratic terms of the O'Brien equation. HEC-RAS also includes the non-linear Herschel–Bulkley (HB) approach, which quantifies shear-thickening or shear-thinning processes. The study used these non-Newtonian models in HEC-RAS to simulate 10 of Parsons et al.’s (2001) flume experiments, which measured the snout and plug velocity of fluids with high solid concentrations (C_v = 68–74%) and a broad range of material gradations (d₅₀ = 0.05–1 mm, d₁₅ = 0.006–0.1 mm). The experiments also measured and back-calculated Bingham and HB parameters of the materials, finding HB powers between 0.45 and 1.25 (i.e. fluids that are dilatant, pseudo-plastic and visco-plastic). The rheological models incorporated into DebrisLib and implemented in HEC-RAS reproduce experimental data well for most experiments. The Bingham model generated a plug velocity root-mean-square error (RMSE) of 0.21 m/s using standard flow parameters and Parsons et al.’s calibrated parameters, a substantial improvement over the unmodified shallow water flow equations (RMSE 0.77 m/s). Experiments with strong snout effects tended to generate higher residuals, especially in the snout velocity. The RMSE associated with the O'Brien equation was larger with the Parsons et al. fit parameters, but similar (0.23 m/s) with measured parameters. The turbulent parameter was the largest (often the dominant) parameter in most O'Brien simulations, with the dispersive stress only proving significant for the coarsest material. DebrisLib had to use a modified version of the dispersive term to simulate these concentrations. Both the 2D depth-averaged shallow water equation (SWE) and diffusion wave equation were used to simulate the experiments. The best results were obtained with the SWE with horizontal mixing.     

Dependence of magnetic field parameters at the subsolar point of the magnetosphere on the interplanetary magnetic field according to the data of the THEMIS experiment

We analyze the dependence of the magnitude of the magnetic field, its three components, and the clock angle in the magnetosheath just in front of the magnetopause on the same values in the solar wind before a shock wave using the data of the THEMIS experiment. We take into account the time delay of the solar wind arrival at the subsolar point of the magnetopause. We obtain dependencies of the components of the magnetic field and the clock angle at the magnetopause on the corresponding quantities in the solar wind for different averaging intervals. We point to the events for which the direction of the magnetic field at the magnetopause is highly different from the direction of the magnetic field in the solar wind up to the sign change.     

地震前后井-含水层系统潮汐参数变化特征分析——以云南弥勒井为例

本文对弥勒井2004～2016年井水位同震响应形态及地震前后井水位潮汐参数变化关系特征进行分析。结果表明,在弥勒井水流运动方向总体是以垂向流为主、径向流为辅的前提下,天然因素(地震)和人工因素(洗井、装置改造)均可改变含水层水流运动方向。当震前潮汐因子和相位差反向变化时,能量较大的地震可以改变其水流运动方向,使震后潮汐因子和相位差呈同向变化,水流运动方向由以垂向为主变为以径向为主;能量较小的地震只是改变其含水层渗透系数,没有改变其水流运动方向,震后仍以垂向流为主。当震前潮汐因子和相位差呈同向变化时,能量较大的地震发生可以使径向含水层渗透系数增大,水流运动方向不变,震后仍以径向为主。通过分析潮汐参数变化,可以进一步得知井-含水层径、垂向流渗透系数的变化,为进一步精确计算含水层参数提供了参考。     

渤海湾盆地新生代拉张中心迁移的热-流变机制

在区域拉张应力作用下, 岩石圈拉张减薄, 减薄的中心地带成为拉张中心, 渤海湾盆地新生代期间拉张中心具有自西向东、自南向北的迁移规律.为了探究渤海湾盆地新生代期间拉张中心的迁移与岩石圈热-流变结构的关系, 本文在构造-热演化模拟的基础上, 对渤海湾盆地各坳陷新生代岩石圈热-流变结构进行了计算.结果表明: 渤海湾盆地新生代岩石圈拉张中心的迁移主要受岩石圈流变强度的非均匀性所控制.拉张中心开始位于岩石圈强度最低的济阳坳陷, 拉张及伴随的的热冷却过程使得济阳坳陷岩石圈强度增大, 并且超过了渤中坳陷岩石圈强度, 因此, 沙四段时期, 拉张中心迁移到了渤中坳陷; 沙三段时期, 拉张中心迁移到了黄骅以及辽东湾坳陷.这是由于拉张速度较慢时(对于约90 km厚度的岩石圈而言, 如果初始地壳厚度＞34 km, 拉张速度0.5 cm·a^-1), 拉张初期会造成岩石圈强度的降低, 但随着拉张的持续进行, 岩石圈整体强度反而增加.孔店组时期, 在岩石圈强度较大的冀中及临清坳陷形成拉张中心可能是叠加了其他地质因素(如前期构造、断裂活动、远程应力等)的影响.东营组时期, 在岩石圈强度较大的渤中坳陷形成拉张中心, 这可能是该时期岩石圈强度同时受到了郯庐断裂带的活动或小规模的地幔对流的影响, 导致其实际强度比模拟计算的强度更低, 因而使得应变集中在渤中以及辽东湾坳陷.

     

THE EFFECT OF HARD INTERLAYER THICKNESS ONTHE SITE SEISMIC RESPONSE

Studies on the effect of near-surface overburden soil layers on seismic motion have shown that the overburden soil layers have a significant impact on the seismic effect of the site due to the formation age, genetic type, thickness difference, structure, and dynamic characteristics of the soil layers. In this paper, the one-dimensional seismic response analysis of a nuclear power plant site containing a thick hard interlayer was conducted to discuss the influence of the hard interlayer thickness on the site seismic response, so as to provide a basis for determining the seismic motion parameters for seismic design of similar sites. Based on the engineering geological data of a nuclear power plant site, five models of one-dimensional soil-layer seismic response analysis were built, and the equivalent linear method of the one-dimensional site seismic response was applied to analyze the effect of the interlayer thickness on the peak acceleration and the acceleration response spectra of the site seismic response. The seismic response characteristics of the site and influence rules of the hard interlayer thickness are summarized as follows:1)Under different input seismic motion levels, the peak acceleration at the top of the hard interlayer was less than the input peak acceleration, and the peak acceleration at the ground surface of site was greater than the input peak acceleration. 2)Under the same input seismic motion, the ratios of the peak accelerations at the top of hard interlayer to the input peak accelerations were smaller than the ratios of the peak accelerations at the ground surface to the input peak acceleration, and these ratios first decreased and then increased gradually with the increase of the hard interlayer thickness; while for the same hard interlayer thickness, these ratios gradually decreased as the input peak acceleration increasing. 3)For the same input seismic motion, the ratios of the peak accelerations at the ground surface of site to those at the top of the hard interlayer increased gradually as the hard interlayer thickness increased; however, corresponding to different hard interlayer thicknesses, the variation characteristics of ratios which are the peak accelerations at the ground surface of site to those at the top of the hard interlayer were inconsistent with the increase of the input peak acceleration. 4)The hard interlayer had a significant influence on the short-period acceleration response spectrum and the thicker the hard interlayer was, the wider the influence frequency band would be; while for a special hard interlayer thickness, the influence frequency band is certain, and the hard interlayer had little effect on the acceleration response spectrum coordinates outside this frequency band, the longer the period is, the less the influence of the hard interlayer on the acceleration response spectrum coordinates. The seismic response characteristics of the site and influence rules of the hard interlayer thickness indicate that the hard interlayer thickness has a significant impact on the peak acceleration and the acceleration response spectra of the site seismic response, and the hard interlayer has obvious isolation effect at the seismic motion, and the increase of its thickness reduces the nonlinear effect of the site and leads to the wider influence frequency band. Meanwhile, the higher the input peak acceleration is, the stronger the nonlinear effect of the site, and it's remarkable that the soft layer overlying the hard interlayer has a significant amplification effect on the seismic motion.     

陆地下垫面与环境气流对锋面气旋发展过程影响的理想试验

地面摩擦和大尺度流场是影响锋面气旋结构演变的重要因子,本文使用WRF模式并采用湿物理方案,通过理想化试验,综合考虑陆面摩擦、气旋式扰动相对于急流位置和大尺度流场对锋面气旋结构变化的影响.结果表明：当仅考虑单一因子时,气旋式扰动位于急流南侧和辐合流流场有利于气旋形成Shapiro-Keyser（S-K）模型结构.当同时考虑地面摩擦和大尺度辐合时,气旋式扰动位于急流北侧的气旋发展整体向经典的挪威气旋模型转变;扰动位于急流南侧的的气旋发展则整体呈现S-K模型结构,此时辐合流场有利于S-K模型结构出现.当同时考虑地面摩擦和大尺度辐散时,扰动处于急流北侧的气旋呈现挪威气旋模型结构;由于气旋式扰动穿越急流和辐散流场同时有利于暖锋后弯发展以及冷暖锋距离加大和锋消,扰动处于急流南侧的气旋呈现典型S-K气旋模型结构.这个结果解释了在东亚大陆地区辐散场形势下出现的S-K模型结构气旋个例.

     

利用最小二乘拟合法进行多站地磁日变基值归算

多站日变改正技术是解决远海区大范围磁力测量地磁日变改正的关键技术,然而多站日变改正中,为了减小磁场水平差异的影响,分站日变基值必须向主站归算.基于地磁日变化的时空特点,尝试将最小二乘拟合法引入到多站日变基值归算中,并采用多站同步观测数据对方法的有效性进行了验证.结论表明:当主站和分站日变性质相似时,最小二乘拟合法可以取得良好的效果,且对同步观测时段低.但是当日变性质差异较大时,则应当采用传统的同步比对法.     

空间尘埃等离子体电磁散射特性研究

基于尘埃等离子体Mie-Debye散射模型,研究平衡态下尘埃等离子体中带电尘埃的电磁散射特性.根据输运理论计算高层大气环境中尘埃等离子体层的电磁散射场,给出电磁波垂直入射时尘埃等离子体层的反射函数,同时分析不同浓度、不同粒径分布状态下尘埃等离子体层对电磁波传播的影响.结果表明:尘埃粒子尺度对其电磁散射特性影响非常大;在对数正态分布情况下,当入射波长远大于Debye半径时,尘埃等离子体的散射主要表现为Debye散射,而入射波长远小于Debye半径时,Mie散射占主要部分;并且尘埃等离子体层对电磁波的衰减随尘埃粒子浓度以及尘埃粒子半径的增大而明显增强.     

Late stage thermal history of the Songliao Basin and its tectonic implications: Evidence from apatite fission track (AFT) analyses

Apatite Fission Track (AFT) data from the Songliao Basin indicates that the late stage tectonic movements in the Songliao Basin have zoning in space and episodes in time. The late stage tectonic movements started from the east part of the basin and migrated westward. AFT ages in the east part of the basin are older than those in the west part of the basin, suggesting that the uplift occurred earlier in the east than in the west. The denudation thickness in the east part of the basin is significantly greater than that in the centre and west. The thermal history evolved two episodes of rapid cooling and subsequent slow cooling processes. Age-depth relationship derived from the AFT data indicates a four-episode denudation history. Further Monte Carlo random simulation of the AFT data reveals the four changing points of the thermal evolution at 65 Ma, 43.5 Ma, 28 Ma and 15 Ma, respectively. The uplifting and denudation rates from different episodes of evolution are proportional to the plate convergence rate. Based on the above analyses and the regional geologic background, it is concluded that the late stage thermal events in the Songliao Basin are the far field response to the subduction of the Pacific Plate under the Eurasian Plate. The first episode of the rapid cooling probably started at the end of the Nenjiang Formation, climaxed at the end of the Cretaceous and ceased at the Late Eocene. The subsequent slow cooling lasts another 15 Ma. The first episode of the evolution is the far field response to the major episode of the Yanshan Movement and subsequent series of the tectonic reorganization, especially the directional change of the Pacific Movement and also the subduction of the Indian Plate underneath the Eurasian Plate. While the second episode of the evolution is the far field response to the extension and closure of the Sea of Japan. Extension led to the migration and converging of the mantle heat flow to the Sea of Japan and resulted in the rapid cooling of the Songliao Basin.     

曲线桥梁摩擦滑移隔振试验及隔振机理研究

针对曲线桥梁在地震中因固定支座处及桥墩底部受力过大而破坏严重的现象,提出考虑支座摩擦滑移隔震的曲线桥梁抗震设计方法,通过大型振动台模型试验,研究曲线桥梁考虑支座摩擦滑移时的隔震机理。试验结果表明:支座摩擦滑移可有效减小地震作用下曲线桥梁的桥墩加速度响应,但同时会增加桥墩与梁体间的相对位移;曲线桥梁进行隔震设计后,桥墩的应变响应明显减小,证明该设计方法可有效避免地震作用下桥墩底部产生较大的损伤;最后,探讨了考虑支座摩擦滑移时曲线桥梁的隔震机理并与震害进行了对比验证。研究成果可为曲线桥梁的减震隔设计提供新的视角。     

Variation of the intensity of the Earth’s magnetic field in Portugal in the 1st millennium BC

The study of magnetization of the ceramic material from 21 archeological monuments of Portugal (the Evora province), dated archeologically from the Bronze Age to the end of the Iron Age has been carried out. For the purpose of more detailed timing of the material from the monuments the method of ceramic age dating on the basis of its porosity has been used. In order to take into account the distorted factors in the determination of the parameters of the ancient geomagnetic field with the aim of the maximal approximation to the actual values the diagnostic features of magnetite weathering have been considered and the level of weathering of the magnetic fraction in the ceramics from archeological monuments has been determined. The data of geomagnetic field-strength variation in the time interval of the 12th century BC to the beginning of the Common Era have been obtained. The field-strength at this time interval varied in the range of 60–90 micro Tesla with the maximal values in the 9th, 8th, and the second half of the 5th to the beginning of the 4th century BC. In addition, the timing of the ceramic material from the urns of the megalithic complex Monte de Tera of the Evora province has been clarified.     

土层对地震的随机反应分析

利用随机振动理论,本文研究了非均匀土层对地震的动力反应问题,假定土层的剪切模量随深度呈指数函数形式分布,关于基岩输入地震动加速度的功率谱密度函数,考虑了两种形式：白噪声谱和有色谱,为了比较起见,文中考虑了均匀土层的情形,数值计算表明;（１）在基岩输入地震动加速度的功率谱密度为白谱的情况下,土层的最大期望反应（相对位移,绝对加速度,剪应变）大于基岩输入地震动的功率谱为有色谱的情形;（２）在土层的上部     

The structure of convection in the spherical mantle with internal heating

Thermal convection in the mantle is caused by the heat transported upwards from the core and by the heat produced by the internal radioactive sources. According to the data on the heat transfer by the mantle plumes and geochemical evidence, only 20% of the total heat of the Earth is supplied to the mantle from the core, whereas most of the heat is generated by the internal sources. Along with the models that correctly allow for the internal heat sources, there are also many publications (including monographs) on the models of mantle convection that completely ignore the internal heating or the heat flux from below. In this study, we analyze to what extent these approximations could be correct. The analytical distributions of temperature and heat flux in the case of internal heating without convection and the results of the numerical modeling for convection with different intensity are presented. It is shown that the structure of thermal convection is governed by the distribution of the heat flux in the mantle but not by the heat balance, as it is typically implicitly assumed in most works. Heat production by the internal sources causes the growth of the heat flux as a function of radius. However, in the spherical mantle of the Earth, the heat flux decreases with radius due to the geometry. It turned out that with the parameters of the present Earth, both these effects compensate each other to a considerable extent, and the resulting heat flux turns out to be nearly constant as a function of radius. Since the structure of the convective flows in the mantle is determined by the distributions of heat flux and total heat flux, in the Cartesian models of the mantle convection the effective contribution of internal heating is small, and ignoring the heat flux from the core significantly distorts the structure of the convective currents and temperature distributions in the mantle.     

风的垂直切变对滞弹性大气中静力适应过程的影响

本文基于描写滞弹性大气静力适应过程的线性方程组,从波动频散关系、气团运动规律和能量转换的角度出发,研究了水平基流及其垂直切变对该模式大气静力适应过程的影响.构造四种水平基流垂直分布模型进行比较,分别为常数型、线性切变型、反气旋切变型和气旋切变型,得到结论：（1）具有重力波性质的波动是滞弹性大气静力适应过程中扰动能量传播的方式,当垂直折射指数大于零时,基本气流及其垂直切变的存在,不仅改变了波动频率的大小,而且改变波动传播的方向;（2）在静力适应过程中气块的运动轨迹呈椭圆形,水平基流及其垂直切变影响椭圆的扁率,同时也影响扰动物理量之间的偏振关系;（3）水平基流的垂直切变是扰动能量和水平基流能量发生转换的媒介,当存在垂直向上的动量输送时,正的垂直风切变对应扰动能量的衰减,水平基流能量的增加,负的垂直风切变对应扰动能量的增加,水平基流能量的衰减;（4）不同的风的垂直切变模型,对静力适应过程的影响不同;对于水平基流呈反气旋切变型和气旋切变型,扰动发展的波动垂直结构为,急流轴上方波动等相位线自下而上向西倾斜,急流轴下方波动等相位线自下而上向东倾斜,反之亦然.

     

General inverse of Stokes, Vening-Meinesz and Molodensky formulae

The undulation of the geoid, the gravity anomaly and the deflection of the vertical are the three basic observations describing the shape and the gravity field of the earth. The Stokes’ formula that computes the undulation of the geoid using the gravity anomaly on the geoid under spherical approximate conditions was first put forward by Stokes[1]. According to Stokes’ theory, The Vening-Meinesz formula that computes the meridian and the prime vertical components of the deflection of the ve…     

水平层状介质中深埋装置系统地电阻率影响系数特征研究

讨论了地电阻率观测装置系统深埋条件下的影响系数特征. 采用这类装置的目的在于减小甚至消除来自表层产生的干扰,以提取深部可能的地震信息. 采用水平层状介质模型,应用边界积分方程法计算了不同装置埋深、不同装置参数下的影响系数及其与结构参数间的关系. 结果表明,不同地层影响系数的大小与电阻率结构、电极埋深、供电极距有密切且复杂的关系. 这表明为了有效地压制表层干扰并观测到底层变化,首先需要精细地探查测区电性结构,然后在此基础上,通过理论分析确认装置埋深及选定装置参数.