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.     

Modelling the effect of changing precipitation inputs on deep soil water utilization

Forests in the Southeastern United States are predicted to experience future changes in seasonal patterns of precipitation inputs as well as more variable precipitation events. These climate change‐induced alterations could increase drought and lower soil water availability. Drought could alter rooting patterns and increase the importance of deep roots that access subsurface water resources. To address plant response to drought in both deep rooting and soil water utilization as well as soil drainage, we utilize a throughfall reduction experiment in a loblolly pine plantation of the Southeastern United States to calibrate and validate a hydrological model. The model was accurately calibrated against field measured soil moisture data under ambient rainfall and validated using 30% throughfall reduction data. Using this model, we then tested these scenarios: (a) evenly reduced precipitation; (b) less precipitation in summer, more in winter; (c) same total amount of precipitation with less frequent but heavier storms; and (d) shallower rooting depth under the above 3 scenarios. When less precipitation was received, drainage decreased proportionally much faster than evapotranspiration implying plants will acquire water first to the detriment of drainage. When precipitation was reduced by more than 30%, plants relied on stored soil water to satisfy evapotranspiration suggesting 30% may be a threshold that if sustained over the long term would deplete plant available soil water. Under the third scenario, evapotranspiration and drainage decreased, whereas surface run‐off increased. Changes in root biomass measured before and 4 years after the throughfall reduction experiment were not detected among treatments. Model simulations, however, indicated gains in evapotranspiration with deeper roots under evenly reduced precipitation and seasonal precipitation redistribution scenarios but not when precipitation frequency was adjusted. Deep soil and deep rooting can provide an important buffer capacity when precipitation alone cannot satisfy the evapotranspirational demand of forests. How this buffering capacity will persist in the face of changing precipitation inputs, however, will depend less on seasonal redistribution than on the magnitude of reductions and changes in rainfall frequency.     

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.     

Evaluating soil water routing approaches in watershed-scale,ecohydrologic modelling

Soil water dynamics are central in linking and regulating natural cycles in ecohydrology, however, mathematical representation of soil water processes in models is challenging given the complexity of these interactions. To assess the impacts of soil water simulation approaches on various model outputs, the Soil and Water Assessment Tool was modified to accommodate an alternative soil water percolation method and tested at two geographically and climatically distinct, instrumented watersheds in the United States. Soil water was evaluated at the site scale via measured observations, and hydrologic and biophysical outputs were analysed at the watershed scale. Results demonstrated an improved Kling–Gupta Efficiency of up to 0.3 and a reduction in percent bias from 5 to 25% at the site scale, when soil water percolation was changed from a threshold, bucket-based approach to an alternative approach based on variable hydraulic conductivity. The primary difference between the approaches was attributed to the ability to simulate soil water content above field capacity for successive days; however, regardless of the approach, a lack of site-specific characterization of soil properties by the soils database at the site scale was found to severely limit the analysis. Differences in approach led to a regime shift in percolation from a few, high magnitude events to frequent, low magnitude events. At the watershed scale, the variable hydraulic conductivity-based approach reduced average annual percolation by 20–50 mm, directly impacting the water balance and subsequently biophysical predictions. For instance, annual denitrification increased by 14–24 kg/ha for the new approach. Overall, the study demonstrates the need for continued efforts to enhance soil water model representation for improving biophysical process simulations.     

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.     

Wetland hydroperiod classification in the western prairies using multitemporal synthetic aperture radar

Wetlands represent one of the world's most biodiverse and threatened ecosystem types and were diminished globally by about two‐thirds in the 20th century. There is continuing decline in wetland quantity and function due to infilling and other human activities. In addition, with climate change, warmer temperatures and changes in precipitation and evapotranspiration are reducing wetland surface and groundwater supplies, further altering wetland hydrology and vegetation. There is a need to automate inventory and monitoring of wetlands, and as a study system, we investigated the Shepard Slough wetlands complex, which includes numerous wetlands in urban, suburban, and agricultural zones in the prairie pothole region of southern Alberta, Canada. Here, wetlands are generally confined to depressions in the undulating terrain, challenging wetlands inventory and monitoring. This study applied threshold and frequency analysis routines for high‐resolution, single‐polarization (HH) RADARSAT‐2, synthetic aperture radar mapping. This enabled a growing season surface water extent hyroperiod‐based wetland classification, which can support water and wetland resource monitoring. This 3‐year study demonstrated synthetic aperture radar‐derived multitemporal open‐water masks provided an effective index of wetland permanence class, with overall accuracies of 89% to 95% compared with optical validation data, and RMSE between 0.2 and 0.7 m between model and field validation data. This allowed for characterizing the distribution and dynamics of 4 marsh wetlands hydroperiod classes, temporary, seasonal, semipermanent, and permanent, and mapping of the sequential vegetation bands that included emergent, obligate wetland, facultative wetland, and upland plant communities. Hydroperiod variation and surface water extent were found to be influenced by short‐term rainfall events in both wet and dry years. Seasonal hydroperiods in wetlands were particularly variable if there was a decrease in the temporary or semipermanent hydroperiod classes. In years with extreme rain events, the temporary wetlands especially increased relative to longer lasting wetlands (84% in 2015 with significant rainfall events, compared with 42% otherwise).     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.     

黑龙江哈尔滨白渔泡国家湿地公园沼泽、林地和农田土壤物理、化学和生物性质差异

为了揭示黑龙江哈尔滨白渔泡国家湿地公园沼泽、林地和农田土壤物理、化学和生物性质的差异,于2018年7月25日~8月2日,在湿地公园内,在天然芦苇(Phragmites australis)沼泽、林地、旱田和水田中设置采样地,采集不同深度(0~10 cm、10~20 cm和20~30 cm)的土壤样品,测定土壤样品的物理、化学和生物指标。研究结果表明,白渔泡国家湿地公园不同采样地土壤指标存在差异;与天然芦苇沼泽土壤相比,其它采样地土壤的含水量明显偏低,土壤全氮、全磷、碱解氮和有机质含量都明显偏小,水田土壤速效磷含量偏大;天然芦苇沼泽土壤脲酶、硝酸还原酶、纤维素酶、蛋白酶和β-葡萄糖苷酶活性都高于林地和农田土壤,水田0~10 cm和10~20 cm深度土壤的硝酸还原酶活性显著高于旱田和林地;与天然芦苇沼泽土壤相比,旱田土壤小于0.25 mm的小团聚体含量偏大,而其它采样地土壤的各粒级团聚体的比例变化较小,水田土壤团聚体平均重量直径比天然芦苇沼泽和旱田土壤低。     

建鲤基因组DNA的RAPD分析

用 1 5个随机引物对建鲤基因组DNA进行RAPD分析 ,结果显示 ,建鲤群体内的遗传相似系数为 0 91 81± 0 0 73 8,多态位点比例为 0 41 67,平均杂合度为 0 1 884。表明建鲤种质较纯 ,群体内的遗传变异程度较小。建鲤优良性状产生的分子学基础是等位基因的杂合程度较高