Influence of a permeable sand layer on the mechanism of backward erosion piping using 3D pipe depth measurements

Backward erosion piping (BEP) poses a threat to the stability of water-retaining structures. This can lead to severe erosion and collapse of embankments. A novel economically appealing measure against BEP is the coarse sand barrier (CSB). The CSB is a trench filled with coarse sand that is placed below the blanket layer on the landward side of the embankment, which prevents the pipe from developing upstream when it encounters the CSB. Inclusion of a CSB creates a vertically layered sand, which is the situation that can also exist in practice but is different from traditional BEP tests with one homogeneous sand. This paper presents new observations and measurements in medium-scale laboratory tests. 3D measurements of the pipe depth and dimensions are presented and analysed. This analysis indicates how the pipe dimensions evolve during the piping process and shows the erosion mechanism for BEP in vertically layered sands. The findings demonstrate the significance of three-dimensional study of the pipe rather than two dimensions. The pipe depth, width and depth-to-width ratios at the pipe tip in critical erosion stages are measured and presented. In the presented tests, two different erosion behaviours (stepwise pipe progression until failure and straight failure) are found and analysed with respect to possible influential parameters. Higher head drops and flow rates are found in tests with straight failure at the stage before progression. A linear relationship between the hydraulic conductivity contrast (k_c) and the critical head drops (h_c) is found and observations are used to investigate deviations from the line.

     

Deciphering petrogenic processes using Pb isotope ratios from time-series samples at Bezymianny and Klyuchevskoy volcanoes,Central Kamchatka Depression

The Klyuchevskoy group of volcanoes in the Kamchatka arc erupts compositionally diverse magmas (high-Mg basalts to dacites) over small spatial scales. New high-precision Pb isotope data from modern juvenile (1956–present) erupted products and hosted enclaves and xenoliths from Bezymianny volcano reveal that Bezymianny and Klyuchevskoy volcanoes, separated by only 9 km, undergo varying degrees of crustal processing through independent crustal columns. Lead isotope compositions of Klyuchevskoy basalts–basaltic andesites are more radiogenic than Bezymianny andesites (²⁰⁸Pb/²⁰⁴Pb = 37.850–37.903, ²⁰⁷Pb/²⁰⁴Pb = 15.468–15.480, and ²⁰⁶Pb/²⁰⁴Pb = 18.249–18.278 at Bezymianny; ²⁰⁸Pb/²⁰⁴Pb = 37.907–37.949, ²⁰⁷Pb/²⁰⁴Pb = 15.478–15.487, and ²⁰⁶Pb/²⁰⁴Pb = 18.289–18.305 at Klyuchevskoy). A mid-crustal xenolith with a crystallization pressure of 5.2 ± 0.6 kbars inferred from two-pyroxene geobarometry and basaltic andesite enclaves from Bezymianny record less radiogenic Pb isotope compositions than their host magmas. Hence, assimilation of such lithologies in the middle or lower crust can explain the Pb isotope data in Bezymianny andesites, although a component of magma mixing with less radiogenic mafic recharge magmas and possible mantle heterogeneity cannot be excluded. Lead isotope compositions for the Klyuchevskoy Group are less radiogenic than other arc segments (Karymsky—Eastern Volcanic Zone; Shiveluch—Northern Central Kamchatka Depression), which indicate increased lower-crustal assimilation beneath the Klyuchevskoy Group. Decadal timescale Pb isotope variations at Klyuchevskoy demonstrate rapid changes in the magnitude of assimilation at a volcanic center. Lead isotope data coupled with trace element data reflect the influence of crustal processes on magma compositions even in thin mafic volcanic arcs.     

Long-term soil moisture dynamics derived from GNSS interferometric reflectometry: a case study for Sutherland,South Africa

Soil moisture is a geophysical key observable for predicting floods and droughts, modeling weather and climate and optimizing agricultural management. Currently available in situ observations are limited to small sampling volumes and restricted number of sites, whereas measurements from satellites lack spatial resolution. Global navigation satellite system (GNSS) receivers can be used to estimate soil moisture time series at an intermediate scale of about 1000 m². In this study, GNSS signal-to-noise ratio (SNR) data at the station Sutherland, South Africa, are used to estimate soil moisture variations during 2008–2014. The results capture the wetting and drying cycles in response to rainfall. The GNSS Volumetric Water Content (VWC) is highly correlated (r ² = 0.8) with in situ observations by time-domain reflectometry sensors and is accurate to 0.05 m³/m³. The soil moisture estimates derived from the SNR of the L1 and L2P signals compared to the L2C show small differences with a RMSE of 0.03 m³/m³. A reduction in the SNR sampling rate from 1 to 30 s has very little impact on the accuracy of the soil moisture estimates (RMSE of the VWC difference 1–30 s is 0.01 m³/m³). The results show that the existing data of the global tracking network with continuous observations of the L1 and L2P signals with a 30-s sampling rate over the last two decades can provide valuable complementary soil moisture observations worldwide.     

Incentives for field hydrology and data sharing: collaboration and compensation: reply to “A need for incentivizing field hydrology,especially in an era of open data”*

ABSTRACT

We thank Allen and Berghuijs for continuing the discussion on field hydrology and data sharing and discuss two incentives to promote data collection and sharing in hydrological sciences: a collaborative attitude and additional funding to make data publicly available.     

Coarsening of three-dimensional structured and unstructured grids for subsurface flow

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.     

Late Quaternary climatic history of the Chilean Channels based on fossil pollen and beetle analyses,with an analysis of the modern vegetation and pollen rain

Fossil beetles and pollen were examined from an intermorainal bog at Puerto Edén, Isla Wellington, Chile (latitude 49°08'S, longitude 74°25'W). Wood from near the base of the section has an age of 12 960 ± 150 yr BP. Occurrence of flightless beetle species in the basal peat sample is evidence that some members of the biota survived the last glacial maximum in refugia. The assumption that the Chilean Channels were entirely ice-covered is incorrect. Plants and insects that invaded the deglaciated terrain were those of an Empetrum heathland in which patches of Nothofagus forest were restricted to sheltered locations. The climate supporting the heathland is inferred to have been windier and probably drier than that of the present day. From 13 000 yr BP to 9500 yr BP Nothofagus forest expanded, possibly in response to less windiness and more available moisture. Neither the fossil beetle nor pollen data support a return to significantly colder conditions between 11 000 and 10 000 yr BP at the time of the Younger Dryas Stade. From 9500 to 5500 yr BP the climate was as wet as that of the present day, based on an increased representation of the pollen of moorland plants and of aquatic beetle species. From 5500 to 3000 yr BP the climate was drier, as indicated by the expansion of Empetrum heath and the reduction in mesic habitats. From 3000 yr BP to the present-day mesic habitats dominated as the climate returned to a wetter mode. The alternatively wetter and drier episodes are attributed to latitudinal shifts in the position of storm tracks in the belt of Southern Westerlies.     

Surveillance for previously unmonitored organic contaminants in the San Francisco Estuary

The San Francisco Estuary Regional Monitoring Program initiated surveillance monitoring to identify previously unmonitored synthetic organic contaminants in the San Francisco Estuary. Organic extracts of water samples were analyzed using gas chromatography-mass spectrometry in full scan mode. The major contaminant classes identified in the samples were fire retardants, pesticides, personal care product ingredients, and plasticizers. Evidence from the literature suggests that some of these contaminants can persist in the environment, induce toxicity, and accumulate in marine biota and in higher food chain consumers. The major sources of these contaminants into the marine environment are the discharge of municipal and industrial wastewater effluents, urban stormwater, and agricultural runoff. As a proactive effort, it is suggested that surveillance studies be used routinely in monitoring programs to identify and prevent potential problem contaminants from harming the marine environment.     

Kriging Regionalized Positive Variables Revisited: Sample Space and Scale Considerations

Frequently, regionalized positive variables are treated by preliminarily applying a logarithm, and kriging estimates are back-transformed using classical formulae for the expectation of a lognormal random variable. This practice has several problems (lack of robustness, non-optimal confidence intervals, etc.), particularly when estimating block averages. Therefore, many practitioners take exponentials of the kriging estimates, although the final estimations are deemed as non-optimal. Another approach arises when the nature of the sample space and the scale of the data are considered. Since these concepts can be suitably captured by an Euclidean space structure, we may define an optimal kriging estimator for positive variables, with all properties analogous to those of linear geostatistical techniques, even for the estimation of block averages. In this particular case, no assumption on preservation of lognormality is needed. From a practical point of view, the proposed method coincides with the median estimator and offers theoretical ground to this extended practice. Thus, existing software and routines remain fully applicable.