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.     

Geology and the Allied liberation of Normandy: highlights to help mark the 75th anniversary of D‐Day, 6 June 1944

During the Second World War, the Allied invasion of the French coast of Normandy on D‐Day, 6 June 1944, was the greatest amphibious assault in world history. An article in Geology Today (v.11, for 1995, pp.58–63) marked the 50th anniversary of the end of the war in Europe, on 8 May 1945, by describing how British military geologists had participated in planning for D‐Day and in the NW Europe campaign that followed it. The work of these geologists provides a classic case history, revealing that ‘military geology’ has many potential applications. Geological factors influenced site selection for temporary airfields, predictions of trafficability for the Normandy beaches, the development of potable water supplies, and quarrying for road metal—and more besides. This new article helps to mark the 75th anniversary of D‐Day by further details of how geologists and geology contributed to Allied victory.     

Crinoid localities of the Devonian of the British Isles

Who collects the Devonian crinoids of south-west England? Since the nineteenth century, almost nobody. Few palaeontologists pursue these fossils, undoubtedly handsome in life, but indifferently preserved. Herein, we make a survey of four of the most important crinoid sites, all of which have the potential to produce new specimens.     

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.     

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.     

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.     

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).     

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.