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.     

Carbon dioxide and radon gas hazard in the Alban Hills area (central Italy)

The sudden and catastrophic, or slow and continuous, release at surface of naturally occurring toxic gases like CO₂, H₂S and Rn poses a serious health risk to people living in geologically active regions. In general this problem receives little attention from local governments, although public concern is raised periodically when anomalous toxic-gas concentrations suddenly kill humans or livestock. For example, elevated CO₂ concentrations have been linked to the death of at least 10 people in the central Italian region of Lazio over the last 20 years, while it was the CO₂ asphyxiation of 30 cows in a heavily populated area near Rome in 1999 which prompted the present soil-gas study into the distribution of the local health risk. A detailed geochemical survey was carried out in an area of about 4 km² in the Ciampino and Marino districts, whereby a total of 274 soil-gas samples were collected and analysed for more than 10 major and trace gas species. Data were then processed using both statistical and geostatistical methods, and the resulting maps were examined in order to highlight areas of elevated risk. General trends of elevated CO₂ and Rn concentrations imply the presence of preferential pathways (i.e. faults and fractures) along which deep gases are able to migrate towards the surface. The CO₂ and Rn anomalous trends often correspond to and are usually elongated parallel to the Apennine mountain range, the controlling structural feature in central Italy. Because of this fundamental anisotropy in the factors controlling the soil-gas distribution, it was found that a geostatistical approach using variogram analysis allowed for a better interpretation of the data. With regard to the health risk to local inhabitants, it was found that although some high risk areas had been zoned as parkland, others had been heavily developed for residential purposes. For example, many new houses were found to have been built on ground which has soil-gas CO₂ concentrations of more than 70% and radon values of more than 250 kBq m⁻³. It is recommended that land-use planners incorporate soil-gas and/or gas flux measurements in environmental assessments in areas of possible risk (i.e. volcanic or structurally active areas).     

An Mw based earthquake Catalogue for central, northern and northwestern Europe using a hierarchy of magnitude conversions

Data from 25 local catalogues and 30special studies of earthquakes in central,northern and northwestern Europe have beenincorporated into a Databank. The dataprocessing includes discriminating eventtypes, eliminating fake events and dupletsand converting different magnitudes andintensities to M_w if this is not givenby the original source. The magnitudeconversion is a key task of the study andimplies establishment of regressionequations where no local relations exist.The Catalogue contains tectonic events fromthe Databank within the area44°N–72°N,25°W–32°E and the time period1300–1993. The lower magnitude level forthe Catalogue entries is setat M_w == 3.50. The area covered by thedifferent catalogues are associated withpolygons. Within each polygon only datafrom one or a small number of the localcatalogues, supplemented by data fromspecial studies, enter the Catalogue. Ifthere are two or more such catalogues orstudies providing a solution for an event,a priority algorithm selects one entry forthe Catalogue. Then M_w is calculatedfrom one of the magnitude types, or frommacroseismic data, given by the selectedentry according to another priority scheme.The origin time, location, M_w magnitude and reference are specified for eachentry of the Catalogue. So is theepicentral intensity, I₀, if providedby the original source. Following thesecriteria, a total of about 5,000earthquakes constitute the Catalogue.Although originally derived for the purposeof seismic hazard calculation within GSHAP,the Catalogue provides a data base for manytypes of seismicity and seismic hazardstudies.     

Regional Stress Field and Seismic Dynamics Along the Border Zone Between Fujian, Guangdong and Jiangxi Provinces

The regional stress field and seismic dynamics along the border zone between Fujian, Guangdong and Jiangxi Provinces are studied based on the seismo-geological data, GPS measurement, and seismicity. The results show that: (1) the principal compressional stress of the stress field is oriented in NW-SE direction and the principal extensional stress is in NE-SW direction; (2) the WNW-ward compression and collision of the Philippine Sea Plate to the eastern coast of Taiwan Island are the most direct and most important dynamic source for preparation and occurrence of strong earthquakes in the Taiwan area and along the border zone between Fujian, Guangdong and Jiangxi Provinces.     

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

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