Intercomparison of measurements of bulk snow density and water equivalent of snow cover with snow core samplers: Instrumental bias and variability induced by observers

Manually collected snow data are often considered as ground truth for many applications such as climatological or hydrological studies. However, there are many sources of uncertainty that are not quantified in detail. For the determination of water equivalent of snow cover (SWE), different snow core samplers and scales are used, but they are all based on the same measurement principle. We conducted two field campaigns with 9 samplers commonly used in observational measurements and research in Europe and northern America to better quantify uncertainties when measuring depth, density and SWE with core samplers. During the first campaign, as a first approach to distinguish snow variability measured at the plot and at the point scale, repeated measurements were taken along two 20 m long snow pits. The results revealed a much higher variability of SWE at the plot scale (resulting from both natural variability and instrumental bias) compared to repeated measurements at the same spot (resulting mostly from error induced by observers or very small scale variability of snow depth). The exceptionally homogeneous snowpack found in the second campaign permitted to almost neglect the natural variability of the snowpack properties and focus on the separation between instrumental bias and error induced by observers. Reported uncertainties refer to a shallow, homogeneous tundra-taiga snowpack less than 1 m deep (loose, mostly recrystallised snow and no wind impact). Under such measurement conditions, the uncertainty in bulk snow density estimation is about 5% for an individual instrument and is close to 10% among different instruments. Results confirmed that instrumental bias exceeded both the natural variability and the error induced by observers, even in the case when observers were not familiar with a given snow core sampler.     

Characterization of complex groundwater flows in the environment of singular buildings by combining hydrogeological and non-destructive geophysical (ground-penetrating radar) techniques: Punta Begoña Galleries (Getxo,Spain)

Locating and quantifying groundwater flow in many built-up areas are a priority with regard to its complete restoration. In this work, a hydrogeological survey of the surroundings of the Punta Begoña Galleries (Getxo, Bizkaia), built on a coastal cliff, was completed by using ground penetrating radar (GPR) testing. Thus, the preliminary characterization of soils and rocks in accessible areas of the cliff was first improved by hydrogeological information gathered from a single survey borehole, including permeability measurements by low pressure injection tests (LPTs) and continuous water level monitoring. As a complementary method, the non-destructive GPR technique was performed during both dry and wet hydrological periods and in tandem with the injection tests, providing more complete spatial and temporal images of water flows. Specifically, GPR allows mapping of flow paths in soils and assessing the continuity of fractures in rock masses. Altogether, this complementary approach provides greater knowledge of complex underground flow dynamics in built environments, thus making it easier to make decisions for their management.     

Climatic vs. tectonic signals in a continental extensional basin (Teruel,NE Spain) from stable isotope (δ18O) and sequence stratigraphical evolution

Alluvial fans and shallow carbonate lakes interfered in the Teruel half‐graben during the Late Miocene–Pliocene. Tectonic influence is recorded in alluvial and lacustrine–palustrine successions, with long‐term climate changes being recorded in detail in the isotopic signatures of carbonates. Episodes of tectonic activity induced alluvial fan progradation and lake retraction in the whole basin. Three lacustrine stages have been identified, which support the idea that climate also exerted an important control on sedimentation. The transition between stages 1 and 2 occurred during a tectonically calm episode due to an increase in aridity in the Early Turolian; small fans with source areas next to the lake margin prograded, inducing lake‐shore retraction. The transition from stage 2 to 3 was caused by the superimposition of increasing tectonic activity and aridity effects. Our study demonstrates that discrimination of allogenic factors controlling sedimentation in continental closed basins is possible using sequence stratigraphy in combination with other techniques such as geochemistry of carbonates.     

Grain‐scale pressure variations and chemical equilibrium in high‐grade metamorphic rocks

In the classical view of metamorphic microstructures, fast viscous relaxation (and so constant pressure) is assumed, with diffusion being the limiting factor in equilibration. This contribution is focused on the only other possible scenario – fast diffusion and slow viscous relaxation – and brings an alternative interpretation of microstructures typical of high‐grade metamorphic rocks. In contrast to the pressure vessel mechanical model applied to pressure variation associated with coesite inclusions in various host minerals, a multi‐anvil mechanical model is proposed in which strong single crystals and weak grain boundaries can maintain pressure variation at geological time‐scales in a polycrystalline material. In such a mechanical context, exsolution lamellae in feldspar are used to show that feldspar can sustain large differential stresses (>10 kbar) at geological time‐scales. Furthermore, it is argued that the existence of grain‐scale pressure gradients combined with diffusional equilibrium may explain chemical zoning preserved in reaction rims. Assuming zero net flux across the microstructure, an equilibrium thermodynamic method is introduced for inferring pressure variation corresponding to the chemical zoning. This new barometric method is applied to plagioclase rims around kyanite in felsic granulite (Bohemian Massif, Czech Republic), yielding a grain‐scale pressure variation of 8 kbar. In this approach, kinetic factors are not invoked to account for mineral composition zoning preserved in rocks metamorphosed at high grade.     

Correction to: Spatial distribution and estimation of rainfall trends and erosivity in the Epitácio Pessoa reservoir catchment,Paraíba,Brazil

Natural Hazards - The article “Spatial distribution and estimation of rainfall trends and erosivity in the Epitácio Pessoa reservoir catchment, Paraíba, Brazil”, written by da...     

Mercury cycling during acid rain recovery at the forested Lesní potok catchment,Czech Republic

From 2011 to 2019, mercury (Hg) stores and fluxes were studied in the small forested catchment Lesní potok (LES) in the central Czech Republic using the watershed mass balance approach together with internal measurements. Mean input fluxes of Hg via open bulk deposition, beech throughfall and spruce throughfall during the periodwere 2.9, 3.9 and 7.6 μg m⁻² year⁻¹, respectively. These values were considerably lower than corresponding deposition Hg fluxes reported in the early years of the 21st century from catchments in Germany. Current bulk precipitation inputs at unimpacted Czech mountainous sites were lower than those in Germany. The largest Hg inputs to the catchment were via litterfall, averaging 22.6 and 17.8 μg m⁻² year⁻¹ for beech and spruce stands. The average Hg input, based on the sum of mean litterfall and throughfall deposition, was 23.0 μg m⁻² year⁻¹, compared to the estimated Hg output in runoff of 0.5 μg m⁻² year⁻¹, which is low compared to other reported values. Thus, only ~2% of Hg input is exported in stream runoff. Stream water Hg was only weakly related to dissolved organic carbon (DOC) but both concentrations were positively correlated with water temperature. The estimated total soil Hg pool averaged 47.5 mg m⁻², only 4% of which was in the O-horizon. Thus Hg in the O-horizon pool represents 72 years of deposition at the current input flux and 3800 years of export at the current runoff flux. Age-dating by ¹⁴C suggested that organic soil contains Hg from recent deposition, while mineral soil at 40–80 cm depth contained 4400-year old carbon, suggesting the soil had accumulated atmospheric Hg inputs through millennia to reach the highest soil Hg pool of the soil profile. These findings suggest that industrial era intensification of the Hg cycle is superimposed on a slower-paced Hg cycle during most of the Holocene.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

Late Quaternary (Marine Isotope Stage 3 to Recent) sedimentary evolution of the Basque shelf (southern Bay of Biscay)

Late Quaternary (MIS 3 to Recent) oceanographic evolution of the Basque shelf has been analysed for the first time based on the sedimentological analysis of three cores obtained from the middle and outer shelves. The cores are located in two interfluves separated by the San Sebastian canyon. The variability of the percentage of the planktonic foraminifera species Neogloboquadrina pachyderma sin. and of δ¹⁸O_bull allowed us to identify the influence of colder and warmer waters in the Basque shelf during the late Quaternary. From ～56 cal. ka BP to the end of the Last Glacial Maximum (19 cal. ka BP) the sedimentary record shows a decreasing trend in the mean grain size that correlates with the eustatic sea‐level fall. The last Deglaciation (19–11.5 cal. ka BP) is characterized by a sea‐level rise that produced an important hiatus in the western outer shelf. During the Holocene, the middle and outer shelves present different behaviours. From 11.5 to 6.7 cal. ka BP, in the outer shelf the sea‐level rise that started during the Deglaciation produced a hiatus, whereas in the middle shelf the sedimentary succession records the presence of warm to temperate waters. Between 6.7–4.9 cal. ka BP, the entrance of cold surface water‐masses that only affected the middle shelf has been identified, and temperate to warm waters occurred in the outer shelf. The cold surface water‐masses retreated during 4.9–4.3 cal. ka BP in the middle shelf. Finally, from 4.3 cal. ka BP to Recent, the middle shelf registers a hiatus due to sea‐level stabilization after a generalized transgression, synchronous to a decrease in the energy of the water‐masses in the outer shelf. In conclusion, the environmental changes detected in the Basque shelf are attributed to both regional and eustatic sea‐level changes.