Electrical resistivity imaging in Shibaozhai area using magnetotelluric sounding:Enlightenment of hydrocarbon reservoirs on the eastern edge of the Sichuan Basin

Shibaozhai area is located in the northeast Sichuan Basin in Southwest China.In recent years,oil and gas reservoirs with abundant reserves have been discovered in nearby areas.Six magnetotelluric(MT)survey lines were arranged to investigate the deep geoelectrical structures and hidden faults that controlled the hydrocarbon migration along the eastern Kaijiang–Liangping trough and Dachiganjing tectonic belt in this area.The MT sounding data for 181 sites,which covered an area of 15×6 km,were collected.Twodimensional transverse-magnetic-mode inversion of the data revealed four geoelectric layers with a total thickness of approximately 10 km.A low resistivity layer with a thickness of 4 km was discovered in the middle of the exploration area.The combination of faults within the study area revealed a northeast-direction thrust-fault system,which was displaced by northwest-direction faults.The regional geological structures played a significant role in controlling hydrocarbon and gas accumulation.The results of this study could be beneficial in analyzing the geological structure of hydrocarbon reservoirs in this area and can be extended to the exploration of areas in the boundary of a basin.

Reverse time migration imaging of tunnels via the finite element method using an unstructured mesh

Wavefield extrapolation is critical in reverse time migration(RTM).The finite diff erence method is primarily used to achieve wavefi eld extrapolation in case of the RTM imaging of tunnels.However,complex tunnel models,including those for karsts and fault fracture zones,are constructed using regular grids with straight curves,which can cause numerical dispersion and reduce the imaging accuracy.In this study,wavefi eld extrapolation was conducted for tunnel RTM using the finite element method,wherein an unstructured mesh was considered to be the body-fi tted partition in a complex model.Further,a Poynting vector calculation equation suitable for the unstructured mesh considered in the fi nite element method was established to suppress the interference owing to low-frequency noise.The tunnel space was considered during wavefi eld extrapolation to suppress the mirror artifacts based on the fl exibility of mesh generation.Finally,the infl uence of the survey layouts(one and two sidewalls)on the tunnel imaging results was investigated.The RTM results obtained for a simple tunnel model with an inclined interface demonstrate that the method based on unstructured meshes can effectively suppress the low-frequency noise and mirror artifacts,obtaining clear imaging results.Furthermore,the two-sidewall tunnel survey layout can be used to accurately obtain the real position of the inclined interface ahead of the tunnel face.The complex tunnel numerical modeling and actual data migration results denote the eff ectiveness of the fi nite element method in which an unstructured mesh is used.

An alternative 2D ICEEMDAN-based denoising method and its application in processing magnetic anomaly data

Due to environmental noise and human factors,magnetic data collected in the field often contain various noises and interferences that significantly affect the subsequent data processing and interpretation.Empirical Mode Decomposition(EMD),an adaptive multiscale analysis method for nonlinear and non-stationary signals,is widely used in geophysical and geodetic data processing.Compared with traditional EMD,Improved Complete Ensemble EMD with Adaptive Noise(ICEEMDAN)is more effective in addressing the problem of mode mixing.Based on the principles of 1D ICEEMDAN,this paper presents an alternative algorithm for 2D ICEEMDAN,extending its application to two-dimensional scenarios.The effectiveness of the proposed approach is demonstrated through synthetic signal experiments,which show that the 2D ICEEMDAN exhibits a weaker mode mixing effect compared to the traditional bidimensional EMD(BEMD)method.Furthermore,to improve the performance of the denoising method based on 2D ICEEMDAN and preserve useful signals in high-frequency components,an improved soft thresholding technique is introduced.Synthetic magnetic anomaly data testing indicates that our denoising method effectively preserves signal continuity and outperforms traditional soft thresholding methods.To validate the practical application of this improved threshold denoising method based on 2D ICEEMDAN,it is applied to ground magnetic survey data in the Yandun area of Xinjiang.The results demonstrate the effectiveness of the method in removing noise while retaining essential information from practical magnetic anomaly data.In particular,practical applications suggest that 2D ICEEMDAN can extract trend signals more accurately than the BEMD.In conclusion,as a potential tool for multi-scale decomposition,the 2D ICEEMDAN is versatile in processing and analyzing 2D geophysical and geodetic data.

Co-seismic deformation for the 2015 M_(W)7.8 Gorkha earthquake(Nepal)using near-field GPS data

Seasonal variations and common mode errors affect the precision of the Global Positioning System(GPS)time series.In this paper,we explore to improve the precision of coordinate time series,thereby providing a better detection of weak or transient deformation signals,particularly co-seismic signals.Based on 97 GPS stations,including the campaign and continuous GPS stations in Nepal and southern Tibet,we first consider seasonal variations and common errors,then obtain co-seismic deformation of the 2015 Gorkha earthquake in Nepal and southern Tibet.Our co-seismic rupture model is characterized by a shallow ramp and a deeper detachment fault,in agreement with the relocated aftershock sequence.Our results indicate that the earthquake rupture is mainly distributed in the upper-crustal fault,and the maximum slip is up to 8.0 m at~15.0 km depth located in the approximate-80 km east of the epicenter.The average slip is more than 5 m,and the total modelled magnitude is M_(W)7.84,consistent with the observed seismic moment.Our rupture model for the 2015 Gorkha earthquake suggests that the rupture zone is not only in the upper crustal Main Himalayan Thrust(MHT),but also spreads to the northern segment of the MHT.

Rastall gravity:accretion disk image in the context of radiation fields and visual transformations compared to Reissner-Nordström black hole

In this study,we investigated the astronomical implications of Rastall gravity,particularly its behavior amidst a radiation field compared to Reissner-Nordström(RN)black holes.We found a crucial correlation between the dynamics of the accretion disk and the parameters Q and Nr,which properly reflect the influence of spacetime metrics on the disk’s appearance.Elevated electric charge Q causes contraction in the disk’s orbit due to enhanced gravitational effects,while higher Nr values lead to outward expansion,influenced by the attributes of the radiation field.Interestingly,the charged black holes surrounded by radiation fields exhibit distinct visual disparities from RN black holes.Brightness decreases and expansion occurs within the innermost stable circular orbit of the accretion disk with rising Nr values.Our study also reveals the process by which the accretion disk transitions from a conventional disk-like structure to a hat-like form at different observation angles,with the redshift effect gradually intensifying.Moreover,the results of the considered Rastall gravity radiation field are consistent with the constraints of the gravitational lensing of the host galaxy on Rastall gravity parameters,thereby enhancing the consistency between theoretical predictions and actual observations.

Age of the Marwar Supergroup,NW India:A note on the U–Pb geochronology of Jodhpur Group felsic volcanics

The Marwar Supergroup(NW Peninsular India)is thought to be of Ediacaran-Cambrian age,based on previous paleontological and geochronological studies.However,direct constraints on the onset of sedimentation within the Marwar basin are still scarce.In this study,we report U–Pb zircon,LA-ICP-MS,and SIMS ages from the Chhoti Khatu felsic volcanic rocks,interlayered with the Jodhpur Group sandstones(Lower Marwar Supergroup).The cathodoluminescence images of the zircons indicate complex morphologies,and core-rim textures coupled with the wide range of ages indicate that they are likely inherited or in the case of thin poorly indurated ash-beds,detrital in origin.The age spectra of 68 zircon analyses from our sampling display a dominant 800–900 Ma age peak corresponding to the age of basement"Erinpura granite"rocks in the region.The youngest inherited zircon from a felsic ash layer yielded a U–Pb age of651 Ma±18 Ma that,together with previous studies and paleontological evidence,indicates a postCryogenian age for the initiation of Marwar sedimentation following a~125 Ma hiatus between the end of Malani magmatism and Marwar deposition.

Coseismic and postseismic slip ruptures for 2015Mw 6.4 Pishan earthquake constrained by static GPS solutions

On 3 July 2015, a Mw 6.4 earthquake occurred on a blind fault struck Pishan, Xinjiang,China. By combining Crustal Movement Observation Network of China（CMONOC） and other Static Global Positioning System（GPS） sites surrounding Pishan region, it provides a rare chance for us to constrain the slip rupture for such a moderate event. The maximum displacement is up to 12 cm, 2 cm for coseismic and postseismic deformation, respectively,and both the deformation patterns show a same direction moving northeastward. With rectangular dislocation model, a magnitude of Mw6.48, Mw6.3 is calculated based on coseismic, postseismic deformation respectively. Our result indicates the western Kunlun range is still moving toward Tarim Basin followed by an obvious postseismic slip associated with this earthquake. To determine a more reasonable model for postseismic deformation, a longer GPS dataset will be needed.

2.0 Ga orogenic graphite deposits and associated 13C-enriched meta-carbonate rocks from South China Craton:Implications for global Lomagundi event

The Lomagundi(-Jatuli)event,characterized by extremely high positive global inorganic carbon isotope excursion at about 2.2 billion years ago,is pivotal in investigating the causes and consequences of great oxygenation event,inventory and sequestration of carbon on the Earth’s surface,evolution of life,and more profoundly tectonic control on Earth’s environment.However,the reasons that caused the isotopic excursion are not resolved yet.Herein,we report the discovery of meta-carbonate rocks with distinct positive carbon isotopic excursion from the Paleoproterozoic continental collision zone of the Kongling Complex,South China Craton.The δ^(13)C_(V-PDB) values for meta-carbonate rocks show positive values in the range from+5.5‰to+11.6‰,whereas the δ^(13)C_(V-PDB) values of associated graphite deposits range from-25.8‰to-9.5‰.Zircon U-Pb-Hf isotopes from zircon-bearing meta-carbonate sample yielded weighted average _(207)Pb/_(206)Pb age of 2001.3±9.5 Ma,with correspondingε_(Hf)(t)range from-7.05 to-3.16,comparable to the values of local 2.9–2.6 Ga basement rocks.Geochemical characteristics of meta-carbonate rocks,such as their rare earth element patterns and the trace element parameters of La,Ce,Eu,and Gd anomalies and Y/Ho ratio,suggest that the carbonate deposition took place in passive continental margin in association with large volumes of organic carbon.The extensive graphite deposits from Kongling Complex in South China Craton,their equivalents in the North China Craton and elsewhere across the globe prove that the burial of ^(12)C-enriched organic carbon has eventually resulted in the global enrichment of ^(13)C in the atmospheric CO_(2),which is recorded in the marine carbonate rocks.Isotopic mass balance estimates indicate that more than half of the organic carbon was buried during the oceanic closure.Hence,the observed global shift could be directly related to the continent collision event in greater China,thus resolving the long-standing paradox of the Lomagundi global positive carbon isotope excursion.Moreover,the present results suggest that orogenesis play a significant role in sequestration of carbon into the continental crust.

Seismic wave modeling in viscoelastic VTI media using spectral element method

Spectral element method（SEM） for elastic media is well known for its great flexibility and high accuracy in solving problems with complex geometries.It is an advanced choice for wave simulations.Due to anelasticity of earth media,SEM for elastic media is no longer appropriate.On fundamental of the second-order elastic SEM,this work takes the viscoelastic wave equations and the vertical transversely isotropic（VTI） media into consideration,and establishes the second-order SEM for wave modeling in viscoelastic VTI media.The second-order perfectly matched layer for viscoelastic VTI media is also introduced.The problem of handling the overlapped absorbed corners is solved.A comparison with the analytical solution in a twodimensional viscoelastic homogeneous medium shows that the method is accurate in the wave-field modeling.Furtherly,numerical validation also presents its great flexibility in solving wave propagation problems in complex heterogeneous media.This second-order SEM with perfectly matched layer for viscoelastic VTI media can be easily applied in wave modeling in a limited region.

Removal of the airwave effect by main-part decomposition of the anomalous field of MCSEM data

The airwave effect greatly influences the observational data from controlledsource electromagnetic exploration in shallow seas, which obscures the abnormal effects generated by exploration targets and, hence, affects the accuracy of the late exploration data interpretation. In this study, we propose a method to separate the main part from the anomalous field of marine controlled-source electromagnetic method （MCSEM） data based on Stratton-Chu integral transforms to eliminate the airwave effect, which dominates observed electromagnetic （EM） response in shallow seawater. This method of separating the main part from the anomalous field is a type of finite impulse response filter based on a discrete data set. Theoretical analysis proved that the method is stable and able to effectively depress noise. A numerical test indicated that the method could successfully eliminate the airwave effect from the observed EM signals generated by an air water interface and a seawater layer. This technique is applicable for seawater models with either flat or rough seabeds.

Near-field triggering of microearthquakes along the Longitudinal Valley fault in eastern Taiwan

Dynamic earthquake triggering as a well-documented phenomenon can provide valuable information for studying the stress loading cycle from failure on faults.In this study,seismicity rate changes were investigated in the Longitudinal Valley fault(LVF)following the 2019 M_(L)5.2 Hualien earthquake,which occurred offshore in eastern Taiwan.After the matched filter technique was applied to continuous waveform data,twice as many microearthquakes were newly detected in the vicinity of the LVF compared with the number listed in the Taiwan Weather Bureau catalog.Seismicity rates in the northern segment of the LVF increased immediately following the Hualien mainshock,which indicated dynamic triggering during the passage of seismic waves.Statistical analysis suggested that following seismic events might be attributed to fault slipping or creeping.These findings show that the observation of earthquake triggering can provide valuable assistance in monitoring the stress perturbations of active faults.

A new backtracking-based sparsity adaptive algorithm for distributed compressed sensing

A new iterative greedy algorithm based on the backtracking technique was proposed for distributed compressed sensing(DCS) problem. The algorithm applies two mechanisms for precise recovery soft thresholding and cutting. It can reconstruct several compressed signals simultaneously even without any prior information of the sparsity, which makes it a potential candidate for many practical applications, but the numbers of non-zero(significant) coefficients of signals are not available. Numerical experiments are conducted to demonstrate the validity and high performance of the proposed algorithm, as compared to other existing strong DCS algorithms.

Magnetic Signature of Serpentinization at Zedang in the South Tibetan Ophiolite Belt

Magnetic signature of serpentinized mantle peridotite has crucial importance in understanding the serpentinization process and interpreting the origin of strong magnetization anomalies at ultramafic-hosted hydrothermal settings. However, different groups of serpentinized peridotites from both ocean drillings and ophiolite complexes have shown considerable variations in the abundance of magnetite(Oufi et al., 2002;Bonnemains et al., 2016;Li et al., 2017). We examined the magnetic properties, petrography and mineral chemistry of variably serpentinized peridotites from Zedang ophiolite in the eastern Yarlung-Zangbo suture in south Tibet to evaluate the conditions of serpentinization and magnetite formation as well as magnetic sources in suture zones. The studied samples were 0–90% serpentinized with densities from 3.316 to 2.593 g cm–3 and show typical mesh textures of olivine replaced by serpentine on thin sections of core specimen. Serpentines were divided into type-1 Fe-poor serpentine mesh(1.84–2.88 wt% FeO) associated with magnetite in the early stage and type-2 Fe-rich serpentine cores(3.92–5.12 wt% FeO) with no formation of magnetite in the late serpentinization. Brucite vein appeared in central serpentine veins and show Mg/(Mg+Fe) values of 0.74–0.87 at ~50–70% of serpentinization. Pure magnetite was identified as the main magnetic carrier by thermomagnetic analyses, but minor Cr-magnetite(~0.8 mole fractions of Fe3O4) was also detected due to oxidation of early spinel. All the peridotite samples show a rapid increase of magnetic susceptibility from ~0.001 to ~0.03 SI before 40–50% of serpentinization and a following flat trend in values 0.02–0.03 SI at > 50% of serpentinization. This density-susceptibility relationship differs from the rapid production of magnetite above 60-70% of serpentinization for many abyssal peridotites(Oufi et al., 2002;Bach et al., 2006) and suggests that magnetite formation was coupled with hydration of olivine in the early serpentinization but the two decoupled at ~ 40–50% of serpentinization. This transition is consistent with the petrographic observation that magnetite-free serpentinization was developed in higher degrees(> 50%) of serpentinization. Prior studies suggested that serpentinization of < 200℃ would generate Fe-rich brucite, serpentine and little magnetite, whereas magnetite-rich serpentinization was associated with Fe-poor brucite and occurred at higher temperatures of 200–300℃(Klein et al., 2014). The petromagnetic features of serpentinized peridotites from the Zedang ophiolite indicate that the serpentinization process took place initially above 250℃(estimate from brucite composition) and continued to lower temperatures of < 200℃, probably during the mantle lithosphere cooling down in forearc settings(Xiong et al., 2017). These serpentinized peridotites have higher magnetization intensities(average 2.26 Am-1) than mafic dolerite dykes and basaltic volcanic rocks(mostly < 1 Am-1) and should be significant sources of aeromagnetic highs in the Yarlung-Zangbo suture.

Theoretical dispersion curves for borehole real-valued wave modes in vertically transverse isotropic formations

The dispersion curves of real-valued modes in a fluid-filled borehole are widely used in acoustic well logging.The accurate dispersion curves are the precondition of theoretical analysis and inversion process.Generally,these curves can be obtained by solving the conventional dispersion equation for isotropic formations and most vertically transverse isotropy(VTI)formations.However,if the real-valued solutions exist when the radial wavenumbers for the formation quasi-P and quasi-S equals to each other,the existed methods based on the conventional dispersion equation could lead to incorrect results for some VTI formations.Few studies have focused on the influence of these real-valued solutions on dispersion curve extraction.To remove these real-valued solutions,we have proposed a modified dispersion equation and its corresponding solving process.When solving the dispersion equation,the Scholte wave velocity of VTI formation at high frequency is used as the initial guess.The two synthetic examples including fast and slow VTI formations validate that these real-valued solutions do not contribute to the wavefield,and the new dispersion curve extraction method is suitable for all kinds of VTI formations.Consequently,the method can provide reliable dispersion curves for both theoretical analysis and anisotropic parameters inversion in VTI formations.

Fast forward modeling of gravity anomalies for two-dimensional bodies of arbitrary shape and density distribution

A fast and high precision spatial domain algorithm is presented for forward modeling of two-dimensional(2D)body gravity anomalies of arbitrary shape and density distribution.The new algorithm takes advantage of the convolution properties of the expression for 2D gravity anomalies,uses a rectangular cell as a grid subdivision unit,and then 2D bodies with irregular cross-sections are approximated by a combination of 2D bodies with a rectangular cross section.The closed-form expression is used to calculate the gravitational anomalies of the combination of 2D bodies with a rectangular cross section.To improve computing effi ciency,the new algorithm uses a fast algorithm for the implementation of the Toeplitz matrix and vector multiplication.The synthetic 2D models with rectangular and circular cross-sections and constant and variable densities are designed to evaluate the computational accuracy and speed of the new algorithm.The experiment results show that the computation costs less than 6 s for a grid subdivision with 10000×10000 elements.Compared to the traditional forward modeling methods,the proposed method significantly improved computational effi ciency while guaranteeing computational accuracy.

The Present Status and Prospect of Earthquake Focal Depth Locating

Locating an earthquakes focal depth is always a key project in seismology. Precise focal depth is of critical importance for evaluating seismic hazards, deciphering dynamic mechanisms of earthquake generating,estimating aftershock evolutions and risk,as well as monitoring nuclear tests. However,how we determine an accurate focal depth is always a challenge in seismological studies. Aiming to solve these problems, we analyzed and summarized the present status and the future development of earthquake focal depth locating. In this paper we first reviewed the present status of focal depth locating in the world,and summarized the frequently-used relocating methods and ideas at present,and introduced two types of focal depth relocating ideas: arrival time relocating and waveform modeling methods. For these ideas,we systematically described the S-P and the Pn-Pg methods that belong to arrival time method,and polarization focal depth locating and amplitude focal depth locating that belongs to waveform modeling,and further analyzed the advantages and limitations of these methods. Since the depth phase methods are highly sensitive to focal depth,and are relatively free from the uncertainties of crustal models,we mainly reviewed the depth phases of s Pm P,s PL,s Pn,and s Sn,and quantitatively evaluated their availabilities and characteristics. Second,we also discussed the effects of crustal velocity models on the reliability of focal depth locating,and reviewed the advancements of seismic tomography techniques over recent years. Finally,based on the present status of the progress on the focal depth locating,and studies of seismic velocity structures,we proposed an idea of combining multiple datasets and relocating methods,jointly utilizing seismologic and geodetic techniques to relocate focal depth,which should be the major research field in investigating focal depth and source parameters in the near future.

Crust and Upper Mantle Density Structures beneath the Eastern Tianshan Region and Its Tectonic Implications

The deformation mechanisms of the Tianshan orogenic belt(TOB)are one of the most important unresolved issues in the collision of the Indian and Eurasian plates.To better understand the lithospheric deformation of the eastern Tianshan orogenic belt,we combined the S-wave tomography and gravity data to develop a three-dimensional(3D)density model of the crust and upper mantle beneath the eastern Tianshan area.Results show that the crust of the eastern Tianshan is mainly characterized by positive density anomalies,revealing widespread subduction-related magmatism during the Paleozoic.We however have also observed extensive low-density anomalies beneath the eastern Tianshan at depths deeper than~100 km,which is likely linked to a relatively hot mantle.The most fundamental differences of the lithosphere within the eastern Tianshan occur in the uppermost mantle.The uppermost mantle layers in the Bogda Shan and Harlik Shan are relatively dense.This is likely associated with an eclogite body in the uppermost mantle.The most significant negative anomaly of the uppermost mantle is however found in the Jueluotage tectonic belt and the central Tianshan Block and is possibly associated with depleted mantle material.We suggest that these differences related to compositional changes may control the strength of the lithospheric mantle and have affected the uplift of the northern and southern segments of the eastern Tianshan after the Permian.

Probing the interseismic locking state of the Xianshuihe fault based on a viscoelastic deformation model

The interseismic locking state of tectonic faults is essential for regional seismic hazard assessments.However,it is challenging to obtain this parameter reliably due to the weak deformation and complex model configurations.To better probe the fault locking state,more reliable physical models and well-covered observations are required.Here we investigate the locking state of the Xianshuihe fault based on a new-developed viscoelastic deformation model.Meanwhile,we combine GPS velocities from 13 new near-field stations and existing stations in this region to improve the spatial resolution.Similar to the theoretical predictions,our results indicate that the elastic model will clearly overestimate the fault locking depth and seismic moment accumulation rate,and the fault slip rate inferred from the elastic model is slightly lower than that from the viscoelastic model.Relying on the locking distribution inferred from the viscoelastic model,we identify four potential asperities on the Xianshuihe fault.More importantly,we find a clear spatial correlation between the fault locking distribution and the rupture extent of historical earthquakes,which indicates that the fault locking state may control the rupture extent and thus the magnitude of earthquakes.In addition,our results show that the 2022 M6.8 Luding earthquake only ruptured the south part of a potential asperity,and the accumulated energy in the northern unruptured area is equivalent to an Mw6.9 earthquake,where the seismic hazard deserves special attention.

Seismogenic Structure of the 1605 Qiongshan M7½ Earthquake and Its Holocene Activity History in Northern Hainan Island, China: Evidence from Cross-Section Drilling and Shallow Seismic Profile

The 1605 M7½ Earthquake is the only earthquake in the history of China that has caused large-scale land subsidence into the sea,with the total area of land subsidence exceeding 100 km2.The disaster has led to the sinking of 72 villages.There is still no clear understanding of the source seismogenic fault of this earthquake.In this work,we conducted a detailed study of the middle segment of the Maniao-Puqian fault(MPF),which is the epicenter area,through geomorphological survey,data collection,shallow seismic exploration,cross-section drilling,and chronological dating.The results showed that the middle segment of the MPF zone is composed of three nearly parallel normal faults with a dextral strike-slip:“Macun-Luodou fault(F2-1),Haixiu-Dongyuan fault(F2-2),and ChangliuZhuxihe fault(F2-3)”.And F2-2 is composed of two secondary faults,namely F2-2′and F2-2″,with a flower-shaped structure buried under the ground.It is distributed nearly east-west,dipping to the north and has experienced at least five stages of activities since the Miocene.The vertical activity rates of F2-2′and F2-2″are~2.32 and~2.5 mm/a,since the Holocene,respectively.There were eight cycles of transgression and regression since the Miocene.The fault activity resulted in the thickening of the Holocene strata with a slight dip to the south,on the hanging wall,showing V-shaped characteristics.The MPF is likely the source seismogenic fault of the M7½ earthquake that hit Qiongshan in 1605.

A novel CGBoost deep learning algorithm for coseismic landslide susceptibility prediction

The accurate prediction of landslide susceptibility shortly after a violent earthquake is quite vital to the emergency rescue in the 72-h‘‘golden window”.However,the limited quantity of interpreted landslides shortly after a massive earthquake makes landslide susceptibility prediction become a challenge.To address this gap,this work suggests an integrated method of Crossing Graph attention network and xgBoost(CGBoost).This method contains three branches,which extract the interrelations among pixels within a slope unit,the interrelations among various slope units,and the relevance between influencing factors and landslide probability,respectively,and obtain rich and discriminative features by an adaptive fusion mechanism.Thus,the difficulty of susceptibility modeling under a small number of coseismic landslides can be reduced.As a basic module of CGBoost,the proposed Crossing graph attention network(Crossgat)could characterize the spatial heterogeneity within and among slope units to reduce the false alarm in the susceptibility results.Moreover,the rainfall dynamic factors are utilized as prediction indices to improve the susceptibility performance,and the prediction index set is established by terrain,geology,human activity,environment,meteorology,and earthquake factors.CGBoost is applied to predict landslide susceptibility in the Gorkha meizoseismal area.3.43%of coseismic landslides are randomly selected,of which 70%are used for training,and the others for testing.In the testing set,the values of Overall Accuracy,Precision,Recall,F1-score,and Kappa coefficient of CGBoost attain 0.9800,0.9577,0.9999,0.9784,and 0.9598,respectively.Validated by all the coseismic landslides,CGBoost outperforms the current major landslide susceptibility assessment methods.The suggested CGBoost can be also applied to landslide susceptibility prediction in new earthquakes in the future.