Forest operations,tree species composition and decline in rainfall explain runoff changes in the Nacimiento experimental catchments,south central Chile

Few long-term studies have explored how intensively managed short rotation forest plantations interact with climate variability. We examine how prolonged severe drought and forest operations affect runoff in 11 experimental catchments on private corporate forest land near Nacimiento in south central Chile over the period 2008–2019. The catchments (7.7–414 ha) contain forest plantations of exotic fast-growing species (Pinus radiata, Eucalyptus spp.) at various stages of growth in a Mediterranean climate (mean long-term annual rainfall = 1381 mm). Since 2010, a drought, unprecedented in recent history, has reduced rainfall at Nacimiento by 20%, relative to the long-term mean. Pre-drought runoff ratios were <0.2 under 8-year-old Eucalyptus; >0.4 under 21-year-old Radiata pine and >0.8 where herbicide treatments had controlled vegetation for 2 years in 38% of the catchment area. Early in the study period, clearcutting of Radiata pine (85%–95% of catchment area) increased streamflow by 150 mm as compared with the year before harvest, while clearcutting and partial cuts of Eucalyptus did not increase streamflow. During 2008–2019, the combination of emerging drought and forestry treatments (replanting with Eucalyptus after clearcutting of Radiata pine and Eucalyptus) reduced streamflow by 400–500 mm, and regeneration of previously herbicide-treated vegetation combined with growth of Eucalyptus plantations reduced streamflow by 1125 mm (87% of mean annual precipitation 2010–2019). These results from one of the most comprehensive forest catchment studies in the world on private industrial forest land indicate that multiple decades of forest management have reduced deep soil moisture reservoirs. This effect has been exacerbated by drought and conversion from Radiata pine to Eucalyptus, apparently largely eliminating subsurface supply to streamflow. The findings reveal tradeoffs between wood production and water supply, provide lessons for adapting forest management to the projected future drier climate in Chile, and underscore the need for continued experimental work in managed forest plantations.     

Local ice caps in Finderup Land,North Greenland,survived the Holocene Thermal Maximum

Local glaciers and ice caps (GICs) comprise only ~5.4% of the total ice volume, but account for ~14–20% of the current ice loss in Greenland. The glacial history of GICs is not well constrained, however, and little is known about how they reacted to Holocene climate changes. Specifically, in North Greenland, there is limited knowledge about past GIC fluctuations and whether they survived the Holocene Thermal Maximum (HTM, ~8 to 5 ka). In this study, we use proglacial lake records to constrain the ice‐marginal fluctuations of three local ice caps in North Greenland including Flade Isblink, the largest ice cap in Greenland. Additionally, we have radiocarbon dated reworked marine molluscs in Little Ice Age (LIA) moraines adjacent to the Flade Isblink, which reveal when the ice cap was smaller than present. We found that outlet glaciers from Flade Isblink retreated inland of their present extent from ~9.4 to 0.2 cal. ka BP. The proglacial lake records, however, demonstrate that the lakes continued to receive glacial meltwater throughout the entire Holocene. This implies that GICs in Finderup Land survived the HTM. Our results are consistent with other observations from North Greenland but differ from locations in southern Greenland where all records show that the local ice caps at low and intermediate elevations disappeared completely during the HTM. We explain the north–south gradient in glacier response as a result of sensitivity to increased temperature and precipitation. While the increased temperatures during the HTM led to a complete melting of GICs in southern Greenland, GICs remained in North Greenland probably because the melting was counterbalanced by increased precipitation due to a reduction in Arctic sea‐ice extent and/or increased poleward moisture transport.     

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.     

The isotopic signature of monsoon conditions,cloud modes,and rainfall type

This work provides a comprehensive physically based framework for the interpretation of the north Australian rainfall stable isotope record (δ¹⁸O and δ²H). Until now, interpretations mainly relied on statistical relationships between rainfall amount and isotopic values on monthly timescales. Here, we use multiseason daily rainfall stable isotope and high resolution (10 min) ground‐based C‐band polarimetric radar data and show that the five weather types (monsoon regimes) that constitute the Australian wet season each have a characteristic isotope ratio. The data suggest that this is not only due to changes in regional rainfall amount during these regimes but, more importantly, is due to different rain and cloud types that are associated with the large scale circulation regimes. Negative (positive) isotope anomalies occurred when stratiform rainfall fractions were large (small) and the horizontal extent of raining areas were largest (smallest). Intense, yet isolated, convective conditions were associated with enriched isotope values whereas more depleted isotope values were observed when convection was widespread but less intense. This means that isotopic proxy records may record the frequency of which these typical wet season regimes occur. Positive anomalies in paleoclimatic records are most likely associated with periods where continental convection dominates and convection is sea‐breeze forced. Negative anomalies may be interpreted as periods when the monsoon trough is active, convection is of the oceanic type, less electric, and stratiform areas are wide spread. This connection between variability of rainfall isotope anomalies and the intrinsic properties of convection and its large‐scale environment has important implications for all fields of research that use rainfall stable isotopes.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

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.     

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.     

Changes in hillslope hydrology in a perched,shallow soil system due to clearcutting and residual biomass removal

Sustainable fuels legislation and volatility in energy prices have put additional pressures on the forestry sector to intensify the harvesting of biomass for “advanced biofuel” production. To better understand how residual biomass removal after harvest affects forest hydrology in relatively low slope terrain, a Before-After-Control-Impact (BACI) study was conducted in the USDA Forest Service's Marcell Experimental Forest, Minnesota, USA. Hydrological measurements were made from 2010–2013 on a forested hillslope that was divided into three treatment blocks, where one block was harvested and residual biomass removed (Biomass Removed), the second was harvested and residual biomass left (Biomass Left), and the last block was left as an Unharvested Control. The pre-harvest period (2 years) was 2010–11 and post-harvest (2 years) was 2012–13. Water table elevation at the upslope and downslope position, subsurface runoff, and soil moisture were measured between May–November. Mixed effect statistical models were used to compare both the before-after and “control” treatment ratios (ratios between harvested hillslopes and the Unharvested Control hillslope). Subsurface runoff significantly increased (p < .05) at both harvested hillslopes but to a greater degree on the Biomass Left hillslope. Greater subsurface runoff volumes at both harvested hillslopes were driven by substantial increases during fall, with additional significant increases during summer on the Biomass Left hillslope. The hydrological connectivity, inferred from event runoff ratios, increased due to harvesting at both hillslopes but only significantly on the Biomass Left hillslope. The winter harvest minimized soil disturbance, resulting in no change to the effective hydraulic conductivity distribution with depth. Thus, the observed hydrological changes were driven by increased effective precipitation and decreased evapotranspiration, increasing the duration that both harvested hillslopes were hydrologically active. The harvesting of residual biomass appears to lessen hydrological connectivity relative to leaving residual biomass on the hillslope, potentially decreasing downstream hydrological impacts of similar forestry operations.