Effect of anisotropy on solute transport in degraded fen peat soils

Peat soils are heterogeneous, anisotropic porous media. Compared to mineral soils, there is still limited understanding of physical and solute transport properties of fen peat soils. In this study, we aimed to explore the effect of soil anisotropy on solute transport in degraded fen peat. Undisturbed soil cores, taken in vertical and horizontal direction, were collected from one drained and one restored fen peatland both in a comparable state of soil degradation. Saturated hydraulic conductivity (K _s) and chemical properties of peat were determined for all soil cores. Miscible displacement experiments were conducted under saturated steady state conditions using potassium bromide as a conservative tracer. The results showed that (1) the K _s in vertical direction (K _sv) was significantly higher than that in horizontal direction (K_sh), indicating that K _s of degraded fen peat behaves anisotropically; (2) pronounced preferential flow occurred in vertical direction with a higher immobile water fraction and a higher pore water velocity; (3) the 5% arrival time (a proxy for the strength of preferential flow) was affected by soil anisotropy as well as study site. A strong correlation was found between 5% arrival time and dispersivity, K _s and mobile water fraction; (4) phosphate release was observed from drained peat only. The impact of soil heterogeneity on phosphate leaching was more pronounced than soil anisotropy. The soil core with the strongest preferential flow released the highest amount of phosphate. We conclude that soil anisotropy is crucial in peatland hydrology but additional research is required to fully understand anisotropy effects on solute transport.     

天山构造带地震各向异性及动力学机制研究进展

天山构造带及邻区的深部动力学机制是地球动力学研究的热点,而地震各向异性是区域构造深部动力学机制的一个重要性质。研究表明,天山构造带上地壳各向异性结果呈现区域性分区,受到构造带与断裂走向和区域应力影响;上地幔各向异性的结果认为快波偏振方向和构造带走向基本平行,但在伊塞克湖附近、塔里木盆地和准噶尔盆地挤压区域各向异性快波方向变化复杂,垂直方向上的变化可能由区域性双层各向异性引起,但局部复杂性原因有待进一步探讨。诸多研究支持天山构造带的地壳与上地幔垂直连贯变形机制。此外,地幔柱、软流圈变形、小尺度地幔对流等概念均被用来解释天山构造带的动力学背景,表明该地区的深部动力学机制非常复杂,需要更深入的探讨。     

Bayesian Markov chain Monte Carlo inversion for anisotropy of PP- and PS-wave in weakly anisotropic and heterogeneous media

A single set of vertically aligned cracks embedded in a purely isotropic background may be considered as a long-wavelength effective transversely isotropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly anisotropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Markov chain Monte Carlo inversion method of PP- and PS-wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inversion scheme taking only the PP-wave into account.     

Determining flow patterns and emplacement dynamics from tsunami deposits with no visible sedimentary structure

In the absence of eyewitness reports or clear sedimentary structures, it can be difficult to interpret tsunami deposits or reconstruct tsunami inundation patterns. The emplacement dynamics of two historical tsunami deposits were investigated at seven transects in Okains Bay, New Zealand, using a combined geospatial, geomagnetic and sedimentological approach. The tsunami deposits are present as layers of sand and silt intercalated between soils and become finer and thinner with distance inland. The deposits are attributed to the 1960 and possibly the 1868 tsunamis, based on radiometric dating and correlation with historical records. Measurements of Magnetic Fabric (MF: Anisotropy of Magnetic Susceptibility) and particle size were used to reconstruct the evolution of flow dynamics laterally and vertically. A combination of statistical methods, including spatial autocorrelation testing, Spearman's rank order correlation, Principal Component Analysis (PCA) and K‐means cluster analysis, was applied to examine relationships between MF parameters and sediment texture, and infer depositional hydrodynamics. Flow patterns deduced from MF show the estuary channel acted as a conduit for inundation, with flow commonly aligned sub‐perpendicular to the estuary bed. MF and sediment data suggest deposition occurred from settling during laminar flow. Evidence of both uprush and backwash deposition, as well as wave reflection from infrastructure, was found. Statistical analysis of data showed significant relationships between grain size parameters and MF parameters associated with flow speed and magnetic fabric type. PCA and cluster analysis differentiated samples into two primary hydrodynamic groups: (1) samples deposited from laminar flow; and (2) samples deposited close to the limit of inundation, which includes samples deposited further inland, those affected by flow convergence, and those in the upper part of tsunami deposits. This approach has potential as a tool for reconstructing hydrodynamic conditions for palaeotsunamis and by combining spatial and statistical analyses, large‐scale investigations can be more easily performed. Copyright © 2016 John Wiley & Sons, Ltd.     

基于电流密度连续性条件的直流电阻率各向异性问题自适应有限元模拟

直流电阻率法被广泛应用在工程和环境及水文地球物理、野外采矿、地热探测等领域.地下岩石常具有层理面和裂缝等具有方向依赖性的结构,岩石电导率常常呈现各向异性特征,因此研究复杂直流电阻率各向异性问题的高精度正演算法具有迫切的理论和学术需求.本文利用面向目标的自适应有限元算法和非结构化网格相结合的方式,解决了带地形任意复杂直流电阻率各向异性问题的高精度正演这一难题.有别于前人的研究成果,本文提出了一种特别的二次虚拟场算法实现带源的任意起伏地形问题模拟;另外,本文第一次基于电流密度连续性条件构建适合直流电阻率各向异性问题的后验误差估计算法,有效地驱动面向目标有限元网格自适应加密过程.最后,通过三组电阻率各向异性模型验证本文提出算法的正确性和适应性,测试结果表明：对于任意复杂直流电阻率各向异性问题,本文提出的算法具有精度高、适应性强等特点;另外,我们还发现电流密度连续性条件可用于设计直流电阻率问题的有效后验误差估计算法.     

中国东北诺敏河火山岩石圈变形:来自剪切波分裂的证据

中国东北新生代板内火山广泛发育,其中诺敏河火山由于上地幔结构研究的匮乏,火山成因尚不明确.利用布设在诺敏河火山周围的40个流动台站所记录到的远震剪切波数据,测量得到82对各向异性参数和219个无效分裂结果.结果表明,研究区快慢波延迟时间变化范围为0.4～1.4s,平均0.78±0.21s;各向异性快波方向范围为N77°W—N18°E,绝大多数快波方向集中在N6.9°W±9.87°,平行于中生代晚期岩石圈伸展变形方向,推测由残留在岩石圈中的化石各向异性所引起.同时,在火山中心及周边部分台站,只观测到无效分裂而没有观测到有效分裂结果,可能是由于残存在岩石圈内的古老形变被上涌的热地幔物质所侵蚀.     

大地电磁相位超象限现象的各向异性模型研究——以上下结构为例

大地电磁观测数据中的相位超象限现象可以由不同的电性结构产生.本文在已实现大地电磁三维任意各向异性有限差分正演的基础上,以具有上下结构关系的三维各向异性模型为例,分析各向异性体的规模以及参数变化对阻抗相位的影响.在下部各向异性体规模明显大于上部各向异性体或表现为层状特征的情况下,上部各向异性体在两个水平方向上的尺度差异较大,可以看作准二维体时容易发生相位超象限;当上部各向异性体在两个水平方向上尺度相近表现为规则三维体时,要产生相位超象限的现象则需要各向异性体具有更高的各向异性比.在同等条件下,增加各向异性体的各向异性比更容易发生相位超象限;而各向异性走向方位角的变化将直接影响到发生相位超象限的范围.对于准二维模型引起相位超象限的条件,沿用二维模型的近似解析分析方法,进一步构建了基于各向异性体的电导率、背景电导率以及各向异性走向角的相位超象限指标函数,从而更加直观地解释在二维或准二维条件下发生相位超象限现象的模型参数特征.     

基于背景噪声成像方法研究郯庐断裂带中南段及邻区地壳速度结构与变形特征

郯庐断裂带是贯穿我国东部北北东走向的一条深大断裂,其中南段及其邻区(115°E-121°E,29.5°N-35°N)穿过了大别造山带、苏鲁造山带、长江中下游成矿带及合肥盆地.为了研究该区域地壳速度结构及变形特征,我们使用安徽省和江苏省及其周边地区105个台站(固定台站98个,流动台站7个)的垂向连续波形数据,时间范围从2014年5月到2015年7月,共计14个月.利用背景噪声互相关方法,从垂直分量互相关函数中最终提取了2590条瑞利面波相速度频散曲线,反演得到周期范围为5~30s的瑞利波方位各向异性相速度分布图,再反演每个网格点瑞利面波相速度频散得到一维层状横波速度模型,然后拼合起来组成三维横波速度模型.根据本文反演结果并综合已有资料,我们得出如下结论:(1)在北大别、蚌埠隆起、长江中下游成矿带、合肥盆地北部大桥凹陷区域存在中地壳横波高速体,可能与岩石圈和下地壳拆沉以及中生代中国东部大范围岩浆活动有直接关系,更深层原因可能与古太平洋俯冲相关;(2)苏鲁造山带南缘,垂直于嘉山响水断裂,从南向北中上地壳低速体深度变浅,这个低速体可能是高压/超高压变质岩与扬子板块接触处的破碎带,是扬子板块与华北板块接触的边界;(3)郯庐断裂合肥-嘉山段两侧以及大别造山带东缘短周期瑞利面波相速度快轴方向与郯庐断裂带走向基本一致,可能是三叠纪碰撞期与白垩纪时期的大规模左旋走滑活动的结果;(4)合肥盆地南部15~20s周期的瑞利波相速度快轴方向为北西-南东向,反应该区域中下地壳快波方向为北西-南东向,推测是大别造山带折返的痕迹;(5)郯庐断裂带的结构和地震活动性存在明显的分段性,嘉山-郯城段郯庐断裂带现今地震活动性弱,但发生过较强的古地震,推断现今郯庐断裂带宿迁段可能处于闭锁状态,从长远来看要注意该地区发生大震的可能.     

Characteristics of Seismic Wave Propagation in the Binchuan Region of Yunnan Using a Dense Seismic Array and Large Volume Airgun Shots

The Binchuan region of Yunnan is a structurally complex region with mountains, basins, and active faults. In this situation,seismic wave propagation exhibits complex characteristics due to strong heterogeneity of underground media instead of following the great-circle path. In order to obtain a high-resolution shallow crustal structure, a dense seismic array was deployed during March 21 to May 30, 2017 in this area. To better understand the complexities of seismic wave propagation in this region, we perform array-based frequency-domain beamforming analysis and single-station based polarization analysis to investigate the characteristics of seismic wave propagation, using airgun-generated P-wave signals recorded by dense array stations in this experiment. The results from these two methods both reveal similar but complex characteristics of seismic wave propagation in the Binchuan basin. The azimuth anomalies off the great-circle path are quite large with values up to 30°, which is caused by strong structural heterogeneity in the very shallow crust. Our research provide a better understanding of the complex geologic structures in this area and provide guidance for detecting concealed faults and distribution of velocity anomalies.     

采用剪切波分裂方向约束的板块绝对运动模型研究

板块绝对运动(即岩石圈板块相对于深部地幔的运动)导致软流圈深度存在较强的剪切作用,为软流圈呈现较强的地震波速各向异性提供了一种物理解释.相应地,软流圈地震各向异性的实测数据为反演板块绝对运动提供了一种定量的约束.本文利用前人发表的由474个剪切波分裂数据组成的全球软流圈地震各向异性方向数据集,结合板块相对运动模型MORVEL,通过加权最小二乘法反演板块绝对运动.计算结果表明,由实测数据约束的板块绝对运动模型不能将观测数据拟合到原始数据测量误差要求的统计水平上,反映出以地震各向异性方向指示板块绝对运动方向存在原始数据测量误差之外的系统性误差.该误差或可归因于板块绝对运动控制软流圈地震各向异性的物理机制复杂性以及小尺度地幔流动的区域复杂性.因此,地震各向异性方向只能从统计平均的意义上约束板块绝对运动方向,两者间的差别除了测量误差外还应包含复杂物理成因的模型误差.在原始数据测量误差的基础上增加20°的模型误差并且剔除一个离群数据后,本文得到了最优拟合剪切波分裂数据集的SKS473模型,并与由热点数据反演得到的板块绝对运动模型进行对比.由于约束板块绝对运动的地震各向异性数据和热点数据都存在误差较大、地理分布不均的局限性,结合两类数据的联合反演或可成为未来建立更高精度板块绝对运动模型的有效途径.