Stability Analysis and Numerical Simulation of Differential Frost Heave

Differential frost heave is often implicated in the formation of patterned ground in regions subject to recurrent freezing and thawing. A linear stability analysis (LSA) indicates that a continuum model of frost heave is linearly unstable under typical natural freezing conditions of silty-clay soils. A two-dimensional non-linear numerical analysis corroborates the frozen time LSA results, and also indicates the importance of non-linear and time-dependent terms that ultimately lead to a preferred mode, which the LSA fails to predict. Instability of the one-dimensional solution occurs at shallow freezing depths and near-zero surface loads when positive perturbations in the ice content at the freezing front lead to a concomitant increase in thermomolecular pressure and upward ice velocity. Differential frost heave can then occur because of the increased heat flux from the perturbed surfaces. A three-dimensional model using random initial surface perturbations indicates that regular surface patterns will evolve with a length scale in the order of 2–4 meters, which corresponds quite closely with naturally-occurring non-sorted patterned ground.     

Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.     

研究“源兆”的若干问题及途径分析

从位场理论的观点讨论了场和源的概念,建议将“场兆”改称为构造前兆,将“源兆”称为震源前兆。文中还分析了震源体与各种前兆观测方法探测的深度,并根据地磁与重力观测方法不同的物理基础和地质基础,指出震磁效应对于研究震源前兆的特殊意义,文中还介绍了１９７６唐山地震观测到震源前兆的部分实际地磁资料。     

形变,应变短临前兆标志体系

建立各手段的前兆异常标志,从而组合各类前兆标志体系是“八五”地震短临预报攻关研究项目的重要目标之一。本文以形变,应变手段的主要研究对象,初步建立了一个形变应变短临前兆标志体系,并总结出了各种异常的判别标志,为利用系统优化理论和专家系统理论进一步研究形变应前兆系统和地震之间的关系,从而为提高中短期地震预报的水平打下了一定的基础。     

暖切变型江淮梅雨锋结构及其形成和维持机制

对在1991年7月上旬特别稳定的梅雨形势下持续10天的特大暴雨期作了合成和典型过程研究,并进一步探讨了暖切变型江淮梅雨锋天气系统模型。它在自由大气层表现为结构竖立的相当正压切变线和一种不同于经典地转动量近似的带状准二维运动场。边界层的浅层锋生倾斜环流以及偏于切变线轴以南的Ekman抽吸与自由大气层中基本上由积云对流加热驱动的竖直环流相互依存和耦合,造成持续的大暴雨。这种系统的初生和中断后复苏的机制可能是由斜压性极弱的副热带高空东移的短波扰动在低层“湿度锋”带紧南侧触发导致的一种带状CISK。     

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.     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

The emplacement of intermediate volume ignimbrites: A case study from Roccamonfina Volcano,Southern Italy

A model is presented for the emplacement of intermediate volume ignimbrites based on a study of two 6 km³ volume ignimbrites on Roccamonfina Volcano, Italy. The model considers that the flows were slow moving, and quickly deflated from turbulent to non-turbulent conditions. Yield strength and density increased whereas fluidisation decreased with time and runout of the pyroclastic flows. In proximal locations, on the caldera rim, heterogeneous exposures including discontinuous lithic breccias, stratified and cross-stratified units interbedded with massive ignimbrite suggest deposition from turbulent flows. In medial locations thick, massive ignimbrite occurs associated with three types of co-ignimbrite lithic breccia which we interpret as being emplaced by non-turbulent flows. Multiple grading of different breccia/lithic concentration types within single flow units indicates that internal shear occurred producing overriding or overlapping of the rear of the flow onto the slower-moving front part. This overriding of different parts of non-turbulent pyroclastic flows could be caused by at least two different mechanisms: (1) changes in flow regime, such as hydraulic jumps that may occur at breaks in slope; and (2) periods of increased discharge rate, possibly associated with caldera collapse, producing fresh pulses of lithic-rich material that sheared onto the slower-moving part of the flow in front.We propose that ground surge deposits enriched in pumice compared with their associated ignimbrite probably formed by a flow separation mechanism from the top and front of the pyroclastic flow. These turbulent clouds moved ahead of the non-turbulent lower part of the flow to form stratified pumice-rich deposits. In distal regions well-developed coarse, often clast-supported, pumice concentrations zones and coarse intra-flow-unit lithic concentrations occur within the massive ignimbrite. We suggest that the flows were non-turbulent, possessed a relatively high yield strength and may have moved by plug flow prior to emplacement.     

施工顺序对软土地坪复合地基的影响探讨

场地回填土与搅拌桩施工顺序将直接影响软土地坪复合地基的处理效果,其中沉降控制是影响地坪安全运行的关键。以近海软土地区工业厂房水泥土搅拌桩复合地基工程为例,探讨了复合地基水泥土搅拌桩、回填土两种不同施工顺序的影响因素及相应处理效果;通过数值模拟试验,计算分析了不同填土厚度下两种工况的地基固结总沉降、施工沉降、工后沉降。结果表明,采用“先土后桩”的施工顺序的施工期沉降大于“先桩后土”,但工后沉降大大减小,提高了软土地坪地基处理效果,有利于处理后地坪的安全运行,验证了前述分析结论,为合理设计方案的确定提供了重要的参考依据。     

Possible creep-related tilt precursors obtained in the Central Apennines (Italy) and in the Southern Caucasus (Georgia)

Daily averaged tilt component data from two sites of the Central Apennines (Italy) and of the Southern Caucasus (Georgia), respectively, revealed intermediate-term tilts as possible precursors to earthquakes (M=3.0÷4.7) which occurred in the above-mentioned seismic areas within a distance of 50 km from the sites. A good temporal correlation as well as a fair spatial correspondence between these residual tilts (with amplitude and duration of some microradians and months, respectively) and main shocks were pointed out, by removing both secular trends and seasonal thermoelastic effects from the raw tilts. An attempt was made to justify the above-mentioned results, based on the assumption that the observed intermediate-term preseismic tilts are the manifestation of aseismic creep episodes of comparable duration in the fault materials of thrust faults close to the tilt sites. The mechanism refers to a strain field slowly propagating from the preparation focal area to the tilt site, through crustal blocks separated by weak transition zones. This propagation is thought to be the cause of the local aseismic fault slip recorded by the tiltmeters. Previously, both discrete structures and strain propagation effects were revealed in the Central Apennines and are thought also to exist in the Southern Caucasus. As in the past, the rheological properties of fault materials are revealed as viscoelastic ones. In fact, creep equations obtained by applying several viscoelastic models on our data, proved to fit quite well some of the observed tilt precursors, producing viscosity and rigidity values very similar to those reported in literature.Professor Petr Viktorovich Manjgaladze died during the writing of this paper