Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important...Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with smallscale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference （FD） method of Biot＇s poroelastic equations, with an arbitrary even-order （2L） accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer （CPML） absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultra- sonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thick- nesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Qsc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments.展开更多
During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects o...During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggered-grid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh （R1） waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method.展开更多
The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the M...The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.展开更多
Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectiv...Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectively,thereby limiting time step selection.Based on the definition of temporal and spatial FD operators,we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy.First,to simulate wave propagation beyond the conventional CFL stability limit,the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones.Second,to preserve the modeling accuracy,the velocity-dependent lengths of the temporal and spatial operators are adaptively varied.The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples,respectively.We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.展开更多
The paper presents a staggered-grid any even-order accurate finite-difference scheme for two-dimensional (2D), three-component (3C), first-order stress-velocity elastic wave equation and its stability condition in the...The paper presents a staggered-grid any even-order accurate finite-difference scheme for two-dimensional (2D), three-component (3C), first-order stress-velocity elastic wave equation and its stability condition in the arbitrary tilt anisotropic media; and derives a perfectly matched absorbing layer (PML) boundary condition and its stag- gered-grid any even-order accurate difference scheme in the 2D arbitrary tilt anisotropic media. The results of nu- merical modeling indicate that the modeling precision is high, the calculation efficiency is satisfactory and the absorbing boundary condition is better. The wave-front shapes of elastic waves are complex in the anisotropic media, and the velocity of qP wave is not always faster than that of qS wave. The wave-front triplication of qS wave and its events in both reflected domain and propagated domain, which are not commonly hyperbola, is a common phenomenon. When the symmetry axis is tilted in the TI media, the phenomenon of S-wave splitting is clearly observed in the snaps of three components and synthetic seismograms, and the events of all kinds of waves are asymmetric.展开更多
Forward modeling of elastic wave propagation in porous media has great importance for understanding and interpreting the influences of rock properties on characteristics of seismic wavefield. However,the finite-differ...Forward modeling of elastic wave propagation in porous media has great importance for understanding and interpreting the influences of rock properties on characteristics of seismic wavefield. However,the finite-difference forward-modeling method is usually implemented with global spatial grid-size and time-step; it consumes large amounts of computational cost when small-scaled oil/gas-bearing structures or large velocity-contrast exist underground. To overcome this handicap,combined with variable grid-size and time-step,this paper developed a staggered-grid finite-difference scheme for elastic wave modeling in porous media. Variable finite-difference coefficients and wavefield interpolation were used to realize the transition of wave propagation between regions of different grid-size. The accuracy and efficiency of the algorithm were shown by numerical examples. The proposed method is advanced with low computational cost in elastic wave simulation for heterogeneous oil/gas reservoirs.展开更多
In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.Ho...In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.However,convolution operator is intractable in finite-difference time-domain(FDTD)method.A great deal of progress has been made in using time stepping instead of convolution in FDTD.To incorporate PML into viscoelastic media,more memory variables need to be introduced,which increases the code complexity and computation costs.By modifying the nonsplitting PML formulation,I propose a viscoelastic model,which can be used as a viscoelastic material and/or a PML just by adjusting the parameters.The proposed viscoelastic model is essentially equivalent to a Maxwell model.Compared with existing PML methods,the proposed method requires less memory and its implementation in existing finite-difference codes is much easier.The attenuation and phase velocity of P-and S-waves are frequency independent in the viscoelastic model if the related quality factors(Q)are greater than 10.The numerical examples show that the method is stable for materials with high absorption(Q=1),and for heterogeneous media with large contrast of acoustic impedance and large contrast of viscosity.展开更多
Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a...Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.展开更多
To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of ...To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
In this study,copper extraction from low-grade oxide-sulfide ores was investigated using a leaching method combined with response surface methodology(RSM)to optimize operational conditions and assess leaching kinetics...In this study,copper extraction from low-grade oxide-sulfide ores was investigated using a leaching method combined with response surface methodology(RSM)to optimize operational conditions and assess leaching kinetics.Given copper's extensive industrial applications,sustainable recovery from low-grade ores is critical.Five key parameters-acid concentration,leaching time,particle size,temperature,and solids percentage-were identified as major influences on copper recovery.The results revealed that leaching time and solids percentage,along with interactions between temperature-time and temperature-solids percentage,had the most significant effects.Optimal conditions for 80% copper recovery while minimizing iron recovery below 3% included an acid concentration of 1.21 mol L^(-1),a leaching time of 108 min,a particle size of 438μm,a temperature of 45℃,and a solids percentage of 18.2%.Leaching kinetics were analyzed using shrinking core models,with the Dickinson model best describing the process,showing an activation energy of 32.63 kJ mol^(-1),indicative of mixed diffusion and chemical reaction control.The final kinetic model effectively predicted the influence of key parameters.These findings highlight the importance of optimizing process variables and selecting suitable kinetic models to enhance extraction efficiency,reduce costs,and improve sustainability in copper recovery.展开更多
Background Multibreed genomic prediction(MBGP)is crucial for improving prediction accuracy for breeds with small populations,for which limited data are often available.Recent studies have demonstrated that partitionin...Background Multibreed genomic prediction(MBGP)is crucial for improving prediction accuracy for breeds with small populations,for which limited data are often available.Recent studies have demonstrated that partitioning the genome into nonoverlapping blocks to model heterogeneous genetic(co)variance in multitrait models can achieve higher joint prediction accuracy.However,the block partitioning method,a key factor influencing model performance,has not been extensively explored.Results We introduce mbBayesABLD,a novel Bayesian MBGP model that partitions each chromosome into nonoverlapping blocks on the basis of linkage disequilibrium(LD)patterns.In this model,marker effects within each block are assumed to follow normal distributions with block-specific parameters.We employ simulated data as well as empirical datasets from pigs and beans to assess genomic prediction accuracy across different models using cross-validation.The results demonstrate that mbBayesABLD significantly outperforms conventional MBGP models,such as GBLUP and BayesR.For the meat marbling score trait in pigs,compared with GBLUP,which does not account for heterogeneous genetic(co)variance,mbBayesABLD improves the prediction accuracy for the small-population breed Landrace by 15.6%.Furthermore,our findings indicate that a moderate level of similarity in LD patterns between breeds(with an average correlation of 0.6)is sufficient to improve the prediction accuracy of the target breed.Conclusions This study presents a novel LD block-based approach for multibreed genomic prediction.Our work provides a practical tool for livestock breeding programs and offers new insights into leveraging genetic diversity across breeds for improved genomic prediction.展开更多
基金supported by the National Natural Science Foundation of China (40925013)the Strategic Leading Science and Technology Programme (Class B) of the Chinese Academy of Sciences (Grant No.XDB10010400)the China National Major Science and Technology Project (2011ZX05023-005004)
文摘Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with smallscale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference （FD） method of Biot＇s poroelastic equations, with an arbitrary even-order （2L） accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer （CPML） absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultra- sonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thick- nesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Qsc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments.
基金sponsed by National Natural Science Foundation of China(NSFC,Grant No.41304077)the Natural Basic Research Program of China(the“973 Project,”Grant No.2013CB733303)Postdoctoral Science Foundation of China(Grant No.2014T70740)
文摘During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggered-grid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh （R1） waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method.
基金supported by the Natural Science Foundation of China(No.41374078)Geological Survey Projects of Ministry of Land and Resources of China(No.12120113086100 and 12120113101300)
文摘The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.
基金the National Natural Science Foundation of China(No.41874144)the Research Foundation of China University of PetroleumBeijing at Karamay(RCYJ2018A-01-001).
文摘Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectively,thereby limiting time step selection.Based on the definition of temporal and spatial FD operators,we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy.First,to simulate wave propagation beyond the conventional CFL stability limit,the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones.Second,to preserve the modeling accuracy,the velocity-dependent lengths of the temporal and spatial operators are adaptively varied.The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples,respectively.We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.
基金Fund Project of Key Lab of Geophysical Exploration of China National Petroleum Corporation (GPR0408).
文摘The paper presents a staggered-grid any even-order accurate finite-difference scheme for two-dimensional (2D), three-component (3C), first-order stress-velocity elastic wave equation and its stability condition in the arbitrary tilt anisotropic media; and derives a perfectly matched absorbing layer (PML) boundary condition and its stag- gered-grid any even-order accurate difference scheme in the 2D arbitrary tilt anisotropic media. The results of nu- merical modeling indicate that the modeling precision is high, the calculation efficiency is satisfactory and the absorbing boundary condition is better. The wave-front shapes of elastic waves are complex in the anisotropic media, and the velocity of qP wave is not always faster than that of qS wave. The wave-front triplication of qS wave and its events in both reflected domain and propagated domain, which are not commonly hyperbola, is a common phenomenon. When the symmetry axis is tilted in the TI media, the phenomenon of S-wave splitting is clearly observed in the snaps of three components and synthetic seismograms, and the events of all kinds of waves are asymmetric.
基金supported by the National Basic Research Program of China (No. 2013CB228604)the National Science and Technology Major Project (No. 2011ZX05030-004-002,2011ZX05019-003)the National Natural Science Foundation (No. 41004050)
文摘Forward modeling of elastic wave propagation in porous media has great importance for understanding and interpreting the influences of rock properties on characteristics of seismic wavefield. However,the finite-difference forward-modeling method is usually implemented with global spatial grid-size and time-step; it consumes large amounts of computational cost when small-scaled oil/gas-bearing structures or large velocity-contrast exist underground. To overcome this handicap,combined with variable grid-size and time-step,this paper developed a staggered-grid finite-difference scheme for elastic wave modeling in porous media. Variable finite-difference coefficients and wavefield interpolation were used to realize the transition of wave propagation between regions of different grid-size. The accuracy and efficiency of the algorithm were shown by numerical examples. The proposed method is advanced with low computational cost in elastic wave simulation for heterogeneous oil/gas reservoirs.
文摘In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.However,convolution operator is intractable in finite-difference time-domain(FDTD)method.A great deal of progress has been made in using time stepping instead of convolution in FDTD.To incorporate PML into viscoelastic media,more memory variables need to be introduced,which increases the code complexity and computation costs.By modifying the nonsplitting PML formulation,I propose a viscoelastic model,which can be used as a viscoelastic material and/or a PML just by adjusting the parameters.The proposed viscoelastic model is essentially equivalent to a Maxwell model.Compared with existing PML methods,the proposed method requires less memory and its implementation in existing finite-difference codes is much easier.The attenuation and phase velocity of P-and S-waves are frequency independent in the viscoelastic model if the related quality factors(Q)are greater than 10.The numerical examples show that the method is stable for materials with high absorption(Q=1),and for heterogeneous media with large contrast of acoustic impedance and large contrast of viscosity.
基金supported by the National Natural Science Foundation of China(No.41474110)Shell Ph.D. Scholarship to support excellence in geophysical research
文摘Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.
文摘To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.
基金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.
基金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 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.
基金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.
基金Open Access funding enabled and organized by Projekt DEAL.
文摘In this study,copper extraction from low-grade oxide-sulfide ores was investigated using a leaching method combined with response surface methodology(RSM)to optimize operational conditions and assess leaching kinetics.Given copper's extensive industrial applications,sustainable recovery from low-grade ores is critical.Five key parameters-acid concentration,leaching time,particle size,temperature,and solids percentage-were identified as major influences on copper recovery.The results revealed that leaching time and solids percentage,along with interactions between temperature-time and temperature-solids percentage,had the most significant effects.Optimal conditions for 80% copper recovery while minimizing iron recovery below 3% included an acid concentration of 1.21 mol L^(-1),a leaching time of 108 min,a particle size of 438μm,a temperature of 45℃,and a solids percentage of 18.2%.Leaching kinetics were analyzed using shrinking core models,with the Dickinson model best describing the process,showing an activation energy of 32.63 kJ mol^(-1),indicative of mixed diffusion and chemical reaction control.The final kinetic model effectively predicted the influence of key parameters.These findings highlight the importance of optimizing process variables and selecting suitable kinetic models to enhance extraction efficiency,reduce costs,and improve sustainability in copper recovery.
基金supported by the Biological Breeding-Major Projects in National Science and Technology(No.2023ZD0404405)the Earmarked Fund for China Agriculture Research System(No.CARS-pig-35)+2 种基金the National Natural Science Foundation of China(No.3227284,32302708)the 2115 Talent Development Program of China Agricultural University,the Chinese Universities Scientific Fund(No.2023TC196)the Seed Industry Revitalization Action Project of Guangdong Province(No.2024-XPY-06-001)。
文摘Background Multibreed genomic prediction(MBGP)is crucial for improving prediction accuracy for breeds with small populations,for which limited data are often available.Recent studies have demonstrated that partitioning the genome into nonoverlapping blocks to model heterogeneous genetic(co)variance in multitrait models can achieve higher joint prediction accuracy.However,the block partitioning method,a key factor influencing model performance,has not been extensively explored.Results We introduce mbBayesABLD,a novel Bayesian MBGP model that partitions each chromosome into nonoverlapping blocks on the basis of linkage disequilibrium(LD)patterns.In this model,marker effects within each block are assumed to follow normal distributions with block-specific parameters.We employ simulated data as well as empirical datasets from pigs and beans to assess genomic prediction accuracy across different models using cross-validation.The results demonstrate that mbBayesABLD significantly outperforms conventional MBGP models,such as GBLUP and BayesR.For the meat marbling score trait in pigs,compared with GBLUP,which does not account for heterogeneous genetic(co)variance,mbBayesABLD improves the prediction accuracy for the small-population breed Landrace by 15.6%.Furthermore,our findings indicate that a moderate level of similarity in LD patterns between breeds(with an average correlation of 0.6)is sufficient to improve the prediction accuracy of the target breed.Conclusions This study presents a novel LD block-based approach for multibreed genomic prediction.Our work provides a practical tool for livestock breeding programs and offers new insights into leveraging genetic diversity across breeds for improved genomic prediction.
