Assessment of heavy metal pollution in Laizhou Bay(China)using the ecological risk index and the integrated biomarker response of the goby Acanthogobius ommaturus

We used the Integrated Biological Responses version 2(IBRv2)method to evaluate the biological eff ects of heavy metals in the sediments in Laizhou Bay,China on the benthic goby Acanthogobius ommaturus.In December 2018,gobies and sediments were collected from 15 stations.We measured the activities of defense enzymes and the contents of malondialdehyde(MDA)and metallothionein(MT)in the goby liver as well as the levels of heavy metals in the sediments and goby muscle tissue.Most of the heavy metal concentrations in sediment at each station were below the Class I criteria set by Chinese Standards for Marine Sediment Quality,and the Håkanson ecological risk index suggested low risk for the heavy metals.We found that A.ommaturus could eff ectively accumulate mercury,cadmium,arsenic,and zinc and that the contents of MT and MDA and the activities of glutathione peroxidase and glutathione reductase were suitable biomarkers of heavy metal pollution in this species.The IBRv2 method integrated these four biomarkers and discriminated stations according to heavy metal pollution.Higher IBRv2 values suggested more adverse eff ects in gobies,corroborating more serious heavy metal contamination.The stations with high IBRv2 values and high contents of heavy metals were mainly distributed in the west and northeast parts of the bay.These results show that the IBRv2 approach is a feasible strategy for assessing heavy metal pollution through biological response and biological status and that it can be implemented for environmental monitoring in Laizhou Bay.

Study of Caspian Goby Neogobius sp. Karyotype Flexibility from Several Biotops

The karyotype of Caspian goby was studied, which has been identified as Caspian bighead goby-Neogobius gorlap before. The results of cytogenetic analyses have shown that the diploid set of goby varies from 38 to 40 chromosomes, which is different from bighead goby （2n = 43-46）. We assumed that Caspian goby has independence as a species.

Application of otolith shape analysis for stock discrimination and species identification of five goby species(Perciformes: Gobiidae) in the northern Chinese coastal waters

We tested the use of otolith shape analysis to discriminate between species and stocks of five goby species( Ctenotrypauchen chinensis, Odontamblyopus lacepedii, Amblychaeturichthys hexanema, Chaeturichthys stigmatias, and Acanthogobius hasta) found in northern Chinese coastal waters. The five species were well differentiated with high overall classification success using shape indices(83.7%), elliptic Fourier coefficients(98.6%), or the combination of both methods(94.9%). However, shape analysis alone was only moderately successful at discriminating among the four stocks(Liaodong Bay, LD; Bohai Bay, BH; Huanghe(Yellow) River estuary HRE, and Jiaozhou Bay, JZ stocks) of A. hasta(50%–54%) and C. stigmatias(65.7%–75.8%). For these two species, shape analysis was moderately successful at discriminating the HRE or JZ stocks from other stocks, but failed to effectively identify the LD and BH stocks. A large number of otoliths were misclassified between the HRE and JZ stocks, which are geographically well separated. The classification success for stock discrimination was higher using elliptic Fourier coefficients alone(70.2%) or in combination with shape indices(75.8%) than using only shape indices(65.7%) in C. stigmatias whereas there was little difference among the three methods for A. hasta. Our results supported the common belief that otolith shape analysis is generally more effective for interspecific identification than intraspecific discrimination. Moreover, compared with shape indices analysis, Fourier analysis improves classification success during inter- and intra-species discrimination by otolith shape analysis, although this did not necessarily always occur in all fish species.

Protective benefit of folded linear HDPE board sand fences along the Golmud-Korla Railway,China:Field observation and wind tunnel study

The Milan Gobi area of the Golmud-Korla Railway in northwest China is located in the lower dispersal area of the mountain pass and has strong winds with evident double wind direction characteristics.This study introduced a novel sand fence deployment technique,termed‘folded linear deployment',designed to position the sand fence orthogonally to the two predominant wind directions for optimal protection.This study used wind tunnel and field tests to evaluate the wind and sand flow characteristics,as well as the windproof and sandresistant performance of folded linear HDPE(Highdensity polyethylene)board sand fences.The results suggest that the airflow around the fence creates clear zoning characteristics.The deceleration area on the BSF(backwind side of the sand fence)is much larger than that on the DSF(downwind side of the sand fence).Thus,sand particles are primarily deposited on the BSF.At different wind speeds,the airflow at 2 and 5 h on the DSF is not disturbed.The WSP(wind speed profile)presents a logarithmic distribution.The airflow is disturbed at 1-20 h on the BSF,and the WSP gradually deviates from the logarithmic law.However,as the airflow moves away from the fence,the WSP gradually approaches a logarithmic distribution.Meanwhile,the WPE(windproof efficiency)and SRE(sand-resistant efficiency)of the sand fence exceed 80%.In addition,the results of wind tunnel tests are compared with those of field tests.The overall dispersion is good,and the best dispersion is found at z/H=2.00,indicating good agreement between the two test results.This study provides a scientific basis for the design of sand hazard control measures,similar to the railway project in the Gobi Gale area.

Three-dimensional(3D)parametric measurements of individual gravels in the Gobi region using point cloud technique

Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materials constituting the Gobi result in notable differences in saltation processes across various Gobi surfaces.It is challenging to describe these processes according to a uniform morphology.Therefore,it becomes imperative to articulate surface characteristics through parameters such as the three-dimensional(3D)size and shape of gravel.Collecting morphology information for Gobi gravels is essential for studying its genesis and sand saltation.To enhance the efficiency and information yield of gravel parameter measurements,this study conducted field experiments in the Gobi region across Dunhuang City,Guazhou County,and Yumen City(administrated by Jiuquan City),Gansu Province,China in March 2023.A research framework and methodology for measuring 3D parameters of gravel using point cloud were developed,alongside improved calculation formulas for 3D parameters including gravel grain size,volume,flatness,roundness,sphericity,and equivalent grain size.Leveraging multi-view geometry technology for 3D reconstruction allowed for establishing an optimal data acquisition scheme characterized by high point cloud reconstruction efficiency and clear quality.Additionally,the proposed methodology incorporated point cloud clustering,segmentation,and filtering techniques to isolate individual gravel point clouds.Advanced point cloud algorithms,including the Oriented Bounding Box(OBB),point cloud slicing method,and point cloud triangulation,were then deployed to calculate the 3D parameters of individual gravels.These systematic processes allow precise and detailed characterization of individual gravels.For gravel grain size and volume,the correlation coefficients between point cloud and manual measurements all exceeded 0.9000,confirming the feasibility of the proposed methodology for measuring 3D parameters of individual gravels.The proposed workflow yields accurate calculations of relevant parameters for Gobi gravels,providing essential data support for subsequent studies on Gobi environments.

Behavioral dependent dispersal in the invasive round goby Neogobius melanostomus depends on population age

Biological invasions cause major ecological and economic costs in invaded habitats. The round goby Neogobius me- lanostomus is a successful invasive species and a major threat to the biodiversity and ecological function of the Baltic Sea. It is native to the Ponto-Caspian region and has, via ballast water transport of ships, invaded the Gulf of Gdansk in Poland. Since 1990, it has spread as far north as Raahe in Northern Finland （64~41 ＂04＂N, 24~28＂44＂E）. Over the past decade, consistent indi- vidual differences of behavioral expressions have been shown to explain various ecological processes such as dispersal, survival or reproduction. We have previously shown that new and old populations differ in personality trait expression. Individuals in new populations are bolder, less sociable and more active than in old populations. Here we investigate if the behavioral differentiation can be explained by phenotype-dependent dispersal. This was investigated by measuring activity, boldness and sociability of in- dividually marked gobies, and subsequently allowing them to disperse in a system composed of five consecutive tanks connected by tubes. Individual dispersal tendency and distance was measured. Our results revealed that in newly established populations, more active individuals disperse sooner and that latency of a group to disperse depends on the mean sociability of the group. This indicates the presence of personality dependent dispersal in this species and that it is maintained at the invasion front but lost as the populations get older

Species Identification of a Commonly Believed Sinarapan Fish Mystichtys luzonensis （Smith, 1902） in Lakes Buhi and Bato of Bicol Region

The paper aimed to determine the true identity of a presumed or commonly believed Sinarapan fish and compare their phylogenetic relationships by using DNA barcoding. The fish samples were collected by researchers from a fisheries research and development institution in April 2015 from four different lakes of Bicol Region, i.e., Lakes Buhi, Danao, Manapao and Bato. DNA was extracted using cetyl trimethyl ammonium bromide （CTAB） extraction buffer. The polymerase chain reaction （PCR） profile for the reaction was 94 ℃ for 10 min, followed by 35 cycles of 1 min at 94℃, 1 min at 48℃ and 1.5 min at 72℃, and a final extension of 10 min at 72℃. The CO 1 region with approximately 650 bp was amplified because of its capability to differentiate taxa. Sequencing was done by Macrogen while phylogenetic analysis was produced using a Molecular Evolutionary Genetics Analysis （MEGA） software version 6.0. The findings showed that CO1 can be used as a DNA marker in molecular identification of the fish samples. Samples from two of the four lakes were successfully sequenced. From basic local alignment search tool （BLAST） results, the maximum parsimony （MP） and neighbor-joining （N J） trees revealed that samples from Lakes Buhi and Bato are not species of Sinarapan but rather species of Leiopotherapon plumbeus and Rhinogobius giurinus, respectively. Furthermore, DNA barcoding is very useful in proving the true identity of unknown samples.

The Expanding Earth: Indisputable Evidences of the Gobi Desert

The most striking contrasts that are found on the continents in paleogeographic reconstructions of the end of the Mesozoic era are the occurrence on the place of the disappeared humid subtropics of the largest Gobi Desert in Eurasia with air temperatures falling below 50°from the freezing point and annual precipitation totals at the level of 100 mm. Science does not know the processes that can lead to a cooling of the atmosphere at 70°and other equally radical changes in nature with a stable position of the blocks of the earth’s crust in space. Changes in the environment of this magnitude can only be the result of land moving northward for a distance equal to about half the radius of the Earth. Titanosaurs, described by the remains in the Gobi deposits, had a body volume, which at modern gravity corresponds to a mass of 10 to 30 ton. However, animals with such a mass and such growth could not exist now. To create the necessary pressure in the vessels and provide energy to the body, Mongolian sauropods would need a heart of 2 - 3 m in diameter. Known types of muscle tissue are unable to maintain an elongated neck and head with a mass of more than a ton. The femur bones of four-legged dinosaurs had strength sufficient to move on land only animals weighing no more than 5 - 7 tons. The bones of giant bipedal dinosaurs at a constant gravitational field would have to be subjected to specific loads, several times greater than the bones of modern elephants, which is excluded. The natural conclusion about the action of a lesser gravity in the Mesozoic provides a solution to the mystery of the truly global spread in that era of bipedal mode of movement as the most energy efficient.

Spatial pattern of plant species diversity and the influencing factors in a Gobi Desert within the Heihe River Basin, Northwest China

Understanding the spatial pattern of plant species diversity and the influencing factors has important implications for the conservation and management of ecosystem biodiversity. The transitional zone between biomes in desert ecosystems, however, has received little attention in that regard. In this study, we conducted a quantitative field survey （including 187 sampling plots） in a 40-km2 study area to determine the spatial pattern of plant species diversity and analyze the influencing factors in a Gobi Desert within the Heihe River Basin, Northwest China. A total of 42 plant species belonging to 16 families and 39 genera were recorded. Shrub and semi-shrub species generally represented the major part of the plant communities （covering 90% of the land surface）, while annual and perennial herbaceous species occupied a large proportion of the total recorded species （71%）. Patrick richness index （R）, Shannon-Wiener diversity index （H）, Simpson＇s dominance index （D）, and Pielou＇s evenness index （I） were all moderately spadaUy variable, and the variability increased with increasing sampling area. The semivariograms for R and H＇ were best fitted with Gaussian models while the semivariograms for D andJ were best fitted with exponential models. Nugget-to-still ratios indicated a moderate spatial autocorrelation for R, H＇, and D while a strong spatial autocorrelation was observed for J. The spatial patterns of R and H＇ were closely related to the geographic location within the study area, with lower values near the oasis and higher values near the mountains. However, there was an opposite trend for D. R, H＇, and D were significantly correlated with elevation, soil texture, bulk density, saturated hydraulic conductivity, and total porosity （P〈0.05）. Generally speaking, locations at higher elevations tended to have higher species richness and diversity and the higher elevations were characterized by higher values in sand and gravel contents, bulk density, and saturated hydraulic conductivity and also by lower values in total porosity. Furthermore, spatial variability of plant species diversity was dependent on the sampling area.

Observation and Study of Land Surface Parameters over Gobi in Typical Arid Region

According to the need of popular land surface process models, characteristics and rules of some key land surface process and soil parameters over Gobi in typical arid region of Northwest China are analyzed by using the data observed during the intensive observation period of the Dunhuang Land–Surface Process Field Experiment (DLSPFE) (May–June 2000). Using the relative reflection as weighting factor, the weighted mean of the surface albedo over Dunhuang Gobi in typical arid region is calculated and its values are 0.255 ± 0.021. After canceling the interference of the buildings, the mean values of the roughness length averaged with logarithm is 0.0019 ± 0.00071 m. After removing the influence of the oasis, the soil wetness factor computed with data under condition of no precipitation is 0.0045. After removing the influence of the precipitation , the mean values of the soil heat capacity over Dunhuang Gobi in typical arid region is 1.12 × 106 J m−3K−1, a bit smaller than the values observed in HEIFE. But the soil heat diffusivity and conductivity are about one of those observed in HEIFE. The soil water content over Dunhuang Gobi in typical synoptic condition is very little and does not exceed 1% basically.

Numerical Simulation of the Critical Scale of Oasis Maintenance and Development in the Arid Regions of Northwest China

Oasis is a special geographic landscape among the vast desert/Gobi in Northwest China (NWC). The surface sensitive heat flux and latent heat flux at Zhangye Oasis during 1 to 11 August 1991 are simulated using the NCAR nonhydrostatic mesoscale model MM5 Version 3. The horizontal grid resolution is set as 1km. By comparing the simulation results with HEIFE observations, it is proved that the model can be used to simulate the surface energy and water mass exchange of arid and semiarid regions in NWC. Based on the above results, the influence of different oasis scales on the local atmospheric field near the ground surface, and the critical scale of oasis maintenance, in NWC are studied dynamically. The following conclusion is obtained: the local thermal circulation between the oasis and the desert/Gobi is formed in the oasis downstream if the oasis scale is larger than 4 km. This local thermal circulation between the oasis and the desert adjacent to the oasis helps to conserve water vapor over the oasis. At the same time, it transfers the abundant water vapor from the oasis into the desert/Gobi near to the oasis to supply relatively plentiful water vapor for desert crops to grow on the fringe of the oasis. So, it is advantageous for oasis extension. However, if the scale of the oasis is smaller than 4 km, it is not easy for the local thermal circulation between the oasis and the desert/Gobi to take shape. This study provides a new standpoint for oasis maintenance and development.

A Parameterization Scheme of Surface Turbulent Momentum and Sensible Heat over the Gobi Underlying Surface

The parameterization of surface turbulent fluxes over the Gobi Desert in arid regions is studied by using rationally screened observational data. First, the characteristics of Monin-Obukhov similarity functions are analyzed and their empirical formulae are fitted. The results show that fitted curves of changes of similarity functions of wind speed and temperature with stability parameter differ little from the typical empirical curves and are within the ranges of scatter of the empirical curves, but their values in the neutral condition arc different from the typical values to some extent. Furthermore, average values of momentum and scalar (sensible heat) roughness lengths as well as changes of scalar roughness length with friction velocity are determined by utilizing the data. It is found that the average values of scalar roughness length are about one order smaller than that of the momentum roughness length and decrease with increasing friction velocity, but they are evidently larger than their theoretically forecasted values.

Characteristics of Hydrologic Transfer between Soil and Atmosphere over Gobi near Oasis at the End of Summer

By utilizing the data observed at Dunhuang during August and September 2000 in the 'Field Experiment on Interaction between Land and Atmosphere in the Arid Region of Northwest China (FEILARNC)', the characteristics of the soil moisture, temperature, and atmospheric humidity are analyzed. It is found that the thickness of the soil temperature active layer is about 5 cm and much thinner than is typical, that not only the atmospheric humidity gradient is often inverted but also the soil moisture gradient in the shallow layer in the Gobi near oasis, that the diurnal variation of soil moisture can be divided into the four stages that are called the wet stage, the losing-water stage, the dry stage, and the attaining-water stage. It is shown that in soil moisture profiles, the depth of the soil moisture active layer is about 10 cm and soil moisture inversion is the main feature in the shallow layer during the wet stage. Such a feature as soil moisture inversion indicates that soil in the shallow layer can inhale moisture from the air through condensation in the nighttime and exhale moisture to the air through evaporation in the daytime. The condensation and evaporation constitute together the full respiration process of moisture on the ground. The formation of soil moisture inversion is related with the state of soil temperature and moisture, the intensity of atmospheric humidity inversion, and the atmospheric thermodynamic stability.

Vegetation patterns and nature reserve construction in an extremely-arid desert in Anxi, NW China's Gansu Province

Anxi County is located in the northwestern part of the Hexi Corridor in Gansu Province and has the sole national level nature reserve of extremely arid desert in China. Phytosociological methods (Braun Blanquet, 1964) are used to classify plant community types in this area. Eleven are distinguished, including six of deserts, four of oases and one transitional type between deserts and oases. Direct gradient analysis(DCA) is employed to correlate the distribution of plant communities to physiogeographic conditions. This study makes clear that water is the most important ecological factor for the distribution of plant species and communities in this area. The effects of water have been demonstrated in different ways. A vegetation gradient from lower altitude to higher altitude in the southern part of the reserve is driven by a precipitation gradient. The effects of the depth of ground water table contribute to the differentiation of vegetation from desert to oasis in the flat area. In a finer scale, the washed gullies have obviously higher species richness and also higher vegetation cover than the surround gobi surfaces, possibly caused by the effects of floods. The vegetation patterns demonstrate that the area of Anxi County is a complete landscape unit. The range of the current nature reserve is not large enough for the purpose of conserving the unique biodiversity in this area.

Folk nomenclature of plants in Cistanche deserticola-associated community in South Gobi, Mongolia

Cistanche deserticola is an important medicinal plant in Mongolia.Despite its significant role in local healing systems,little traditional knowledge had been reported.The present study investigated folk names of C.deserticola and other species of the same community in Umnugobi Province,South Gobi region of Mongolia,based on ethnobotanical approaches.The high correspondence between folk names and scientific names of plant species occurring in Cistanche-associated community shows the scientific meaning of folk nomenclature and classification in Mongolia.The Mongolian and folk names of plants were formed on the basis of observations and understanding of wild plants including their morphology,phenology and traditional uses as well.Results from this study will support the conservation of C.deserticola itself,a rare and endangered plant species listed in the Monglian Red Data Book.Our documentation of folk nomenclature based on 96 plant species in the Cistanche community,as a part of traditional knowledge associated with biodiversity,will be very helpful for making strategy of plant biodiversity conservation in Mongolia.

The slope creep law for a soft rock in an open-pit mine in the Gobi region of Xinjiang, China

The lithology of the strata in the Gobi region of the Xinjiang autonomous region of China is mainly composed of mudstone,silty mudstone,and other soft rocks.Because of the low strength of the rock mass and the serious effects of physical weathering in this area,the slope stability in open-pit mines is poor,and creep deformation and instability can readily occur.Taking the Dananhu No.2 open-pit mine as a typical example,the creep test of a mudstone sample under different stress levels was studied.Then,based on a bottom friction experiment and a FLAC3D numerical simulation,the deformation and failure processes of the slope were analyzed.The stress–displacement curve and the displacement–time curve for the monitoring points were plotted to obtain the relationship between the stress and displacement for the slope of the soft rock.The results showed that the long-term strength of the mudstone was between 8.0 and 8.8 MPa,and that stable creep occurred when the slope was under low stress.The potential failure mode for this type of slope is that the front edge creeps along the weak layer and then a crack is formed at the trailing edge of the slope.When the crack penetrates the weak layer,cutting bedding and bedding sliding occur.The deformation process of the stable creep slope includes an initial deformation stage,an initial creep stage,a constant velocity creep stage and a deceleration creep stage.

Wind tunnel test on the effect of metal net fences on sand flux in a Gobi Desert, China

The Lanzhou-Xinjiang High-speed Railway runs through an expansive windy area in a Gobi Desert, and sand-blocking fences were built to protect the railway from destruction by wind-blown sand. However, the shielding effect of the sand-blocking fence is below the expectation. In this study, effects of metal net fences with porosities of 0.5 and 0.7 were tested in a wind tunnel to determine the effectiveness of the employed two kinds of fences in reducing wind velocity and restraining wind-blown sand. Specifically, the horizontal wind velocities and sediment flux densities above the gravel surface were measured under different free-stream wind velocities for the following conditions： no fence at all, single fence with a porosity of 0.5, single fence with a porosity of 0.7, double fences with a porosity of 0.5, and double fences with a porosity of 0.7. Experimental results showed that the horizontal wind velocity was more significantly decreased by the fence with a porosity of 0.5, especially for the double fences. The horizontal wind velocity decreased approximately 65% at a distance of 3.25 m（i.e., 13 H, where H denotes the fence height） downwind the double fences, and no reverse flow or vortex was observed on the leeward side. The sediment flux density decreased exponentially with height above the gravel surface downwind in all tested fences. The reduction percentage of total sediment flux density was higher for the fence with a porosity of 0.5 than for the fence with a porosity of 0.7, especially for the double fences. Furthermore, the decreasing percentage of total sediment flux density decreased with increasing free-stream wind velocity. The results suggest that compared with metal net fence with a porosity of 0.7, the metal net fence with a porosity of 0.5 is more effective for controlling wind-blown sand in the expansive windy area where the Lanzhou-Xinjiang High-speed Railway runs through.

Effect of stones on the sand saltation threshold during natural sand and dust storms in a stony desert in Tsogt-Ovoo in the Gobi Desert,Mongolia

Non-erodible elements such as stones and vegetation are key to controlling wind erosion and dust emission in drylands.Stony deserts are widely distributed in the Gobi Desert,but the effect of stones on wind erosion and dust emission have not been well studied,except under artificial conditions.In this study,we evaluated the effect of stones on wind erosion and dust emission by measuring the sand saltation threshold in a stony desert in Tsogt-Ovoo in the Gobi Desert,Mongolia,under natural surface conditions during sand and dust storms.We quantified the amount of stones by measuring the roughness density,and determined the threshold friction velocity for sand saltation by measuring wind speed and sand saltation count.Our results showed that the threshold friction velocity increased with the roughness density of stones.In the northern part of the study area,where neither a surface crust nor vegetation was observed,the roughness density of stones was 0.000 in a topographic depression(TD),0.050 on a northern slope(N.SL),and 0.160 on the northern mountain(N.MT).The mean threshold friction velocity values were 0.23,0.41,and 0.57 m/s at the TD,N.SL,and N.MT sites,respectively.In the southern part of the study area,the roughness density values of stones were 0.000 and 0.070-0.320 at the TD and southern slope sites,respectively,and the mean threshold friction velocities were 0.23 and 0.45-0.71 m/s,respectively.We further compared the observed threshold friction velocities with simulated threshold friction velocities using Raupach's theoretical roughness correction and the measured roughness density values,and found that Raupach's roughness correction worked very well in the simulation of threshold friction velocity in the stony desert.This means that the results of our stone measurement can be applied to a numerical dust model.

Overall efficiency of a V-shaped nylon net fence in preventing sand damage to the Mogao Grottoes

A V-shaped nylon net fence installed in 1990 on top of the Mogao Grottoes is shown to be effectively resisting aeolian sand damage to the grottoes. The structure guides and causes deposition of sand from westerly wind （the primary hard wind）, but to some extent hinders the inverse function of easterly wind carrying sand away from the grottoes toward Mount Mingsha. The gobi side by the fence experiences higher wind speed, so that on which are easily formed undercutting pits, and the deposited sands on it generally form double-peak structures due to abundant sand sources. If the earth surface characteristics in gobi areas by both sides of the fence are quite varied, the erosion and deposition features of the accumulating sand section are similar in different seasons; however, if the earth surface characteristics are similar, the features become irregular in different seasons. Sand depositions with long slope feet are formed along the south and north sides of V-shaped nylon net fence. Disrupted by strong westerly wind and northwesterly wind, sand accumulations by north of the fence are in form of single peaks. Although the operation duration of the V-shaped nylon net fence has exceeded its design life （10 years）, our observations indicate that it is still effective in reducing wind-driven sand damages to the Mogao Grottoes, so it should not be withdrawn.

Impacts of late Quaternary environmental change on the long-tailed ground squirrel(Urocitellus undulatus)in Mongolia

Impacts of Quaternary environmental changes on mammal faunas of central Asia remain poorly understood due to a lack of comprehensive phylogeographic sampling for most species.To help address this knowledge gap,we conducted the most extensive molecular analysis to date of the long-tailed ground squirrel （Urocitellus undulatus Pallas 1778） in Mongolia,a country that comprises the southern core of this species＇ range.Drawing on material from recent collaborative field expeditions,we genotyped 128 individuals at two mitochondrial genes （cytochrome b and cytochrome oxidase Ⅰ;1797 bp total）.Phylogenetic inference supports the existence of two deeply divergent infraspecific lineages （corresponding to subspecies U.u.undulatus and U.u.eversmanni）,a result in agreement with previous molecular investigations but discordant with patterns of range-wide craniometric and external phenotypic variation.In the widespread western eversmanni lineage,we recovered geographically-associated clades from the：（a） Khangai,（b） Mongolian Altai,and （c） Govi Altai mountain ranges.Phylogeographic structure in U.u.eversmanni is consistent with an isolation-by-distance model;however,genetic distances are significantly lower than among subspecies,and intra-clade relationships are largely unresolved.The latter patterns,as well as the relatively higher nucleotide polymorphism of populations from the Great Lakes Depression of northwestern Mongolia,suggest a history of range shifts into these lowland areas in response to Pleistocene glaciation and environmental change,followed by upslope movements and mitochondrial lineage sorting with Holocene aridification.Our study illuminates possible historical mechanisms responsible for U.undulatus genetic structure and contributes to a framework for ongoing exploration of mammalian response to past and present climate change in central Asia.