Stability of a circular tunnel in cohesive-frictional soil subjected to surcharge loading

The stability of circular tunnels in cohesive-frictional soils subjected to surcharge loading has been investigated theoretically and numerically assuming plane strain conditions. Despite the importance of this problem, previous research on the subject is very limited. At present, no generally accepted design or analysis method is available to evaluate the stability of tunnels/openings in cohesive-frictional soils. In this study, continuous loading is applied to the ground surface, and both smooth and rough interface conditions are modelled. For a series of tunnel diameter-to-depth ratios and material properties, rigorous lower- and upper-bound solutions for the ultimate surcharge loading are obtained by applying finite element limit analysis techniques. For practical use, the results are presented in the form of dimensionless stability charts with the actual tunnel stability numbers being closely bracketed from above and below. As an additional check on the solutions, upper-bound rigid-block mechanisms have been developed and the predicted collapse loads from these are compared with those from finite element limit analysis. Finally, an expression that approximates the ultimate surcharge load has been devised which is convenient for use by practising engineers.     

Vertical Screening Distance Criteria to Evaluate Vapor Intrusion Risk from 1,2-Dichloroethane (1,2-DCA)

Vapor intrusion (VI) involves migration of volatile contaminants from subsurface through unsaturated soil into overlying buildings. In 2015, the US EPA recommended an approach for screening VI risks associated with gasoline releases from underground storage tank (UST) sites. Additional assessment of the VI risk from petroleum hydrocarbons was deemed unnecessary for buildings separated from vapor sources by more than recommended vertical screening distances. However, these vertical screening distances did not apply to potential VI risks associated with releases of former leaded gasoline containing 1,2-dichloroethane (1,2-DCA), because of a lack of empirical data on the attenuation of 1,2-DCA in soil gas. This study empirically evaluated 144 paired measurements of 1,2-DCA concentrations in soil gas and groundwater collected at 47 petroleum UST sites combined with BioVapor modeling. This included (1) assessing the frequency of 1,2-DCA detections in soil gas below 10⁻⁶ risk-based screening levels at different vertical separation distances and (2) comparing the US EPA recommended vertical screening distances with those predicted by BioVapor modeling. Vertical screening distances were predicted for different soil types using aerobic biodegradation rate constants estimated from the measured soil-gas data combined with conservative estimates of source concentrations. The modeling indicates that the vertical screening distance of 6 feet (1.8 m) recommended for dissolved-phase sources is applicable for 1,2-DCA below certain threshold concentrations in groundwater, while 15 feet (4.6 m) recommended for light nonaqueous phase liquid (LNAPL) sources is applicable for sites with clay and loam soils in the vadose zone, but not sand, if 1,2-DCA concentrations in groundwater exceed 150 μg/L. This dependence of the predicted vertical screening distances on soil type places added emphasis on proper soil characterization for VI screening at sites with 1,2-DCA sources. The soil-gas data suggests that a vertical screening distance of 15 feet (4.6 m) is necessary for both dissolved-phase and LNAPL sources.     

Rapid open-system assembly of a large silicic magma body: time-resolved evidence from cored plagioclase crystals in the Oruanui eruption deposits, New Zealand

Rhyolitic pumices in the 26.5 ka Oruanui eruption (Taupo volcano, New Zealand) contain an average of 10 wt% crystals. About 2 wt% of the crystal population is feldspar crystals that display bluish–grey cloudy cores, the colour being imparted by exsolved needles of rutile. The volume of cloudy-cored feldspars thus amounts to ~1.0 km³ in a total magma volume of ~530 km³. The cored feldspars show great variability in detail, but in general have a rounded cloudy core bounded by a zone rich in glass and mineral inclusions, that was then overgrown by a euhedral clear rim. Sr-isotopic variations in eight representative crystals were measured on micromilled samples of selected growth zones in the cores and rims, and linked to feldspar compositions through microprobe traverses. The cloudy cores range from ⁸⁷Sr/⁸⁶Sr = 0.70547 to 0.70657, with compositions of An₄₃ to An₇₈. The overgrowth rims display wider variations: inner parts show extreme ranges in composition (maxima ⁸⁷Sr/⁸⁶Sr = 0.70764 and An₇₈), while outer parts in seven of eight crystals are zoned, with outward-decreasing Sr-isotopic and An values to figures that are in accord with the bulk pumice and other, clear-feldspar values, respectively. The three parts of the crystals represent distinct regimes. The cloudy cores are inherited from an intermediate plutonic protolith that has been subjected to melting. The inner overgrowth rims were crystallised from a high temperature, relatively radiogenic melt derived from Mesozoic-Palaeozoic metasedimentary rocks (“greywacke”). The outer euhedral rims reflect mixing into and continued growth within the growing Oruanui magma body. The cloudy-cored feldspars also contain rare zircon inclusions. Twenty one zircons were recovered by HF digestion of a bulk sample of cloudy feldspars and analysed by SHRIMP for U–Th isotopes with which to calculate model ages. Eighteen of 21 crystals returned finite ages, the model-age spectrum of which is similar to the age spectra from free zircons in Oruanui pumices. Assembly of the Oruanui magma body was not only rapid (over ~40 kyr, as shown by other data) but involved a wide open system, with significant contributions from partly-melted intermediate-composition igneous intrusions (cloudy cores) and greywacke melts (inner overgrowths) being introduced into the magma body up to the point of eruption. Such open system behaviour contrasts with that proposed in models for comparably voluminous silicic magmas derived dominantly by fractionation (such as the Bishop Tuff) where the magma and its crystal cargo were better insulated thermally and chemically from country-rock interaction.     

Environmental impact of the 1.8 ka Taupo eruption, New Zealand: Landscape responses to a large-scale explosive rhyolite eruption

Large-scale ignimbrite eruptions from rhyolitic caldera volcanoes can trigger geologically instantaneous changes in sedimentary systems over huge areas by either burying existing environments or overloading them with vast quantities of unconsolidated particulate material. The post-eruption readjustment of the landscape to such perturbations is one of the most dramatic processes in physical sedimentology, exemplified here by the 1.8 ka Taupo eruption in the central North Island of New Zealand. This eruption generated voluminous fall deposits, then climaxed with emplacement of a c. 30 km³ non-welded ignimbrite over a near-circular area of c. 20 000 km². Approximately 90% of the area, but < 50% of the ignimbrite volume, is represented by a landscape-mantling unit that covered the pre-eruption topography to a depth varying from c. 10 m in proximal areas to less than 15–30 cm distally. The remainder of the ignimbrite deposit is represented by landscape-modifying material that ponded in valley bottoms and depressions to thicknesses of up to 70 m, with no systematic variation in thickness with distance from source.The headwaters of many of the North Island's largest rivers were impacted by both the primary pyroclastic fall and flow material. Large-scale post-eruption remobilisation of this material, coupled with the re-establishment of fluvial systems, occurred in a distinct sequence as recorded by the evolution of sedimentary facies in different sub-environments. Following an initial period dominated by mass flows, re-establishment of fluvial systems began with the headward erosion of box canyons through the ponded ignimbrite deposits, a process often associated with the break-out of temporary lakes. Aggradational streams developed in these channels rapidly evolved from shallow, ephemeral, sediment-laden outbursts associated with flash flood events to deeper, permanent braided rivers, before declining sediment yields led to retrenchment of single thread rivers and a return to pre-eruption gradients and bedloads years to decades later. Typically the modern profile of many streams and rivers follow closely their pre-eruption profiles, and incision and erosion is overwhelmingly confined to the deposits of the eruption itself.Although the general remobilisation pattern is similar for all impacted river systems, detailed studies of the Waikato, Rangitaiki, Mohaka, Ngaruroro and Whanganui catchments show that the relative timing and scale of each eruption response phase differs between each catchment. These reflect differences in catchment physiography and hydrology, and the volume and type of pyroclastic material deposited in each. Ultimately, the landscape response reflects the relative spatial distributions of, and the volumetric ratios between, the volumes of pyroclastic debris, water, and accommodation space in the basin (cf. Kataoka and Manville, this volume).     

Assessment of shallow groundwater in Lake Nyos catchment (Cameroon,Central-Africa): implications for hydrogeochemical controls and uses

In 1986, carbon dioxide gas exploded from Lake Nyos and killed about 1,800 people. After that disaster, various administrative and research activities have been conducted to mitigate subsequent disasters. However, none of those endeavors have characterized the groundwater chemistry to identify hydrogeochemical processes that control the water chemistry, and the quality of the water for domestic and agricultural uses that support the lives of un-official resettlers around Lake Nyos. Conventional hydrochemical techniques coupled with statistical and graphical analysis were therefore employed to establish the baseline hydrochemical conditions, assess processes controlling solutes distribution in shallow groundwater in the Lake Nyos catchment and explore its usability. Groundwater samples were analyzed for their physical and chemical properties. The wide ranges of electrical conductivity and total dissolved solid values reveal the heterogeneous distribution of groundwater within the watershed. The relative abundance of major dissolved species was Ca > Mg > Na > K for cations and HCO₃ >>> Cl > SO₄ > NO₃ for anions. Piper diagram classified almost all water samples into mixed CaMg–HCO₃ water type. Major ion geochemistry reveals that, in addition to silicates weathering (water–rock interaction), ion exchange processes regulate the groundwater chemistry. Principal component analysis supports the occurrence of water rock interaction. Hierarchical cluster analysis showed that the chemistry of groundwater in the study area is controlled by three main factors, and suggests no hydraulic connectivity between deep lake water and groundwater in the catchment. The quality assessment of the groundwater showed that groundwater parameters are within the acceptable limit of the World Health Organization and Nigeria guidelines for drinking and domestic uses, and water found to be good for irrigation.     

Mapping and analysing virtual outcrops

Laser scanning is a very efficient way to generate realistic, high-resolution digital models of 3-D geological outcrops. This paper discusses the methodologies involved in the creation and analysis of virtual outcrops, based on laser scanner data. The visualisation of the laser scanner data as a photorealistic 3-D object is described. Geological features picked out on the virtual outcrop (e.g. fractures, faults or bedding planes) can be extrapolated outward, into space, and inward, into the subsurface, using tension surfaces.Electronic Supplementary Material Supplementary material is available in the online version of this article at Reviewed by: J.D. Clemens, D. Yuen     

Structural vibration reduction using self-tuning fuzzy control of magnetorheological dampers

Among the control devices considered for dissipating seismic energy and reducing structural vibrations is the magnetorheological (MR) damper which consists of a hydraulic cylinder filled with a suspension of micron-sized, magnetically polarizable iron particles capable of reversibly changing from free-flowing, linear viscous fluid to semi-solid with the application of a magnetic field. Several algorithms have been proposed for regulating the amount of damping provided by MR dampers. An attractive option is the use of fuzzy controllers because they are simple, intrinsically robust, and they do not depend on a model of the system. Tuning of these controllers, however, has shown to be a difficult task because of the large number of parameters involved. This paper proposes a self-tuning fuzzy controller to regulate MR dampers’ properties and reduce structural responses of single degree-of-freedom seismically excited structures. Robustness to changes in seismic motions and structural characteristics was assessed by subjecting a rigid and a flexible building to different earthquake records. Results show that the self-tuning controller proposed effectively reduced responses of both structures to all earthquakes considered. In addition, results were compared to those of a fuzzy controller with constant scaling factors and to those of two passive strategies: “passive on” and “passive off”, where the current to the MR dampers was set to its maximum allowable value, and zero, respectively.     

Factors influencing the height of Hawaiian lava fountains: implications for the use of fountain height as an indicator of magma gas content

The heights of lava fountains formed in Hawaiian-style eruptions are controlled by magma gas content, volume flux and the amounts of lava re-entrainment and gas bubble coalescence. Theoretical models of lava fountaining are used to analyse data on lava fountain height variations collected during the 1983–1986 Pu'u 'O'o vent of Kilauea volcano, Hawaii. The results show that the variable fountain heights can be largely explained by the impact of variations in volume flux and amount of lava re-entrainment on erupting magmas with a constant gas content of 0.32 wt.% H₂O. However, the gas content of the magma apparently declined by 0.05 wt.% during the last 10 episodes of the eruption series and this decline is attributed to more extensive pre-eruption degassing due to a shallowing of the sub-vent feeder dike. It is concluded that variations in lava fountain height cannot be simply interpreted as variations in gas content, as has previously been suggested, but that fountain height can still be a useful guide to minimum gas contents. Where sufficient data are available on eruptive volume fluxes and extent of lava entrainment, greatly improved estimates can be made of magma gas content from lava fountain height.     

Internal erosion of soil pipes: Sediment rating curves for soil pipes

The collapse of soil pipes due to internal erosion can result in fully mature gullies. Few studies have measured the rates of sediment detachment and transport through soil pipes in situ. The objectives of this work were to determine suspended sediment concentration (SSC) in soil pipes as a function of pipeflow rate to develop sediment rating curves (SRC) and measure the bedload transport as a function of cumulative flow per storm event. H-flumes were installed in seven discontinuous gullies formed by pipe collapse and instrumented for pipe discharge measurements and suspended sediment sampling. The typical response to pipeflow was an initial flush of high concentration of suspended sediment followed by a decrease as pipeflow increased (rising limb of hydrograph). Pipeflows were often so dynamic that it was difficult to consistently capture the initial flush of sediment, resulting in weak to non-existent SRCs. The falling limb of the hydrograph tended to have a relatively low SSC. Thus, soil pipe SRCs tended to be better represented by hysteretic SRCs, although relationships between SSC and flow rate were poorly represented by SRCs. A power law equation given by SSC = aQ^b was adopted to represent the SRC relationships. Fitting this equation to data showed a correlation between the offset, a, and the slope, b, with the slope decreasing as the offset increases. Both SRC parameters (a and b) were correlated to the contributing area of the individual pipe. Bedload appeared to be an important contributor to sediment transport, with bedload – expressed as an average event sediment concentration (mg l⁻¹) – decreasing as the volume of the event discharge (m³) increased. A significant portion (11–31%) of the bedload material was gravel and aggregates (>2 mm diameter material). While this work was the first to determine SRCs for soil pipes, refined sampling and measurement techniques are needed. © 2020 John Wiley & Sons, Ltd.