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.     

Plant diversity and community structure of Brazilian <Emphasis Type="Italic">Páramos</Emphasis>

In Eastern South America, high altitude grasslands represent a mountain system that has a high number of endemic species. However, studies on the ecology of plant communities in these environments remain scarce. We aimed to evaluate the patterns of biodiversity and structure of plant communities from rocky outcrops in high altitude grasslands of three areas at the Caparaó National Park, southeastern Brazil, by sampling 300 randomly distributed plots. Then, we compared the floristic composition, relative abundance, and biological and vegetation spectra among areas. We classified species as endemic and non-endemic and verified the occurrence of endangered species. Species richness was evaluated by rarefaction analysis on the sampling units. The importance value and species abundance distribution (SAD) models were assessed. We also performed an indicator species analysis. We sampled 58 species belonging to 49 genera and 32 families. The number of species decreased with increasing altitude, with significant differences being observed among areas regarding richness, abundance, and cover. Of the total number of species, 10 are endemic to the Caparaó National Park and 17 are listed on the Brazilian Red List of endangered species. The dominant families on all peaks were Asteraceae and Poaceae. The SAD models showed lognormal and geometric distributions, corroborating the fact that 10 species that were common to all three areas were also the most dominant ones in the communities and showed the highest importance values, which ranged between 35% and 60%. Indicator species analysis revealed that 28 species (48.27%) were indicators. Of these, 42.85% had maximum specificity, meaning that they occurred only in one area. Thus, the number of species per life form ratio was similar among areas, yet vegetation spectra differed, especially for hemicryptophytes. The altimetric difference among the areas showed to be a very important driver in the community assembly, influencing the evaluated variables, however, other drivers as soil depth, slope and water could also influence the community structure on a smaller and local spatial scale.     

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.     

Creating an isotopically similar Earth–Moon system with correct angular momentum from a giant impact

The giant impact hypothesis is the dominant theory explaining the formation of our Moon. However, the inability to produce an isotopically similar Earth–Moon system with correct angular momentum has cast a shadow on its validity. Computer-generated impacts have been successful in producing virtual systems that possess many of the observed physical properties. However, addressing the isotopic similarities between the Earth and Moon coupled with correct angular momentum has proven to be challenging. Equilibration and evection resonance have been proposed as means of reconciling the models. In the summer of 2013, the Royal Society called a meeting solely to discuss the formation of the Moon. In this meeting, evection resonance and equilibration were both questioned as viable means of removing the deficiencies from giant impact models. The main concerns were that models were multi-staged and too complex. We present here initial impact conditions that produce an isotopically similar Earth–Moon system with correct angular momentum. This is done in a single-staged simulation. The initial parameters are straightforward and the results evolve solely from the impact. This was accomplished by colliding two roughly half-Earth-sized impactors, rotating in approximately the same plane in a high-energy, off-centered impact, where both impactors spin into the collision.     

Simulating the Propagation of Infrasonic Waves and Estimating the Energy of the Chelyabinsk Meteoroid Explosion Observed on February 15, 2013

Results obtained from simulating the propagation of infrasonic waves from the Chelyabinsk meteoroid explosion observed on February 15, 2013, are given. The pseudodifferential parabolic equation (PDPE) method has been used for calculations. Data on infrasonic waves recorded at the IS31 station (Aktyubinsk, Kazakhstan), located 542.7 km from the likely location of the explosion, have been analyzed. Six infrasonic arrivals (isolated clearly defined pulse signals) were recorded. It is shown that the first “fast” arrival (F) corresponds to the propagation of infrasound in a surface acoustic waveguide. The rest of the arrivals (T1–T5) are thermospheric. The agreement between the results of calculations based on the PDPE method and experimental data is satisfactory. The energy E of the explosion has been estimated using two methods. One of these methods is based on the law of conservation of the acoustic pulse I, which is a product of the wave profile area S/2 of the signal under analysis and the distance to its source E _I [kt] = 1.38 × 10^–10 (I [kg/s])^1.482. The other method is based on the relation between the energy of explosion and the dominant period T of recorded signal E _T [kt] = 1.02 × (T [s]²/σ)^3/2, where σ is the dimensionless distance determining the degree of nonlinear effects during the propagation of sound along ray trajectories. According to the data, the explosion energy E _I,T ranges from 1.87 to 32 kt TNT.     

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.