Simulating the 1999 Capbreton canyon turbidity current with a Cellular Automata model     

T. Salles  T. Mulder  M. Gaudin  M.C. Cacas  S. Lopez  P. Cirac 《Geomorphology》2008,97(3-4):516-537

A numerical model has been developed for the simulation of turbidity currents driven by nonuniform, non cohesive sediment and flowing over a complex three dimensional submarine topography. The model is based on an alternative approach known as Cellular Automata paradigm. The model is validated by comparing a simulation with a reported field-scale event. The chosen case is a turbidity current which occurred in Capbreton Canyon and was initiated by a storm in December 1999. Using data from recent oceanographic cruises, the deposit of the event has been precisely described, which constrain values of model parameters. The model simulates the 1999 turbidity current over the actual canyon topography and related turbidite using three different types of particle. The model successfully simulates areas of erosion and deposition in the canyon. It predicts the vertical and longitudinal grain size evolution, and shows that the fining-up sequence can be deposited by several phases of deposition and erosion related to the current energetic variation during its evolution. This result could explain the presence of intrabed contacts or the frequent lack of facies in Bouma sequences.  相似文献   

Determination and analysis of long-wavelength transition zone structure using SS precursors     

C. Houser †  G. Masters  M. Flanagan  P. Shearer 《Geophysical Journal International》2008,174(1):178-194

Global mapping of 410 and 660 km discontinuity topography and transition zone thickness has proven to be a powerful tool for constraining mantle chemistry, dynamics and mineralogy. Numerous seismic and mineral physics studies suggest that the 410 km discontinuity results from the phase change of olivine to wadsleyite and the 660 km discontinuity results from the phase change of ringwoodite to perovskite and magnesiowustite. Underside reflections of the 410 and 660 km discontinuities arrive as precursors to SS . With the recent development of a semi-automated method of determining SS arrivals, we have more than tripled the Flanagan and Shearer (1998a) data set of handpicked SS waveforms. We are able to increase resolution by stacking waveforms in 5° rather than 10° radius bins as well as increasing data coverage significantly in the southern hemisphere. The resulting SS-S410S and SS-S660S times are heavily influenced by upper-mantle velocity structure. We perform a joint inversion for discontinuity topography and velocity heterogeneity as well as performing a simple velocity correction to the precursor differential times and find little difference between the two methods. The 660 km discontinuity topography and transition zone thickness are correlated with velocities in the transition zone whereas the 410 km discontinuity topography is not. In addition, the 410 km discontinuity topography is not correlated with the 660 km discontinuity topography, rather anticorrelated, as expected due to the opposite signs of the Clapeyron slopes of their respective phase changes. These results suggest that, whereas the topography of 660 km discontinuity could be dominated by thermal effects, the topography of the 410 km discontinuity is likely dominated by compositional effects. In addition, unlike previous studies which find less topography on the 410 km discontinuity than on the 660 km discontinuity, our 410 and 660 km topography have similar amplitudes.  相似文献   

基于紧凑存储的任意四边形网格有限元法地震波场数值模拟     

刘瑞合  周建科  周学锋  印兴耀  曹丹平 《CT理论与应用研究》2015,24(5):667-680

采用有限元法精确模拟复杂介质条件下的地震波场,一般采用三角单元,但在节点数相同的情况下,三角单元的计算精度不如矩形单元高,采用三角单元模拟复杂界面时,编制确定结构刚度矩阵非零元素位置的程序也较为麻烦。采用矩形单元离散含有倾斜或起伏界面的地质模型时,无法避免绕射噪声,若加密网格又会增加计算量。为此,本文基于任意四边形单元模拟声波的传播,在倾斜或起伏界面条件下,可以有效避免因离散引起的“阶梯状”界面,在不增加计算量以及内存占用的前提下,有效地消除离散绕射噪声。采用对角的集中质量矩阵代替一致质量矩阵,避免矩阵的求逆运算,从而提高显式有限元法的计算效率;对结构刚度矩阵采用紧凑存储格式,每一行需要存储的元素最多为5个,同时零元素不参与运算,既减少内存的占用,又极大地提高计算效率。   相似文献   

Climate Change Impacts on Flow,Sediment and Nutrient Export in a Great Lakes Watershed Using SWAT          下载免费PDF全文

Siddhartha Verma  Rabin Bhattarai  Nathan S. Bosch  Richard C. Cooke  Prasanta K. Kalita  Momcilo Markus 《洁净——土壤、空气、水》2015,43(11):1464-1474

  相似文献   

青藏高原东北缘玛雅雪山晚第四纪冰川发育的气候和构造耦合   总被引：2，自引：1，他引：1  

刘蓓蓓  张威  崔之久  刘亮 《冰川冻土》2015,37(3):701-710

青藏高原东北缘的玛雅雪山(海拔4 447 m)保存着确切的第四纪冰川遗迹. 野外地貌调查与光释光测年方法相结合, 确认玛雅雪山晚第四纪主要经历3次冰川作用: 第Ⅰ组冰碛时代为新冰期; 第Ⅱ组冰碛物年龄为(23.2±1.0)ka, 其上覆泥石流年龄为(2.9±0.3)~(2.3±0.1)ka, 上层土壤年龄为(3.6±0.2)ka, 对应于深海氧同位素2阶段(MIS 2)的末次冰盛期(LGM); 第Ⅲ组冰碛年龄为(42.6±1.9)~(45.7±3.0) ka, 属于末次冰期中冰阶, 对应MIS 3中期. 采用最新综合因子法计算玛雅雪山现代冰川物质平衡线为海拔4 605 m. 依据冰川地貌形态, 计算末次冰期平衡线为海拔3 800 m. 通过庄浪河阶地的拔河高度及各级阶地的年代, 以河流的下切速率代表玛雅雪山的抬升速率, 计算得到末次冰期中期以来玛雅雪山抬升了50~60 m. 利用玛雅雪山周边的达里加山和太白山冰川漂砾的¹⁰Be 数据近似代表流域侵蚀速率, 推算出玛雅雪山剥蚀速率大约为29 mm·ka^-1, 推断MIS 3以来流域的剥蚀量为1~2 m. 综合末次冰期中期以来的构造抬升量和剥蚀量, 恢复末次冰期中期时的流域高度为海拔4 200 m, 平衡线高度为海拔3 750 m. 研究结果显示: 研究区在MIS 3时, 流域平均高度已经在平衡线之上, 在流域平均高度到主峰之间冰川开始积累, 发育冰川. 结合其他环境指标综合推断, 玛雅雪山晚第四纪冰川的发育是气候和构造耦合的产物.  相似文献   

基于点稳定系数法的斜坡稳定性分析     

马建全  王念秦  张新社 《吉林大学学报(地球科学版)》2015,45(2):533-540

基于莫尔-库伦强度理论构架,界定了点稳定系数的概念,并推导其计算公式。利用Geostudio软件建立了均质斜坡模型及计算其应力分布,并在此基础上结合MATLAB软件计算斜坡模型中各点的点稳定系数,勾绘出斜坡体内不同稳定度区域,探析了斜坡稳定性,并与传统极限平衡法进行了对比。对比结果表明:对直立斜坡,两种方法的计算结果均为不稳定,但点稳定性系数法勾绘出坡脚及坡脚底部存在两处不稳定区域;对60°斜坡,点稳定系数法的计算结果表明坡脚处存在潜在不稳定区域,而极限平衡法的计算结果表明坡体处于稳定状态;对45°斜坡,两种方法的计算结果均为稳定,计算结果一致。进一步分析得到结论:点稳定系数法不需要假设或指定某一形状滑面进行斜坡稳定性评价,且可考虑应力集中对坡体稳定性的影响;极限平衡法以稳定系数表达计算结果,而点稳定系数法以不稳定区域表达计算结果。在分析了应力和岩土体力学参数因素对点稳定系数法计算结果的敏感性后发现:相对于极限平衡法,岩土体力学参数对点稳定系数法影响更为敏感,存在黏聚力界限点和内摩擦角界限点,且对均质斜坡破坏形式(局部滑动变形破坏或整体压缩变形破坏)起着非常重要的作用。  相似文献   

Comparison of 1D linear,equivalent-linear,and nonlinear site response models at six KiK-net validation sites     

《Soil Dynamics and Earthquake Engineering》2015

Vertical seismometer arrays represent a unique interaction between observed and predicted ground motions, and they are especially helpful for validating and comparing site response models. In this study, we perform comprehensive linear, equivalent-linear, and nonlinear site response analyses of 191 ground motions recorded at six validation sites in the Kiban–Kyoshin network (KiK-net) of vertical seismometer arrays in Japan. These sites, which span a range of geologic conditions, are selected because they meet the basic assumptions of one-dimensional (1D) wave propagation, and are therefore ideal for validating and calibrating 1D nonlinear soil models. We employ the equivalent-linear site response program SHAKE, the nonlinear site response program DEEPSOIL, and a nonlinear site response overlay model within the general finite element program Abaqus/Explicit. Using the results from this broad range of ground motions, we quantify the uncertainties of the alternative site response models, measure the strain levels at which the models break down, and provide general recommendations for performing site response analyses. Specifically, we find that at peak shear strains from 0.01% to 0.1%, linear site response models fail to accurately predict short-period ground motions; equivalent-linear and nonlinear models offer a significant improvement at strains beyond this level, with nonlinear models exhibiting a slight improvement over equivalent-linear models at strains greater than approximately 0.05%.  相似文献   

Evaluation of active and passive seismic earth pressures considering internal friction and cohesion     

《Soil Dynamics and Earthquake Engineering》2015

The Log-Spiral-Rankine (LSR) model, which is a generalized formulation for assessing the active and passive seismic earth pressures considering the internal friction and cohesion of backfill soil, is reviewed and improved in this study. System inconsistencies in the LSR model are identified, which result from an inaccurate assumption on the vertical normal stress field (σ_z=γz) in a general c–ϕ soil medium, and from omitting the effect of soil cohesion when solving for the stress states along the failure surface. The remedies to the said inconsistencies are presented, and local and global iteration schemes are introduced to solve the resulting highly coupled multivariate nonlinear system of equations. The modified LSR model provides a more complete and accurate solution for earth retaining systems, including the geometry of the mobilized soil body, the stress state along the failure surface, as well as the magnitude and the point of application of the resultant earth thrust.  相似文献   

Newmark sliding block model for pile-reinforced slopes under earthquake loading     

《Soil Dynamics and Earthquake Engineering》2015

Recent studies have demonstrated that the use of a discretely-spaced row of piles can be effective in reducing the deformations of slopes in earthquakes. In this paper, an approximate strain-dependant Newmark sliding-block procedure for pile-reinforced slopes has been developed, for use in analysis and design of the piling scheme, and the model is validated against centrifuge test data. The interaction of the pile within the slipping soil was idealised using a non-linear elasto-plastic (P–y) model, while the interaction within the underlying stable soil was modelled using an elastic response model in which (degraded) soil stiffness is selected for an appropriate amount of shear strain. This combined soil–pile interaction model was incorporated into the improved Newmark methodology for unreinforced slopes presented by Al-defae et al. [1], so that the final method additionally incorporates strain-dependent geometric hardening (slope re-grading). When combined with the strain-dependent pile resistance, the method is therefore applicable to analysis of both the mainshock and subsequent aftershocks acting on the deformed slope. It was observed that the single pile resistance is mobilised rapidly at the start of a strong earthquake and that this and the permanent slope deformation are therefore strongly influenced by pile stiffness properties, pile spacing and the depth of the slip surface. The model shows good agreement with the centrifuge test data in terms of the prediction of permanent deformation at the crest of the slope (important in design for selecting an appropriate pile layout/spacing i.e. S/B) and in terms of the maximum permanent bending moments induced in the piles (important for appropriate structural detailing of the piles), so long as the slip surface depth can be accurately predicted. A method for doing this, based on limit analysis, is also presented and validated.  相似文献   

CityEngine CGA支持下的苗族特色民居三维建模     

李昕娟  王加胜 《测绘通报》2017,(12):112-116

少数民族民居是中华民族文化的重要组成部分,而在现今民族大融合的背景下,传统少数民族建筑的传承和发展正呈现衰退之势,采用三维建模技术对少数民族民居进行数字化保存意义重大。苗族民居造型复杂,具有独特的吊脚楼结构,给规则建模带来困难。为此,本文以CityEngine软件为建模平台,巧妙地组合与运用挤出、拆分、平移等基础函数,编写了苗族民居底层支柱、二层主屋和歇山型屋顶的规则,完整地实现了苗族典型民居的规则建模,并实现了民居样式的动态调整与民族村寨的快速生成。研究结果表明,基于CGA规则对复杂的民族民居三维建模方法可行,为数字化保存少数民族民居提供了便捷的方式。  相似文献