Sediment deposition from flow at low gradients into a buffer strip—a critical test of re-entrainment theory

The spatial and size distribution of sediment deposited from short periods of overland flow due to the effect of a simulated grass buffer strip was measured for low slopes of 1.6, 3.4 and 5.1%. These data were analysed so as to critically evaluate two alternative models of the process of re-entrainment of recently deposited sediment. A model of re-entrainment, previously thought to be appropriate only for a steady-state or equilibrium situation, was found to give better agreement with experiments than did a model previously used in the literature on this subject.     

Drainage in finite-sized unsaturated zones

After the initiation of gravity drainage, water is often assumed to be either (a) draining under unit gradient, or (b) at capillary/gravity equilibrium. Both of these simplifications can be useful, but the regimes of validity of each assumption must be delineated. Water pressures are measured versus time and distance as water drains out of a 1.6 m long sand column to determine the relative effects of capillary and gravitational forces during drainage. For medium sized sands (0.15–0.3 mm in diameter), the capillary pressure is constant in space in a large region of the column for over 12 days, and the water continues to flow under unit gradient for relatively long time scales. Similar results are seen for finer sands, but with a much faster approach to equilibrium. Numerical simulations and analytical estimates are presented and compare favorably to the measurements. Together, the experimental, theoretical and analytical results are used to calculate when capillary/gravity equilibrium is reached as a function of porous media properties and length of the unsaturated zone. The ratio of the length of the unsaturated zone to the bubbling pressure is a key parameter in determining the drainage regime, and that even for relatively short unsaturated zones the equilibrium time scale can be on the order of years.     

Parameters affecting maximum fluid transport in large aperture fractures

We present results of laboratory experiments to study the behavior of liquids moving in unsaturated wide-aperture fractures. A 5-mm-thick glass plate cut with a 1.7-mm aperture was used as a fractured rock analog to study behavior of film and capillary droplet flow modes. Flow rates ranged between 0.6 and 6.0 ml/min. Variability in the ambient barometric pressure, resulting from weather conditions, seemed to play a role in the natural selection of flow mode. For droplet mode, constant input conditions resulted in highly variable transport properties within the fracture. The advancing meniscus exhibited a dynamic contact angle that was a function of the droplet speed and much larger than the static contact angle. Flow rate was reduced due to the larger contact angle. Analytical expressions for droplet velocity and flow capacity are presented as a function of the dynamic rather than the static contact angle.     

Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.     

Flow and fracture maps for basaltic rock deformation at high temperatures

Understanding how the strength of basaltic rock varies with the extrinsic conditions of stress state, pressure and temperature, and the intrinsic rock physical properties is fundamental to understanding the dynamics of volcanic systems. In particular it is essential to understand how rock strength at high temperatures is limited by fracture. We have collated and analysed laboratory data for basaltic rocks from over 500 rock deformation experiments and plotted these on principal stress failure maps. We have fitted an empirical flow law (Norton’s law) and a theoretical fracture criterion to these data. The principal stress failure map is a graphical representation of ductile and brittle experimental data together with flow and fracture envelopes under varying strain rate, temperature and pressure. We have used these maps to re-interpret the ductile–brittle transition in basaltic rocks at high temperatures and show, conceptually, how these failure maps can be applied to volcanic systems, using lava flows as an example.     

Exploration and discovery in Yellowstone Lake: results from high-resolution sonar imaging, seismic reflection profiling, and submersible studies

‘No portion of the American continent is perhaps so rich in wonders as the Yellow Stone’ (F.V. Hayden, September 2, 1874)Discoveries from multi-beam sonar mapping and seismic reflection surveys of the northern, central, and West Thumb basins of Yellowstone Lake provide new insight into the extent of post-collapse volcanism and active hydrothermal processes occurring in a large lake environment above a large magma chamber. Yellowstone Lake has an irregular bottom covered with dozens of features directly related to hydrothermal, tectonic, volcanic, and sedimentary processes. Detailed bathymetric, seismic reflection, and magnetic evidence reveals that rhyolitic lava flows underlie much of Yellowstone Lake and exert fundamental control on lake bathymetry and localization of hydrothermal activity. Many previously unknown features have been identified and include over 250 hydrothermal vents, several very large (>500 m diameter) hydrothermal explosion craters, many small hydrothermal vent craters (1–200 m diameter), domed lacustrine sediments related to hydrothermal activity, elongate fissures cutting post-glacial sediments, siliceous hydrothermal spire structures, sublacustrine landslide deposits, submerged former shorelines, and a recently active graben. Sampling and observations with a submersible remotely operated vehicle confirm and extend our understanding of the identified features. Faults, fissures, hydrothermally inflated domal structures, hydrothermal explosion craters, and sublacustrine landslides constitute potentially significant geologic hazards. Toxic elements derived from hydrothermal processes also may significantly affect the Yellowstone ecosystem.     

Earthen barriers to control lava flows in the 2001 eruption of Mt. Etna

Preceded by four days of intense seismicity and marked ground deformation, a new eruption of Mt. Etna started on 17 July and lasted until 9 August 2001. It produced lava emission and strombolian and phreatomagmatic activity from four different main vents located on a complex fracture system extending from the southeast summit cone for about 4.5 km southwards, from 3000 to 2100 m elevation (a.s.l.). The lava emitted from the lowest vent cut up an important road on the volcano and destroyed other rural roads and a few isolated country houses. Its front descended southwards to about 4 km distance from the villages of Nicolosi and Belpasso. A plan of intervention, including diversion and retaining barriers and possibly lava flow interruption, was prepared but not activated because the flow front stopped as a consequence of a decrease in the effusion rate. Extensive interventions were carried out in order to protect some important tourist facilities of the Sapienza and Mts. Silvestri zones (1900 m elevation) from being destroyed by the lava emitted from vents located at 2700 m and 2550 m elevation. Thirteen earthen barriers (with a maximum length of 370 m, height of 10–12 m, base width of 15 m and volume of 25 000 m³) were built to divert the lava flow away from the facilities towards a path implying considerably less damage. Most of the barriers were oriented diagonally (110–135°) to the direction of the flow. They were made of loose material excavated nearby and worked very nicely, resisting the thrust of the lava without any difficulty. After the interventions carried out on Mt. Etna in 1983 and in 1991–1992, those of 2001 confirm that earthen barriers can be very effective in controlling lava flows.     

Tensor time domain electromagnetic resistivity measurements at Ngatamariki geothermal field, New Zealand

Experimental measurements in the Ngatamariki geothermal field, North Island, New Zealand were made to test the applicability of the time domain electromagnetic method for detailed investigation of the resistivity structure within a geothermal field. Low-frequency square wave signals were transmitted through three grounded bipole current sources sited about 8 km from the measurement lines. Despite high levels of electrical noise, transient electric field vectors could be determined reliably for times between 0.02 and 3.3 s after each step in the source current. Instantaneous apparent resistivity tensors were then calculated. Apparent resistivity pseudosections along the two measurement lines show smooth variations of resistivity from site to site. Over most of the field the images consistently show a three-layer resistivity structure with a conductive middle layer (3–10 Ωm) representing the conductive upper part of the thermal reservoir. A deep-seated region of low resistivity in the northwest of the field may indicate a conductive structure at about 1 km associated with a deeper diorite intrusion. Measurements sited closer than about 100 m to drillholes appear to have been disturbed by metallic casing in the holes. A change in resistivity structure in the east of the field may indicate a major geological or hydrothermal boundary.     

哈尔滨市磨盘山水利枢纽区三维渗流数值模拟模型研究

哈尔滨市磨盘山水库防渗型式拟采用垂直砼防渗墙及帷幕灌浆方案。为了分析左岸不同长度防渗对水库渗漏量及坝基渗透稳定性的影响。对磨盘山水利枢纽区渗流场进行了三维数值模拟研究。结果表明，随着左岸防渗长度的增大。渗漏量及下游渗流出口渗透比降基本呈线性趋势减小。当左岸防渗长度为200m时。渗漏量及各渗透比降均小于允许值，为推荐的渗控方案。     

The submarine eruption of the Bombarda volcano,Milos Island,Cyclades, Greece

The Milos volcanic field includes a well-exposed volcaniclastic succession which records a long history of submarine explosive volcanism. The Bombarda volcano, a rhyolitic monogenetic center, erupted ∼1.7 Ma at a depth <200 m below sea level. The aphyric products are represented by a volcaniclastic apron (up to 50 m thick) and a lava dome. The apron is composed of pale gray juvenile fragments and accessory lithic clasts ranging from ash to blocks. The juvenile clasts are highly vesicular to non-vesicular; the vesicles are dominantly tube vesicles. The volcaniclastic apron is made up of three fades: massive to normally graded pumice-lithic breccia, stratified pumice-lithic breccia, and laminated ash with pumice blocks. We interpret the apron beds to be the result of water-supported, volcaniclastic mass-How emplacement, derived directly from the collapse of a small-volume, subaqueous eruption column and from syn-eruptive, down-slope resedimentation of volcaniclastic debris. During this eruptive phase, the activity could have involved a complex combination of phreatomagmatic explosions and minor submarine effusion. The lava dome, emplaced later in the source area, is made up of flow-banded lava and separated from the apron by an obsidian carapace a few meters thick. The near-vertical orientation of the carapace suggests that the dome was intruded within the apron. Remobilization of pyroclastic debris could have been triggered by seismic activity and the lava dome emplacement. Published online: 30 January 2003 Editorial responsibility: J. McPhie