A laboratory model of surface crust formation and disruption on lava flows through non-uniform channels

Crust formation on basaltic lava flows dictates conditions of both flow cooling and emplacement. For this reason, flow histories are dramatically different depending on whether lava is transported through enclosed lava tubes or through open channels. Recent analog experiments in straight uniform channels (Griffiths et al. J Fluid Mech 496:33–62, 2003) have demonstrated that tube flow, dictated by a stationary surface crust, can be distinguished from a mobile crust regime, where a central solid crust is separated from channel walls by crust-free shear zones, by a simple dimensionless parameter ϑ, such that ϑ<25 produces tube flow and ϑ>25 describes the mobile crust regime. ϑ combines a previously determined parameter ψ, which describes the balance between the formation rate of surface solid and the shear strain that disrupts the solid crust, with the effects of thermal convection (described by the Rayleigh number Ra).Here we explore ways in which ϑ can be used to describe the behavior of basaltic lava channels. To do this we have extended the experimental approach to examine the effects of channel irregularities (expansions, contractions, sinuosity, and bottom roughness) on crust formation and disruption. We find that such changes affect local flow behavior and can thus change channel values of ϑ. For example, gradual widening of a channel results in a decrease in flow velocity that causes a decrease in ϑ and may allow a down-flow transition from the mobile crust to the tube regime. In contrast, narrowing of the channel causes an increase in flow velocity (increasing ϑ), thus inhibiting tube formation.We also quantify the fraction of surface covered by crust in the mobile crust regime. In shallow channels, variations in crust width (d _c) with channel width (W) are predicted to follow d _c∼W ^5/3. Analysis of channelized lava flows in Hawaii shows crustal coverage consistent with this theoretical result along gradually widening or narrowing channel reaches. An additional control on crustal coverage in both laboratory and basaltic flows is disruption of surface crust because of flow acceleration through constrictions, around bends, and over breaks in slope. Crustal breakage increases local rates of cooling and may cause local blockage of the channel, if crusts rotate and jam in narrow channel reaches. Together these observations illustrate the importance of both flow conditions and channel geometry on surface crust development and thus, by extension, on rates and mechanisms of flow cooling. Moreover, we note that this type of analysis could be easily extended through combined use of FLIR and LiDAR imaging to measure crustal coverage and channel geometry directly.Editorial responsibility: A. Harris     

Viscous Newtonian laminar flow in a rectangular channel: application to Etna lava flows

 We introduce a 3D model for near-vent channelized lava flows. We assume the lava to be an isothermal Newtonian liquid flowing in a rectangular channel down a constant slope. The flow velocity is calculated with an analytical steady-state solution of the Navier-Stokes equation. The surface velocity and the flow rate are calculated as functions of the flow thickness for different flow widths, and the results are compared with those of a 2D model. For typical Etna lava flow parameters, the influence of levees on the flow dynamics is significant when the flow width is less than 25 m. The model predicts the volume flow rate corresponding to the surface velocity, taking into account that both depend on flow thickness. The effusion rate is a critical parameter to evaluate lava flow hazard. We propose a model to calculate the effusion rate given the lava flow width, the topograhic slope, the lava density, the surface flow velocity, and either the lava viscosity or the flow thickness. Received: 20 January 1998 / Accepted: 8 January 1999     

Simulations of the 2004 lava flow at Etna volcano using the magflow cellular automata model

Since the mechanical properties of lava change over time, lava flows represent a challenge for physically based modeling. This change is ruled by a temperature field which needs to be modeled. MAGFLOW Cellular Automata (CA) model was developed for physically based simulations of lava flows in near real-time. We introduced an algorithm based on the Monte Carlo approach to solve the anisotropic problem. As transition rule of CA, a steady-state solution of Navier-Stokes equations was adopted in the case of isothermal laminar pressure-driven Bingham fluid. For the cooling mechanism, we consider only the radiative heat loss from the surface of the flow and the change of the temperature due to mixture of lavas between cells with different temperatures. The model was applied to reproduce a real lava flow that occurred during the 2004–2005 Etna eruption. The simulations were computed using three different empirical relationships between viscosity and temperature.     

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.     

Voluminous lava flows at Oldoinyo Lengai in 2006: chronology of events and insights into the shallow magmatic system

The largest natrocarbonatite lava flow eruption ever documented at Oldoinyo Lengai, NW Tanzania, occurred from March 25 to April 5, 2006, in two main phases. It was associated with hornito collapse, rapid extrusion of lava covering a third of the crater and emplacement of a 3-km long compound rubbly pahoehoe to blocky aa-like flow on the W flank. The eruption was followed by rapid enlargement of a pit crater. The erupted natrocarbonatite lava has high silica content (3% SiO₂). The eruption chronology is reconstructed from eyewitness and news media reports and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data, which provide the most reliable evidence to constrain the eruption’s onset and variations in activity. The eruption products were mapped in the field and the total erupted lava volume estimated at 9.2 ± 3.0 × 10⁵ m³. The event chronology and field evidence are consistent with vent construct instability causing magma mixing and rapid extrusion from shallow reservoirs. It provides new insights into and highlights the evolution of the shallow magmatic system at this unique natrocarbonatite volcano.     

Kinematic significance of mingling-rolling structures in lava flows: a case study from Porri Volcano (Salina, Southern Tyrrhenian Sea)

 A basaltic andesite lava flow from Porri Volcano (Salina, Southern Tyrrhenian Sea) is composed of two different magmas. Magma A (51 vol.% of crystals) has a dacitic glass composition, and magma B (18 vol.% of crystals), a basaltic glass composition. Magma B is hosted in A and consists of sub-spherical enclaves and boudin-like, banding and rolling structures (RS). Four types of RS have been recognized: σ–type;δ–type; complex σ-δ–types and transitional structures between sub-spherical enclaves and rolling structures. An analysis of the RS has been performed in order to reconstruct the flow kinematics and the mechanism of flow emplacement. Rolling structures have been selected in three sites located at different distances from the vent. In all sites most RS show the same sense of shear. Kinematic analysis of RS allows the degree of flow non-coaxiality to be determined. The non-coaxiality is expressed by the kinematic vorticity number Wk, a measure of the ratio Sr between pure shear strain rate and simple shear strain rate. The values of Wk calculated from the measured shapes of microscopic RS increase with increasing distance from the vent, from approximately 0.5 to 0.9. Results of the structural analysis reveal that the RS formed during the early–intermediate stage of flow emplacement. They represent originally sub-spherical enclaves deformed at low shear strain. At higher strain, RS deformed to give boudin-like and stretched banding structures. Results of the kinematic analysis suggest that high viscosity lava flows are heterogeneous non-ideal shear flows in which the degree of non-coaxiality increases with the distance from the vent. In the vent area, deformation is intermediate between simple shear and pure shear. Farther from the vent, deformation approaches ideal simple shear. Lateral extension processes occur only in the near-vent zone, where they develop in response to the lateral push of magma extruded from the vent. Lateral shortening processes develop in the distal zone and record the gravity-driven movement of the lava. The lava flow advanced by two main mechanisms, lateral translation and rolling motion. Lateral translation equals rolling near the vent, while rolling motion prevailed in the distal zones. Received: 6 November 1997 / Accepted: 29 November 1997     

Prediction and assessment of the disturbances of the coal mining in Kailuan to karst groundwater system

Coal resources and water resources play an essential and strategic role in the development of China's social and economic development, being the priority for China's medium and long technological development. As the mining of the coal extraction is increasingly deep, the mine water inrush of high-pressure confined karst water becomes much more a problem. This paper carried out research on the hundred-year old Kailuan coal mine's karst groundwater system. With the help of advanced Visual Modflow software and numerical simulation method, the paper assessed the flow field of karst water area under large-scale exploitation. It also predicted the evolution ofgroundwaterflow field under different mining schemes of Kailuan Corp. The result shows that two cones of depression are formed in the karst flow field of Zhaogezhuang mining area and Tangshan mining area, and the water levels in two cone centers are −270 m and −31 m respectively, and the groundwater generally flows from the northeast to the southwest. Given some potential closed mines in the future, the mine discharge will decrease and the water level of Ordovician limestone will increase slightly. Conversely, given increase of coal yield, the mine drainage will increase, falling depression cone of Ordovician limestone flow field will enlarge. And in Tangshan's urban district, central water level of the depression cone will move slightly towards north due to pumping of a few mines in the north.     

Application of grey correlation method to evaluate potential groundwater recharge sites

Artificial recharge is a practical tool available for increasing the groundwater storage capacity. The efficiency of artificial recharge is related to various hydrogeological factors of the target area. In this study, a variable saturated groundwater flow model, FEMWATER, was used to evaluate the arrival times of recharged water that infiltrates from an artificial recharge pond to the groundwater table under various hydrogeological conditions. Forty-five arrival times were generated by FEMWATER. The relationships between the arrival times and hydrogeological factors used in the simulation of FEMWATER were analyzed by the grey correlation method. The results show the order of importance of the factors as they influence the arrival time. In order from high to low importance, they are α, D _g, θ _e, D _p, K _S and β. D _g and D _p are interpreted as the potential for movement of the recharge water; θ_e is the water storage capacity of soil, and K _S represents the ability of soil to transport water. α and β describe the characteristic curve of the unsaturated soil. The method was applied to evaluate a suitable site for artificial recharge in the Yun-Lin area. Grey correlation analysis was performed to obtain the grey correlation grade using the minimum arrival time as a reference sequence. An index is proposed herein to determine the recharge efficiency of 20 sampling sites. A contour mapping of index values at the 20 sampling sites identified three areas for artificial aquifer recharge in Yun-Lin. Area A in the upper plain is considered more appropriate for groundwater recharge than areas B and C in the coast.     

井孔水温对远场巨震同震响应及其机制的数模研究

热量的流动必然伴随着温度的改变, 地震活动期间的井孔水温同震响应也服从热力学的基本规律。 本文以汶川8.0级地震和日本9.0级地震为例, 以海口ZK26井水温同震响应的观测数据为基础, 运用热力学传导方程的数值模拟方法, 研究了同震响应过程中井孔水温度变化与热量传导之间的关系及相关的热力学机制问题.通过正演的方法, 得出沿井深方向不同时刻水温同震响应变化的数值模拟曲线, 水温数值模拟曲线与实际测量数据一致. 研究结果表明, 井孔中水温同震响应方式(上升、 下降或不变), 与水温传感器的位置、 热源的位置和分布、 传感器与热源之间的相对空间等因素有关。     

井孔水温对远场巨震同震响应及其机制的数模研究

热量的流动必然伴随着温度的改变,地震活动期间的井孔水温同震响应也服从热力学的基本规律。本文以汶川8.0级地震和日本福岛9.0级地震为例,以海口ZK26井水温同震响应的观测数据为基础,运用热力学传导方程的数值模拟方法,研究了同震响应过程中井孔水温度变化与热量传导之间的关系及相关的热力学机制问题。通过正演的方法,得出沿井深方向不同时刻水温同震响应变化的数值模拟曲线,水温数值模拟曲线与实际测量数据一致。研究结果表明,井孔中水温同震响应方式（上升、下降或不变）,与水温传感器的位置、热源的位置和分布、传感器与热源之间的相对空间等因素有关。     

Infiltration through series of soil aggregates: Neutron radiography and modeling

Soils are often structured as fine-porous aggregates separated by large inter-aggregate pores. Under unsaturated conditions, water is mostly stored in the aggregates and water flow depends on the properties of the aggregates as well as on those of the contacts between aggregates. The goal of this study is to model and evaluate the hydraulic properties of the contacts. We used neutron radiography to monitor the infiltration of water through series of aggregates. The flow process was numerically simulated by considering the hydraulically conducting contact area between aggregates as a variable that depends on the capillary pressure. This contact area was evaluated by matching the observed and simulated water flow across aggregates. We determined the conductivity of the contacts assuming that it scales with the contact area. We also measured the equivalent conductivity of series of aggregates. We found that during drainage the hydraulically conducting contact area drastically decreases and the conductivity of the contacts becomes much smaller than that of the aggregates. We also found that the equivalent conductivity of the aggregate series decreases as the conductivity of the contacts. We concluded that the contacts control the flow: they are highly conductive when wet, but act as bottle-necks under drained conditions. The abrupt transition between these two limiting cases indicates that the contact region is more rapidly drained than the aggregate interior. Our findings might be extended to describe water transfer to an evaporating boundary, infiltration and storage of water in unsaturated aggregated soils.     

Changes in the epipelic diatom assemblage in nutrient rich streams due to the variations of simultaneous stressors

Benthic diatoms are often used for assessing environmental conditions, such as water quality and habitat conditions in stream and river systems. Although laboratory experiments have shown that each diatom species have different levels of tolerance to different stressors, few studies have been conducted in laboratory settings that analyze the responses of the diatom assemblage to the effects of multiple simultaneous variables. The aim of this study was to evaluate some structural responses (such as species composition and diversity) of the diatom assemblage on a short time scale to the effects of the simultaneous increase in four variables that are directly linked to the environmental changes affecting the Pampean streams: turbidity, nutrients (phosphorous and nitrogen), water velocity and temperature. To this end we conducted a five-week laboratory experiment using artificial channels where we simulated two environmental conditions (LOW and HIGH) employing epipelic biofilm from a mesotrophic stream. The results obtained in the experiment show that the structure of the diatom assemblage in the epipelic biofilm is affected by the simultaneous modification of temperature, water velocity, nutrient concentration and turbidity. These modifications in the assemblage included moderate decreases in diversity, small decreases in the proportion of species sensitive to eutrophication and saprobity, moderate increases in the IDP (Pampean Diatom Index) values and moderate changes in the percentages of the stalked growth-forms. The relative abundance of species such as Luticola mutica, Navicula cryptocephala and Navicula lanceolata were negatively affected by both treatments; other species such as Planothidium lanceolatum, Caloneis bacillum, Encyonema minutum, Humidophila contenta, Luticola kotschyi, Nitzschia amphibia, Navicula veneta, Pinnularia subcapitata var. subcapitata were positively affected by the HIGH treatment; and Nitzschia fonticola was positively affected by both treatments. The results suggest that, in the very short term of the bioassay conducted, the diatom assemblage can modify its structure to respond in a sensitive manner to the abrupt changes in multiple physical–chemical variables.     

太平洋板块俯冲对中国东北深浅震影响机理的数值模拟

本文采用与深度有关的不同分层结构模型,并考虑太平洋板块俯冲角度差异等特征,建立太平洋板块向中国东北地区俯冲的2D纵向静力学模型.以太平洋板片俯冲速度为约束条件,通过变化俯冲板块的俯冲角度,数值模拟太平洋板块向中国东北的俯冲过程,探讨太平洋板块俯冲作用对于我国东北深浅震的影响,得到不同俯冲角度模型的深浅部应力场分布,揭示区域构造应力场的总体特征和断裂带局部特性,并讨论了活动断裂带以及邻近区域对东北地区深浅部构造应力场的响应.结果表明,太平洋板块俯冲角度的变化对于中国东北地区深浅震的地震活动格局具有重要影响,断裂带构造环境是浅源地震孕育的基本条件.     

电力系统地震功能失效数值模拟与控制分析

本文研究了电力供应系统地震功能失效评估方法,尤其对电力供应系统功能潮流分析方法与负荷的控制问题作了深入的探讨.不仅分析了供电系统的潮流(包括电压和功率)变化对供电系统的影响问题,而且建立了电流控制模型,利用DSEA控制算法对供电系统进行负荷控制,并在此基础上进行了系统过负荷时的数值模拟研究;另外,通过一实例分析,对本文的研究方法进行了验证,从而为供电系统抗震分析、功能评估及灾害预防提供理论基础.     

From the pore scale to the lab scale: 3-D lab experiment and numerical simulation of drainage in heterogeneous porous media

A well-controlled 3-D experiment with pre-defined block heterogeneities is conducted, where neutron tomography is used to map 3-D water distribution after two successive drainage steps. The material and hydraulic properties of the two sands are first measured in the laboratory with multistep outflow experiments. Additionally, the pore structure of the sands is acquired by means of image analysis of synchrotron tomography data and the structure is used for pore-scale simulation of one- and two-phase flow with Lattice-Boltzmann methods. This gives us another set of material and hydraulic parameters of the sands. The two sets of hydraulic properties (from the lab scale and from the pore scale) are then used in numerical simulations of the 3-D experiment.     

A sensitivity analysis of lake water level response to changes in climate and river regimes

Lake water level regimes are influenced by climate, hydrology and land use. Intensive land use has led to a decline in lake levels in many regions, with direct impacts on lake hydrology, ecology and ecosystem services. This study examined the role of climate and river flow regime in controlling lake regimes using three different lakes with different hydraulic characteristics (volume-inflow ratio, CIR). The regime changes in the lakes were determined for five different river inflows and five different climate patterns (hot-arid, tropical, moderate, cold-arid, cold-wet), giving 75 different combinations of governing factors in lake hydrology. The input data were scaled to unify them for lake comparisons. By considering the historical lake volume fluctuations, the duration (number of months) of lake volume in different ‘wetness’ regimes from ‘dry’ to ‘wet’ was used to develop a new index for lake regime characterisation, ‘Degree of Lake Wetness’ (DLW). DLW is presented as two indices: DLW₁, providing a measure of lake filling percentage based on observed values and lake geometry, and DLW₂, providing an index for lake regimes based on historical fluctuation patterns. These indices were used to classify lake types based on their historical time series for variable climate and river inflow. The lake response time to changes in hydrology or climate was evaluated. Both DLW₁ and DLW₂ were sensitive to climate and hydrological changes. The results showed that lake level in high CIR systems depends on climate, whereas in systems with low CIR it depends more on river regime.