The 1987–88 Alaskan Bight tsunamis: Deep ocean data and model comparisons

Excellent deep ocean records have been obtained of two tsunamis recently generated in the Alaskan Bight on 30 November 1987 and 6 March 1988, providing the best available data set to date for comparison with tsunami generation/propagation models. Simulations have been performed with SWAN, a nonlinear shallow water numerical model, using source terms estimated by a seafloor deformation model based on the rectangular fault plane formalism. The tsunami waveform obtained from the model is quite sensitive to the specific source assumed. Significant differences were found between the computations and observations of the 30 November 1987 tsunami, suggesting inadequate knowledge of the source characteristics. Fair agreement was found between the data and the model for the first few waves of the 6 March 1988 tsunami. Model estimates of the seismic moment and total slip along the fault plane are also in fair agreement with those derived from the published Harvard centroid solution for the 6 March 1988 event, implying that the computed seafloor deformation does bear some similarity to the actual source.     

Insights on the small tsunami from January 28, 2020, Caribbean Sea MW7.7 earthquake by numerical simulation and spectral analysis

A major left-lateral strike-slip Mw7.7 earthquake occurred in the vicinity of the Caribbean Sea on January 28, 2020. As a result, a small-scale tsunami was generated. The properties of the seismogenic source were described using observational data gathered for the earthquake and tsunami, as well as information on the regional tectonic setting. The tsunami was simulated with the COMCOT model and Okada’s dislocation model from finite fault solutions for M_W7.7 Caribbean Sea earthquakes published by the United States Geological Survey. The simulation results were compared to tide gauge records to validate whether the seafloor’s vertical displacement generated by the strike-slip fault caused a small-scale tsunami. We conducted a spectral analysis of the tsunami to better understand the characteristics of tsunami records. The tsunami simulation results showed that the co-seismic vertical displacement caused by a strike-slip M_W7.7 earthquake could have contributed to the small-scale tsunami, but the anomalously large high-frequency tsunami waves recorded by the George Town tide gauge 11 min after the earthquake were unrelated to the earthquake-generated tsunami. According to the spectrum analysis, the predominant period of noticeable high-frequency tsunami waves recorded by the George Town tide gauge occurred only two minutes after the earthquake. This indicates that the source of the small-scale tsunami was close to the George Town station and the possible tsunami source was 150 km away from George Town station. These facts suggest that a submarine landslide was caused by the strike-slip earthquake. The comprehensive analysis showed that the small-scale tsunami was not caused solely by co-seismic seafloor deformation from the strike-slip event but that an earthquake-triggered submarine landslide was the primary cause. Therefore, the combined impact of two sources led to the small-scale tsunami.

     

基岩岛屿地下水流场数值模拟研究

岛屿是我国领土的重要组成部分,对国家安全和国防军事意义非凡。岛屿的地下水资源尤为珍贵,地下水流场模拟是研究地下水分布规律的重要手段。水文地质条件复杂、可利用的观测井较少等原因,造成了基岩岛屿水文地质模型概化精度不高、初始条件难以获取等问题。为克服基岩岛屿地下水流场模拟的诸多困难,将珠海外伶仃岛作为研究区,利用数字高程模型数据开展地质地貌建模; 利用探地雷达法、直流电法与地质分析对岛屿进行探测,获取了地层数据; 采用地下水遥感评估法,利用实测井位数据,确定了地下水的初始水位,进而对基岩岛屿地下水流场进行建模; 最终,通过对外伶仃岛地下水流场的数值模拟得出地下水模拟流场图。岛上多个测点的探测水位值与模拟水位相关性较好,其拟合优度R²为0.872 2。由此可见,综合遥感、物探、水文地质手段等技术方法获取的数据,采用地下水模拟软件或程序实现基岩岛屿地下水流场的数值模拟,是基岩岛屿地下水资源研究的一个有效方法。     

地面核磁共振模型约束反演含水层参数

导电模型的地面核磁共振感应电动势是含水量的非线性函数.引入模型约束的迭代反演方法求解该非线性问题的反演.在反演过程中,根据理论公式计算迭代过程中灵敏度矩阵,并采用平坦模型和光滑模型两种约束.对均匀半空间、层状导电模型和实际数据进行了反演模拟,结果表明,模型约束迭代反演方法能从地面核磁共振感应电动势获得含水层较为合理的含水量及分布,且该结果可以从作为初始模型的均匀含水量分布反演得到.对无噪音数据,平坦模型和光滑模型约束对反演的含水量分布影响不大;但当数据存在一定的噪音时,平坦模型约束将比光滑模型约束获得更为精确的含水层参数.     

Tsunami simulation due to seaquake at Manila Trench and Sulu Trench

Seaquake is a phenomenon where there are water disturbance at the sea, caused by earthquake or submarine eruption. The scope of this study focuses on tsunami simulation due to Manila Trench and Sulu Trench seaquake which is prone to harm Malaysia offshore areas. Manila Trench is a highly potential earthquake source that can generate tsunami in South China Sea. Meanwhile, Sulu Trench could be a threat to east of Sabah offshore areas. In this study, TUNA-M2 model was utilized to perform tsunami simulation at South China Sea and Sulu Sea. TUNA-M2 model applied Okada source model to create tsunami generation due to earthquake. It utilized linear shallow water equation during tsunami propagation with its radiant boundary condition. Five simulations performed at each study region. Forecast points at South China Sea areas were divided into three separate locations which are at the Peninsular Malaysia, west of Sabah and Sarawak offshore areas. Forecast points at Sulu Sea were focused at the east of Sabah offshore areas. This paper will present the simulation results of tsunami wave height and arrival time at various forecast points. The findings of this study show that the range of tsunami wave height at Sulu Sea is higher than that of South China Sea. The tsunami arrival time at Sulu Sea is less than South China Sea. It can be concluded that Sulu Sea poses worse tsunami threat than South China Sea to the Malaysian offshore areas.     

Simulation of the Arabian Sea Tsunami propagation generated due to 1945 Makran Earthquake and its effect on western parts of Gujarat (India)

The 1945 Tsunami generated due to Makran Earthquake in the Arabian Sea was the most devastating tsunami in the history of the Arabian Sea and caused severe damage to property and loss of life. It occurred on 28th November 1945, 21:56 UTC (03:26 IST) with a magnitude of 8.0 (M _w), originating off the Makran Coast of Pakistan in the Arabian Sea. It has impacted as far as Mumbai in India and was noticed up to Karvar Coast, Karnataka. More than 4,000 people were killed as a result of the earthquake and the tsunami. In this paper an attempt is made for a numerical simulation of the tsunami generation from the source, its propagation into the Arabian Sea and its effect on the western coast of India through the use of a numerical model, referred to as Tunami-N2. The present simulation is carried out for a duration of 300 min. It is observed from the results that the simulated arrival time of tsunami waves at the western coast of India is in good agreement with the available data sources. The paper also presents run-up elevation maps prepared using Shuttle Radar Topographic Mission (SRTM) data, showing the possible area of inundation due to various wave heights along different parts of the Gujarat Coast. Thus, these results will be useful in planning the protection measures against inundation due to tsunami and in the implementation of a warning system.     

Modelling of Tsunamis Generated by Pyroclastic Flows (Ignimbrites)

Pyroclastic flows entering the sea played a major role in generating the largest tsunamiwaves, arising from the 1883 eruption of Krakatau, Indonesia, which caused a considerabledeath toll, most deaths resulting from the tsunamis. The potential exists for similar eventsto occur in Indonesia and New Zealand.Processes leading to tsunami generation by pyroclastic flows, especially those associatedwith Krakatau-type eruptions, are reviewed. The major processes include:1. Deposition at the shoreline causing a lateral displacement as the zone of depositionmoves offshore.2. Upward and lateral displacement of water caused by the propagation of a watersupported mass-flow.3. Downward and lateral displacement of water caused by the sinking of debris from a segregated flow travelling over the water surface.4. Upward displacement of a large volume of water due to the deposition of acaldera-infill ignimbrite or pyroclastic flow deposit.The pyroclastic flow is modelled as a horizontal piston forcingwater displacement. The flow behaves as a wedge of material displacingseawater horizontally and vertically as it moves outwards from the source.Individual pyroclastic flows are treated as linear features that travel alonga specific direction from the volcano, exhibiting limited lateral spreading.The event duration for the formation of a large pyroclastic flow and thedeposition of the ignimbrite is taken as 200–400 s, with flow velocitiesdependent on the volume of material erupted.For simulations it is assumed that the ignimbrite deposit is elliptical with relativelyuniform thickness and the principal axis orientated along the flow direction. Therefore the tsunami is generated by defining an elliptical source region and defining an effective displacement behaviour at each node within that region. The effective displacement is defined by a start time, a finish time and a vertical velocity. These three parameters determine when the seafloor starts to rise and how far it travels during a model time step. The result is a seafloor disturbance that propagates away from the source.The major difficulty with this approach is determination of the appropriate verticalvelocity. With a real pyroclastic flow the effective vertical velocity at any point isvery high. However the model needs to average the displacement spatially andtemporally. Accordingly we apply the model to pyroclastic flows from Mayor Island, New Zealand to examine the influence of model parameters. To further calibrate the numerical model this study is being undertaken in conjunction with physical modelling of the Krakatau 1883 eruption at the Indonesian Tsunami Research Center, BPPT, Jakarta. Historical data will also be used to refine and calibrate the pyroclastic flow model.     

北京房山世界地质公园中元古界雾迷山组地震-海啸序列及地质特征——以野三坡园区为例

野三坡园区中元古界(1800~1200Ma B.P)位于燕辽裂陷槽南部中轴断层南支(古紫荆关断裂)之东盘。野三坡雾迷山组笫二段白云岩中的地震序列自下而上分为：A.阶梯状断层与震碎角砾岩变形单元;B.微褶皱-微断层变形单元;C.液化均一变形单元;D.液化卷曲变形单元。这4个变形单元分别代表海底之下深度不足10m内已固结成岩—半固结—弱固结沉积层的变形特征,都经历了前震及主震期（P波和S波）不同程度的影响。海啸系列（海啸岩）自下而上为：E.津浪丘状层单元;F.丘-槽构造层单元和G.粒序均一层单元。海啸是在主震期数十秒后发生的。海底瞬时大幅度抬升,然后突然大幅度下降,使外海海水涌入,引发海啸,形成丘状层及有关变形。余震阶段的震荡流与一次沉降事件分别形成F和G单元。本区中轴断层（古紫荆关断层）是一条海底直立的断层带,对雾迷山组二段中的地震-海啸及沉积过程均起到激发与控制作用。此次震中位于中轴断层带西缘的白石山,震级为7.0~7.5里氏级,是震源浅、裂度大的海底地震。     

龚家方4号斜坡涌浪数值模拟分析

对于库区滑坡来说,不能只考虑滑坡体本身造成的灾害,还要考虑滑坡体引起的涌浪灾害,为了研究滑坡涌浪的传播、衰减规律,在Geo-wave软件的技术上,二次开发形成FAST软件。以三峡库区龚家方4号斜坡为研究对象,分别在175、156、145 m的库水位条件下,在长约23 km、宽约10.4 km的区域内进行涌浪数值模拟,获得涌浪传播模拟数据。经过模拟软件数据处理模块的计算分析,形成了分析涌浪传播规律的一系列图件。对不同水位下涌浪模拟的计算结果进行对比分析发现,随着库水位的下降,滑坡产生的最大涌浪值和在对岸的爬高值都有增长的趋势,但其对航道存在威胁的时间逐渐变短。模拟区各位置的最大波高空间分布形态具有中间内凹、两翼沿岸坡延伸的特征。涌浪传播的急剧衰减区基本分布在涌浪源附近1 km的范围内,涌浪源处的波高越大,单位距离内的涌浪下降高度也越大。由于涌浪在岸边有叠加、壅高现象,建议航道内船只经过地质灾害点附近时应沿江中心快速通行。     

Hydraulic and numerical study on the generation of a subaqueous landslide-induced tsunami along the coast

By carrying out the hydraulic experiments in a one-dimensional open channel and two-dimensional basin, we clarified the process of how a landslide on a uniform slope causes the generation of a tsunami. The effect of the interactive force that occurs between the debris flow layer and the tsunami is significant in the generation of a tsunami. The continuous flow of the debris into the water makes the wave period of the tsunami short. The present experiments apply numerical simulation using the two-layer model with shear stress models on the bottom and interface, and the results are compared. The simulated debris flow shows good agreement with the measured results and ensures the rushing process into the water. We propose that the model use a Manning coefficient of 0.01 for the smooth slope and 0.015 for the rough slope, and a horizontal viscosity of 0.01 m²/s for the landslide; an interactive force of 0.2 for each layer is recommended. The dispersion effect should be included in the numerical model for the propagation from the shore.     

A new landslide-induced tsunami simulation model and its application to the 1792 Unzen-Mayuyama landslide-and-tsunami disaster

By combining landslide dynamics research and tsunami research, we present an integrated series of numerical models quantitatively simulating the complete evolution of a landslide-induced tsunami. The integrated model simulating the landslide initiation and motion uses measured landslide dynamic parameters from a high-stress undrained dynamic-loading ring shear apparatus. It provides the numerical data of a landslide mass entering and moving under water to the tsunami simulation model as the trigger of tsunami. The series of landslide and tsunami simulation models were applied to the 1792 Unzen-Mayuyama megaslide and the ensuing tsunami disaster, which is the largest landslide disaster, the largest volcanic disaster, and the largest landslide-induced tsunami disaster to have occurred in Japan. Both the 1792 megaslide and the tsunami portions of the disaster are well documented, making this an excellent test of the reliability and precision of the new simulation model. The simulated tsunami heights at the coasts well match the historical tsunami heights recorded by “Tsunami-Dome-Ishi” (a stone showing the tsunami reaching point) and memorial stone pillars.     

Far-Field Tsunami Potential and a Real-Time Forecast System for the Pacific Using the Inversion Method

Estimating tsunami potential is anessential part of mitigating tsunami disasters. Weproposed a new method to estimate the far-fieldtsunami potential by assuming faultmodels on the Pacific Rim. We find thata tsunami that generates in the areas wherethere is no tsunami in the history can damagethe Japanese coast. This shows that it isimportant to estimate tsunami potential byassuming fault models other than the pastearthquake data.Another important activity to mitigate tsunamidisasters is to provide appropriatewarnings to coastal communities when dangerfrom a tsunami is imminent. We applied anew inversion method using wavelet transformto a part of the real-time tsunami forecastsystem for the Pacific. Because this inversionmethod does not require fault location, it ispossible to analyze a tsunami in real timewithout all seismic information. In order tocheck the usability of the system, anumerical simulation was executed assuming anearthquake at sea off Taiwan. The correlationcoefficient for the estimated initialwaveform to the assumed one was calculatedto be 0.78. It takes 90 min to capturetime-series waveform data from tsunamigauges and 5 sec to estimate the 2-D initialwaveform using the inversion method. After that,it takes 2 minutes to forecast thetsunami heights at the Japanese coast. Since thesum of these times is less than the 105minutes transit time of the tsunami fromTaiwan to Japan, it is possible to give a warningto the residents before the tsunami attacksthe Japanese coast. Comparing the tsunamiheights forecasted by this system with thosecalculated by the fault model, the averageerror was 0.39 m. The average error ofthe arrival time was 0.007 min.     

Distribution Functions of Tsunami Wave Heights

The problem of describing the distribution functions of tsunami wave heights is discussed. Data on runup heights obtained in field surveys of several tsunamis for the last decade are used to calculate the empirical distribution functions. It is shown that the log-normal distribution describes the observed data well. This means that the irregular topography and coastline are major factors which influence the height distribution. The power distribution related with the geometric decay of the propagated wave is a good approximation for one event (Sulawesi, January 1, 1996) only. Results of a numerical simulation of the tsunami event in the Japan (East) Sea on July 12, 1993 are presented. It is shown that the computed wave height distribution, obtained by using the runup correction in the framework of nonlinear shallow-water theory, is in good agreement with the observed height distribution. Simulations are used to study the transformation of the distribution function on different distances from the source.     

海洋可控源电磁法多参数正演响应特征分析

鉴于介质的色散现象,充分考虑介质电磁参数(ε、μ、σ),获得各层介质中的电磁场和势函数。为了考察岩、矿石电磁参数、海水深度和海底地层厚度等变化时海洋电磁响应特征规律,给出典型地电模型和观测系统参数,借助高精度快速汉克尔滤波系数,采用30点高斯勒让德数值求积方法,计算水平电性源频率域可控源电磁法(CSEM)在海底各观测点电场Ez和磁场Bx的正演响应。结果表明,随海水深度变浅,电场Ez和磁场Bx振幅曲线变化幅度增大,空气波逐渐占据主导地位,当水深超过3km时,观测区内可以忽略空气波的影响;电场Ez分量对海底高阻层引起的异常大,而中低阻层引起的异常小,磁场Bx分量幅值对目标层电导率变化与Ez相同,因此,水平电性源频率域CSEM法不适宜探测海底低阻目标。并且对厚层目标层的探测能力高于薄层。针对岩层中磁导率变化的情况,海底沉积物电性均匀或者呈层状分布时,电场Ez和磁场Bx幅值均受到影响。此外,电场Ez和磁场Bx的幅值基本不受介电常数变化而影响。海洋可控源的探测效果与偏移距的选择有密切关系,有利偏移距范围为3000m^12000m。     

Constitutive modeling of the destructuration and anisotropy of natural soft clay

The mechanical behavior of natural clays is affected by their inherent anisotropy and metastable soil structure. A simple hierarchical model that considers initial anisotropy and destructuration was formulated within the framework of critical state soil mechanics. In the proposed model, stress sensitivity and a destructuration index were introduced to account for the degree of bonding and the rate of destructuration, respectively. An inclined yield surface was used to incorporate the effect of the initial anisotropy. The proposed model can be degenerated to the Modified Cam Clay model by setting the initial stress sensitivity equal to unity and using a horizontal yield surface. Reasonable agreement between the model simulations and the experimental results on a variety of stress paths demonstrated that the proposed model can capture well the deformation behavior of natural clay and reconstituted soil. The model was implemented into the finite element program for the numerical analysis of an embankment on soft clay improved with prefabricated vertical drains. The numerical predictions were compared with the field-measured data in terms of embankment settlement. Additionally, the numerical simulations were analyzed in terms of horizontal displacements, excess pore water pressure, mean effective stress and volumetric strain. All of the simulations and comparisons indicate the importance of considering the effects of plastic anisotropy, interparticle bonding and destructuration caused by loading beyond yield stress and field disturbance in analyzing the behavior of an embankment on natural soft clay.     

黄河三峡库区滑坡诱发的涌浪预测方法

黄河三峡库区的涌浪灾害风险不容忽视,经验公式是宜优先考虑的涌浪快捷评价方法.对黄河三峡焦家崖头2012年2月7日的黄土滑坡和涌浪进行调查,分析了黄土滑坡及涌浪的特征.采用9种涌浪经典计算公式,计算了涌浪的初始浪高、对岸爬高等特征参数.与调查结果对比表明,采用美国土木工程师协会推荐法、水科院算法、Huber and Hager模型和潘家铮算法获取的焦家崖头黄土滑坡诱发的涌浪特征参数均接近实际,其确定的校正系数分别为2.14、1.92、0.6和0.66,对比考虑安全性和经济性后推荐采用潘家铮算法预测黄河三峡的涌浪.     

Discontinuous deformation analysis of super section tunnel surrounding rock stability based on joint distribution simulation

Based on binocular stereo photography measurement technology, the actual distribution information of joints and fractures in the tunnel face were obtained via image processing and feature extraction, and the preliminary evaluation of surrounding rock stability of Laohushan tunnel was conducted according to surrounding rock classification method. Since all available surrounding rock classification approaches didn’t consider the influences of the size effect of tunnel excavation span and unfavorable geologic bodies such as weak-fracture zones, an improved discontinuous deformation analysis (DDARF) method was adopted to conduct a numerical simulation of the deformation and fracture processes of the surrounding rock. Using the traveling wave method, triangle DDA blocks were automatically generated in the calculation zone, and the block boundaries were divided into real joints and virtual joints. Based on the real joint information obtained via aforementioned photographic measurement, the real joints in the tunnel face were dynamically modified, in order to achieve the simulation of joint distribution. The results revealed that the fracture evolution regularity, deformation failure mechanism, and block dropping phenomenon that the DDARF joint simulation model calculated are in good agreement with actual conditions, while those obtained based on conventional models present differences from field situation. Additionally, focusing on the localised rockfall phenomenon of Laohushan tunnel, the crack extension rate was introduced to conduct a quantitative comparison of the rock crack evolution process with or without anchor supporting. The research results offer practical guidance for field construction and anchorage support scheme optimization.     

Patterns of stress at midocean ridges and their offsets due to seafloor subsidence

The effect of the seafloor subsidence on the horizontal stress field is investigated by combining the finite element method with a formulation that allows us to compute the two-dimensional (2D) horizontal stresses arising from isostatically compensated vertical loads. The topographic load created by the elevation of midocean ridges relative to old ocean floor is shown to be a significant source of ridge-parallel tensile stresses. These may predominate over the ridge-perpendicular stresses and explain observations at midocean ridge offsets such as (1) oblique normal faulting at ridge-transform intersections trending up to 60° relative to the ridge axis, and (2) nontransform offsets consisting of structures oriented at 45° relative to the ridge trend. At midocean ridge overlaps, rotation of the ridge-parallel tensile stresses favours rift propagation at more than 45° relative to the ridge trend. It is suggested that propagating rift tips that bend abruptly lead to partially unlocked offsets, and as a result large overlaps may eventually start to rotate and evolve into a microplate.     

Pacific Basin tsunami hazards associated with mass flows in the Aleutian arc of Alaska

We analyze mass-flow tsunami generation for selected areas within the Aleutian arc of Alaska using results from numerical simulation of hypothetical but plausible mass-flow sources such as submarine landslides and volcanic debris avalanches. The Aleutian arc consists of a chain of volcanic mountains, volcanic islands, and submarine canyons, surrounded by a low-relief continental shelf above about 1000–2000 m water depth. Parts of the arc are fragmented into a series of fault-bounded blocks, tens to hundreds of kilometers in length, and separated from one another by distinctive fault-controlled canyons that are roughly normal to the arc axis. The canyons are natural regions for the accumulation and conveyance of sediment derived from glacial and volcanic processes. The volcanic islands in the region include a number of historically active volcanoes and some possess geological evidence for large-scale sector collapse into the sea. Large scale mass-flow deposits have not been mapped on the seafloor south of the Aleutian Islands, in part because most of the area has never been examined at the resolution required to identify such features, and in part because of the complex nature of erosional and depositional processes. Extensive submarine landslide deposits and debris flows are known on the north side of the arc and are common in similar settings elsewhere and thus they likely exist on the trench slope south of the Aleutian Islands. Because the Aleutian arc is surrounded by deep, open ocean, mass flows of unconsolidated debris that originate either as submarine landslides or as volcanic debris avalanches entering the sea may be potential tsunami sources.To test this hypothesis we present a series of numerical simulations of submarine mass-flow initiated tsunamis from eight different source areas. We consider four submarine mass flows originating in submarine canyons and four flows that evolve from submarine landslides on the trench slope. The flows have lengths that range from 40 to 80 km, maximum thicknesses of 400–800 m, and maximum widths of 10–40 km. We also evaluate tsunami generation by volcanic debris avalanches associated with flank collapse, at four locations (Makushin, Cleveland, Seguam and Yunaska SW volcanoes), which represent large to moderate sized events in this region. We calculate tsunami sources using the numerical model TOPICS and simulate wave propagation across the Pacific using a spherical Boussinesq model, which is a modified version of the public domain code FUNWAVE. Our numerical simulations indicate that geologically plausible mass flows originating in the North Pacific near the Aleutian Islands can indeed generate large local tsunamis as well as large transoceanic tsunamis. These waves may be several meters in elevation at distal locations, such as Japan, Hawaii, and along the North and South American coastlines where they would constitute significant hazards.     

Detection and evaluation of horizontal fractures in earth dams using the self-potential method

The self-potential method can be applied to evaluate the degree of water seepage into earth or rockfill dams. Spatial distribution of measured self-potential data can indicate possible anomalous water flow. Phenomena, such as piping, can be modelled as cylindrical bodies. Internal erosion can develop structures, which can be represented by a sphere or point source. Differential settlement in the dam structure occasionally creates horizontal fractures, which require equivalent geometrical body formula. An equation, which allows the calculation of a self-potential profile over a horizontal contact, was developed. This equation can also be applied to the inspection of water flow into horizontal drainage filters. When inverse interpretation returns bodies of relatively small amplitudes, then their probability can be tested statistically. A test, based on cross-correlation between a modelled curve and field data, can be used to evaluate their existence at a given probability level. An acceptance criterion is computed, using the concept of likelihood coefficient. The equation was applied to two case histories. The first is an example of water flow evaluation into a horizontal drainage filter. Several anomalous bodies were interpreted from residual self-potential data, i.e. from the difference between the theoretical response of the filter and the measured self-potential values, then, most of the interpreted bodies were statistically tested. The second case deals with detection and evaluation of a horizontal fissure. The interpreted parameters of the detected body corresponded well to a horizontal fracture found when the water level in the reservoir was lowered.