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.     

Impact of the sampling duration on the uncertainty of averaged velocity measurements with acoustic instruments

Average velocity in streams is a key variable for the analysis and modelling of hydrological and hydraulic processes underpinning water resources science and practice. The present study evaluates the impact of the sampling duration on the quality of average velocity measurements acquired with contemporary instruments such as Acoustic Doppler Velocimeters (ADV) an Acoustic Doppler Current Profilers (ADCP). The evaluation combines considerations on turbulent flows and principles and configurations of acoustic instruments with practical experience in conducting customized analysis for uncertainty analysis purposes. The study sheds new insights on the spatial and temporal variability of the uncertainty in the measurement of average velocities due to variable sampling durations acting in isolation from other sources of uncertainties. Sampling durations of 90 and 150 s are found sufficient for ADV and ADCP, respectively, to obtain reliable average velocities in a flow affected only by natural turbulence and instrument noise. Larger sampling durations are needed for measurements in most of the natural streams exposed to additional sources of data variability.     

Effects of changes in canopy interception on stream runoff response and recovery following clear-cutting of a Japanese coniferous forest in Fukuroyamasawa Experimental Watershed in Japan

Understanding changes in evapotranspiration during forest regrowth is essential to predict changes of stream runoff and recovery after forest cutting. Canopy interception (Ic) is an important component of evapotranspiration, however Ic changes and the impact on stream runoff during regrowth after cutting remains unclear due to limited observations. The objective of this study was to examine the effects of Ic changes on long-term stream runoff in a regrowth Japanese cedar and Japanese cypress forest following clear-cutting. This study was conducted in two 1-ha paired headwater catchments at Fukuroyamasawa Experimental Watershed in Japan. The catchments were 100% covered by Japanese coniferous plantation forest, one of which was 100% clear-cut in 1999 when the forest was 70 years old. In the treated catchment, annual runoff increased by 301 mm/year (14% of precipitation) the year following clear-cutting, and remained 185 mm/year (7.9% of precipitation) higher in the young regrowth forest for 12–14 years compared to the estimated runoff assuming no clear-cutting. The Ic change was −358 mm/year (17% of precipitation) after cutting and was −168 mm/year (6.7% of precipitation) in the 12–14 years old regrowth forest compared to the observed Ic during the pre-cutting period. Stream runoff increased in all seasons, and the Ic change was the main fraction of evapotranspiration change in all seasons throughout the observation period. These results suggest that the change in Ic accounted for most of the runoff response following forest cutting and the subsequent runoff recovery in this coniferous forest.     

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.     

A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Subsurface flow measurements using passive flux meters in variably-saturated cold-regions landscapes

To date, passive flux meters have predominantly been applied in temperate environments for tracking the movement of contaminants in groundwater. This study applies these instruments to reduce uncertainty in (typically instantaneous) flux measurements made in a low-gradient, wetland dominated, discontinuous permafrost environment. This method supports improved estimation of unsaturated and over-winter subsurface flows which are very difficult to quantify using hydraulic gradient-based approaches. Improved subsurface flow estimates can play a key role in understanding the water budget of this landscape.     

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.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

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

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

Muon density measurements with the KASCADE central detector

Frequency distributions of local muon densities in high-energy extensive air showers (EAS) are presented as signature of the primary cosmic ray energy spectrum in the knee region. Together with the gross shower variables like shower core position, angle of incidence, and the shower sizes, the KASCADE experiment is able to measure local muon densities for two different muon energy thresholds. The spectra have been reconstructed for various core distances, as well as for particular subsamples, classified on the basis of the shower size ratio N_μ/N_e. The measured density spectra of the total sample exhibit clear kinks reflecting the knee of the primary energy spectrum. While relatively sharp changes of the slopes are observed in the spectrum of EAS with small values of the shower size ratio, no such feature is detected at EAS of large N_μ/N_e ratio in the energy range of 1–10 PeV. Comparing the spectra for various thresholds and core distances with detailed Monte Carlo simulations the validity of EAS simulations is discussed.