The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

福建省第三产业发展的区域比较分析

运用可拓工程方法和主成分分析方法建立了福建省第三产业发展水平的综合评判物元模型。模型的输出结果表明：福建省各地区第三产业发展水平差距悬殊；综合关联度将第三产业发展水平划分成4个等级，单指标关联度则揭示了各地区第三产业发展的优、劣势。     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Une approche innovante pour modéliser la biodégradation des composés organochlorés volatils en aquifères poreux

The goal of this study was to develop an innovative chloroethene biodegradation module based on biological, thermodynamical and mechanistic concepts. The biodegradation scheme was based on the postulate that in each part of an aquifer only one degradation mechanism is dominant: the one involving the most energetic electron acceptor. Thus, the selection of the active degradation mechanism was a function of the concentration of different electron acceptors. Modified Monod-type kinetics was used in order to take into account the possible influence of some compounds on the biodegradation of a given organic compound. The numerical model developed was applied to a simple test case, whose results are presented here. To cite this article: F. Nex et al., C. R. Geoscience 338 (2006).     

Assessment of sand encroachment in Kuwait using GIS

Assessment of sand encroachment in Kuwait using Geographical Information System (GIS) technology has been formulated as a Multi-Criteria Decision Making problem. The Delphi method and Analytical Hierarchy Process were adopted as evaluating techniques, in which experts’ judgments were analyzed for objectively estimating and weighting control factors. Seven triggering factors, depicted in the form of maps, were identified and ordered according to their priority. These factors are (1) wind energy; (2) surface sediment; (3) vegetation density; (4) land use; (5) drainage density; (6) topographic change and (7) vegetation type. The factor maps were digitized, converted to raster data and overlaid to determine their possible spatial relationships. Applying a susceptibility model, a map of sand encroachment susceptibility in Kuwait was developed. The map showed that the areas of very high and high sand encroachment susceptibility are located within the main corridor of sand pathway that coincides with the northwesterly dominant wind direction.     

龙永煤田危机矿山资源接替的几种勘探模式

通过龙永煤田现有煤矿状况调查分析发现：许多矿井已步入“花甲之年”，接近矿井服务年限。为了保障我省能源安全，应当采用有效勘探方法，方能使它们“返老还童”。研究发现危机矿山延长服务年限的几种勘探模式是：其一为扩大（延伸）勘探模式；其二为小型井田（块段）勘探模式；其三为资料分析研究模式。经探采对比发现，采用这三种勘探模式后，产生了巨大的经济效益，使一大批老矿井青春再度焕发。     

单一滑面边坡的非稳定蠕变模型分析

对单一滑面边坡的非稳定蠕变分析是边坡稳定性研究的基础。文中对岩质边坡单一滑面的流变变形机制进行了分析，提出了一种能较好反映这一变形机制的非稳定蠕变模型。特别是描述单一滑面边坡的非稳定蠕变破坏阶段，并对其稳定性进行了讨论。结合一些工程实例进行了对比验证，为滑坡灾害的预测和防治提供依据。