体波谱振幅相关系数计算程序的开发与实现

序列地震的体波谱振幅相关系数是地震序列震后快速判定的重要依据.体波谱振幅相关系数计算程序是在地震序列发生后,读取其地震波形数据,快速计算其相关系数,绘制滑动相关系数的时序图或聚类图,判定震源机制一致性程度,为地震序列的震后判定提供依据.     

Monitoring and numerical modelling of riverbank erosion processes: a case study along the Cecina River (central Italy)

Riverbank retreat along a bend of the Cecina River, Tuscany (central Italy) was monitored across a near annual cycle (autumn 2003 to summer 2004) with the aim of better understanding the factors influencing bank changes and processes at a seasonal scale. Seven flow events occurred during the period of investigation, with the largest having an estimated return period of about 1·5 years. Bank simulations were performed by linking hydrodynamic, fluvial erosion, groundwater flow and bank stability models, for the seven flow events, which are representative of the typical range of hydrographs that normally occur during an annual cycle. The simulations allowed identification of (i) the time of onset and cessation of mass failure and fluvial erosion episodes, (ii) the contributions to total bank retreat made by specific fluvial erosion and mass‐wasting processes, and (iii) the causes of retreat. The results show that the occurrence of bank erosion processes (fluvial erosion, slide failure, cantilever failure) and their relative dominance differ significantly for each event, depending on seasonal hydrological conditions and initial bank geometry. Due to the specific planimetric configuration of the study bend, which steers the core of high velocity fluid away from the bank at higher flow discharges, fluvial erosion tends to occur during particular phases of the hydrograph. As a result fluvial erosion is ineffective at higher peak discharges, and depends more on the duration of more moderate discharges. Slide failures appear to be closely related to the magnitude of peak river stages, typically occurring in close proximity to the peak phase (preferentially during the falling limb, but in some cases even before the peak), while cantilever failures more typically occur in the late phase of the flow hydrograph, when they may be induced by the cumulative effects of any fluvial erosion. Copyright © 2008 John Wiley & Sons, Ltd.     

A finite volume upwind scheme for the solution of the linear advection–diffusion equation with sharp gradients in multiple dimensions

A finite volume upwind numerical scheme for the solution of the linear advection equation in multiple dimensions on Cartesian grids is presented. The small-stencil, Modified Discontinuous Profile Method (MDPM) uses a sub-cell piecewise constant reconstruction and additional information at the cell interfaces, rather than a spatial extension of the stencil as in usual methods. This paper presents the MDPM profile reconstruction method in one dimension and its generalization and algorithm to two- and three-dimensional problems. The method is extended to the advection–diffusion equation in multiple dimensions. The MDPM is tested against the MUSCL scheme on two- and three-dimensional test cases. It is shown to give high-quality results for sharp gradients problems, although some scattering appears. For smooth gradients, extreme values are best preserved with the MDPM than with the MUSCL scheme, while the MDPM does not maintain the smoothness of the original shape as well as the MUSCL scheme. However the MDPM is proved to be more efficient on coarse grids in terms of error and CPU time, while on fine grids the MUSCL scheme provides a better accuracy at a lower CPU.     

Physically based modelling of sheet erosion (detachment and deposition processes) in complex hillslopes

A one‐dimensional uncoupled model governed by this research is a physics‐based modelling of the rainfall‐runoff induced erosion process. The presented model is composed of three parts of a three‐dimensional (3D) hillslope geometry, a nonlinear storage (kinematic wave) model for hillslope hydrological response, and an unsteady physically based surface erosion model. The 3D hillslope geometry model allows describing of the hillslope morphology by defining their plan shape and profile curvature. By changing these two topographic parameters, nine basic hillslope types are derived. The modelling of hillslope hydrological response is based on a flow continuity equation as the relation of discharge and flow depth is passed on kinematic wave approximation. The erosion model is based on a mass conservation equation for unsteady flow. The model assumes that suspended sediment does not affect flow dynamics. The model also accounts for the effect of flow depth plus loose soil depth on soil detachment. The presented model was run for two different precipitations, slope content, and length, and results were plotted for sediment detachment/deposition rate. Based on the obtained results, in hillslopes with convex and straight profile curvatures, sediment detachment only occurred in the whole length of the hillslope. However, in concave ones, sediment detachment and deposition only occurred together in hillslope. The hillslopes with straight profiles and convergent plans have the highest rate of detachment. Also, results show that most detachment rates occur in convex profile curvatures, which are about 15 times more than in straight profiles. Copyright © 2015 John Wiley & Sons, Ltd.     

Magnitude and processes of bank erosion at a small stream in Denmark

River banks are important sources of sediment and phosphorus to fluvial systems, and the erosion processes operating on the banks are complex and change over time. This study explores the magnitude of bank erosion on a cohesive streambank within a small channelized stream and studies the various types of erosion processes taking place. Repeat field surveys of erosion pin plots were carried out during a 4‐year period and observations were supplemented by continuous monitoring of volumetric soil water content, soil temperature, ground water level and exposure of a PEEP sensor. Bank erosion rates (17·6–30·1 mm year^?1) and total P content on the banks were relatively high, which makes the bank an important source of sediment and phosphorus to the stream, and it was estimated that 0·27 kg P_tot year^?1 ha^?1 may potentially be supplied to the stream from the banks. Yearly pin erosion rates exceeding 5 cm year^?1 were mainly found at the lower parts of the bank and were associated with fluvial erosion. Negative erosion pin readings were widespread with a net advance of the bank during the monitoring period mainly attributed to subaerial processes and bank failure. It was found that dry periods characterized by low soil water content and freeze–thaw cycles during winter triggered bank failures. The great spatial variability, in combination with the temporal interaction of processes operating at different scales, requires new tools such as 3‐D topographical surveying to better capture bank erosion rates. An understanding of the processes governing bank erosion is required for riparian management using vegetational measures as root size and structure play different roles when it comes to controlling bank erosion processes. Copyright © 2010 John Wiley & Sons, Ltd.     

Critical caving erosion width for cantilever failures of river bank

The cantilever failure is one of the typical bank failures, in which the lateral caving erosion at the bottom of the bank plays an important role. When the caving erosion width is larger than a certain value, the cantilever failures such as shear, toppling and stress failures may occur. In order to understand the condition of the cantilever failure, the collapse mechanisms of the cantilever failures are studied based on the bank stability theory and flume experiment. According to the bank stability equation with the lateral erosion, the critical caving erosion width (CCEW) formulas for the shear and toppling failures of simple slope bank were derived in this paper. The formulas show that the CCEW increases as the overhanging soil thickness and soil cohesion increase, and decreases as the crack depth on the bank surface and the slope angle of the bank increase. And these formulas were tested with experimental data, which shows the predicted values are good agreement with experimental data. The paper reveals a quantitative expression on the process of the river cantilever failure.     

Subaerial river bank erosion processes and their interaction with other bank erosion mechanisms on the River Arrow,Warwickshire, UK

River bank erosion occurs primarily through a combination of three mechanisms: mass failure, fluvial entrainment, and subaerial weakening and weathering. Subaerial processes are often viewed as ‘preparatory’ processes, weakening the bank face prior to fluvial erosion. Within a river basin downstream process ‘domains’ occur, with subaerial processes dominating the upper reaches, fluvial erosion the middle, and mass failure the lower reaches of a river. The aim of this paper is to demonstrate that (a) subaerial processes may be underestimated as an erosive agent, and (b) process dominance has a temporal, as well as spatial, aspect. Bank erosion on the River Arrow, Warwickshire, UK, was monitored for 16 months (December 1996 to March 1998) using erosion pins. Variations in the rate and aerial extent of erosion are considered with reference to meteorological data. Throughout the first 15 months all erosion recorded was subaerial, resulting in up to 181 mm a^?1 of bank retreat, compared with 13 to 27 mm a^?1 reported by previous researchers. While the role of subaerial processes as ‘preparatory’ is not contended, it is suggested that such processes can also be erosive. The three bank erosion mechanisms operate at different levels of magnitude and frequency, and the River Arrow data demonstrate this. Thus the concept of process dominance has a temporal, as well as spatial aspect, particularly over the short time‐periods often used for studying processes in the field. Perception of the relative efficacy of each erosive mechanism will therefore be influenced by the temporal scale at which the bank is considered. With the advent of global climate change, both these magnitude–frequency characteristics and the consequent interaction of bank erosion mechanisms may alter. It is therefore likely that recognition of this temporal aspect of process dominance will become increasingly important to studies of bank erosion processes. Copyright © 2001 John Wiley & Sons, Ltd.     

Quantifying hillslope erosion rates and processes for a coastal California landscape over varying timescales

The Earth's surface erodes by processes that occur over different spatial and temporal scales. Both continuous, low‐magnitude processes as well as infrequent, high‐magnitude events drive erosion of hilly soil‐mantled landscapes. To determine the potential variability of erosion rates we applied three independent, field‐based methods to a well‐studied catchment in the Marin Headlands of northern California. We present short‐term, basin‐wide erosion rates determined by measuring pond sediment volume (40 years) and measured activities of the fallout nuclides ¹³⁷Cs and ²¹⁰Pb (40–50 years) for comparison with long‐term (>10 ka) rates previously determined from in situ‐produced cosmogenic ¹⁰Be and ²⁶Al analyses. In addition to determining basin‐averaged rates, ¹³⁷Cs and ²¹⁰Pb enable us to calculate point‐specific erosion rates and use these rates to infer dominant erosion processes across the landscape. When examined in the context of established geomorphic transport laws, the correlations between point rates of soil loss from ¹³⁷Cs and ²¹⁰Pb inventories and landscape morphometry (i.e. topographic curvature and upslope drainage area) demonstrate that slope‐driven processes dominate on convex areas while overland flow processes dominate in concave hollows and channels. We show a good agreement in erosion rates determined by three independent methods: equivalent denudation rates of 143 ± 41 m Ma^?1 from pond sediment volume, 136 ± 36 m Ma^?1 from the combination of ¹³⁷Cs and ²¹⁰Pb, and 102 ± 25 m Ma^?1 from ¹⁰Be and ²⁶Al. Such agreement suggests that erosion of this landscape is not dominated by extreme events; rather, the rates and processes observed today are indicative of those operating for at least the past 10 000 years. Copyright © 2006 John Wiley & Sons, Ltd.     

A simple method for the determination of ionic diffusion coefficients in flooded soils

Soil cores from river marginal wetlands from the Torridge and Severn catchments in the UK were collected to study rates of soil denitrification at different sites and at two stations (levee and backplain depression) at the river margin. Half the cores were sterilized prior to flooding to destroy the denitrifying bacteria. After flooding and equilibration, monitoring the concentration of amended nitrate in the supernatant of the sterile cores over a period of 7 days provided a simple procedure for the estimation of the diffusion coefficient of the nitrate ion in the flooded soils. An expression was developed that permitted this diffusion coefficient to be extracted from the slope of a plot of supernatant concentration versus (time)^1/2. The values obtained, at 15 °C, varied from 2·4 to 6·8 × 10^?10m²s^?1. Sterile cores are usually treated as controls in denitrification experiments; this work develops a procedure whereby they may yield useful soil process information. Copyright © 2001 John Wiley & Sons, Ltd.     

Investigation on the stochastic simulation of strong ground motions for Bucharest area

This short article evaluates the stochastic method of ground motion simulation for Bucharest area using both the single-corner frequency model and recently introduced double-corner frequency models. A dedicated Q model is derived using ground motions obtained during the largest Vrancea earthquakes from the past 30 years. The simulated ground motions are tested against the observed data from the Vrancea earthquakes of August 1986 and May 1990. Moreover, the observed data are also compared against simulations obtained using the Q model derived by Oth et al. (2008). Finally, the results of the simulations show that the derived Q model is better suited for larger magnitude events, while the Q model of Oth et al. (2008) provides better results for smaller earthquakes.     

Wind tunnel tests of the dynamic processes that control wind erosion of a sand bed

Aeolian sand transport is a complicated process that is affected by many factors (e.g. wind velocity, sand particle size, surface microtopography). Under different experimental conditions, erosion processes will therefore produce different results. In this study, we conducted a series of wind tunnel experiments across a range of wind velocities capable of entraining sand particles (8.0, 10.0, 12.0, and 14.0 m s^-1) to study the dynamic changes of the shear velocity, aerodynamic roughness length, and sand transport. We found that the shear velocity and aerodynamic roughness length are not constant; rather, they change dynamically over time, and the rules that describe their changes depend on the free-stream air velocity. For wind tunnel experiments without feeding sand into the airflow, the sand bed elevation decreases with increasing erosion time, and this change significantly affected the values of shear velocity and aerodynamic roughness length. A Gaussian distribution function described the relationships between the sand transport rate (q_T) and the duration of wind erosion (T). It is therefore necessary for modelers to consider both deflation of the bed and the time scale used when calculating sand transport or erosion rates. © 2018 John Wiley & Sons, Ltd.     

Numerical modelling of aeolian erosion over a surface with non‐uniformly distributed roughness elements

The present study is focused on the analysis of the mean wall friction velocity on a surface including roughness elements exposed to a turbulent boundary layer. These roughness elements represent non‐erodible particles over an erodible surface of an agglomeration of granular material on industrial sites. A first study has proposed a formulation that describes the evolution of the friction velocity as a function of geometrical parameters and cover rate with different uniform roughness distributions. The present simulations deal with non‐uniform distributions of particles with a random sampling of diameters, heights, positions and arrangements. The evolution (relative to geometrical parameters of the roughness elements) of the friction velocity for several non‐uniform distributions of roughness elements was analysed by the equation proposed in the literature and compared to the results obtained with the numerical simulations. This comparison showed very good agreement. Thus, the formulation developed for uniform particles was found also to be valid for a larger spectrum of particles noted on industrial sites. The present work aims also to investigate in detail the fluid mechanics over several roughness particles. Copyright © 2013 John Wiley & Sons, Ltd.     

利用不同倾斜仪和应变仪检测地球自由振荡的对比与分析

中国地震局地壳形变观测台网布设有垂直摆倾斜仪、钻孔倾斜仪、洞体应变仪、分量式钻孔应变仪与体应变仪等地形变观测仪器。这些观测仪器均记录到了2011年日本M_W9.0大地震激发的自由振荡信号。本文分别利用单台数据和多台数据叠积,检测到_0S_3~_0S_(30)全部的球型自由振荡基频振型和_0T_3~_0T_(20)全部的环型自由振荡基频振型及部分谐频振型。此外,通过对这些检测结果的对比,分析了它们对不同自由振荡类型、不同频段的振型检测能力。分析发现垂直摆倾斜仪对球型自由振荡的检测结果最佳,且由于在低频段有较高的噪声干扰,钻孔倾斜仪无法检测到低阶的球型自由振荡。对于环型自由振荡的检测,分量式钻孔应变仪检测结果最佳。     

小北湖火山台迁建前后观测数据对比分析

小北湖火山台原址位于原始森林保护区,按照林区防火要求,2019年选址迁建.选取迁建前后小北湖火山台观测数据,通过计算不同频段范围内背景噪声记录的加速度功率谱密度,分析了迁建前后地震观测台站的噪声特征.结果 表明,在1-20 Hz频段内,迁建后功率谱噪声值普遍降低10-20 dB;从RMS值来看,UD向降低22.0％,E...     

Revisiting the Phasor-Walkout method for detailed investigation of harmonic signals in time series

A well-known but rarely used powerful method to investigate the presence of harmonic signals in time series is the Phasor Walkout method (other names are: Graphical Fourier Transform, Summation Dial, Complex Demodulation). At a given test frequency the complex contributions (phasors) to the Fourier Transform of each sample in an equidistantly sampled series are added vectorially in the complex plane. The resulting pattern, the walkout, reveals information about the properties of the signal which is not easy to obtain by other methods. Synthetic examples are used to demonstrate the resolving power of the method. The following geophysical examples for the application of this method are shown: determination of the frequency of the breathing mode₀ S ₀ of the earth after a large earthquake; the study of superconducting gravimeter records after a large deep earthquake used in the core mode interpretation of a spectral peak, the study of the residualS ₃ (8 h period) signal in a tidal record and the bichromatic Rayleigh-waves from Mount Pinatubo on June 15, 1991.     

地震动反应谱变异系数分析

利用随机振动理论,分别考虑3种地震动随机过程模型,本文对反应谱变异系数及相关问题进行了深入分析。研究结果表明：在一定范围内,地震动模型、场地条件和结构体系阻尼比对结构加速度反应越零率有明显影响,但这三者对反应谱变异系数的量值影响不大。在工程设计中,可以采用白噪声模型按本文给出公式计算反应谱变异系数。     

Runoff dependent erosivity and slope length factors suitable for modelling annual erosion using the Universal Soil Loss Equation

Despite revisions and refinements, the Revised Universal Soil Loss Equation (RUSLE), which is the revised version of the Universal Soil Loss Equation (USLE), over predicts small annual soil losses and under predicts large annual soil losses. To some large extent, this results from the equation over estimating small event soil losses and under estimating large event soil losses. Replacing the USLE/RUSLE event erosivity index (EI₃₀) by the product of EI₃₀ and the runoff ratio (Q_R) significantly reduces the errors in estimating event erosion when runoff is measured, but the USLE‐M, the USLE variant that uses the Q_REI₃₀ index, requires crop and support practice factors that differ from those used in the RUSLE. The theory which enables the Q_REI₃₀ index to be used in association with the RUSLE crop and support practice factors is presented. In addition, the USLE/RUSLE approach was developed for conditions where runoff is produced uniformly over a hill slope. A runoff dependent slope length factor that takes account of runoff variations over a hill slope is presented and demonstrated for the situation where runoff from a low runoff producing area passes onto an area where runoff is produced more readily. Copyright © 2007 John Wiley & Sons, Ltd.