P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

The response of an elastic half‐space under a momentary shear line impulse

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.     

Low-cost volcano surveillance from space: case studies from Etna, Krafla, Cerro Negro, Fogo, Lascar and Erebus

 Satellite data offer a means of supplementing ground-based monitoring during volcanic eruptions, especially at times or locations where ground-based monitoring is difficult. Being directly and freely available several times a day, data from the advanced very high resolution radiometer (AVHRR) offers great potential for near real-time monitoring of all volcanoes across large (3000×3000 km) areas. Herein we describe techniques to detect and locate activity; estimate lava area, thermal flux, effusion rates and cumulative volume; and distinguish types of activity. Application is demonstrated using data for active lavas at Krafla, Etna, Fogo, Cerro Negro and Erebus; a pyroclastic flow at Lascar; and open vent systems at Etna and Stromboli. Automated near real-time analysis of AVHRR data could be achieved at existing, or cheap to install, receiving stations, offering a supplement to conventional monitoring methods. Received: 21 January 1997 / Accepted: 3 April 1997     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Infiltration at the Tai rain forest (Ivory Coast): Measurements and modelling

Infiltration experiments have been performed at three sites along a well-known catena under virgin tropical rain forest using a portable sprinkling infiltrometer. Experimentally determined infiltration curves are presented. Infiltration curves are also simulated on the basis of the Mein-Larson equation. The parameters for this model have been obtained from the infiltration curves (saturated conductivity) and simple soil moisture determinations (fillable porosity). The agreement between experimentally determined and modelled infiltration is reasonable, provided (a) saturated conductivity as derived from the experimental data is corrected, (b) a storage parameter, also derived from the experimental data, is added to the Mein-Larson model, and (c) the decline in soil porosity with depth is either small or occurs abruptly at shallow depth. Comparison of observed infiltration rates with rainfall intensity shows that Horton Overland Flow has to occur naturally at least on the middle and lower section of the catena. Despite the fact that most parameters can be estimated in principle from basic soil data, it remains advisable to obtain sprinkling infiltrometer field measurements, because of soil variability due to dynamic surface conditions, macroporosity, air entrapment, and irregularity of the wetting front.     

North European last glacial–interglacial transition (LGIT; 15–9 ka) tephrochronology: extended limits and new events

High‐precision correlation of palaeoclimatic and palaeoenvironmental records is crucial for testing hypotheses of synchronous change. Although radiocarbon is the traditional method for dating late Quaternary sedimentary sequences, particularly during the last glacial–interglacial transition (LGIT; 15–9 ka), there are inherent problems with the method, particularly during periods of climate change which are often accompanied by major perturbations in atmospheric radiocarbon content. An alternative method is the use of tephras that act as time‐parallel marker horizons. Within Europe, numerous volcanic centres are known to have erupted during the LGIT, providing considerable potential for high‐precision correlation independent of past radiocarbon fluctuations. Here we report the first identification of the Vedde Ash and Askja Tephra in Ireland, significantly extending the known provenance of these events. We have also identified two new horizons (the Roddans Port Tephras A and B) and tentatively recognise an additional horizon from Vallensgård Mose (Denmark) that provide crucial additional chronological control for the LGIT. Two phases of the Laacher See Tephra (LST) are reported, the lower Laacher See Tephra (LLST) and probably the C2 phase of the Middle Laacher See Tephra (MLST‐C2) indicating a more northeasterly distribution of this fan than reported previously. Copyright © 2006 John Wiley & Sons, Ltd.     

淇、卫河流域“96·8”致洪暴雨过程的数值模拟

应用LASGREMM模式 ,对淇、卫河流域“96·8”致洪暴雨过程进行了模拟预报分析。结果表明 ,LAS GREM模式不仅能较准确地预报淇、卫河流域暴雨落区 ,而且可以进一步分析天气过程的特点及演变情况     

Soil respiration in tropical seasonal rain forest in Xishuangbanna, SW China

With the static opaque chamber and gas chromatography technique, from January 2003 to January 2004 soil respiration was investigated in a tropical seasonal rain forest in Xishuangbanna, SW China. In this study three treatments were applied, each with three replicates: A (bare soil), B (soil+litter), and C (soil+litter+seedling). The results showed that soil respiration varied seasonally, low from December 2003 to February 2004, and high from June to July 2004. The annual average values of CO2 efflux from soil respiration differed among the treatments at 1% level, with the rank of C (14642 mgCO2· m-2. h-1)＞B (12807 mgCO2· m-2. h-1)＞A (9532 mgCO2· m-2. h-1). Diurnal variation in soil respiration was not apparent due to little diurnal temperate change in Xishuangbanna. There was a parabola relationship between soil respiration and soil moisture at 1% level. Soil respiration rates were higher when soil moisture ranged from 35% to 45%. There was an exponential relationship between soil respiration and soil temperature (at a depth of 5cm in mineral soil) at 1% level. The calculated Q1o values in this study,ranging from 2.03 to 2.36, were very near to those of tropical soil reported. The CO2 efflux in 2003was 5.34 kgCO2· m-2. a-1 from soil plus litter plus seedling, of them 3.48 kgCO2· m-2. a-1 from soil (accounting for 62.5%), 1.19 kgCO2· m-2. a-1 from litter (22.3%) and 0.67 kgCO2·m-2. a-1 from seedling (12.5%).     

广西钦江河口枯季水动力特征及其对冷空气的响应

河口作为陆海相互作用的关键带,是河流入海泥沙及污染物的主要归宿地。钦江作为广西第二大河流,在茅尾海资源开发利用和生态环境保护与修复中发挥着重要作用。为进一步认识钦江河口地区的沉积动力过程,2021年10月27日至12月8日在此展开水动力和水体环境要素的连续观测。现场观测与分析结果表明,观测期间钦江河口潮汐类型为正规全日潮,浅水分潮显著,平均潮差为3.07 m,具有落潮优势;潮流以全日分潮流为主导,平均流速为0.12 m/s,运动形式主要为往复流;余流流向主要为西南方向。钦江河口潮汐和潮流性质较外湾和茅尾海中部海域全日潮特征更显著。枯季期间河口表层沉积物随潮汐运动表现为侵蚀-沉降交替的变化规律,垂向上侵蚀通量大于沉降通量,水平方向上悬沙净向海输运。冷空气过境带来的降雨使得钦江河口水体盐度降低、浊度增大,悬沙浓度及输沙量增加,同时冷空气南下时北风增强引起钦江河口减水效应并使潮差略有增加,水动力强度增大,增加表层沉积物的活动性,从而引起底部沉积物再悬浮强度和频率增加。     

预测有效波高的深度学习模型研究

研究基于RNN、LSTM、GRU深度学习模型,针对NOAA浮标数据集中的44013、44014、44017浮标的数据,通过斯皮尔曼相关性分析提高模型预测效果。实验结果表明,在进行相关性分析后,S-RNN、S-LSTM、 S-GRU的预测效果均比原始RNN、LSTM、GRU模型预测效果好。此外,提出一种基于LSTM的LSTM-Attention 波高预测模型,并进行相关实验,量化LSTM-Attention模型的预测效果,实验结果表明LSTM-Attention模型有更好的预测效果。为评估模型的泛化能力,研究还提出了一种采用邻近浮标数据进行学习,预测浮标缺失数据的方 法。实验结果表明,该方法的预测精度可以达到97.93%。本研究为海浪预测提供了新的方法和思路,也为未来深 度学习模型在海浪预测中的应用提供了参考。