Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.     

对青藏东北缘现今块体划分、运动及变形的初步研究

利用2维非连续变形分析方法（DDA），以位移代替围压作为边界约束力，研究青藏东北缘现今块体划分及其运动变形。根据该地区地质构造及地震活动，以GPS点测量位移作为模拟结果约束点，得出了较合理的块体划分模型和随时间演化的主应变分布图，并把应变高值区与近几年来发生的5级以上地震作对比，得出了研究区内地震危险性可能较大的区域。另外，对模拟的甘青块体与阿拉善块体的边缘带断裂左旋运动做了大概计算。     

低轨道人造卫星(CHAMP、GRACE、GOCE)与高精度地球重力场——卫星重力大地测量的最新发展及其对地球科学的重大影响

评述了卫星重力大地测量的最新发展及其对地球科学的重大影响。为了更好地理解地球内部物理构造与海洋动力学，以及大陆，冰川和海洋的相互作用，改善现有地球重力场模型（包括精度和空间解析度）是非常重要的。IUGG等国际组织对此已经强调了很多年。最近，由德国的GFZ（GeoForschungsZentrum)，美国的NASA（National Aeronautics and Space Adminitration)以及欧洲宇航局ESA（European Space Agency)开发研制了最先进的地球监测技术－SST（Satellite-to-Sateilite Tracking)。其主要特点是利用现有的GPS连续追踪新发射低轨道卫星，并由低轨道卫星对地球重力场作精密观测。已经发射和即将发射的卫星有3颗：GHAMP（Challenging Mini-Satellite Payload for Geophysical Research an Application)已经于2000年发射；GRACE（Gravity Recovery and Climate Experimert)定于2002年发射；GOCE（Gravity Field and Steady-state Ocean Cirulation Explorer)计划2004年发射，它们可以统称为重力卫星。载有SST技术的人造卫星的主要目的是获得具有前所未有的高精度和高空间解析度的全球重力场和大地水准面模型，加强人们对地球内部构造的理解并为海洋和气象研究提供更好地参考。上述3个重力卫星工作在有明显区别的不同波谱内，它们有不同的科学应用，仅有一小部分重合。所以，就应用而言它们是完全互补的。它们在地球科学中的应用将是广泛的，特别对于固体地球物理学，海洋学以及大地测量学等领域，它们将会带来革命性的变化，其意义不亚于GPS。     

THE POSITION AND ROLES OF PORTS AND RAILWAYS IN THE EURASIA LAND BRIDGE TRANSPORTATION IN NORTHEAST CHINA

This paper systematically analyzes and proves the favorable factors of utilizing the ports and the railways in Northeast China for the Asia-Europe land bridge transportation. It will be more beneficial to Dalian Port and HaDa (Harbin-Dalian) and Binzhou (Harbin-Manzhouli) railways in Northeast China to the Siberia railway in Russia than 1) to other ports and the corresponding railways in China, 2) to the Nakhodka Port in Russia and the corresponding railways and 3) to the Chongjin Port in North Korea and the corresponding railways. This paper also puts forward the reform measures to adopt the ports and the railways in the northeast region in China for the land bridge transportation and the problems in transportation policies and management systems.     

沂水崔家峪玻璃用石英砂岩矿地质特征

沂水崔家峪玻璃用石英砂岩矿床赋存于早寒武世李官组砂岩段中．呈近水平的层状产出。因其岩石坚硬．矿体呈环山的平台状分布。矿体厚度大，矿石品级高，特级品矿石二氧化硅平均含量98．47％，铁杂质平均含量0．043％(选矿后，铁杂质含量可降至0．02％以下)，为一优质玻璃硅质原料矿床。矿石为细一中粒石英砂岩，粒度以中粒为主，矿石由碎屑颗粒和胶结物组成，碎屑成分含量为97％～98％．其中绝大部分是石英颗粒．具典型的砂屑结掏。该矿床属滨海陆源沉积矿床。     

Fluid Inclusion Feature of Gold-Bearing Quartz Veins from the Jinniushan Gold Deposit, Shandong, China

The auriferous veins at Jinniushan occurs within the Jinniushan faulted zone in the Kunyushan Granite. Optical observation reveals that gold ore body formed during the main stage of hydrothermal activity. Detailed geothermometric studies of fluid inclusions from the veins show that the forming temperature ranges between 130℃ and 370℃ and the salinity is from 4.01 to 15.21 wvt percent NaCl. The ore-forming fluid is featured by low to moderate salinity, and low to moderate temperature. According to investigations of the values of vapor/liquid and temperatures of the ore-forming fluids, we propose that the boiling fluid inclusions exist in the main mineralization stages. Fluid boiling is suggested as a mechanism for the precipitation of gold from the hydrothermal fluid in the Jinniushan gold deposit.