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     

Accelerating ground warming in the Canadian Prairie provinces: Is it a result of global warming?

The results of precision temperature logs made to depths of several hundred meters in some 80 wells in Western Canada, most of which are located in the Prairie Provinces, show evidence of warming at the ground surface in the 0.5 K to 3.5 K range (average=2.2±0.7 K, for 80 unevenly distributed sites). Modeling shows that this warming mostly pertains to this century and it has been most substantal in the last four decades if the ramp function of the linear increase of surface temperature is assumed. Using the step function model's increase of surface temperature (land clearing, forest fires, etc.) the calculated onset of warming would pertain mostly to the last two decades. Contour maps of ground temperatures currently and previously and a contour map of the ground warming magnitude dilineate a large regional character of the ground temperature change at the southern marigin of permafrost for the large area of the Prairie Provinces. In many cases however, the magnitude of ground warming is much larger than the magnitude of air warming. This is especially evident for the northern areas of Alberta in the boreal forest ecoprovince. The magnitude of ground warming is equal to the magnitude of surface air warming in southern Alberta in the grassland and aspen parkland ecoprovinces. The analysis of the temperature depth response to the surface warming from well data shows the integrated effect of surface air warming together with the increases in ground temperature due to natural terrain effects and other anthropogenical changes to the surface of the earth.     

Evolution of cool skin and direct air-sea gas transfer coefficient during daytime

Absorption of solar radiation within the thermal molecular sublayer of the ocean can modify the temperature difference across the cool skin as well as the air-sea gas transfer. Our model of renewal type is based on the assumption that the thermal and diffusive molecular sublayers below the ocean surface undergo cyclic growth and destruction, the heat and gas transfer between the successive burst events are performed by molecular diffusion. The model has been upgraded to include heating due to solar radiation. The renewal time is parameterized as a function of the surface Richardson number and the Keulegan number. A Rayleigh number criterion characterizes the convective instability of the cool skin under solar heating. Under low wind speed conditions, the solar heating can damp the convective instability, strongly increasing the renewal time and correspondingly decreasing the interfacial gas exchange. In the ocean, an additional convective instability caused by salinity flux due to evaporation becomes of importance in such cases. The new parameterization is compared with the cool skin data obtained in the western equatorial Pacific during the Tropical Ocean Global Atmosphere Coupled Ocean Atmosphere Response Experiment in February 1993. In combination with a model of the diurnal thermocline it describes main features of the field data both in nighttime and daytime. Under low wind speed conditions (< 5 m s^-1) diurnal variations of the sea surface temperature due to the formation of a diurnal thermocline were substantially larger than those across the cool skin. Under wind speeds > 5 m s^-1, diurnal variations of the surface temperature due to the variations of the thermal molecular sublayer become more important.     

Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

地下水非稳定流的LTFAM算法

渗流域内应用拉普拉斯变换(LT)建立相应的有限分析(FAM)方程,顾及渗流域内地下水流的初始条件和边界条件,可在LT空间构成一个封闭的以水头像函数为变量的线性方程组。将此方程组所得的解,通过Stehfest数值反演公式,可归化为时间域的解(水头)。由于时间t被隐含在数值方程内,从而克服了传统数值法按时段(△t)逐步迭代的缺陷,提高了计算效率,也为用嵌入法建立地下水流管理模型提供了一条捷径。     

胶东海阳所超镁铁岩中强亲铁元素的地球化学特征及其构造环境探讨

强亲铁元素与亲石元素具有不同的地球化学行为，因此能够从不同的角度为造山带中超镁铁岩的成因及演化提供重要信息。位于苏鲁造山带东北端的胶东海阳所超镁铁岩主要由橄榄岩和辉石岩组成，它们常以团块状赋存于花岗质片麻岩中。虽然前人对这些超镁铁岩已经开展大量岩石学研究，但关于其成因及构造属性仍存在较大争议。本文开展了海阳所超镁铁岩的全岩主微量元素、强亲铁元素及Re-Os同位素的分析工作，结果显示蛇纹石化橄榄岩具有较高的MgO和Fe₂O₃^T含量，较低的Al₂O₃、TiO₂和CaO含量，明显富集流体迁移元素（U、Pb），亏损高场强元素（Zr、Hf），强亲铁元素没有发生明显分异，但Ru显示正异常，表明海阳所蛇纹石化橄榄岩是经历了低-中等程度部分熔融及熔/流体交代作用影响的残余地幔橄榄岩。海阳所辉石岩的主量元素表现出明显的结晶分异特征，稀土元素较原始地幔富集，铂族元素（PGEs）含量较低且发生了明显的分异，表明辉石岩的地幔源区经历过高程度的部分熔融和硫化物的分离。海阳所蛇纹石化橄榄岩的Os同位素地球化学特征表现出大洋亲和性，与辉石岩不具有熔体-残留体的关系。由于该地区发育较深层次的韧性剪切带，蛇纹石化橄榄岩中的橄榄石与辉石表现出韧性变形的特征，同时有辉石岩侵入到橄榄岩的现象，表明该地区的蛇纹石化地幔橄榄岩与辉石岩既不同时，也不同源，因此，暗示了该套岩石组合可能形成于大洋核杂岩（OCC）与洋脊型蛇绿岩（MOR）堆晶岩交互发育环境。     

地铁施工对济南白泉泉群流量的影响

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。     

漳州热田的对流热流和传导热流的研究

漳州地热系统属对流型地热系统.漳州热田是我国东南沿海地区目前所见热田中温度最高（121.5℃）的一个.地表热流值的研究表明,热田中心具有最大的实测热流值（359mw/m²）.本文根据热田内152个钻孔的测温资料和56块岩石样品的热导率数据,采用三种不同方法计算出漳州热田及其邻近地区的大地热流值,并讨论了热流值分布的特点.     

云南大地热流及地热地质问题

本文报道了云南20多个大地热流测试数据。据热流分布特点并参考前人有关资料,结合地球物理测深资料、岩石热物性与放射性生热率测试数据,对各区地壳深度范围的温度、热流配分组成作了粗略分析,划分出各有特征的若干地质-热流区。云南普遍地热偏高(滇东南小范围除外)的主要原因是新构造运动强烈;怒江以西近代岩浆活动亦强烈,具有典型的板块交汇带的基本地热特点