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.

     

Large debris avalanches and associated eruptions in the Holocene eruptive history of Shiveluch Volcano, Kamchatka, Russia

 Shiveluch Volcano, located in the Central Kamchatka Depression, has experienced multiple flank failures during its lifetime, most recently in 1964. The overlapping deposits of at least 13 large Holocene debris avalanches cover an area of approximately 200 km² of the southern sector of the volcano. Deposits of two debris avalanches associated with flank extrusive domes are, in addition, located on its western slope. The maximum travel distance of individual Holocene avalanches exceeds 20 km, and their volumes reach ∼3 km³. The deposits of most avalanches typically have a hummocky surface, are poorly sorted and graded, and contain angular heterogeneous rock fragments of various sizes surrounded by coarse to fine matrix. The deposits differ in color, indicating different sources on the edifice. Tephrochronological and radiocarbon dating of the avalanches shows that the first large Holocene avalanches were emplaced approximately 4530–4350 BC. From ∼2490 BC at least 13 avalanches occurred after intervals of 30–900 years. Six large avalanches were emplaced between 120 and 970 AD, with recurrence intervals of 30–340 years. All the debris avalanches were followed by eruptions that produced various types of pyroclastic deposits. Features of some surge deposits suggest that they might have originated as a result of directed blasts triggered by rockslides. Most avalanche deposits are composed of fresh andesitic rocks of extrusive domes, so the avalanches might have resulted from the high magma supply rate and the repetitive formation of the domes. No trace of the 1854 summit failure mentioned in historical records has been found beyond 8 km from the crater; perhaps witnesses exaggerated or misinterpreted the events. Received: 18 August 1997 / Accepted: 19 December 1997     

Process length variation of the cyst of the dinoflagellate Protoceratium reticulatum in the North Pacific and Baltic‐Skagerrak region: calibration as an annual density proxy and first evidence of pseudo‐cryptic speciation

Process length variation of cysts of the dinoflagellate Protoceratium reticulatum (Claparède et Lachmann) Bütschli in surface sediments from the North Pacific was investigated. The average process length showed a significant inverse relation to annual seawater density: σ_{t annual} = ?0.8674 × average process length + 1029.3 (R² = 0.84), with a standard error of 0.78 kg m^?3. A sediment trap study from Effingham Inlet in British Columbia revealed the same relationship between average process length and local seawater density variations. In the Baltic–Skagerrak region, the average process length variation was related significantly to annual seawater density: σ_{t annual} = 3.5457 × average process length ? 993.28 (R² = 0.86), with a standard error of 3.09 kg m^?3. These calibrations cannot be reconciled, which accentuates the regional character of the calibrations. This can be related to variations in molecular data (small subunit, long subunit and internal transcribed spacer sequences), which show the presence of several genotypes and the occurrence of pseudo‐cryptic speciation within this species. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Paleomagnetic studies of Paleolithic site deposits in the Matzuka Cave (Northern Caucasus)

As a result of detailed paleomagnetic and magnetic studies of Paleolithic site deposits in the Matuzka Cave, a record of the Matuzka geomagnetic excursion in lithologic layer 7 has been discovered and studied. Such characteristic features as the geomagnetic field direction, position of the virtual geomagnetic pole, geomagnetic field intensity (roughly estimated) after and during the excursion, and climatic conditions coeval with its existence make the Matuzka excursion similar to the ～130-ka Blake excursion. This dates at ～130 ka the formation of layer 7 with ancient archaeological findings.     

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.