Species’ geographical range, environmental range and traits lead to specimen collection preference of dominant plant species of grasslands in Northern China

Many different factors,such as species traits,socio-economic factors,geographical and environmental factors,can lead to specimen collection preference.This study aims to determine whether grassland specimen collection in China is preferred by species traits(i.e.,plant height,flowering and fruiting period),environmental range(i.e.,the temperature and precipitation range)and geographical range(i.e.,distribution range and altitudinal range).Ordinary least squares models and phylogenetic generalized linear mixed models were used to analyze the relationships between specimen number and the explanatory variables.Random Forest models were then used to find the most parsimonious multivariate model.The results showed that interannual variation in specimen number between 1900 and 2020 was considerable.Specimen number of these species in southeast China was notably lower than that in northwest China.Environmental range and geographical range of species had significant positive correlations with specimen number.In addition,there were relatively weak but significant associations between specimen number and species trait(i.e.,plant height and flowering and fruiting period).Random Forest models indicated that distribution range was the most important variable,followed by flowering and fruiting period,and altitudinal range.These findings suggest that future floristic surveys should pay more attention to species with small geographical range,narrow environmental range,short plant height,and short flowering and fruiting period.The correction of specimen collection preference will also make the results of species distribution model,species evolution and other works based on specimen data more accurate.

Comparison of Microwave Imaging Algorithms for Short-Range Scenarios

Three dimensional(3-D)imaging algorithms with irregular planar multiple-input-multiple-output(MIMO)arrays are discussed and compared with each other.Based on the same MIMO array,a modified back projection algorithm(MBPA)is accordingly proposed and four imaging algorithms are used for comparison,back-projection method(BP),back-projection one in time domain(BP-TD),modified back-projection one and fast Fourier transform(FFT)-based MIMO range migration algorithm(FFT-based MIMO RMA).All of the algorithms have been implemented in practical application scenarios by use of the proposed imaging system.Back to the practical applications,MIMO array-based imaging system with wide-bandwidth properties provides an efficient tool to detect objects hidden behind a wall.An MIMO imaging radar system,composed of a vector network analyzer(VNA),a set of switches,and an array of Vivaldi antennas,have been designed,fabricated,and tested.Then,these algorithms have been applied to measured data collected in different scenarios constituted by five metallic spheres in the absence and in the presence of a wall between the antennas and the targets in simulation and pliers in free space for experimental test.Finally,the focusing properties and time consumption of the above algorithms are compared.

The interaction between temperature and precipitation on the potential distribution range of Betula ermanii in the alpine treeline ecotone on the Changbai Mountain

Alpine treeline ecotones are highly sensitive to climate warming.The low temperature-determined alpine treeline is expected to shift upwards in response to global warming.However,little is known about how temperature interacts with other important factors to influence the distribution range of tree species within and beyond the alpine treeline ecotone.Hence,we used a GF-2 satellite image,along with bioclimatic and topographic variables,to develop an ensemble suitable habitat model based on the species distribution modeling algorithms in Biomod2.We investigated the distribution of suitable habitats for B.ermanii under three climate change scenarios(i.e.,low(SSP126),moderate(SSP370)and extreme(SSP585)future emission trajectories)between two consecutive time periods(i.e.,current-2055,and 2055-2085).By 2055,the potential distribution range of B.ermanii will expand under all three climate scenarios.The medium and high suitable areas will decline under SSP370 and SSP585scenarios from 2055 to 2085.Moreover,under the three climate scenarios,the uppermost altitudes of low suitable habitat will rise to 2,329 m a.s.l.,while the altitudes of medium and high suitable habitats will fall to 2,201 and2,051 m a.s.l.by 2085,respectively.Warming promotes the expansion of B.ermanii distribution range in Changbai Mountain,and this expansion will be modified by precipitation as climate warming continues.This interaction between temperature and precipitation plays a significant role in shaping the potential distribution range of B.ermanii in the alpine treeline ecotone.This study reveals the link between environmental factors,habitat distribution,and species distribution in the alpine treeline ecotone,providing valuable insights into the impacts of climate change on high-elevation vegetation,and contributing to mountain biodiversity conservation and sustainable development.

Spatiotemporal Evaluation and Future Projection of Diurnal Temperature Range over the Tibetan Plateau in CMIP6 Models

The diurnal temperature range(DTR) serves as a vital indicator reflecting both natural climate variability and anthropogenic climate change. This study investigates the historical and projected multitemporal DTR variations over the Tibetan Plateau. It assesses 23 climate models from phase 6 of the Coupled Model Intercomparison Project(CMIP6) using CN05.1 observational data as validation, evaluating their ability to simulate DTR over the Tibetan Plateau. Then, the evolution of DTR over the Tibetan Plateau under different shared socioeconomic pathway(SSP) scenarios for the near,middle, and long term of future projection are analyzed using 11 selected robustly performing models. Key findings reveal:(1) Among the models examined, BCC-CSM2-MR, EC-Earth3, EC-Earth3-CC, EC-Earth3-Veg, EC-Earth3-Veg-LR,FGOALS-g3, FIO-ESM-2-0, GFDL-ESM4, MPI-ESM1-2-HR, MPI-ESM1-2-LR, and INM-CM5-0 exhibit superior integrated simulation capability for capturing the spatiotemporal variability of DTR over the Tibetan Plateau.(2) Projection indicates a slightly increasing trend in DTR on the Tibetan Plateau in the SSP1-2.6 scenario, and decreasing trends in the SSP2-4.5, SSP3-7.0, and SPP5-8.5 scenarios. In certain areas, such as the southeastern edge of the Tibetan Plateau, western hinterland of the Tibetan Plateau, southern Kunlun, and the Qaidam basins, the changes in DTR are relatively large.(3) Notably, the warming rate of maximum temperature under SSP2-4.5, SSP3-7.0, and SPP5-8.5 is slower compared to that of minimum temperature, and it emerges as the primary contributor to the projected decrease in DTR over the Tibetan Plateau in the future.

Design of a high sensitivity and wide range angular rate sensor based on exceptional surface

It is found that when the parity–time symmetry phenomenon is introduced into the resonant optical gyro system and it works near the exceptional point,the sensitivity can in theory be significantly amplified at low angular rate.However,in fact,the exceptional point is easily disturbed by external environmental variables,which means that it depends on harsh experimental environment and strong control ability,so it is difficult to move towards practical application.Here,we propose a new angular rate sensor structure based on exceptional surface,which has the advantages of high sensitivity and high robustness.The system consists of two fiber-optic ring resonators and two optical loop mirrors,and one of the resonators contains a variable ratio coupler and a variable optical attenuator.We theoretically analyze the system response,and the effects of phase and coupling ratio on the system response.Finally,compared with the conventional resonant gyro,the sensitivity of this exceptional surface angular rate sensor can be improved by about 300 times at low speed.In addition,by changing the loss coefficient in the ring resonator,we can achieve a wide range of 600 rad/s.This scheme provides a new approach for the development of ultra-high sensitivity and wide range angular rate sensors in the future.

Robust Timing Constraints for Granitic Magmatism and Hydrothermal Mineralization in the Tieshanlong W-Sn Ore Field,Eastern Nanling Range,South China

The Tieshanlong ore field is an important part of the Nanling Range,which is famous worldwide for its W-Sn mineralization.Notably,the mineralization age of the Tieshanlong ore field is not well constrained,and our field investigation reveals that granitic emplacement occurred at different stages.However,previous studies have not distinguished these multiple stages of magmatism.The Tieshanlong granite complex is closely related to the Huangsha quartz vein-type W-Sn deposit and Tongling skarn-type Cu-W-Sn deposit in this field.Through field investigations and isotopic age analyses,this work studies the relationship between multistage magmatic activity and mineralization in the Tieshanlong ore field.LA-ICP-MS zircon U-Pb isotope analyses revealed that the first-and second-staged granites formed at 154.2±0.6 Ma(MSDW=1.4)and 151.2±0.4 Ma(MSDW=1.5),with zirconε_(Hf)(t)values ranging from-13.1 to-10.5 and from-14.7 to-11.1,respectively.These data suggest that the Tieshanlong granite complex was derived from the partial melting of ancient crustal material.LA-ICP-MS U-Pb dating of wolframite and cassiterite reveals that W-Sn mineralization occurred at 160-150 Ma,which agrees well with the U-Pb dating results of the second-staged granite within analytical errors.The magmatic activity in this ore field can be divided into three stages:175-154 Ma,154-150 Ma and 150-145 Ma.The quartz vein-and skarn-type W-Sn mineralization is closely related to second-staged fine-grained twomica granite,and formed earlier than skarn-type Cu-mineralization.This study establishes a metallogenic model for the Tieshanlong ore field,and this model has important practical significance for identifying concealed W-Sn(-Cu)deposits around other granitic complexes in the Nanling Range.

Influence of network structure on spreading dynamics via tie range

There are various phenomena of malicious information spreading in the real society, which cause many negative impacts on the society. In order to better control the spreading, it is crucial to reveal the influence of network structure on network spreading. Motifs, as fundamental structures within a network, play a significant role in spreading. Therefore, it is of interest to investigate the influence of the structural characteristics of basic network motifs on spreading dynamics.Considering the edges of the basic network motifs in an undirected network correspond to different tie ranges, two edge removal strategies are proposed, short ties priority removal strategy and long ties priority removal strategy. The tie range represents the second shortest path length between two connected nodes. The study focuses on analyzing how the proposed strategies impact network spreading and network structure, as well as examining the influence of network structure on network spreading. Our findings indicate that the long ties priority removal strategy is most effective in controlling network spreading, especially in terms of spread range and spread velocity. In terms of network structure, the clustering coefficient and the diameter of network also have an effect on the network spreading, and the triangular structure as an important motif structure effectively inhibits the spreading.

Study of the response of 10B-doped MCP to wide-energy range neutrons from eV to MeV

Neutron-sensitive microchannel plates(nMCPs)have applications in neutron detection,including energy spectrum measurements,neutron-induced cross sections,and neutron imaging.10B-doped MCPs(B-MCPs)have attracted significant attention owing to their potential for exhibiting a high neutron detection efficiency over a large neutron energy range.Good spatial and temporal resolutions are useful for neutron energy-resolved imaging.However,their practical applications still face many technical challenges.In this study,a B-MCP with 10 mol%10B was tested for its response to wide-energy neutrons from eV to MeV at the Back-n white neutron source at the China Spallation Neutron Source.The neutron detection efficiency was calibrated at 1 eV,which is approximately 300 times that of an ordinary MCP and indicates the success of 10 B doping.The factors that caused the reduction in the detection efficiency were simulated and discussed.The neutron energy spectrum obtained using B-MCP was compared with that obtained by other measurement methods,and showed very good consistency for neutron energies below tens of keV.The response is more complicated at higher neutron energy,at which point the elastic and nonelastic reactions of all nuclides of B-MCP gradually become dominant.This is beneficial for the detection of neutrons,as it compensates for the detection efficiency of B-MCP for high-energy neutrons.

Image Steganography by Pixel-Value Differencing Using General Quantization Ranges

A new steganographic method by pixel-value differencing(PVD)using general quantization ranges of pixel pairs’difference values is proposed.The objective of this method is to provide a data embedding technique with a range table with range widths not limited to powers of 2,extending PVD-based methods to enhance their flexibility and data-embedding rates without changing their capabilities to resist security attacks.Specifically,the conventional PVD technique partitions a grayscale image into 1×2 non-overlapping blocks.The entire range[0,255]of all possible absolute values of the pixel pairs’grayscale differences in the blocks is divided into multiple quantization ranges.The width of each quantization range is a power of two to facilitate the direct embedding of the bit information with high embedding rates.Without using power-of-two range widths,the embedding rates can drop using conventional embedding techniques.In contrast,the proposed method uses general quantization range widths,and a multiple-based number conversion mechanism is employed skillfully to implement the use of nonpower-of-two range widths,with each pixel pair being employed to embed a digit in the multiple-based number.All the message bits are converted into a big multiple-based number whose digits can be embedded into the pixel pairs with a higher embedding rate.Good experimental results showed the feasibility of the proposed method and its resistance to security attacks.In addition,implementation examples are provided,where the proposed method adopts non-power-of-two range widths and employsmultiple-based number conversion to expand the data-hiding and steganalysis-resisting capabilities of other PVD methods.

Simulation Study on the Migration Range of CO_(2) in the Offshore Saline Aquifer

The geological storage of carbon dioxide(CO_(2)) is a crucial technology for mitigating climate change. Offshore deep saline aquifers have elicited increased attention due to their remarkable potential for storing CO_(2). During long-term storage, CO_(2) migration in a deep saline aquifer needs special attention to prevent it from reaching risk points and leading to security issues. In this paper, a mechanism model is established according to the geological characteristics of saline aquifers in an offshore sedimentary basin in China. The CO_(2) migration over 100 years is simulated considering geological changes such as permeability, dip angle, thickness, and salinity. The effects of injection conditions on the CO_(2) migration range are also investigated. Results reveal that the migration range of CO_(2) in the injection period exceeds 70%, even if the postinjection period's duration is five times longer than that of the injection period. As the values of the above geological parameters increase, the migration range of CO_(2) increases, and permeability has a particularly substantial influence. Moreover, the influences of injection rate and well type are considerable. At high injection rates, CO_(2) has a greater likelihood of displacing brine in a piston-like scheme. CO_(2) injected by long horizontal wells migrates farther compared with that injected by vertical wells. In general, the plane migration range is within 3 000 m, although variations in the reservoir and injection parameters of the studied offshore saline aquifers are considered. This paper can offer references for the site selection and injection well deployment of CO_(2) saline aquifer storage. According to the studied offshore aquifers, a distance of at least 3 000 m from potential leakage points, such as spill points, active faults, and old abandoned wells, must be maintained.

Extraction of Acoustic Normal Mode Depth Functions Using Range-Difference Method with Vertical Linear Array Data

Data-derived normal mode extraction is an effective method for extracting normal mode depth functions in the absence of marine environmental data.However,when the corresponding singular vectors become nonunique when two or more singular values obtained from the cross-spectral density matrix diagonalization are nearly equal,this results in unsatisfactory extraction outcomes for the normal mode depth functions.To address this issue,we introduced in this paper a range-difference singular value decomposition method for the extraction of normal mode depth functions.We performed the mode extraction by conducting singular value decomposition on the individual frequency components of the signal's cross-spectral density matrix.This was achieved by using pressure and its range-difference matrices constructed from vertical line array data.The proposed method was validated using simulated data.In addition,modes were successfully extracted from ambient noise.

Analysis of application range of simplified models for field to thermo-field to thermionic emission processes from the cathode

Electron emission plays a dominant role in plasma-cathode interactions and is a key factor in many plasma phenomena and industrial applications.It is necessary to illustrate the various electron emission mechanisms and the corresponding applicable description models to evaluate their impacts on discharge properties.In this study,detailed expressions of the simplified formulas valid for field emission to thermo-field emission to thermionic emission typically used in the numerical simulation are proposed,and the corresponding application ranges are determined in the framework of the Murphy-Good theory,which is commonly regarded as the general model and to be accurate in the full range of conditions of the validity of the theory.Dimensionless parameterization was used to evaluate the emission current density of the Murphy-Good formula,and a deviation factor was defined to obtain the application ranges for different work functions(2.5‒5 eV),cathode temperatures(300‒6000 K),and emitted electric fields(10^(5) to 10^(10) V·m^(-1)).The deviation factor was shown to be a nonmonotonic function of the three parameters.A comparative study of particle number densities in atmospheric gas discharge with a tungsten cathode was performed based on the one-dimensional implicit particle-in-cell(PIC)with the Monte Carlo collision(MCC)method according to the aforementioned application ranges.It was found that small differences in emission current density can lead to variations in the distributions of particle number density due to changes in the collisional environment.This study provides a theoretical basis for selecting emission models for subsequent numerical simulations.

Design and Verification for Dual⁃mode CDFS and High⁃Load Compressor with a Large Flow Regulation Range

This paper presents the design and verification of the dual-mode core driven fan stage(CDFS)and high-load compressor with a large flow regulation range.In view of the characteristics of large flow regulation range of the two modes and high average stage load coefficient,this paper investigates the design technology of the dual-mode high-efficiency compressor with a large flow regulation range and high-load compressor with an average stage load coefficient of 0.504.Building upon this research,the design of the dual-mode CDFS and four-stage compressor is completed,and three-dimensional numerical simulation of the two modes is carried out.Finally,performance experiment is conducted to verify the result of three-dimensional numerical simulation.The experiment results show that the compressor performance is improved for the whole working conditions by using the new design method,which realizes the complete fusion design of the CDFS and high-pressure compressor(HPC).The matching mechanism of stage characteristics of single and double bypass modes and the variation rule of different adjustment angles on performance are studied comprehensively.Furthermore,it effectively reduces the length and weight of compressor,and breaks through the key technologies such as high-load compressor with the average load factor of 0.504.These findings provide valuable data and a methodological foundation for the development of the next generation aeroengine.

Scale effect removal and range migration correction for hypersonic target coherent detection

The detection of hypersonic targets usually confronts range migration(RM)issue before coherent integration(CI).The traditional methods aiming at correcting RM to obtain CI mainly considers the narrow-band radar condition.However,with the increasing requirement of far-range detection,the time bandwidth product,which is corresponding to radar’s mean power,should be promoted in actual application.Thus,the echo signal generates the scale effect(SE)at large time bandwidth product situation,influencing the intra and inter pulse integration performance.To eliminate SE and correct RM,this paper proposes an effective algorithm,i.e.,scaled location rotation transform(ScLRT).The ScLRT can remove SE to obtain the matching pulse compression(PC)as well as correct RM to complete CI via the location rotation transform,being implemented by seeking the actual rotation angle.Compared to the traditional coherent detection algorithms,Sc LRT can address the SE problem to achieve better detection/estimation capabilities.At last,this paper gives several simulations to assess the viability of ScLRT.

LFM Radar Source Passive Localization Algorithm Based on Range Migration

Traditional single-satellite passive localization algorithms are influenced by frequency and angle measurement accuracies,resulting in error estimation of emitter position on the order of kilometers.Subsequently,a single-satellite localization algorithm based on passive synthetic aper-ture(PSA)was introduced,enabling high-precision positioning.However,its estimation of azimuth and range distance is considerably affected by the residual frequency offset(RFO)of uncoopera-tive system transceivers.Furthermore,it requires data containing a satellite flying over the radia-tion source for RFO search.After estimating the RFO,an accurate estimation of azimuth and range distance can be carried out,which is difficult to achieve in practical situations.An LFM radar source passive localization algorithm based on range migration is proposed to address the dif-ficulty in estimating frequency offset.The algorithm first provides a rough estimate of the pulse repetition time(PRT).It processes intercepted signals through range compression,range interpola-tion,and polynomial fitting to obtain range migration observations.Subsequently,it uses the changing information of range migration and an accurate PRT to formulate a system of nonlinear equations,obtaining the emitter position and a more accurate PRT through a two-step localization algorithm.Frequency offset only induces a fixed offset in range migration,which does not affect the changing information.This algorithm can also achieve high-precision localization in squint scenar-ios.Finally,the effectiveness of this algorithm is verified through simulations.

FINITE RANGE SET WITH TRUNCATED MULTIPLICITY FOR MEROMORPHIC FUNCTIONS ON SOME COMPLEX DISC

In this paper,we consider the truncated multiplicity finite range set problem of meromorphic functions on some complex disc.By using the value distribution theory of meromorphic functions,we establish a second main theorem for meromorphic functions with finite growth index which share meromorphic functions(may not be small functions).As its application,we also extend the result of a finite range set with truncated multiplicity.

Relationship between hemoglobin glycation index and risk of hypoglycemia in type 2 diabetes with time-in-range in target

BACKGROUND In patients with type 2 diabetes mellitus(T2DM),the risk of hypoglycemia also occurs in at a time-in-range(TIR)of>70%.The hemoglobin glycation index(HGI)is considered the best single factor for predicting hypoglycemia,and offers new perspectives for the individualized treatment of patients with well-controlled blood glucose levels that are easily ignored in clinical settings.All participants underwent a 7-days continuous glucose monitoring(CGM)using a retrospective CGM system.We obtained glycemic variability indices using the CGM system.We defined HGI as laboratory hemoglobin A1c minus the glucose management indicator.Patients were categorized into low HGI(HGI<0.5)and high HGI groups(HGI≥0.5)according to HGI median(0.5).Logistic regression and receiver operating characteristic curve analyses were used to determine the risk factors for hypoglycemia.RESULTS We included 129 subjects with T2DM(54.84±12.56 years,46%male)in the study.Median TIR score was 90%.The high HGI group exhibited lower TIR and greater time below range with higher hemoglobin A1c than the low HGI group;this suggests more glycemic excursions and an increased incidence of hypoglycemia in the high HGI group.Multivariate analyses revealed that mean blood glucose,standard deviation of blood glucose and HGI were independent risk factors for hypoglycemia.Receiver operating characteristic curve analysis indicated that the HGI was the best predictor of hypoglycemia.In addition,the optimal cut-off points for HGI,mean blood glucose,and standard deviation of blood glucose in predicting hypoglycemia were 0.5%,7.2 mmol/L and 1.4 mmol/L respectively.CONCLUSION High HGI was significantly associated with greater glycemic excursions and increased hypoglycemia in patients with TIR>70%.Our findings indicate that HGI is a reliable predictor of hypoglycemia in this population.

A Residual Frequency Offset Estimation Method Based on Range Migration Fitting

Recently,researchers have proposed an emitter localization method based on passive synthetic aperture.However,the unknown residual frequency offset(RFO)between the transmit-ter and the receiver causes the received Doppler signal to shift,which affects the localization accu-racy.To solve this issue,this paper proposes a RFO estimation method based on range migration fitting.Due to the high frequency modulation slope of the linear frequency modulation(LFM)-mod-ulation radar signal,it is not affected by RFO in range compression.Therefore,the azimuth time can be estimated by fitting the peak value position of the pulse compression in range direction.Then,the matched filters are designed under different RFOs.When the zero-Doppler time obtained by the matched filters is consistent with the estimated azimuth time,the given RFO is the real RFO between the transceivers.The simulation results show that the estimation error of azimuth distance does not exceed 20 m when the received signal duration is not less than 3 s,the pulse repe-tition frequency(PRF)of the transmitter radar signal is not less than 1 kHz,the range detection is not larger than 1000 km,and the signal noise ratio(SNR)is not less than-5 dB.

Extending the solid solution range of sodium ferric pyrophosphate:Off‐stoichiometric Na_(3)Fe(2.5)(P_(2)O_(7))_(2)as a novel cathode for sodium‐ion batteries

Iron‐based pyrophosphates are attractive cathodes for sodium‐ion batteries due to their large framework,cost‐effectiveness,and high energy density.However,the understanding of the crystal structure is scarce and only a limited candidates have been reported so far.In this work,we found for the first time that a continuous solid solution,Na_(4−α)Fe_(2+α)_(2)(P_(2)O_(7))_(2)(0≤α≤1,could be obtained by mutual substitution of cations at center‐symmetric Na3 and Na4 sites while keeping the crystal building blocks of anionic P_(2)O_(7) unchanged.In particular,a novel off‐stoichiometric Na_(3)Fe(2.5)(P_(2)O_(7))_(2)is thus proposed,and its structure,energy storage mechanism,and electrochemical performance are extensively investigated to unveil the structure–function relationship.The as‐prepared off‐stoichiometric electrode delivers appealing performance with a reversible discharge capacity of 83 mAh g^(−1),a working voltage of 2.9 V(vs.Na^(+)/Na),the retention of 89.2%of the initial capacity after 500 cycles,and enhanced rate capability of 51 mAh g^(−1)at a current density of 1600 mA g^(−1).This research shows that sodium ferric pyrophosphate could form extended solid solution composition and promising phase is concealed in the range of Na_(4−α)Fe_(2+α)_(2)(P_(2)O_(7))_(2),offering more chances for exploration of new cathode materials for the construction of high‐performance SIBs.

Spatial-temporal distribution and geochemistry of highly evolved Mesozoic granites in Great Xing’an Range,NE China:Discriminant criteria and geological significance

Highly evolved granite is an important sign of the mature continent crust and closely associated with deposits of rare metals.In this work,the authors undertake systematically zircon U-Pb ages and whole rock elemental data for highly evolved granitic intrusions from the Great Xing’an Range(GXR),NE China,to elucidate their discriminant criteria,spatial-temporal distribution,differentiation and geodynamic mecha-nism.Geochemical data of these highly evolved granites suggest that high w(SiO_(2))(>70%)and differentiation index(DI>88)could be quantified indicators,while strong Eu depletion,high TE_(1,3),lowΣREE and low Zr/Hf,Nb/Ta,K/Rb could only be qualitative indicators.Zircon U-Pb ages suggest that the highly evolved gran-ites in the GXR were mainly formed in Late Mesozoic,which can be divided into two major stages:Late Ju-rassic-early Early Cretaceous(162-136 Ma,peak at 138 Ma),and late Early Cretaceous(136-106 Ma,peak at 126 Ma).The highly evolved granites are mainly distributed in the central-southern GXR,and display a weakly trend of getting younger from northwest to southeast,meanwhile indicating the metallogenic potential of rare metals within the central GXR.The spatial-temporal distribution,combined with regional geological data,indicates the highly evolved Mesozoic granites in the GXR were emplaced in an extensional environ-ment,of which the Late Jurassic-early Early Cretaceous extension was related to the closure of the Mongol-Okhotsk Ocean and roll-back of the Paleo-Pacific Plate,while the late Early Cretaceous extension was mainly related to the roll-back of the Paleo-Pacific Plate.