The Bear Brook Watershed in Maine: Multi-decadal whole-watershed experimental acidification

The Bear Brook Watershed in Maine (BBWM) is a long-term research site established to study the response of forest ecosystem function to environmental disturbances of chronic acidic deposition and ecosystem nitrogen enrichment. Starting in 1989, the West Bear (treated) watershed received bimonthly applications of ammonium sulfate [(NH₄)₂SO₄] fertilizer from above the canopy, whereas East Bear (reference) received ambient deposition. The treatments were stopped in 2016, marking the beginning of the recovery phase. Research at the site has focused on soils, streams, and vegetation. Here, we describe data collected over three decades at the BBWM—input and stream output nutrient fluxes, quantitative soil pits and soil chemistry, and soil temperature and moisture.     

Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

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.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

The Upper Devonian Kellwasser event recorded in a regressive sequence from inner shelf to lagoonal pond,Catalan Coastal Ranges,Spain

This study, conducted in the Catalan Coastal Ranges, north‐east Spain, describes the Upper Devonian Kellwasser event in a shallowing‐upward sequence of black shales, siltstones and quartz arenites. This sequence was deposited in a progradational and regressive coastal system where the sedimentary environment evolved from the inner shelf to a lagoonal pond located landward of the shoreline. Three anomalous succeeding steps have been identified by geochemical analysis. The first one, detected on the inner shelf, was characterized by oxygen depletion and high organic productivity. The second, detected in the nearshore, was caused by hydrothermal activity occurring under normal oxic conditions. The third and most intense step was identified in the muds of the lagoonal pond and has been linked to strong anoxic conditions, elevated clastic input derived from changes in the weathering regime at the source area and moderate hydrothermal activity. The Kellwasser event is thus defined in the study area as stepwise and multi‐causal. This is the first time that the Kellwasser event has been identified in a sedimentary environment behind the shoreline. It is also the first time that it has been reported in the Catalan Coastal Ranges.     

基于噪声谱比法的内蒙古西部地区台站场地响应研究

利用内蒙古西部12个台站的地脉动噪声数据，采用噪声谱比法研究台站的场地响应情况。分析表明，台站场地响应按曲线形态可分为3类，且可能受地形地貌、局部构造和台基状况等条件影响。对比分析噪声谱比法与Moya方法的场地响应结果发现，场地响应曲线形态基本一致，只有极少数台站存在明显差异。结果表明，内蒙古西部区域大部分台站的台基状况总体较好，场地响应曲线较为平坦，无明显频率放大点。     

一种有发展前景的隧道及井巷超前预报技术——电磁导弹探测技术

对于隧道和煤矿井巷的安全掘进 ,超前预报显得特别重要。但由于掌子面前的特殊环境限制及多种干扰的严重影响 ,给隧道及井巷超前预报技术的研究带来许多困难。在瞬变电磁理论的基础上 ,引入新的电磁传播理论———电磁导弹技术 ,作为对该技术的初步探索 ,文章论证了该技术在隧道超前预报中应用的可能性 ,并提出了研究的主要内容和发展方向以及可能取得的新的技术突破和技术支撑     

Activity and kinematics of members of the TW Hydrae association

We present high-resolution echelle spectroscopy of 20 stars in 16 systems catalogued as members of the TW Hydrae association, and 16 stars identified as possible new members. We have calibrated the range of coronal and chromospheric activity expected for such young stars as a function of spectral type by combining our observations with literature data for field and open cluster stars. We also compute space motions for TWA members and candidate members with proper motion measurements, using two techniques to estimate distances to stars lacking direct trigonometric parallax measurements. The mean space motion of the four TWA members with known parallaxes is  ( U , V , W : −10.0, −17.8, −4.6) km s⁻¹  . 14 of the candidates have properties inconsistent with cluster membership; the remaining two are potential new members, although further observations are required to confirm this possibility.     

First evidence of IMF control of Jovian magnetospheric boundary locations: Cassini and Galileo magnetic field measurements compared

The Cassini spacecraft, en route to Saturn, passed close to Jupiter while the Galileo spacecraft was completing its 28th and 29th orbits of Jupiter, thus offering a unique opportunity for direct study of the solar wind-Jovian interaction. Here evidence is given of response of the Jovian magnetopause and bow shock positions to changes of the north-south component of the solar wind magnetic field, a phenomenon long known to occur in equivalent circumstances at Earth. The period analyzed starts with the passage over Cassini of an interplanetary shock far upstream of Jupiter. The shock's arrival at Galileo on the dusk-flank of the magnetosphere caused Galileo to exit into the solar wind. Using inter-spacecraft timing based on the time delay established from the shock arrival at each spacecraft, we point out that Galileo's position with respect to the Jovian bow shock appears to correlate with changes in the disturbed north-south reversing field seen behind the shock. We specifically rule out the alternative of changes in the shape of the bow shock with rotations of the interplanetary magnetic field as the cause.