统一应力场中基底断裂对盖层复杂断块变形的影响——来自砂箱实验的启示

基底先存断裂的活动会对其盖层岩层的变形起到较大的影响,尤其当应力直接作用于基底时更是对变形起到了控制作用。通过一系列的砂箱实验模拟分析了在统一的构造应力场中,当深部的基底断裂作平移滑动时盖层断块的被动变形情况:当其具有伸展分量时,会形成一个近对称的走滑 伸展裂谷形态,当具有挤压分量时,会形成以逆冲走滑断裂为边界的对称的局部挠曲隆起;同时在剖面上会形成典型的走滑构造特征,变形区域的大小与伸展或挤压分量的大小有关。先存的盖层断块受到基底作用力时,除内部变形本身还会发生旋转,形成局部的拉伸和挤压区;当一个地区的基底断裂多次活动甚至发生反转时,就会使地表形成特别复杂的构造现象;郯庐断裂带中段的埕岛-垦东潜山构造带的变形是一个典型的走滑基底控制的情况,实验结果证明基底的走滑反转变形造成这些断块的旋扭,盖层的非完全反转形成了剖面上的“复式花状构造”。     

Discontinuity of the Mozumi–Sukenobu fault low-velocity zone, central Japan, inferred from 3-D finite-difference simulation of fault zone waves excited by explosive sources

We delineate shallow structures of the Mozumi–Sukenobu fault, central Japan, using fault zone waves generated by near-surface explosions and detected by a seismometer array. Two explosive sources, S1 and S2, were placed at a distance of about 2 km from the array, and the other two, S3 and S4, were at a distance of about 4 km. Fault zone head waves and fault zone trapped waves following direct P wave arrivals were clearly identified in the seismograms recorded by a linear seismometer array deployed across the fault in a research tunnel at a depth of 300 m. Synthetic waveforms generated by a 3-D finite-difference (3-D FD) method were compared with observed fault zone waves up to 25 Hz. The best fitting model indicates a 200-m-wide low-velocity zone extending at least to shot site S1 located 2 km east of the seismic array with a 20% decrease in the P wave velocity relative to the wall rock. The width of the low-velocity zone is consistent with the fault zone defined by direct geological observation in the research tunnel. However, the low-velocity zone should disappear just to the east of the site S1 to explain the observed fault zone waves for shot S3 and S4 located 4 km east of the seismometer array. Yet the observation and the simulation show notable trapped wave excitation even though shots S3 and S4 are outside the fault zone. These results indicate that (1) the effective waveguide for seismic waves along the fault does not exist east of source site S1 although the surface traces of the fault are observed in this region, and (2) considerable trapped waves can be excited by sources well outside the fault zone. These results highlight the along-strike variability in fault zone structure.     

Thessaloniki–Gerakarou Fault Zone (TGFZ): the western extension of the 1978 Thessaloniki earthquake fault (Northern Greece) and seismic hazard assessment

Active faulting and seismic properties are re-investigated in the eastern precinct of the city of Thessaloniki (Northern Greece), which was seriously affected by two large earthquakes during the 20th century and severe damage was done by the 1759 event. It is suggested that the earthquake fault associated with the occurrence of the latest destructive 1978 Thessaloniki earthquake continues westwards to the 20-km-long Thessaloniki–Gerakarou Fault Zone (TGFZ), which extends from the Gerakarou village to the city of Thessaloniki. This fault zone exhibits a constant dip to the N and is characterised by a complicated geometry comprised of inherited 100°-trending faults that form multi-level branching (tree-like fault geometry) along with NNE- to NE-trending faults. The TGFZ is compatible with the contemporary regional N–S extensional stress field that tends to modify the pre-existing NW–SE tectonic fabric prevailing in the mountainous region of Thessaloniki. Both the 1978 earthquake fault and TGFZ belong to a ca. 65-km-long E–W-trending rupture fault system that runs through the southern part of the Mygdonia graben from the Strymonikos gulf to Thessaloniki. This fault system, here called Thessaloniki–Rentina Fault System (TRFS), consists of two 17–20-km-long left-stepping 100°-trending main fault strands that form underlapping steps bridged by 8–10-km-long ENE–WSW faults. The occurrence of large (M6.0) historical earthquakes (in 620, 677 and 700 A.D.) demonstrates repeated activation, and therefore the possible reactivation of the westernmost segment, the TGFZ, could be a major threat to the city of Thessaloniki. Changes in the Coulomb failure function (ΔCFF) due to the occurrence of the 1978 earthquake calculated out in this paper indicate that the TGFZ has been brought closer to failure, a convincing argument for future seismic hazard along the TGFZ.     

Ultrashallow seismic experiment on a trenched section of the Chelunpu fault zone, Taiwan

Ultrashallow P-wave seismic reflection experiments were conducted at a model test site and in a trenched shallow fault zone along the Chelunpu fault line. The field layout was designed to have the shallowest undistorted reflection from about 1 m depth with 0.5 m vertical resolution. The smallest group interval tested in this study was 0.05 m with a 0.25 ms sample interval, which can avoid spatial aliasing of ground roll if the target is very shallow and the velocities are low. Data processing was designed to be simple but consistent. As the ultrashallow reflections may be contaminated with high-amplitude coherent noise in many aspects, first break muting and surgical muting were performed on each file as detailed as possible, and f–k filtering was applied mainly for the purpose of attenuating the aliasing energy and back-scattered noise. Data acquired in this study show that the low P-wave velocities (< 200 m/s) and high dominant frequencies (120–200 Hz) of near-surface layers may have a potential vertical resolution of 0.4 m or even better.Comparing the test profile with the ground-penetrating radar (GPR) control profile of the same test site and correlating the results obtained from the study site with those of the geologic cross-section of the trench, this experiment demonstrates the possibility of using seismic methods in investigating shallow structures at depths of less than a few meters with vertical resolution comparable to the GPR technique.     

Temporal variation in the geometry of a strike–slip fault zone: Examples from the Dead Sea Transform

The location of the active fault strands along the Dead Sea Transform fault zone (DST) changed through time. In the western margins of Dead Sea basin, the early activity began a few kilometers west of the preset shores and moved toward the center of the basin in four stages. Similar centerward migration of faulting is apparent in the Hula Valley north of the Sea of Galilee as well as in the Negev and the Sinai Peninsula. In the Arava Valley, seismic surveys reveal a series of buried inactive basins whereas the current active strand is on their eastern margins. In the central Arava the centerward migration of activity was followed by outward migration with Pleistocene faulting along NNE-trending faults nearly 50 km west of the center. Largely the faulting along the DST, which began in the early–middle Miocene over a wide zone of up to 50 km, became localized by the end of the Miocene. The subsidence of fault-controlled basins, which were active in the early stage, stopped at the end of the Miocene. Later during the Plio-Pleistocene new faults were formed in the Negev west of the main transform. They indicate that another cycle has begun with the widening of the fault zone. It is suggested that the localization of faulting goes on as long as there is no change in the stress field. The stresses change because the geometry of the plates must change as they move, and consequently the localization stage ends. The fault zone is rearranged, becomes wide, and a new localization stage begins as slip accumulates. It is hypothesized that alternating periods of widening and narrowing correlate to changes of the plate boundaries, manifest in different Euler poles.     

A grid-search inversion method looking for the best classification of polyphase fault/slip data

A novel, grid search-based stress inversion method is developed in this paper to find the global minima in the solution region for the classification of fault/slip data into many single-phase subsets. Exhaustively repetitious grid searches are taken to deal with possible local minima, in a departure from existing grid search-based inversion methods. Two stopping rules, to stop at the abrupt change of the objective function or at the least change of the classification, are adopted in the method to look for the best classification. Much calculation time is saved by using a modified version of conventional grid search. The feasibility of this proposed method is demonstrated by applying it to two artificial examples and two real examples. However, enormous time in calculation is still needed in the case of a data set either with a large number of data or for a large number of assigned subsets.     

Unravelling the evolution of an Alpine to post-glacially active fault in the Swiss Alps

The Gemmi fault is a prominent NW–SE striking lineament that crosses the Gemmi Pass in the central Swiss Alps. A multidisciplinary investigation of this structure that included geological mapping, joint profiling, cathodoluminescence and scanning electron microscopy, stable isotope measurements, luminescence- and U-TH-dating, 3D ground penetrating radar (GPR) surveying and trenching reveals a history of fault movements from the Miocene to the Holocene. The main fault zone comprises a 0.5–3 m thick calcite cataclasite formed during several cycles of veining and brittle deformation. Displaced Cretaceous rock layers show an apparent dextral slip of 10 m along the fault.A detailed study of a small sediment-filled depression that crosses the fault provides evidence for a post-glacial reactivation of the fault. A trench excavated across the fault exposed a Late-Glacial-age loess layer and late Holocene colluvial-like slope-wash deposits that showed evidence for fault displacement of a few centimeters, indicating a recent strike-slip reactivation of the fault. Focal mechanisms of recent instrumentally recorded earthquakes are consistent with our findings that show that the fault at the Gemmi Pass, together with other parallel faults in this area, may be reactivated in today's stress field. Taking together all the observations of its ancient and recent activity, the Gemmi fault can be viewed as a window through geological space and time.     

遥感技术在断裂研究中的应用——以张家口—蓬莱断裂带为例

张家口一蓬莱断裂带是一条北西西向活动断裂带。本文选取ETM+光学影像和SRTM高程影像作为主要数据源,结合研究区已有地质资料研究发现该带断裂构造的北西西向线性特征明显。从水系分布和错断地形等地貌标志判断,该断裂具有左阶组合样式和左行走滑活动特征。据遥感影像综合特征,可将该带分为张家口段、延庆-怀来段和北京一天津段,影像特征分段性显著,并与断裂带的分段性一致。研究结果表明,张家口-蓬莱断裂具有左行走滑的运动学特征,限制或错断北北东或北东向断裂,并且控制该区域的左阶雁列式第四纪盆地群和第四纪冲洪积物的分布。该断裂带各段对不同规模的水系分布和形态影响比较大,且北京-天津地区的华北平原段断裂对第四纪冲洪积扇和沿海地区的贝壳堤的形态和分布有一定的控制作用。地球物理深部数据和野外地质考察资料也证实了遥感解泽的结果,证明遥感技术在探查断裂构造空间展布和活动性鉴定中有着广泛的应用前景。     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Supra-salt normal fault growth during the rise and fall of a diapir: Perspectives from 3D seismic reflection data,Norwegian North Sea

Normal faulting and the deep subsurface flow of salt are key processes controlling the structural development of many salt-bearing sedimentary basins. However, our detailed understanding of the spatial and temporal relationship between normal faulting and salt movement is poor due to a lack of natural examples constraining their geometric and kinematic relationship in three-dimensions. To improve our understanding of these processes, we here use 3D seismic reflection and borehole data from the Egersund Basin, offshore Norway, to determine the structure and growth of a normal fault array formed during the birth, growth and decay of an array of salt structures. We show that the fault array and salt structures developed in response to: (i) Late Triassic-to-Middle Jurassic extension, which involved thick-skinned, sub-salt and thin-skinned supra-salt faulting with the latter driving reactive diapirism; (ii) Early Cretaceous extensional collapse of the walls; and (iii) Jurassic-to-Neogene, active and passive diapirism, which was at least partly coeval with and occurred along-strike from areas of reactive diapirism and wall collapse. Our study supports physical model predictions, showcasing a three-dimensional example of how protracted, multiphase salt diapirism can influence the structure and growth of normal fault arrays.     

Time evolution of hydraulic and electrokinetic parameters around the Nojima fault, Japan, estimated by an electrokinetic method

The Nojima Fault Zone Probe was designed to study the properties and healing processes of the Nojima fault, which is the surface fault rupture of the Hyogo-ken Nanbu earthquake (M7.2) of 1995 (1995 Kobe earthquake). In this project, water injection experiments were conducted in a borehole of 1800 m depth at the Nojima fault. We set up electrodes around the borehole and observed self-potential variations to investigate the magnitude of electrokinetic and hydraulic parameters around the Nojima fault zone. In the 1997 experiment, self-potential variations were in the range of a few to about 20 mV across 320–450 m electrode dipoles with hydraulic pressure variations from 3.5 to 4 MPa. In the 2000 experiment, self-potential variations were in the range of a few to about 85 mV across 160–260 m electrode dipoles with the hydraulic pressure variations from 3 to 4.5 MPa. In the 2003 experiment, self-potential variations were in the range of a few to about 30 mV across 20–80 m electrode dipoles with hydraulic pressure of 4 MPa. These observed self-potential variations were explained well with an electrokinetic effect due to the underground flow of the injected water. From the observed results, we estimated that the ratio of hydraulic parameters (permeability, porosity, and tortuosity) to electrokinetic parameters (zeta potential and dielectric constant) decreased approximately 40% during eight years after the earthquake. This result suggests that the healing process around the fault zone progress.     

3D seismic analysis of the structure and evolution of a salt-influenced normal fault zone: A test of competing fault growth models

In this paper we determine the structure and evolution of a normal fault system by applying qualitative and quantitative fault analysis techniques to a 3D seismic reflection dataset from the Suez Rift, Egypt. Our analysis indicates that the October Fault Zone is composed of two fault systems that are locally decoupled across a salt-bearing interval of Late Miocene (Messinian) age. The sub-salt system offsets pre-rift crystalline basement, and was active during the Late Oligocene-early Middle Miocene. It is composed of four, planar, NW–SE-striking segments that are hard- linked by N–S-striking segments, and up to 2 km of displacement occurs at top basement, suggesting that this fault system nucleated at or, more likely, below this structural level. The supra-salt system was active during the Pliocene-Holocene, and is composed of four, NW–SE-striking, listric fault segments, which are soft-linked by unbreached relay zones. Segments in the supra-salt fault system nucleated within Pliocene strata and have maximum throws of up to 482 m. Locally, the segments of the supra-salt fault system breach the Messinian salt to hard-link downwards with the underlying, sub-salt fault system, thus forming the upper part of a fault zone composed of: (i) a single, amalgamated fault system below the salt and (ii) a fault system composed of multiple soft-linked segments above the salt. Analysis of throw-distance (T-x) and throw-depth (T-z) plots for the supra-salt fault system, isopach maps of the associated growth strata and backstripping of intervening relay zones indicates that these faults rapidly established their lengths during the early stages of their slip history. The fault tips were then effectively ‘pinned’ and the faults accumulated displacement via predominantly downward propagation. We interpret that the October Fault Zone had the following evolutionary trend; (i) growth of the sub-salt fault system during the Oligocene-to-early Middle Miocene; (ii) cessation of activity on the sub-salt fault system during the Middle Miocene-to-?Early Pliocene; (iii) stretching of the sub- and supra-salt intervals during Pliocene regional extension, which resulted in mild reactivation of the sub-salt fault system and nucleation of the segmented supra-salt fault system, which at this time was geometrically decoupled from the sub-salt fault system; and (iv) Pliocene-to-Holocene growth of the supra-salt fault system by downwards vertical tip line propagation, which resulted in downward breaching of the salt and dip-linkage with the sub-salt fault system. The structure of the October Fault Zone and the rapid establishment of supra-salt fault lengths are compatible with the predictions of the coherent fault model, although we note that individual segments in the supra-salt array grew in accordance with the isolated fault model. Our study thereby indicates that both coherent and isolated fault models may be applicable to the growth of kilometre-scale, basin-bounding faults. Furthermore, we highlight the role that fault reactivation and dip-linkage in mechanically layered sequences can play in controlling the three-dimensional geometry of normal faults.     

The nature and origin of off-fault damage surrounding strike-slip fault zones with a wide range of displacements: A field study from the Atacama fault system, northern Chile

Damage surrounding the core of faults is represented by deformation on a range of scales from microfracturing of the rock matrix to macroscopic fracture networks. The spatial distribution and geometric characterization of damage at various scales can help to predict fault growth processes, subsequent mechanics, bulk hydraulic and seismological properties of a fault zone. Within the excellently exposed Atacama fault system, northern Chile, micro- and macroscale fracture densities and orientation surrounding strike-slip faults with well-constrained displacements ranging over nearly 5 orders of magnitude (0.12 m–5000 m) have been analyzed. Faults have been studied that cut granodiorite and have been passively exhumed from 6 to 10 km depth. This allows direct comparison of the damage surrounding faults of different displacements. The faults consist of a fault core and associated damage zone. Macrofractures in the damage zone are predominantly shear fractures orientated at high angles to the faults studied. They have a reasonably well-defined exponential decrease with distance from the fault core. Microfractures are a combination of open, healed, partially healed and fluid inclusion planes (FIPs). FIPs are the earliest set of fractures and show an exponential decrease in fracture density with perpendicular distance from the fault core. Later microfractures do not show a clear relationship of microfracture density with perpendicular distance from the fault core. Damage zone widths defined by the density of FIPs scale with fault displacement but appear to reach a maximum at a few km displacement. One fault, where damage was characterized on both sides of the fault core shows no damage asymmetry. All faults appear to have a critical microfracture density at the fault core/damage zone boundary that is independent of displacement. An empirical relationship for microfracture density distribution with displacement is presented. Preferred FIP orientations have a high angle to the fault close to the fault core and become more diffuse with distance. Models that predict off-fault damage such as a migrating process zone during fault formation, wear from geometrical irregularities and dynamic rupture are all consistent with our data. We conclude it is very difficult to distinguish between them on the basis of field data alone, at least within the limits of this study.     

Application of a modified pattern informatics method to forecasting the locations of future large earthquakes in the central Japan

We propose a modification of the Pattern Informatics (PI) method that has been developed for forecasting the locations of future large earthquakes. This forecast is based on analyzing the space–time patterns of past earthquakes to find possible locations where future large earthquakes are expected to occur. A characteristic of our modification is that the effect of errors in the locations of past earthquakes on the output forecast is reduced. We apply the modified and original methods to seismicity in the central part of Japan and compared the forecast performances. We also invoke the Relative Intensity (RI) of seismic activity and randomized catalogs to constitute null hypotheses. We do statistical tests using the Molchan and Relative Operating Characteristic (ROC) diagrams and the log-likelihoods and show that the forecast for using the modified PI method is generally better than the competing original-PI forecast and the forecasts from the null hypotheses. Using the bootstrap technique with Monte-Carlo simulations, we further confirm that earthquake sequences simulated based on the modified-PI forecast can be statistically the same as the real earthquake sequence so that the forecast is acceptable. The main and innovative science in this paper is the modification of the PI method and the demonstration of its applicability, showing a considerable promise as an intermediate-term earthquake forecasting tool.     

四川汶川地震断裂带科学钻探2号孔(WFSD-2)岩性特征和断裂带的结构

以WFSD-2钻孔岩心为研究对象,通过详细的岩心编录和岩石学、构造地质学等研究,识别出该钻孔岩心具有6段岩性,从上向下依次为彭灌杂岩(0~599.31m)、三叠系须家河组二段(599.31~1211.49m)、彭灌杂岩(1211.49~1679.51m)、三叠系须家河组三段(1679.51~1715.48m)、彭灌杂岩(1715.48~2081.47m)、三叠系须家河组四段(2081.47~2283.56m)。彭灌杂岩主要以花岗岩和火山岩为主,三叠系须家河组沉积岩以砂岩、粉砂岩、泥岩、页岩、煤层(线)和砾岩为主。3套彭灌杂岩与三叠系须家河组沉积岩重复出现,时代较老的岩性段逆冲覆盖在新的地层之上,表明龙门山构造带由一系列逆冲岩片叠置而成。岩心中断裂岩较为发育,主要为断层角砾岩、碎裂岩和断层泥,反映出脆性变形作用的特点。通过对断裂岩的统计分析,厘定了20余条产状不同、规模不等的次级断裂带,断裂带宽度和断裂密度峰值显示FZ600、FZ720、FZ782、FZ817、FZ922、FZ951、FZ1449、FZ1681、FZ2082为主要断裂带,其中FZ1681系规模最大的一条断裂。依据断裂岩的组合特征可以将岩心中断裂带的结构以断层泥为核部划分为两大类:对称型断裂带和不对称型断裂带。根据地表破裂带、WFSD-1钻孔岩心中主滑移带位置的几何关系、岩性分层等因素,可推断汶川地震主滑移带应位于FZ1134、FZ1449或FZ1681之中,同时也暗示该地区经常发生类似汶川地震的大地震活动。研究表明,龙门山地区经历了强烈的构造缩短和快速隆升作用,暗示龙门山地区构造活动非常强烈。     

四川汶川地震断裂带科学钻探2号孔(WFSD-2)岩性特征和断裂带的结构

以WFSD-2钻孔岩心为研究对象,通过详细的岩心编录和岩石学、构造地质学等研究,识别出该钻孔岩心具有6段岩性,从上向下依次为彭灌杂岩(0~599.31m)、三叠系须家河组二段(599.31~1211.49m)、彭灌杂岩(1211.49~1679.51m)、三叠系须家河组三段(1679.51~1715.48m)、彭灌杂岩(1715.48~2081.47m)、三叠系须家河组四段(2081.47~2283.56m)。彭灌杂岩主要以花岗岩和火山岩为主,三叠系须家河组沉积岩以砂岩、粉砂岩、泥岩、页岩、煤层(线)和砾岩为主。3套彭灌杂岩与三叠系须家河组沉积岩重复出现,时代较老的岩性段逆冲覆盖在新的地层之上,表明龙门山构造带由一系列逆冲岩片叠置而成。岩心中断裂岩较为发育,主要为断层角砾岩、碎裂岩和断层泥,反映出脆性变形作用的特点。通过对断裂岩的统计分析,厘定了20余条产状不同、规模不等的次级断裂带,断裂带宽度和断裂密度峰值显示FZ600、 FZ720、 FZ782、 FZ817、 FZ922、 FZ951、FZ1449、 FZ1681、FZ2082为主要断裂带,其中FZ1681系规模最大的一条断裂。依据断裂岩的组合特征可以将岩心中断裂带的结构以断层泥为核部划分为两大类：对称型断裂带和不对称型断裂带。根据地表破裂带、WFSD-1钻孔岩心中主滑移带位置的几何关系、岩性分层等因素,可推断汶川地震主滑移带应位于FZ1134、FZ1449或FZ1681之中,同时也暗示该地区经常发生类似汶川地震的大地震活动。研究表明,龙门山地区经历了强烈的构造缩短和快速隆升作用,暗示龙门山地区构造活动非常强烈。     

A forward test of the Decelerating–Accelerating Seismic Strain model to western south and central America

Global observations show that strong mainshocks are preceded by decelerating preshocks which occur in the focal (seismogenic) region of the ensuing mainshock and by accelerating preshocks which occur in a broader (critical) region of the mainshock. Predictive properties of these preshocks have been expressed by empirical relations supported by theory and form the Decelerating–Accelerating Seismic Strain (D–AS) model. A respective algorithm has been developed which is used to identify the critical and seismogenic region and estimate (predict) the corresponding ensuing mainshock. In the present work a forward test of this model is performed by attempting intermediate-term prediction of future big (M ≥ 7.7) mainshocks along the western coast of south and central America. Three regions of decelerating shocks and three corresponding regions of accelerating shocks have been identified. The parameters (origin time, magnitude, epicenter coordinates) as well as their uncertainties have been estimated (predicted) for the corresponding probably ensuing three mainshocks. This forward test allows an objective evaluation of the model's ability for an intermediate-term prediction of strong shallow mainshocks.     

Observations of wind and waves in the central Bay of Bengal during BOBMEX-99 and their effect on mixed layer depth variability due to forced mixing

Time-series wind and wave measurements were carried out onboard INS Sagardhwani in the central Bay of Bengal during BOBMEX-99. Various other marine meteorological and oceanographic measurements relevant to monsoon studies were also collected simultaneously. The observed variations of wind and waves and the associated mixed layer depth (MLD) variability based on both temperature and density criteria at 3 hourly intervals are presented in this paper as a case study. At the time-series location (13‡N, 87‡E) wind varied between 6 and 16m/s and the predominant direction was southwesterly. The significant wave height and period varied from 1.9 to 3.7m and 8 to 13 s respectively. Some of the available statistical predictive methods for the determination of MLD by forced mixing are utilized to test the extent of mechanical mixing within the top layer of water by the local wind and wave activity. The same is extended to formulate a new empirical relation for gross estimation of effective depth within which the sound energy is generally trapped during its transmission in the surface duct. The present analysis aiming for estimation of observed MLD variability (35 to 75 m) using the suggested simple empirical relation reveals that, the mixed layer variability observed during the experiment depends on both local ocean variability as well as remote forcing as reported earlier.     

散粒体滑坡涌浪运动特征与能量转化规律研究

滑坡涌浪灾害造成的伤亡和损失一般远超滑坡灾害本身,已引起广泛关注。针对该类灾害的预测评价一直是研究的难点,精细刻画滑坡涌浪运动特征和能量转化规律是解决问题的关键和前提。因此,以天然碎石模拟散粒体滑坡,建立缩尺的矩形水槽滑坡涌浪三维物理试验模型,研究滑坡体积、速度和水深对堆积体形态和涌浪特征的影响以及滑体与水体能量转化规律。研究结论包括：（1）滑体通过冲击、挤压和抬升水体形成涌浪的特征受滑体规模和形态影响显著;厚度大且速度快的滑体趋向产生非线性过渡波,厚度小且速度低的滑体一般产生非线性震荡波;（2）基于试验提出的体积效应和速度效应揭示了一定水深的涌浪规模与滑坡体积和速度呈正相关规律,水深效应则解释了一定波能的涌浪在波形和速度上的分布差异,结果表明浅水比深水工况涌浪高度平均高出67%、波速平均高出51.17%;（3）基于试验的数值反演总结了滑体势能与波能转化率为1.00%~3.07%,由于三维试验中滑体与水体在水平和环向的扩散造成更多能量耗散,相较二维试验转化率较低。通过探讨散粒体滑坡运动过程和首浪产生、传播及爬高特征,揭示滑体-水体能量转化基本规律,研究成果对滑坡涌浪防灾减灾工作具有一定的理论价值和意义。