Step-over of strike-slip faults and overpressure fluid favor occurrence of foreshocks:Insights from the 1975 Haicheng fore-main-aftershock sequence,China

This study analyzed and summarized in detail the spatial and temporal distributions of earthquakes,tidal responses,focal mechanisms,and stress field characteristics for the M 7.3 Haicheng earthquake sequence in February 1975.The foreshocks are related to the main fault and the conjugate faults surrounding the extension step-over in the middle.The initiation timing of the foreshock clusters and the original time of the mainshock were clearly modulated by the Earth's tidal force and coincided with the peak of dilational volumetric tidal strain.As a plausible and testable hypothesis,we proposed a fluid-driven foreshock model,by which all observed seismicity features can be more reasonably interpreted with respect to the results of existing models.Together with some other known examples,the widely existing step-over along strike-slip faults and associated conjugate faults,especially for extensional ones in the presence of deep fluids,favor the occurrence of short-term foreshocks.Although clustered seismicity with characteristics similar to those of the studied case is not a sufficient and necessary condition for large earthquakes to occur under similar tectonic conditions,it undoubtedly has a warning significance for the criticality of the main fault.Subsequent testing would require quantification of true/false positives/negatives.

Dynamic Vision Enabled Contactless Cross-Domain Machine Fault Diagnosis With Neuromorphic Computing

Dear Editor,This letter presents a novel dynamic vision enabled contactless cross-domain fault diagnosis method with neuromorphic computing.The event-based camera is adopted to capture the machine vibration states in the perspective of vision.

Autonomous Recommendation of Fault Detection Algorithms for Spacecraft

Dear Editor, This letter deals with the problem of algorithm recommendation for online fault detection of spacecraft. By transforming the time series data into distributions and introducing a distribution-aware measure, a principal method is designed for quantifying the detectabilities of fault detection algorithms over special datasets.

The MW5.5 earthquake on August 6,2023,in Pingyuan,Shandong,China:A rupture on a buried fault

On August 6,2023,a magnitude MW5.5 earthquake struck Pingyuan County,Dezhou City,Shandong Province,China.This event was significant as no large earthquakes had been recorded in the region for over a century,and no active fault had been previously identified.This study collects 1309 P-wave arrival times and 866 S-wave arrival times from 74 seismic stations less than 200 km to the epicenter to constrain the spatial distribution of the mainshock and its 125 early aftershocks by the double difference earthquake relocation method,and selects 864 P-waveforms from 288 stations located within 800 km of the epicenter to constrain the focal mechanism solution of the mainshock through centroid moment tensor inversion.The relocation and the inversion indicate,the Pingyuan MW5.5 earthquake was caused by a rupture on a buried fault,likely an extensive segment of the Gaotang fault.This buried fault exhibited a dip of approximately 75°to the northwest,with a strike of 222°,similar to the Gaotang fault.The rupture initiated at the depth of 18.6 km and propagated upward and northeastward.However,the ground surface was not broken.The total duration of the rupture was~6.0 s,releasing the scalar moment of 2.5895×1017 N·m,equivalent to MW5.54.The moment rate reached the maximum only 1.4 seconds after the rupture initiation,and the 90%scalar moment was released in the first 4.6 s.In the first 1.4 seconds of the rupture process,the rupture velocity was estimated to be 2.6 km/s,slower than the local S-wave velocity.As the rupture neared its end,the rupture velocity decreased significantly.This study provides valuable insights into the seismic characteristics of the Pingyuan MW5.5 earthquake,shedding light on the previously unidentified buried fault responsible for the seismic activity in the region.Understanding the behavior of such faults is crucial for assessing seismic hazards and enhancing earthquake preparedness in the future.

Instability mechanism of mining roadway passing through fault at different angles in kilometre-deep mine and control measures of roof cutting and NPR cables

The angle α between the fault strike and the axial direction of the roadway produces different damage characteristics. In this paper, the research methodology includes theoretical analyses, numerical simulations and field experiments in the context of the Daqiang coal mine located in Shenyang, China. The stability control countermeasure of "pre-splitting cutting roof + NPR anchor cable"(PSCR-NPR) is simultaneously proposed. According to the different deformation characteristics of the roadway, the faults are innovatively classified into three types, with α of type I being 0°-30°, α of type II being 30°-60°, and α of type III being 60°-90°. The full-cycle stress evolution paths during mining roadway traverses across different types of faults are investigated by numerical simulation. Different pinch angles α lead to high stress concentration areas at different locations in the surrounding rock. The non-uniform stress field formed in the shallow surrounding rock is an important reason for the instability of the roadway. The pre-cracked cut top shifted the high stress region to the deep rock mass and formed a low stress region in the shallow rock mass. The high prestressing NPR anchor cable transforms the non-uniform stress field of the shallow surrounding rock into a uniform stress field. PSCR-NPR is applied in the fault-through roadway of Daqiang mine. The low stress area of the surrounding rock was enlarged by 3-7 times, and the cumulative convergence was reduced by 45%-50%. It provides a reference for the stability control of the deep fault-through mining roadway.

Expert Experience and Data-Driven Based Hybrid Fault Diagnosis for High-SpeedWire Rod Finishing Mills

The reliable operation of high-speed wire rod finishing mills is crucial in the steel production enterprise.As complex system-level equipment,it is difficult for high-speed wire rod finishing mills to realize fault location and real-time monitoring.To solve the above problems,an expert experience and data-driven-based hybrid fault diagnosis method for high-speed wire rod finishing mills is proposed in this paper.First,based on its mechanical structure,time and frequency domain analysis are improved in fault feature extraction.The approach of combining virtual value,peak value with kurtosis value index,is adopted in time domain analysis.Speed adjustment and side frequency analysis are proposed in frequency domain analysis to obtain accurate component characteristic frequency and its corresponding sideband.Then,according to time and frequency domain characteristics,fault location based on expert experience is proposed to get an accurate fault result.Finally,the proposed method is implemented in the equipment intelligent diagnosis system.By taking an equipment fault on site,for example,the effectiveness of the proposed method is illustrated in the system.

Label Recovery and Trajectory Designable Network for Transfer Fault Diagnosis of Machines With Incorrect Annotation

The success of deep transfer learning in fault diagnosis is attributed to the collection of high-quality labeled data from the source domain.However,in engineering scenarios,achieving such high-quality label annotation is difficult and expensive.The incorrect label annotation produces two negative effects:1)the complex decision boundary of diagnosis models lowers the generalization performance on the target domain,and2)the distribution of target domain samples becomes misaligned with the false-labeled samples.To overcome these negative effects,this article proposes a solution called the label recovery and trajectory designable network(LRTDN).LRTDN consists of three parts.First,a residual network with dual classifiers is to learn features from cross-domain samples.Second,an annotation check module is constructed to generate a label anomaly indicator that could modify the abnormal labels of false-labeled samples in the source domain.With the training of relabeled samples,the complexity of diagnosis model is reduced via semi-supervised learning.Third,the adaptation trajectories are designed for sample distributions across domains.This ensures that the target domain samples are only adapted with the pure-labeled samples.The LRTDN is verified by two case studies,in which the diagnosis knowledge of bearings is transferred across different working conditions as well as different yet related machines.The results show that LRTDN offers a high diagnosis accuracy even in the presence of incorrect annotation.

Analysis of faulting destruction and water supply pipeline damage from the first mainshock of the February 6,2023 Türkiye earthquake doublet

In 2023,two consecutive earthquakes exceeding a magnitude of 7 occurred in Türkiye,causing severe casualties and economic losses.The damage to critical urban infrastructure and building structures,including highways,railroads,and water supply pipelines,was particularly severe in areas where these structures intersected the seismogenic fault.Critical infrastructure projects that traverse active faults are susceptible to the influence of fault movement,pulse velocity,and ground motions.In this study,we used a unique approach to analyze the acceleration records obtained from the seismic station array(9 strong ground motion stations)located along the East Anatolian Fault(the seismogenic fault of the MW7.8 mainshock of the 2023 Türkiye earthquake doublet).The acceleration records were filtered and integrated to obtain the velocity and displacement time histories.We used the results of an on-site investigation,jointly conducted by China Earthquake Administration and Türkiye’s AFAD,to analyze the distribution of PGA,PGV,and PGD recorded by the strong motion array of the East Anatolian Fault.We found that the maximum horizontal PGA in this earthquake was 3.0 g,and the maximum co-seismic surface displacement caused by the East Anatolian Fault rupture was 6.50 m.As the fault rupture propagated southwest,the velocity pulse caused by the directional effect of the rupture increased gradually,with the maximum PGA reaching 162.3 cm/s.We also discussed the seismic safety of critical infrastructure projects traversing active faults,using two case studies of water supply pipelines in Türkiye that were damaged by earthquakes.We used a three-dimensional finite element model of the PE(polyethylene)water pipeline at the Islahiye State Hospital and fault displacement observations obtained through on-site investigation to analyze pipeline failure mechanisms.We further investigated the effect of the fault-crossing angle on seismic safety of a pipeline,based on our analysis and the failure performance of the large-diameter Thames Water pipeline during the 1999 Kocaeli earthquake.The seismic method of buried pipelines crossing the fault was summarized.

Path-Following Control With Obstacle Avoidance of Autonomous Surface Vehicles Subject to Actuator Faults

This paper investigates the path-following control problem with obstacle avoidance of autonomous surface vehicles in the presence of actuator faults,uncertainty and external disturbances.Autonomous surface vehicles inevitably suffer from actuator faults in complex sea environments,which may cause existing obstacle avoidance strategies to fail.To reduce the influence of actuator faults,an improved artificial potential function is constructed by introducing the lower bound of actuator efficiency factors.The nonlinear state observer,which only depends on measurable position information of the autonomous surface vehicle,is used to address uncertainties and external disturbances.By using a backstepping technique and adaptive mechanism,a path-following control strategy with obstacle avoidance and fault tolerance is designed which can ensure that the tracking errors converge to a small neighborhood of zero.Compared with existing results,the proposed control strategy has the capability of obstacle avoidance and fault tolerance simultaneously.Finally,the comparison results through simulations are given to verify the effectiveness of the proposed method.

Deep-large faults controlling on the distribution of the venting gas hydrate system in the middle of the Qiongdongnan Basin, South China Sea

Many locations with concentrated hydrates at vents have confirmed the presence of abundant thermogenic gas in the middle of the Qiongdongnan Basin(QDNB).However,the impact of deep structures on gasbearing fluids migration and gas hydrates distribution in tectonically inactive regions is still unclear.In this study,the authors apply high-resolution 3D seismic and logging while drilling(LWD)data from the middle of the QDNB to investigate the influence of deep-large faults on gas chimneys and preferred gasescape pipes.The findings reveal the following:(1)Two significant deep-large faults,F1 and F2,developed on the edge of the Songnan Low Uplift,control the dominant migration of thermogenic hydrocarbons and determine the initial locations of gas chimneys.(2)The formation of gas chimneys is likely related to fault activation and reactivation.Gas chimney 1 is primarily arises from convergent fluid migration resulting from the intersection of the two faults,while the gas chimney 2 benefits from a steeper fault plane and shorter migration distance of fault F2.(3)Most gas-escape pipes are situated near the apex of the two faults.Their reactivations facilitate free gas flow into the GHSZ and contribute to the formation of fracture‐filling hydrates.

Semi-analytical solution for mechanical analysis of tunnels crossing strike-slip fault zone considering nonuniform fault displacement and uncertain fault plane position

The tunnel subjected to strike-slip fault dislocation exhibits severe and catastrophic damage.The existing analysis models frequently assume uniform fault displacement and fixed fault plane position.In contrast,post-earthquake observations indicate that the displacement near the fault zone is typically nonuniform,and the fault plane position is uncertain.In this study,we first established a series of improved governing equations to analyze the mechanical response of tunnels under strike-slip fault dislocation.The proposed methodology incorporated key factors such as nonuniform fault displacement and uncertain fault plane position into the governing equations,thereby significantly enhancing the applicability range and accuracy of the model.In contrast to previous analytical models,the maximum computational error has decreased from 57.1%to 1.1%.Subsequently,we conducted a rigorous validation of the proposed methodology by undertaking a comparative analysis with a 3D finite element numerical model,and the results from both approaches exhibited a high degree of qualitative and quantitative agreement with a maximum error of 9.9%.Finally,the proposed methodology was utilized to perform a parametric analysis to explore the effects of various parameters,such as fault displacement,fault zone width,fault zone strength,the ratio of maximum fault displacement of the hanging wall to the footwall,and fault plane position,on the response of tunnels subjected to strike-slip fault dislocation.The findings indicate a progressive increase in the peak internal forces of the tunnel with the rise in fault displacement and fault zone strength.Conversely,an augmentation in fault zone width is found to contribute to a decrease in the peak internal forces.For example,for a fault zone width of 10 m,the peak values of bending moment,shear force,and axial force are approximately 46.9%,102.4%,and 28.7% higher,respectively,compared to those observed for a fault zone width of 50 m.Furthermore,the position of the peak internal forces is influenced by variations in the ratio of maximum fault displacement of the hanging wall to footwall and the fault plane location,while the peak values of shear force and axial force always align with the fault plane.The maximum peak internal forces are observed when the footwall exclusively bears the entirety of the fault displacement,corresponding to a ratio of 0:1.The peak values of bending moment,shear force,and axial force for the ratio of 0:1 amount to approximately 123.8%,148.6%,and 111.1% of those for the ratio of 0.5:0.5,respectively.

Formation mechanism of fault accommodation zones under combined stress in graben basin:Implications from geomechanical modeling

A fault accommodation zone is a type of structure that is defined as regulating displacement and strain between faults structure.Increasing numbers of fault accommodation zones are being identified in graben basins,indicating the potential exploration target and petroleum accumulation areas.This study aims to analyze the formation mechanism and development of fault accommodation zones under combined stress by a numerical simulation method considering geomechanical modeling.Using three-dimensional(3-D)seismic interpretation and fractal dimension method,exampled with the Dongxin fault zone,the fault activity and fault combination pattern were conducted to quantitatively characterize the activity difference in fault accommodation zones.Combined with mechanical experiment test,a geomehcanical model was established for fault accommodation zones in a graben basin.Integrating the paleostress numerical simulations and structural physical simulation experiment,the developmental characteristics and genetic mechanism of fault accommodation zones were summarized.Influenced by multi movements and combined stresses,three significant tectonic evolution stages of the Dongxing Fault Zone(DXFZ)were distinguished:During the E_(s)^(3)sedimentary period,the large difference in the stress,strain,and rupture distribution in various faults were significant,and this stage was the key generation period for the prototype of the DXFZ,including the FAZ between large-scale faults.During the E_(s)^(2)sedimentary period,the EW-trending symmetric with opposite dipping normal faults and the NE-SW trending faults with large scale were furtherly developed.The junction area of two secondary normal faults were prone to be ruptured,performing significant period for inheriting and developing characteristics of fault accommodation zones.During the Es1 sedimentary period,the high-order faults in the DXFZ exhibited the obvious fault depressions and strike-slip activity,and the fault accommodation zones were furtherly inherited and developed.This stage was the molded and formative period of the FAZ,the low-order faults,and the depression in the DXFZ.

Advances in seismological methods for characterizing fault zone structure

Large earthquakes frequently occur along complex fault systems.Understanding seismic rupture and long-term fault evolution requires constraining the geometric and material properties of fault zone structures.We provide a comprehensive overview of recent advancements in seismological methods used to study fault zone structures,including seismic tomography,fault zone seismic wave analysis,and seismicity analysis.Observational conditions limit our current ability to fully characterize fault zones,for example,insufficient imaging resolution to discern small-scale anomalies,incomplete capture of crucial fault zone seismic waves,and limited precision in event location accuracy.Dense seismic arrays can overcome these limitations and enable more detailed investigations of fault zone structures.Moreover,we present new insights into the structure of the Anninghe-Xiaojiang fault zone in the southeastern margin of the Qinghai-Xizang Plateau based on data collected from a dense seismic array.We found that utilizing a dense seismic array can identify small-scale features within fault zones,aiding in the interpretation of fault zone geometry and material properties.

Characterization and spatial analysis of coseismic landslides triggered by the Luding Ms 6.8 earthquake in the Xianshuihe fault zone, Southwest China

On September 5, 2022, a magnitude Ms 6.8 earthquake occurred along the Moxi fault in the southern part of the Xianshuihe fault zone located in the southeastern margin of the Tibetan Plateau,resulting in severe damage and substantial economic loss. In this study, we established a coseismic landslide database triggered by Luding Ms 6.8 earthquake, which includes 4794 landslides with a total area of 46.79 km^(2). The coseismic landslides primarily consisted of medium and small-sized landslides, characterized by shallow surface sliding. Some exhibited characteristics of high-position initiation resulted in the obstruction or partial obstruction of rivers, leading to the formation of dammed lakes. Our research found that the coseismic landslides were predominantly observed on slopes ranging from 30° to 50°, occurring at between 1000 m and 2500 m, with slope aspects varying from 90° to 180°. Landslides were also highly developed in granitic bodies that had experienced structural fracturing and strong-tomoderate weathering. Coseismic landslides concentrated within a 6 km range on both sides of the Xianshuihe and Daduhe fault zones. The area and number of coseismic landslides exhibited a negative correlation with the distance to fault lines, road networks, and river systems, as they were influenced by fault activity, road excavation, and river erosion. The coseismic landslides were mainly distributed in the southeastern region of the epicenter, exhibiting relatively concentrated patterns within the IX-degree zones such as Moxi Town, Wandong River basin, Detuo Town to Wanggangping Township. Our research findings provide important data on the coseismic landslides triggered by the Luding Ms 6.8 earthquake and reveal the spatial distribution patterns of these landslides. These findings can serve as important references for risk mitigation, reconstruction planning, and regional earthquake disaster research in the earthquake-affected area.

An Insight Survey on Sensor Errors and Fault Detection Techniques in Smart Spaces

The widespread adoption of the Internet of Things (IoT) has transformed various sectors globally, making themmore intelligent and connected. However, this advancement comes with challenges related to the effectiveness ofIoT devices. These devices, present in offices, homes, industries, and more, need constant monitoring to ensuretheir proper functionality. The success of smart systems relies on their seamless operation and ability to handlefaults. Sensors, crucial components of these systems, gather data and contribute to their functionality. Therefore,sensor faults can compromise the system’s reliability and undermine the trustworthiness of smart environments.To address these concerns, various techniques and algorithms can be employed to enhance the performance ofIoT devices through effective fault detection. This paper conducted a thorough review of the existing literature andconducted a detailed analysis.This analysis effectively links sensor errors with a prominent fault detection techniquecapable of addressing them. This study is innovative because it paves theway for future researchers to explore errorsthat have not yet been tackled by existing fault detection methods. Significant, the paper, also highlights essentialfactors for selecting and adopting fault detection techniques, as well as the characteristics of datasets and theircorresponding recommended techniques. Additionally, the paper presents amethodical overview of fault detectiontechniques employed in smart devices, including themetrics used for evaluation. Furthermore, the paper examinesthe body of academic work related to sensor faults and fault detection techniques within the domain. This reflectsthe growing inclination and scholarly attention of researchers and academicians toward strategies for fault detectionwithin the realm of the Internet of Things.

Estimation of Displacement and Extension due to Reverse Drag of Normal Faults: Forward Method

In the case of reverse drag of normal faulting, the displacement and horizontal extension are determined based on the established equations for the three mechanisms: rigid body, vertical shear and inclined shear. There are three sub-cases of basal detachment for the rigid body model: horizontal detachment, antithetic detachment and synthetic detachment. For the rigid body model, the established equations indicate that the total displacement on the synthetic base (Dt2) is the largest, that on the horizontal base (Dt1) is moderate, and that on the antithetic base (Dt3) is the smallest. On the other hand, the value of (Dt1) is larger than the displacement for the vertical shear (Dt4). The value of (Dt1) is larger than or less than the displacement for the inclined shear (Dt5) depending on the original fault dip δ0, bedding angle θ, and the angle of shear direction β. For all original parameters, the value of Dt5 is less than the value of Dt4. Also, by comparing three rotation mechanisms, we find that the inclined shear produces largest extension, the rigid body model with horizontal detachment produces the smallest extension, and the vertical shear model produces moderate extension.

A Cascading Fault Path Prediction Method for Integrated Energy Distribution Networks Based on the Improved OPA Model under Typhoon Disasters

In recent times,the impact of typhoon disasters on integrated energy active distribution networks(IEADNs)has received increasing attention,particularly,in terms of effective cascading fault path prediction and enhanced fault recovery performance.In this study,we propose a modified ORNL-PSerc-Alaska(OPA)model based on optimal power flow(OPF)calculation to forecast IEADN cascading fault paths.We first established the topology and operational model of the IEADNs,and the typical fault scenario was chosen according to the component fault probability and information entropy.The modified OPA model consisted of two layers:An upper-layer model to determine the cascading fault location and a lower-layer model to calculate the OPF by using Yalmip and CPLEX and provide the data to update the upper-layer model.The approach was validated via the modified IEEE 33-node distribution system and two real IEADNs.Simulation results showed that the fault trend forecasted by the novel OPA model corresponded well with the development and movement of the typhoon above the IEADN.The proposed model also increased the load recovery rate by>24%compared to the traditional OPA model.

Experimental investigation on the origin of carbonaceous materials in the fault zone of the Wenchuan earthquake

Carbonaceous materials in seismic fault zones may considerably influence seismic fault slip;however,the formation mechanism of carbonaceous materials remains unclear.In this study,we proposed a novel hypothesis for the formation of carbonaceous materials in fault gouge.Thus,we conducted a CO_(2) hydrogenation experiment in a high-temperature reactor at a co-seismic temperature,with fault gouge formed during the Wenchuan earthquake as the catalyst.Our experimental results demonstrate that carbonaceous materials in fault zones are formed on the fault gouge during the chemical reaction process,suggesting that the carbonaceous materials are possibly generated from the catalytic hydrogenation of CO_(2),followed by thermal cracking of its products.The results of this study provide a theoretical basis for understanding fault behavior and earthquake physics.

A Fault Detection Method for Electric Vehicle Battery System Based on Bayesian Optimization SVDD Considering a Few Faulty Samples

Accurate and reliable fault detection is essential for the safe operation of electric vehicles.Support vector data description(SVDD)has been widely used in the field of fault detection.However,constructing the hypersphere boundary only describes the distribution of unlabeled samples,while the distribution of faulty samples cannot be effectively described and easilymisses detecting faulty data due to the imbalance of sample distribution.Meanwhile,selecting parameters is critical to the detection performance,and empirical parameterization is generally timeconsuming and laborious and may not result in finding the optimal parameters.Therefore,this paper proposes a semi-supervised data-driven method based on which the SVDD algorithm is improved and achieves excellent fault detection performance.By incorporating faulty samples into the underlying SVDD model,training deals better with the problem of missing detection of faulty samples caused by the imbalance in the distribution of abnormal samples,and the hypersphere boundary ismodified to classify the samplesmore accurately.The Bayesian Optimization NSVDD(BO-NSVDD)model was constructed to quickly and accurately optimize hyperparameter combinations.In the experiments,electric vehicle operation data with four common fault types are used to evaluate the performance with other five models,and the results show that the BO-NSVDD model presents superior detection performance for each type of fault data,especially in the imperceptible early and minor faults,which has seen very obvious advantages.Finally,the strong robustness of the proposed method is verified by adding different intensities of noise in the dataset.

Geometry and 3D seismic characterisation of post-rift normal faults in the Pearl River Mouth Basin,northern South China Sea

Based on high-resolution 3D seismic data acquired in the Pearl(Zhujiang)River Mouth Basin of the northern South China Sea,this study investigated the geometry,spatial extension,and throw distribution of the post-rift normal fault through detailed seismic interpretation and fault modeling.A total of 289 post-rift normal faults were identified in the study area and can be classified into four types:(1)isolated normal faults above the carbonate platform;(2)isolated normal faults cutting through the carbonate platform;(3)conjugate normal faults,and(4)connecting normal faults.Throw distribution analysis on the fault planes show that the vertical throw profiles of most normal fault exhibit flat-topped profiles.Isolated normal faults above the carbonate platform exhibit roughly concentric ellipses with maximum throw zones in the central section whereas the normal faults cutting through the carbonate platform miss the lowermost section due to the chaotic seismic reflections in the interior of the carbonate platform.The vertical throws of conjugate normal faults anomalously decrease toward their intersection region on the fault plane whereas the connecting normal faults present two maximum throw zones in the central section of the fault plane.According to the symmetric elliptical distribution model of fault throw,an estimation was made indicating that normal faults cutting through the carbonate platform extended downward between-1308 s and-1780 s(two-way travel time)in depth and may not penetrate the entire Liuhua carbonate platform.Moreover,it is observed that the distribution of karst caves on the top of the carbonate platform disaccord with those of hydrocarbon reservoirs and the post-rift normal faults cutting through the carbonate platform in the study area.We propose that these karst caves formed most probably by corrosive fluids derived from magmatic activities during the Dongsha event,rather than pore waters or hydrocarbons.