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.     

Seasonal connections between meteoric water and streamflow generation along a mountain headwater stream

In snowmelt-driven mountain watersheds, the hydrologic connectivity between meteoric waters and stream flow generation varies strongly with the season, reflecting variable connection to soil and groundwater storage within the watershed. This variable connectivity regulates how streamflow generation mechanisms transform the seasonal and elevational variation in oxygen and hydrogen isotopic composition (δ¹⁸O and δD) of meteoric precipitation. Thus, water isotopes in stream flow can signal immediate connectivity or more prolonged mixing, especially in high-relief mountainous catchments. We characterized δ¹⁸O and δD values in stream water along an elevational gradient in a mountain headwater catchment in southwestern Montana. Stream water isotopic compositions related most strongly to elevation between February and March, exhibiting higher δ¹⁸O and δD values with decreasing elevation. These elevational isotopic lapse rates likely reflect increased connection between stream flow and proximal snow-derived water sources heavily subject to elevational isotopic effects. These patterns disappeared during summer sampling, when consistently lower δ¹⁸O and δD values of stream water reflected contributions from snowmelt or colder rainfall, despite much higher δ¹⁸O and δD values expected in warmer seasonal rainfall. The consistently low isotopic values and absence of a trend with elevation during summer suggest lower connectivity between summer precipitation and stream flow generation as a consequence of drier soils and greater transpiration. As further evidence of intermittent seasonal connectivity between the stream and adjacent groundwaters, we observed a late-winter flush of nitrate into the stream at higher elevations, consistent with increased connection to accumulating mineralized nitrogen in riparian wetlands. This pattern was distinct from mid-summer patterns of nitrate loading at lower elevations that suggested heightened human recreational activity along the stream corridor. These observations provide insights linking stream flow generation and seasonal water storage in high elevation mountainous watersheds. Greater understanding of the connections between surface water, soil water and groundwater in these environments will help predict how the quality and quantity of mountain runoff will respond to changing climate and allow better informed water management decisions.     

Water quality data for Fall Creek,New York,USA: 1972–1995

The 33 086 ha mixed land use Fall Creek watershed in upstate New York is part of the Great Lakes drainage system. Results from more than 3500 water samples are available in a data set that compiles flow data and measurements of various water quality analytes collected between 1972 and 1995 in all seasons and under all flow regimes in Fall Creek and its tributaries. Data is freely accessible at https://ecommons.cornell.edu/handle/1813/8148 and includes measurements of suspended solids, pH, alkalinity, calcium, magnesium, potassium, sodium, chloride, nitrate nitrogen (NO3-N), sulphate sulphur (SO4-S), phosphorus (P) fractions molybdate reactive P (MRP) and total dissolved P (TDP), percent P in sediment, and ammonium nitrogen (NH4-N). Methods, sub-watershed areas, and coordinates for sampling sites are also included. The work represented in this data set has made important scientific contributions to understanding of hydrological and biogeochemical processes that influence loading in mixed use watersheds and that have an impact on algal productivity in receiving water bodies. In addition, the work has been foundational for important regulatory and management decisions in the region.     

Variations in glacier volume and snow cover and their impact on lake storage in the Paiku Co Basin,in the Central Himalayas

Glaciers and snow cover are important constituents of the surface of the Tibetan Plateau. The responses of these phenomena to global environmental changes are sensitive, rapid and intensive due to the high altitudes and arid cold climate of the Tibetan Plateau. Based on multisource remote sensing data, including Landsat images, MOD10A2 snow product, ICESat, Cryosat-2 altimetry data and long-term ground climate observations, we analysed the dynamic changes of glaciers, snow melting and lake in the Paiku Co basin using extraction methods for glaciers and lake, the degree-day model and the ice and lake volume method. The interaction among the climate, ice-snow and the hydrological elements in Paiku Co is revealed. From 2000 to 2018, the basin tended to be drier, and rainfall decreased at a rate of −3.07 mm/a. The seasonal temperature difference in the basin increased, the maximum temperature increased at a rate of 0.02°C/a and the minimum temperature decreased at a rate of −0.06°C/a, which accelerated the melting from glaciers and snow at rates of 0.55 × 10⁷ m³/a and 0.29 × 10⁷ m³/a, respectively. The rate of contribution to the lake from rainfall, snow and glacier melted water was 55.6, 27.7 and 16.7%, respectively. In the past 18 years, the warmer and drier climate has caused the lake to shrink. The water level of the lake continued to decline at a rate of −0.02 m/a, and the lake water volume decreased by 4.85 × 10⁸ m³ at a rate of −0.27 × 10⁸ m³/a from 2000 to 2018. This evaluation is important for understanding how the snow and ice melting in the central Himalayas affect the regional water cycle.     

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.     

High-frequency monitoring of hydrological and biogeochemical fluxes in forested catchments of southern Chile

The variability of rainfall-dependent streamflow at catchment scale modulates many ecosystem processes in wet temperate forests. Runoff in small mountain catchments is characterized by a quick response to rainfall pulses which affects biogeochemical fluxes to all downstream systems. In wet-temperate climates, water erosion is the most important natural factor driving downstream soil and nutrient losses from upland ecosystems. Most hydrochemical studies have focused on water flux measurements at hourly scales, along with weekly or monthly samples for water chemistry. Here, we assessed how water and element flows from broad-leaved, evergreen forested catchments in southwestern South America, are influenced by different successional stages, quantifying runoff, sediment transport and nutrient fluxes during hourly rainfall events of different intensities. Hydrograph comparisons among different successional stages indicated that forested catchments differed in their responses to high intensity rainfall, with greater runoff in areas covered by secondary forests (SF), compared to old-growth forest cover (OG) and dense scrub vegetation (CH). Further, throughfall water was greatly nutrient enriched for all forest types. Suspended sediment loads varied between successional stages. SF catchments exported 455 kg of sediments per ha, followed by OG with 91 kg/ha and CH with 14 kg/ha, corresponding to 11 rainfall events measured from December 2013 to April 2014. Total nitrogen (TN) and phosphorus (TP) concentrations in stream water also varied with rainfall intensity. In seven rainfall events sampled during the study period, CH catchments exported less nutrients (46 kg/ha TN and 7 kg/ha TP) than SF catchments (718 kg/ha TN and 107 kg/ha TP), while OG catchments exported intermediate sediment loads (201 kg/ha TN and 23 kg/ha TP). Further, we found significant effects of successional stage attributes (vegetation structure and soil physical properties) and catchment morphometry on runoff and sediment concentrations, and greater nutrients retention in OG and CH catchments. We conclude that in these southern hemisphere, broad-leaved evergreen temperate forests, hydrological processes are driven by multiple interacting phenomena, including climate, vegetation, soils, topography, and disturbance history.     

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.     

Jurassic igneous rocks of the central Sanandaj–Sirjan zone (Iran) mark a propagating continental rift,not a magmatic arc

Jurassic igneous bodies of the Sanandaj–Sirjan zone (SaSZ) in SW Iran are generally considered as a magmatic arc but critical evaluation of modern geochronology, geochemistry and radiogenic isotopes challenges this conclusion. There is no evidence for sustained igneous activity along the ~1,200 km long SaSZ, as expected for a convergent plate margin; instead activity was brief at most sites and propagated NW at ~20 mm/a. Jurassic igneous rocks define a bimodal suite of gabbro‐diorite and granite. Chemical and isotopic compositions of mafic rocks indicate subcontinental lithospheric mantle sources that mostly lacked subduction‐related modifications. The arc‐like features of S‐type granites reflect massive involvement of Cadomian crust and younger sediments to generate felsic melts in response to mafic intrusions. We conclude that Jurassic SaSZ igneous activity occurred in a continental rift, not an arc. SaSZ igneous rocks do not indicate that subduction along the SW margin of Eurasia began in Jurassic time.