Seismic Displacement of Gravity Walls by a Two-body Model

Gravity walls retaining dry soil are modeled as a system of two bodies: (a) the gravity wall that slides along the wall-foundation soil boundary and (b) the critical soil wedge in the soil behind the wall. The strength of the system is defined by both the frictional and the cohesional components of resistance. The angle of the prism of the critical soil wedge behind the wall is obtained using the limit equilibrium method. The model accounts for changes in the geometry of the backfill soil behind the wall by considering the displacements at the end of each time step under limit equilibrium. The model shows that the standard (single) block model is over-conservative for the extreme case of critical-to-applied-seismic acceleration ratios less than about 0.30, but works well for cases where this ratio ranges between 0.5 and 0.8. The model is applied to predict the seismic displacement of gravity walls (a) tested in the shaking-table and (b) studied numerically by elaborate elasto–plastic analyses.     

Fracture in Westerly Granite under AE Feedback and Constant Strain Rate Loading: Nucleation, Quasi-static Propagation, and the Transition to Unstable Fracture Propagation

New observations of fracture nucleation are presented from three triaxial compression experiments on intact samples of Westerly granite, using Acoustic Emission (AE) monitoring. By conducting the tests under different loading conditions, the fracture process is demonstrated for quasi-static fracture (under AE Feedback load), a slowly developing unstable fracture (loaded at a `slow' constant strain rate of 2.5 × 10<sup>−6</sup> /s) and an unstable fracture that develops near instantaneously (loaded at a `fast' constant strain rate of 5 × 10⁻⁵ /s). By recording a continuous ultrasonic waveform during the critical period of fracture, the entire AE catalogue can be captured and the exact time of fracture defined. Under constant strain loading, three stages are observed: (1) An initial nucleation or stable growth phase at a rate of ~ 1.3 mm/s, (2) a sudden increase to a constant or slowly accelerating propagation speed of ~ 18 mm/s, and (3) unstable, accelerating propagation. In the ~ 100 ms before rupture, the high level of AE activity (as seen on the continuous record) prevented the location of discrete AE events. A lower bound estimate of the average propagation velocity (using the time-to-rupture and the existing fracture length) suggests values of a few m/s. However from a low gain acoustic record, we infer that in the final few ms, the fracture propagation speed increased to 175 m/s. These results demonstrate similarities between fracture nucleation in intact rock and the nucleation of dynamic instabilities in stick slip experiments. It is suggested that the ability to constrain the size of an evolving fracture provides a crucial tool in further understanding the controls on fracture nucleation.     

长岭山北麓断裂上新发现的地震地表破裂带

通过野外地震地质调查,在甘肃省古浪县、景泰县交界处的长岭山地区新发现了丰富的线状分布的地震地表破裂现象,表现为基岩崩塌、黄土滑坡、地震裂缝、地震土林、地震沟槽、陷落坑以及断层陡坎等,局部可见残留的破裂滑动自由面。其分布严格受长岭山北麓活动断裂控制。通过地震地表破裂带本身的特征分析、年代学分析、地表破裂带长度与震级的拟合关系以及结合区域活动断裂资料研究,认为该地震地表破裂带可能为1927年古浪8.0级地震的产物。     

Evaluation of the Urban Tile in MOSES using Surface Energy Balance Observations

The UK Met Office has introduced a new scheme for its urban tile in MOSES 2.2 (Met Office Surface Exchange Scheme version 2.2), which is currently implemented within the operational Met Office weather forecasting model. Here, the performance of the urban tile is evaluated in two urban areas: the historic core of downtown Mexico City and a light industrial site in Vancouver, Canada. The sites differ in terms of building structures and mean building heights. In both cases vegetation cover is less than 5%. The evaluation is based on surface energy balance flux measurements conducted at approximately the blending height, which is the location where the surface scheme passes flux data into the atmospheric model. At both sites, MOSES 2.2 correctly simulates the net radiation, but there are discrepancies in the partitioning of turbulent and storage heat fluxes between predicted and observed values. Of the turbulent fluxes, latent heat fluxes were underpredicted by about one order of magnitude. Multiple model runs revealed MOSES 2.2 to be sensitive to changes in the canopy heat storage and in the ratio between the aerodynamic roughness length and that for heat transfer (temperature). Model performance was optimum with heat capacity values smaller than those generally considered for these sites. The results suggest that the current scheme is probably too simple, and that improvements may be obtained by increasing the complexity of the model.     

Characteristics of secondary circulations over an inhomogeneous surface simulated with large-eddy simulation

Large-eddy simulation is used to study secondary circulations in the convective boundary layer modulated as a result of horizontally varying surface properties and surface heat fluxes over flat terrain. The presence of heat flux heterogeneity and its alignment with respect to geostrophic wind influences the formation, strength and orientation of organized thermals. Results show boundary-attached roll formation along heat flux maxima in the streamwise direction. The streamwise organization of the updrafts and downdrafts formed downwind of heterogeneities leads to counter-rotating secondary circulations in the crosswind plane. The distribution of resolved-scale pressure deviations shows large pressure gradients in the crosswind plane. Spanwise and vertical velocity variances and heat flux profiles depict considerable spatial variability compared to a homogeneous forest simulation. Secondary circulations are observed for various ambient wind scenarios parallel and perpendicular to heterogeneities. In the presence of increased wind speed, thermals emerging from the heat flux heterogeneity are elongated, and organize along and downwind of large-scale heterogeneity in the streamwise direction. Simulation with a reduced heat flux shows a shallower circulation with a lower aspect ratio. Point measurements of heat flux inside the roll circulation could be overestimated by up to 15–25% compared to a homogeneous case.     

The effect of mountainous topography on moisture exchange between the “surface” and the free atmosphere

Typical numerical weather and climate prediction models apply parameterizations to describe the subgrid-scale exchange of moisture, heat and momentum between the surface and the free atmosphere. To a large degree, the underlying assumptions are based on empirical knowledge obtained from measurements in the atmospheric boundary layer over flat and homogeneous topography. It is, however, still unclear what happens if the topography is complex and steep. Not only is the applicability of classical turbulence schemes questionable in principle over such terrain, but mountains additionally induce vertical fluxes on the meso-γ scale. Examples are thermally or mechanically driven valley winds, which are neither resolved nor parameterized by climate models but nevertheless contribute to vertical exchange. Attempts to quantify these processes and to evaluate their impact on climate simulations have so far been scarce. Here, results from a case study in the Riviera Valley in southern Switzerland are presented. In previous work, measurements from the MAP-Riviera field campaign have been used to evaluate and configure a high-resolution large-eddy simulation code (ARPS). This model is here applied with a horizontal grid spacing of 350 m to detect and quantify the relevant exchange processes between the valley atmosphere (i.e. the ground “surface” in a coarse model) and the free atmosphere aloft. As an example, vertical export of moisture is evaluated for three fair-weather summer days. The simulations show that moisture exchange with the free atmosphere is indeed no longer governed by turbulent motions alone. Other mechanisms become important, such as mass export due to topographic narrowing or the interaction of thermally driven cross-valley circulations. Under certain atmospheric conditions, these topographical-related mechanisms exceed the “classical” turbulent contributions a coarse model would see by several times. The study shows that conventional subgrid-scale parameterizations can indeed be far off from reality if applied over complex topography, and that large-eddy simulations could provide a helpful tool for their improvement.     

Kinetics of bubble collision and attachment to hydrophobic solids: I. Effect of surface roughness

Influence of surface roughness of the Teflon plates on kinetics of the bubble attachment was studied. Phenomena occurring during collisions of the air bubble, rising in clean water, with Teflon plates, differing only in their surface roughness, were recorded and analysed using a high-speed camera. Variations of the local velocity of the bubble during the collisions and the time of the bubble attachment were determined. It was found that the Teflon surface roughness was the parameter of a crucial importance for the attachment time of the colliding bubble. Depending on degree of the surface roughness the time of the attachment varied by over order of magnitude (from 3 to over 80 ms). In the case the Teflon surfaces having roughness below 1 μm there were recorded four to five “approach–bounce” cycles prior to the bubble attachment. Moreover, after the first collision the rapid pulsations of the bubble shape (within fraction of millisecond) were recorded. For surfaces of roughness ca. 50 μm and larger the attachment always occurred during the first collision—there was no bouncing observed and the time of the attachment was below 3 ms. It was documented that presence of a micro-bubble at the surface facilitated attachment of the colliding bubble.     

Statistical features of aftershock distribution size for moderate and large earthquakes in Chinese mainland

Based on data of earthquake sequences with MS≥5.0 in Chinese mainland from 1970 to 2004,for different se-quence types and different rupture modes of the main shock,the relationship between aftershock distribution size Rand the magnitude of the main shock M0 has been studied statistically.Considering the rupture mode of the mainshock,we give the quantitative statistical relationships between R and M0 under 95%confidence level for differentsequence types.Qualitatively,lgR,the logarithm of the aftershock distribution size,is positively correlative to theM0,but the data distribution is dispersed.Viewing from different sequence types,the correlation between R and M0is very weak for isolated earthquake type(IET)sequence,R distributes in the range from 5 to 60 km;For main-shock-aftershock type(MAT),lgR is positively correlative to M0;For multiple main shock type(MMT),the core-lation between lgR and M0 is not very obvious when M0≤6.2 and R distributes in the range from 5 to 70 km,whileit shows a linear correlation when M0≥6.3.The statistical results also show that the occupational ratios of differentsequence types for strike-slip and oblique slip are almost the same.But for dip-slip(mostly are thrust mechanisms),the ratio of MAT is higher than that of IET and MMT.Comparing with previous results,it indicates that,when M0is large enough,R is mainly determined by M0 and there is almost no relationship with the rupture mode of themain shock.     

2005年10月8日巴基斯坦7.8级地震热红外异常

收集了2005年巴基斯坦7.8级地震区Noaa17/AVHRR的2004-2005年的热红外遥感资料,并进行地表温度反演,得到了以震中为中心7°×7°范围的地表温度值。经热红外图像解译与构造关系对比分析,表明热红外解译图像与构造分布具有一致性。震源区地表温度时间序列显示:震前大约3个月内震源区存在明显的热红外异常。     

不同标度震级关系和台基影响问题探讨

利用近年中国数字地震台网观测报告中的中国及邻区地震目录资料,取能够同时计算面波震级MS和近震震级ML的地震,拟合出2种震级之间的关系,发现与过去使用30余年的换算公式存在系统差别。考虑现代化地震数字观测系统有动态范围大、频带宽、台站多、覆盖区域广、测定精度高的特点,由此拟合的关系应当更可靠。考虑不同误差因素统计面波震级MS和近震震级ML的关系表明,2种震级之间换算和不换算的差别并不大,因此建议MS和ML震级之间不换算为宜。分析地震台站的场地响应发现:一些基岩台的放大因子在1附近,场地响应不随频率变化;处于沉积土层的台站,放大因子在低频段大于1,高频段小于1;还有一些台站在某个频段范围大于1,或者小于1。另外,沉积土层台站的单台震级和台网平均震级差也随周期增大,呈现由负到正的增长趋势。看来并不存在固定的台站校正系数,因此不能把台站校正作为提高震级测定精度的方法。