The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was devel...The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was developed to simulate the IBC. The simulation of galvanic coupling IBC and electrostatic coupling IBC were implemented along with different signal transmission paths, and their attenuations were calculated. Our study showed that the position near the signal electrode had higher potential than other positions in the two types of IBC, while the potential generally decreased along the axis of the body parts. Both signal attenuations of the two types IBC increased with increasing signal transmission distance, and the electrostatic coupling IBC had comparatively higher receiving potential than the galvanic coupling IBC. The results indicated that the proposed modeling method could be used for the research of biomedical monitoring based on IBC technology.展开更多
Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and ...Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and with the nucleus removed were developed and validatedusing published experimental and clinical data. Then the models with a stainless steel nucleusprosthesis implanted and with polymer nucleus prostheses of various properties implanted were usedfor the 3D finite-element biomechanical analysis. All the above simulation and analysis were carriedout for the L4/L5 disc under a human worst--daily compression load of 2000 N. The results show thatthe polymer materials with Young's modulus of elasticity E = 0.1-100 MPa and Poisson's ratio v=0.35-0.5 are suitable to produce artificial nucleus prosthesis in view of biomechanicalconsideration.展开更多
The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor...Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.展开更多
The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such...The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.展开更多
Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficien...Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.展开更多
In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and cha...In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and change.Likewise,this period has seen a significant increase in our understanding of the physical processes and mechanisms that drive precipitation and its variability across different regions of Africa.By leveraging a large volume of climate model outputs,numerous studies have investigated the model representation of African precipitation as well as underlying physical processes.These studies have assessed whether the physical processes are well depicted and whether the models are fit for informing mitigation and adaptation strategies.This paper provides a review of the progress in precipitation simulation overAfrica in state-of-the-science climate models and discusses the major issues and challenges that remain.展开更多
(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbi...(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbitrary-elevation one-cylinder model.The derived results include a closed-form expression for the space-time correlation function and some quasi-closed-form ones for the space-Doppler power spectrum density,the level crossing rate,and the average fading duration,which are shown to be the generalizations of those previously obtained from the two-dimensional(2-D)one-ring model and the 3-D low-elevation one-cylinder model for terrestrial mobile-to-mobile channels.The close agreements between the theoretical results and the simulations as well as the measurements validate the utility of the derived channel statistics.Based on the derived expressions,the impacts of some parameters on the channel characteristics are investigated in an effective,efficient,and explicable way,which leads to a general guideline on the manual parameter estimation from the measurement description.展开更多
In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical propert...In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical properties of rocks,the cracking processes of pre-cracked rocks have been extensively studied using numerical modeling methods.The peridynamics(PD)exhibits advantages over other numerical methods due to the absence of the requirements for remeshing and external crack growth criterion.However,for modeling pre-cracked rock cracking processes under impact,current PD implementations lack generally applicable rock constitutive models and impact contact models,which leads to difficulties in determining rock material parameters and efficiently calculating impact loads.This paper proposes a non-ordinary state-based peridynamics(NOSBPD)modeling method integrating the Drucker-Prager(DP)plasticity model and an efficient contact model to address the above problems.In the proposed method,the Drucker-Prager plasticity model is integrated into the NOSBPD,thereby equipping NOSBPD with the capability to accurately characterize the nonlinear stress-strain relationship inherent in rocks.An efficient contact model between particles and meshes is designed to calculate the impact loads,which is essentially a coupling method of PD with the finite element method(FEM).The effectiveness of the proposed NOSBPD modeling method is verified by comparison with other numerical methods and experiments.Experimental results indicate that the proposed method can effectively and accurately predict the 3D cracking processes of pre-cracked cracks under impact loading,and the maximum principal stress is the key driver behind wing crack formation in pre-cracked rocks.展开更多
Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally...Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally expensive and typically requires exhaustive search.We introduce a thermodynamics-inspired framework that treats activation distributions as energy-filtered physical systems and employs the free energy of activations as a principled evaluation metric.Phase-transition-like phenomena in the free-energy profile—such as extrema,inflection points,and curvature changes—yield reliable estimates of the critical pruning threshold,providing a theoretically grounded means of predicting sharp accuracy degradation.To further enhance efficiency,we propose a renormalized free energy technique that approximates full-evaluation free energy using only the activation distribution of the unpruned network.This eliminates repeated forward passes,dramatically reducing computational overhead and achieving speedups of up to 550×for MLPs.Extensive experiments across diverse vision architectures(MLP,CNN,ResNet,MobileNet,Vision Transformer)and text models(LSTM,BERT,ELECTRA,T5,GPT-2)on multiple datasets validate the generality,robustness,and computational efficiency of our approach.Overall,this work establishes a theoretically grounded and practically effective framework for activation pruning,bridging the gap between analytical understanding and efficient deployment of sparse neural networks.展开更多
Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evo...Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evolution,and image synthesis to compare directly with HST,LICIACube,ground-based and Lucy observations of the DART impact.Decomposing ejecta into(1)a highvelocity(~1600 m/s)plume exhibiting Na/K resonance,(2)a low-velocity(~1 m/s)conical component shaped by binary gravity and solar radiation pressure,and(3)meter-scale boulders,we quantify each component’s mass and momentum.Fitting photometric decay curves and morphological evolution yields size-velocity distributions and,via scaling laws,estimates of Dimorphos’bulk density,cratering parameters,and cohesive strength that agree with dynamical constraints.Photometric ejecta modeling therefore provides a robust route to constrain momentum enhancement and target properties,improving predictive capability for kinetic-deflection missions.展开更多
Industrial robot dynamics lay the foundation for high-precision and high-speed control, and accurate identification of dynamic parameters is essential for precise dynamic calculations. The choice of friction models is...Industrial robot dynamics lay the foundation for high-precision and high-speed control, and accurate identification of dynamic parameters is essential for precise dynamic calculations. The choice of friction models is a critical component in the identification of industrial robot dynamics. Traditional static friction models struggle to capture the hysteresis effects caused by robot joint elasticity and clearances, leading to large torque prediction errors when the joint velocity crosses zero. Due to the presence of hysteresis effects, the joint velocity crosses zero in the forward direction, and the reverse direction will have different friction patterns. Although the hysteresis effects can be modeled as an ordinary differential equation(ODE), it is difficult to determine the ODE structure that achieves both generalization and accuracy to describe the hysteresis effects of the friction model. To address this issue, we propose the neural hysteresis friction(NHF), which uses neural ODE to model the hysteresis effects in a data-driven manner, thereby mitigating the current inadequacies in the study of dynamic friction characteristics. The experiments on a real 6-axis industrial robot demonstrate that our proposed method can accurately model the friction dynamics during directional switching and outperform other modeling methods. Velocity tracking control experiments show that NHF can effectively reduce tracking errors when the velocity crosses zero.展开更多
Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately...Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately 80 km along a NW-SE axis.With a summit reaching~3214 m(a.s.l.),Mount Kandil is a stratovolcano that,like many other peaks within the EAHP and the Lesser Caucasus,experienced significant environmental changes during Late Pleistocene.Among these,glacial processes stand out as the most profound,having distinctly shaped the mountains geomorphic landscape.This study presents,for the first time,a comprehensive analysis of the glacial morphology of Mount Kandil based on multiple datasets.Field-based morphological observations indicate that an area of approximately 32.62 km^(2)has been sculpted by glacial activity.Within six glaciated regions on Mount Kandil,25 cirques and 6 glacial valleys have been identified.In addition,moraines in various locations exhibit characteristic morphologies.Furthermore,valley glaciers are inferred to have descended to altitudes as low as~2000 m.The paleoequilibrium line(p ELA)was estimated to use AABR method within GIS,yielding a mean pELA of~2730 m.Ice thickness modelling indicates that the thickness of glaciers in the Kandil Mountain valleys reaches up to~350 m.Due to its orographic extension,Mount Kandil is exposed to humid northwest winds and receives substantial frontal precipitation(~686 mm annually).The compiled geomorphic,cartographic and morphometric parameters suggest that the glaciation dynamics of Mount Kandil—situated between the Southeastern Taurus and the Lesser Caucasus—closely resemble those observed in the Lesser Caucasus.This indicates that glaciation was primarily governed by northern atmospheric systems with additional influences from southerly or westerly winds.The integrated data also underscores the role of multiple atmospheric systems in controlling the glaciation regime around the Lesser Caucasus.Additionally,findings on regional pELA question the common belief that pELA increases eastward in EAHP.展开更多
In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicti...In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicting long-term roadway stability,necessitating the development of a reliable constitutive creep model and numerical simulation approach.In this study,creep experiments were conducted on pre-damaged rock with varying initial damage levels to investigate the time-dependent mechanical properties.Based on the experimental results,an accelerated-creep criterion was proposed,and an elastic-viscoplastic creep damage model(EVPCD)was established that simultaneously considers the effects of time-dependent damage and instantaneous damage caused by stress disturbances on rock creep behavior.Subsequently,the effectiveness of the proposed creep model was verified using experimental data,and the secondary development of the EVPCD model was completed based on the FLAC3D platform.Following this,a long-term stability analysis method of deep surrounding rock that accounts for excavation-and mining-induced disturbances was proposed.Using the main roadway of Xutuan Coal Mine as a case study,numerical simulations were carried out to investigate the time-dependent deformation and failure characteristics of the surrounding rock following excavation and mining disturbance.Combined with on-site monitoring of the surrounding rock damage areas,the results indicate that the EVPCD outperforms the CVISC and Nishihara models in predicting the time-dependent behavior of deep surrounding rock.展开更多
Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and ...Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and finite element analyses.Three tests of the slope deformation under earthquake and post-rainfall earthquakes are first studied using image analysis techniques.Then,based on an elastoplastic constitutive model,numerical simulations are carried out using the finite element method and compared with the centrifuge test results.Finally,a parametric study is performed to clarify the effects of antecedent rainfall on earthquake-induced slope deformation.The results show that slope deformation caused by post-rainfall earthquakes differs from that caused by earthquakes without antecedent rainfall.The seepage flow and soil strength of the slope are affected by previous rainfall conditions,such as intensity and duration,which directly influence the slope deformation caused by the subsequent earthquake.Soil displacement and strain become greater and the slip surface is more noticeable during the post-rainfall earthquake of higher intensity.In addition,the time interval between the rainfall and the earthquake has a considerable impact on the detailed characteristics of the slope deformation,and the significant deformation occurs at the time of lowest soil strength when seepage flow reaches the lower part of the slope.Moreover,the repeated intermittent rainfall greatly affects the subsequent earthquake-induced slope deformation,the main characteristics of which are closely related to the changes in saturation and strength of the slope.However,with the prolonged time gap between each round of rainfall,the earthquake-induced slope deformation becomes insignificant.展开更多
Hepatitis B Virus(HBV)infection and heavy alcohol consumption are the two primary pathogenic causes of liver cirrhosis.In this paper,we proposed a deterministic mathematical model and a logistic equation to investigat...Hepatitis B Virus(HBV)infection and heavy alcohol consumption are the two primary pathogenic causes of liver cirrhosis.In this paper,we proposed a deterministic mathematical model and a logistic equation to investigate the dynamics of liver cirrhosis progression as well as to explain the implications of variations in alcohol consumption on chronic hepatitis B patients,respectively.The intricate interactions between liver cirrhosis,recovery,and treatment dynamics are captured by the model.This study aims to show that alcohol consumption by Hepatitis B-infected individuals accelerates liver cirrhosis progression while treatment of acutely infected individuals reduces it.We proved that a unique solution of the proposed model exists,which is positive and bounded.Using the next-generation matrix approach,two basic reproductive numbers R_(A_(0))and R_(A_(max))are calculated to identify future recurrence.The equilibrium points are calculated,and both equilibria are proved locally and globally asymptotically stable when R_(0)is below and above one,respectively.It is shown that bifurcation exists at R_(0)=1 and a detailed proof for forward bifurcation is given.Furthermore,we performed the sensitivity analysis of the model parameters on R_(0).For the confirmation of analytical work,we performed numerical simulations,and the results indicate that the treatment and the inhibitory effects reduce the risk of developing liver cirrhosis in individuals,while heavy alcohol consumption accelerates markedly the liver cirrhosis progression in patients with chronic hepatitis B.展开更多
Accurate mechanical modeling is essential for robotic belt grinding(RBG), a process characterized by compliant contact mechanisms that make force prediction particularly challenging. However, existing mechanical model...Accurate mechanical modeling is essential for robotic belt grinding(RBG), a process characterized by compliant contact mechanisms that make force prediction particularly challenging. However, existing mechanical models predominantly focus on macroscale compliance while neglecting grain-scale compliant motion. Moreover, abrasive grains are typically idealized as regular shapes, overlooking the inherent stochasticity of real grain geometries. This study proposes a shapeequivalence method for modeling stochastic abrasive grains and develops a multiscale compliant force model for RBG. Specifically, an individual grain is represented as a polygonal pyramid with stochastic edges that is mathematically equivalent to a cone;this method unifies the treatment of grain geometries and streamlines the modeling process. The mathematical equivalence relationship for random grain shapes is further derived based on a grain-compliant contact model. By integrating grain geometric characteristics and progressive grain wear, an analytical mechanical model that captures both the static contact force and dynamic grinding force is established, thereby describing the transition from grain-workpiece compliant interaction to belt-workpiece elastic contact. Grinding experiments were conducted using abrasive belts with different grain shape distributions to validate the model. The results demonstrated reliable predictions of the tangential grinding force and its component characteristics. Additional analyses were performed to reveal how the tangential grinding force varies with wear time and grinding parameters.展开更多
Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical si...Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical simulation plays a significant role in quantitatively evaluating current processes and making targeted improvements,but its limitations lie in the inability to dynamically reflect the formation outcomes of castings under varying process conditions,making real-time adjustments to gating and riser designs challenging.In this study,an automated design model for gating and riser systems based on integrated parametric 3D modeling-simulation framework is proposed,which enhances the flexibility and usability of evaluating the casting process by simulation.Firstly,geometric feature extraction technology is employed to obtain the geometric information of the target casting.Based on this information,an automated design framework for gating and riser systems is established,incorporating multiple structural parameters for real-time process control.Subsequently,the simulation results for various structural parameters are analyzed,and the influence of these parameters on casting formation is thoroughly investigated.Finally,the optimal design scheme is generated and validated through experimental verification.Simulation analysis and experimental results show that using a larger gate neck(24 mm in side length) and external risers promotes a more uniform temperature distribution and a more stable flow state,effectively eliminating shrinkage cavities and enhancing process yield by 15%.展开更多
To address the issues of frequent identity switches(IDs)and degraded identification accuracy in multi object tracking(MOT)under complex occlusion scenarios,this study proposes an occlusion-robust tracking framework ba...To address the issues of frequent identity switches(IDs)and degraded identification accuracy in multi object tracking(MOT)under complex occlusion scenarios,this study proposes an occlusion-robust tracking framework based on face-pedestrian joint feature modeling.By constructing a joint tracking model centered on“intra-class independent tracking+cross-category dynamic binding”,designing a multi-modal matching metric with spatio-temporal and appearance constraints,and innovatively introducing a cross-category feature mutual verification mechanism and a dual matching strategy,this work effectively resolves performance degradation in traditional single-category tracking methods caused by short-term occlusion,cross-camera tracking,and crowded environments.Experiments on the Chokepoint_Face_Pedestrian_Track test set demonstrate that in complex scenes,the proposed method improves Face-Pedestrian Matching F1 area under the curve(F1 AUC)by approximately 4 to 43 percentage points compared to several traditional methods.The joint tracking model achieves overall performance metrics of IDF1:85.1825%and MOTA:86.5956%,representing improvements of 0.91 and 0.06 percentage points,respectively,over the baseline model.Ablation studies confirm the effectiveness of key modules such as the Intersection over Area(IoA)/Intersection over Union(IoU)joint metric and dynamic threshold adjustment,validating the significant role of the cross-category identity matching mechanism in enhancing tracking stability.Our_model shows a 16.7%frame per second(FPS)drop vs.fairness of detection and re-identification in multiple object tracking(FairMOT),with its cross-category binding module adding aboute 10%overhead,yet maintains near-real-time performance for essential face-pedestrian tracking at small resolutions.展开更多
The intracontinental subduction of a>200-km-long section of the Tajik-Tarim lithosphere beneath the Pamir Mountains is proposed to explain nearly 30 km of shortening in the Tajik fold-thrust belt and the Pamir upli...The intracontinental subduction of a>200-km-long section of the Tajik-Tarim lithosphere beneath the Pamir Mountains is proposed to explain nearly 30 km of shortening in the Tajik fold-thrust belt and the Pamir uplift.Seismic imaging revealed that the upper slab was scraped and that the lower slab had subducted to a depth of>150 km.These features constitute the tectonic complexity of the Pamirs,as well as the thermal subduction mechanism involved,which remains poorly understood.Hence,in this study,high-resolution three-dimensional(3D)kinematic modeling is applied to investigate the thermal structure and geometry of the subducting slab beneath the Pamirs.The modeled slab configuration reveals distinct along-strike variations,with a steeply dipping slab beneath the southern Pamirs,a more gently inclined slab beneath the northern Pamirs,and apparent upper slab termination at shallow depths beneath the Pamirs.The thermal field reveals a cold slab core after delamination,with temperatures ranging from 400℃to 800℃,enveloped by a hotter mantle reaching~1400℃.The occurrence of intermediate-depth earthquakes aligns primarily with colder slab regions,particularly near the slab tear-off below the southwestern Pamirs,indicating a strong correlation between slab temperature and seismicity.In contrast,the northern Pamirs exhibit reduced seismicity at depth,which is likely associated with thermal weakening and delamination.The central Pamirs show a significant thermal anomaly caused by a concave slab,where the coldest crust does not descend deeply,further suggesting crustal detachment or mechanical failure.The lateral asymmetry in slab temperature possibly explains the mechanism of lateral tearing and differential slab-mantle coupling.展开更多
基金Supported by the National Natural Science Foundation of China(60801050)the Excellent Talent Fund of Beijing(2011)Excellent Young Scholars Research Fund of Beijing Institute ofTechnology(2012)
文摘The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was developed to simulate the IBC. The simulation of galvanic coupling IBC and electrostatic coupling IBC were implemented along with different signal transmission paths, and their attenuations were calculated. Our study showed that the position near the signal electrode had higher potential than other positions in the two types of IBC, while the potential generally decreased along the axis of the body parts. Both signal attenuations of the two types IBC increased with increasing signal transmission distance, and the electrostatic coupling IBC had comparatively higher receiving potential than the galvanic coupling IBC. The results indicated that the proposed modeling method could be used for the research of biomedical monitoring based on IBC technology.
文摘Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and with the nucleus removed were developed and validatedusing published experimental and clinical data. Then the models with a stainless steel nucleusprosthesis implanted and with polymer nucleus prostheses of various properties implanted were usedfor the 3D finite-element biomechanical analysis. All the above simulation and analysis were carriedout for the L4/L5 disc under a human worst--daily compression load of 2000 N. The results show thatthe polymer materials with Young's modulus of elasticity E = 0.1-100 MPa and Poisson's ratio v=0.35-0.5 are suitable to produce artificial nucleus prosthesis in view of biomechanicalconsideration.
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
基金supported by the National Natural Science Foundation of China(No.41104068)the Deep Exploration in China,Sino Probe-03-05
文摘Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.
基金financially supported by the National Natural Science Foundation of China(No.41204055,41164003,and 41104074)Opening Project(No.SMIL-2014-06) of Hubei Subsurface Multi-scale Imaging Lab(SMIL),China University of Geosciences(Wuhan)
文摘The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.
基金National High Technology Research and Development Program(863 Program)(No.2006AA06Z105,2007AA06Z134)
文摘Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.
基金the World Climate Research Programme(WCRP),Climate Variability and Predictability(CLIVAR),and Global Energy and Water Exchanges(GEWEX)for facilitating the coordination of African monsoon researchsupport from the Center for Earth System Modeling,Analysis,and Data at the Pennsylvania State Universitythe support of the Office of Science of the U.S.Department of Energy Biological and Environmental Research as part of the Regional&Global Model Analysis(RGMA)program area。
文摘In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and change.Likewise,this period has seen a significant increase in our understanding of the physical processes and mechanisms that drive precipitation and its variability across different regions of Africa.By leveraging a large volume of climate model outputs,numerous studies have investigated the model representation of African precipitation as well as underlying physical processes.These studies have assessed whether the physical processes are well depicted and whether the models are fit for informing mitigation and adaptation strategies.This paper provides a review of the progress in precipitation simulation overAfrica in state-of-the-science climate models and discusses the major issues and challenges that remain.
基金supported in part by the National Key Research and Development Program of China(2021YFB2900501)in part by the Shaanxi Science and Technology Innovation Team(2023-CX-TD-03)+3 种基金in part by the Science and Technology Program of Shaanxi Province(2021GXLH-Z-038)in part by the Natural Science Foundation of Hunan Province(2023JJ40607 and 2023JJ50045)in part by the Scientific Research Foundation of Hunan Provincial Education Department(23B0713 and 24B0603)in part by the National Natural Science Foundation of China(62401371,62101275,and 62372070).
文摘(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbitrary-elevation one-cylinder model.The derived results include a closed-form expression for the space-time correlation function and some quasi-closed-form ones for the space-Doppler power spectrum density,the level crossing rate,and the average fading duration,which are shown to be the generalizations of those previously obtained from the two-dimensional(2-D)one-ring model and the 3-D low-elevation one-cylinder model for terrestrial mobile-to-mobile channels.The close agreements between the theoretical results and the simulations as well as the measurements validate the utility of the derived channel statistics.Based on the derived expressions,the impacts of some parameters on the channel characteristics are investigated in an effective,efficient,and explicable way,which leads to a general guideline on the manual parameter estimation from the measurement description.
基金support from the National Natural Science Foundation of China(Grant Nos.42277161 and 42230709).
文摘In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical properties of rocks,the cracking processes of pre-cracked rocks have been extensively studied using numerical modeling methods.The peridynamics(PD)exhibits advantages over other numerical methods due to the absence of the requirements for remeshing and external crack growth criterion.However,for modeling pre-cracked rock cracking processes under impact,current PD implementations lack generally applicable rock constitutive models and impact contact models,which leads to difficulties in determining rock material parameters and efficiently calculating impact loads.This paper proposes a non-ordinary state-based peridynamics(NOSBPD)modeling method integrating the Drucker-Prager(DP)plasticity model and an efficient contact model to address the above problems.In the proposed method,the Drucker-Prager plasticity model is integrated into the NOSBPD,thereby equipping NOSBPD with the capability to accurately characterize the nonlinear stress-strain relationship inherent in rocks.An efficient contact model between particles and meshes is designed to calculate the impact loads,which is essentially a coupling method of PD with the finite element method(FEM).The effectiveness of the proposed NOSBPD modeling method is verified by comparison with other numerical methods and experiments.Experimental results indicate that the proposed method can effectively and accurately predict the 3D cracking processes of pre-cracked cracks under impact loading,and the maximum principal stress is the key driver behind wing crack formation in pre-cracked rocks.
基金output of a research project implemented as part of the Basic Research Program at HSE University。
文摘Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally expensive and typically requires exhaustive search.We introduce a thermodynamics-inspired framework that treats activation distributions as energy-filtered physical systems and employs the free energy of activations as a principled evaluation metric.Phase-transition-like phenomena in the free-energy profile—such as extrema,inflection points,and curvature changes—yield reliable estimates of the critical pruning threshold,providing a theoretically grounded means of predicting sharp accuracy degradation.To further enhance efficiency,we propose a renormalized free energy technique that approximates full-evaluation free energy using only the activation distribution of the unpruned network.This eliminates repeated forward passes,dramatically reducing computational overhead and achieving speedups of up to 550×for MLPs.Extensive experiments across diverse vision architectures(MLP,CNN,ResNet,MobileNet,Vision Transformer)and text models(LSTM,BERT,ELECTRA,T5,GPT-2)on multiple datasets validate the generality,robustness,and computational efficiency of our approach.Overall,this work establishes a theoretically grounded and practically effective framework for activation pruning,bridging the gap between analytical understanding and efficient deployment of sparse neural networks.
基金supported by the National Natural Science Foundation of China(Grant No.12272018)the National Key Basic Research Project(2022JCJQZD20600).
文摘Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evolution,and image synthesis to compare directly with HST,LICIACube,ground-based and Lucy observations of the DART impact.Decomposing ejecta into(1)a highvelocity(~1600 m/s)plume exhibiting Na/K resonance,(2)a low-velocity(~1 m/s)conical component shaped by binary gravity and solar radiation pressure,and(3)meter-scale boulders,we quantify each component’s mass and momentum.Fitting photometric decay curves and morphological evolution yields size-velocity distributions and,via scaling laws,estimates of Dimorphos’bulk density,cratering parameters,and cohesive strength that agree with dynamical constraints.Photometric ejecta modeling therefore provides a robust route to constrain momentum enhancement and target properties,improving predictive capability for kinetic-deflection missions.
基金supported by the National Natural Science Foundation of China (Grant No.52188102)。
文摘Industrial robot dynamics lay the foundation for high-precision and high-speed control, and accurate identification of dynamic parameters is essential for precise dynamic calculations. The choice of friction models is a critical component in the identification of industrial robot dynamics. Traditional static friction models struggle to capture the hysteresis effects caused by robot joint elasticity and clearances, leading to large torque prediction errors when the joint velocity crosses zero. Due to the presence of hysteresis effects, the joint velocity crosses zero in the forward direction, and the reverse direction will have different friction patterns. Although the hysteresis effects can be modeled as an ordinary differential equation(ODE), it is difficult to determine the ODE structure that achieves both generalization and accuracy to describe the hysteresis effects of the friction model. To address this issue, we propose the neural hysteresis friction(NHF), which uses neural ODE to model the hysteresis effects in a data-driven manner, thereby mitigating the current inadequacies in the study of dynamic friction characteristics. The experiments on a real 6-axis industrial robot demonstrate that our proposed method can accurately model the friction dynamics during directional switching and outperform other modeling methods. Velocity tracking control experiments show that NHF can effectively reduce tracking errors when the velocity crosses zero.
基金supported by Van Yüzüncü Yıl University,Scientific Research Projects Coordination Unit(Project No:SDK-2025-11935)Van Yüzüncü Yıl University,Scientific Research Projects Coordination Unit for supporting the study。
文摘Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately 80 km along a NW-SE axis.With a summit reaching~3214 m(a.s.l.),Mount Kandil is a stratovolcano that,like many other peaks within the EAHP and the Lesser Caucasus,experienced significant environmental changes during Late Pleistocene.Among these,glacial processes stand out as the most profound,having distinctly shaped the mountains geomorphic landscape.This study presents,for the first time,a comprehensive analysis of the glacial morphology of Mount Kandil based on multiple datasets.Field-based morphological observations indicate that an area of approximately 32.62 km^(2)has been sculpted by glacial activity.Within six glaciated regions on Mount Kandil,25 cirques and 6 glacial valleys have been identified.In addition,moraines in various locations exhibit characteristic morphologies.Furthermore,valley glaciers are inferred to have descended to altitudes as low as~2000 m.The paleoequilibrium line(p ELA)was estimated to use AABR method within GIS,yielding a mean pELA of~2730 m.Ice thickness modelling indicates that the thickness of glaciers in the Kandil Mountain valleys reaches up to~350 m.Due to its orographic extension,Mount Kandil is exposed to humid northwest winds and receives substantial frontal precipitation(~686 mm annually).The compiled geomorphic,cartographic and morphometric parameters suggest that the glaciation dynamics of Mount Kandil—situated between the Southeastern Taurus and the Lesser Caucasus—closely resemble those observed in the Lesser Caucasus.This indicates that glaciation was primarily governed by northern atmospheric systems with additional influences from southerly or westerly winds.The integrated data also underscores the role of multiple atmospheric systems in controlling the glaciation regime around the Lesser Caucasus.Additionally,findings on regional pELA question the common belief that pELA increases eastward in EAHP.
基金funded by the National Natural Science Foundation of China(Nos.52004098,U24B2041,and 52274079)the Key Research and Development Program of Henan Province(No.251111320400)+1 种基金the Key Research Project Plan for Higher Education Institutions in Henan Province(Nos.24A570006 and 25A570002)the Scientific and Technological Research Project in Henan Province(No.242102320061).
文摘In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicting long-term roadway stability,necessitating the development of a reliable constitutive creep model and numerical simulation approach.In this study,creep experiments were conducted on pre-damaged rock with varying initial damage levels to investigate the time-dependent mechanical properties.Based on the experimental results,an accelerated-creep criterion was proposed,and an elastic-viscoplastic creep damage model(EVPCD)was established that simultaneously considers the effects of time-dependent damage and instantaneous damage caused by stress disturbances on rock creep behavior.Subsequently,the effectiveness of the proposed creep model was verified using experimental data,and the secondary development of the EVPCD model was completed based on the FLAC3D platform.Following this,a long-term stability analysis method of deep surrounding rock that accounts for excavation-and mining-induced disturbances was proposed.Using the main roadway of Xutuan Coal Mine as a case study,numerical simulations were carried out to investigate the time-dependent deformation and failure characteristics of the surrounding rock following excavation and mining disturbance.Combined with on-site monitoring of the surrounding rock damage areas,the results indicate that the EVPCD outperforms the CVISC and Nishihara models in predicting the time-dependent behavior of deep surrounding rock.
基金supported by the China Postdoctoral Science Foundation(CPSF)(Grant No.2024M762769)the Natural Science Basic Research Program of Shaanxi(Grant No.2024JC-YBQN-0333)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20232230).
文摘Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and finite element analyses.Three tests of the slope deformation under earthquake and post-rainfall earthquakes are first studied using image analysis techniques.Then,based on an elastoplastic constitutive model,numerical simulations are carried out using the finite element method and compared with the centrifuge test results.Finally,a parametric study is performed to clarify the effects of antecedent rainfall on earthquake-induced slope deformation.The results show that slope deformation caused by post-rainfall earthquakes differs from that caused by earthquakes without antecedent rainfall.The seepage flow and soil strength of the slope are affected by previous rainfall conditions,such as intensity and duration,which directly influence the slope deformation caused by the subsequent earthquake.Soil displacement and strain become greater and the slip surface is more noticeable during the post-rainfall earthquake of higher intensity.In addition,the time interval between the rainfall and the earthquake has a considerable impact on the detailed characteristics of the slope deformation,and the significant deformation occurs at the time of lowest soil strength when seepage flow reaches the lower part of the slope.Moreover,the repeated intermittent rainfall greatly affects the subsequent earthquake-induced slope deformation,the main characteristics of which are closely related to the changes in saturation and strength of the slope.However,with the prolonged time gap between each round of rainfall,the earthquake-induced slope deformation becomes insignificant.
文摘Hepatitis B Virus(HBV)infection and heavy alcohol consumption are the two primary pathogenic causes of liver cirrhosis.In this paper,we proposed a deterministic mathematical model and a logistic equation to investigate the dynamics of liver cirrhosis progression as well as to explain the implications of variations in alcohol consumption on chronic hepatitis B patients,respectively.The intricate interactions between liver cirrhosis,recovery,and treatment dynamics are captured by the model.This study aims to show that alcohol consumption by Hepatitis B-infected individuals accelerates liver cirrhosis progression while treatment of acutely infected individuals reduces it.We proved that a unique solution of the proposed model exists,which is positive and bounded.Using the next-generation matrix approach,two basic reproductive numbers R_(A_(0))and R_(A_(max))are calculated to identify future recurrence.The equilibrium points are calculated,and both equilibria are proved locally and globally asymptotically stable when R_(0)is below and above one,respectively.It is shown that bifurcation exists at R_(0)=1 and a detailed proof for forward bifurcation is given.Furthermore,we performed the sensitivity analysis of the model parameters on R_(0).For the confirmation of analytical work,we performed numerical simulations,and the results indicate that the treatment and the inhibitory effects reduce the risk of developing liver cirrhosis in individuals,while heavy alcohol consumption accelerates markedly the liver cirrhosis progression in patients with chronic hepatitis B.
基金supported by the National Natural Science Foundation of China (Grant Nos.52505554,52575571)the Postdoctoral Fellowship Program of CPSF (Grant No.GZB20250348)。
文摘Accurate mechanical modeling is essential for robotic belt grinding(RBG), a process characterized by compliant contact mechanisms that make force prediction particularly challenging. However, existing mechanical models predominantly focus on macroscale compliance while neglecting grain-scale compliant motion. Moreover, abrasive grains are typically idealized as regular shapes, overlooking the inherent stochasticity of real grain geometries. This study proposes a shapeequivalence method for modeling stochastic abrasive grains and develops a multiscale compliant force model for RBG. Specifically, an individual grain is represented as a polygonal pyramid with stochastic edges that is mathematically equivalent to a cone;this method unifies the treatment of grain geometries and streamlines the modeling process. The mathematical equivalence relationship for random grain shapes is further derived based on a grain-compliant contact model. By integrating grain geometric characteristics and progressive grain wear, an analytical mechanical model that captures both the static contact force and dynamic grinding force is established, thereby describing the transition from grain-workpiece compliant interaction to belt-workpiece elastic contact. Grinding experiments were conducted using abrasive belts with different grain shape distributions to validate the model. The results demonstrated reliable predictions of the tangential grinding force and its component characteristics. Additional analyses were performed to reveal how the tangential grinding force varies with wear time and grinding parameters.
基金financially supported by the National Key Research and Development Program of China (2022YFB3706802)。
文摘Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical simulation plays a significant role in quantitatively evaluating current processes and making targeted improvements,but its limitations lie in the inability to dynamically reflect the formation outcomes of castings under varying process conditions,making real-time adjustments to gating and riser designs challenging.In this study,an automated design model for gating and riser systems based on integrated parametric 3D modeling-simulation framework is proposed,which enhances the flexibility and usability of evaluating the casting process by simulation.Firstly,geometric feature extraction technology is employed to obtain the geometric information of the target casting.Based on this information,an automated design framework for gating and riser systems is established,incorporating multiple structural parameters for real-time process control.Subsequently,the simulation results for various structural parameters are analyzed,and the influence of these parameters on casting formation is thoroughly investigated.Finally,the optimal design scheme is generated and validated through experimental verification.Simulation analysis and experimental results show that using a larger gate neck(24 mm in side length) and external risers promotes a more uniform temperature distribution and a more stable flow state,effectively eliminating shrinkage cavities and enhancing process yield by 15%.
基金supported by the confidential research grant No.a8317。
文摘To address the issues of frequent identity switches(IDs)and degraded identification accuracy in multi object tracking(MOT)under complex occlusion scenarios,this study proposes an occlusion-robust tracking framework based on face-pedestrian joint feature modeling.By constructing a joint tracking model centered on“intra-class independent tracking+cross-category dynamic binding”,designing a multi-modal matching metric with spatio-temporal and appearance constraints,and innovatively introducing a cross-category feature mutual verification mechanism and a dual matching strategy,this work effectively resolves performance degradation in traditional single-category tracking methods caused by short-term occlusion,cross-camera tracking,and crowded environments.Experiments on the Chokepoint_Face_Pedestrian_Track test set demonstrate that in complex scenes,the proposed method improves Face-Pedestrian Matching F1 area under the curve(F1 AUC)by approximately 4 to 43 percentage points compared to several traditional methods.The joint tracking model achieves overall performance metrics of IDF1:85.1825%and MOTA:86.5956%,representing improvements of 0.91 and 0.06 percentage points,respectively,over the baseline model.Ablation studies confirm the effectiveness of key modules such as the Intersection over Area(IoA)/Intersection over Union(IoU)joint metric and dynamic threshold adjustment,validating the significant role of the cross-category identity matching mechanism in enhancing tracking stability.Our_model shows a 16.7%frame per second(FPS)drop vs.fairness of detection and re-identification in multiple object tracking(FairMOT),with its cross-category binding module adding aboute 10%overhead,yet maintains near-real-time performance for essential face-pedestrian tracking at small resolutions.
基金the Chinese Academy of Sciences Pioneer Hundred Talents Program and the Second Tibetan Plateau Scientific Expedition and Research Program(Grant No.2019QZKK0708)supported by a MEXT(Ministry of Education,Culture,Sports,Science and Technology)KAKENHI(Grants-in-Aid for Scientific Research)grant(Grant No.21H05203)Kobe University Strategic International Collaborative Research Grant(Type B Fostering Joint Research).
文摘The intracontinental subduction of a>200-km-long section of the Tajik-Tarim lithosphere beneath the Pamir Mountains is proposed to explain nearly 30 km of shortening in the Tajik fold-thrust belt and the Pamir uplift.Seismic imaging revealed that the upper slab was scraped and that the lower slab had subducted to a depth of>150 km.These features constitute the tectonic complexity of the Pamirs,as well as the thermal subduction mechanism involved,which remains poorly understood.Hence,in this study,high-resolution three-dimensional(3D)kinematic modeling is applied to investigate the thermal structure and geometry of the subducting slab beneath the Pamirs.The modeled slab configuration reveals distinct along-strike variations,with a steeply dipping slab beneath the southern Pamirs,a more gently inclined slab beneath the northern Pamirs,and apparent upper slab termination at shallow depths beneath the Pamirs.The thermal field reveals a cold slab core after delamination,with temperatures ranging from 400℃to 800℃,enveloped by a hotter mantle reaching~1400℃.The occurrence of intermediate-depth earthquakes aligns primarily with colder slab regions,particularly near the slab tear-off below the southwestern Pamirs,indicating a strong correlation between slab temperature and seismicity.In contrast,the northern Pamirs exhibit reduced seismicity at depth,which is likely associated with thermal weakening and delamination.The central Pamirs show a significant thermal anomaly caused by a concave slab,where the coldest crust does not descend deeply,further suggesting crustal detachment or mechanical failure.The lateral asymmetry in slab temperature possibly explains the mechanism of lateral tearing and differential slab-mantle coupling.
