Advances in Grassland Ecosystem Research in Northern Tibetan Plateau

With the aggravation of global change, the response and adaptation of the unique ecosystem in Qinghai-Tibet Plateau to global change have been increasingly concerned by scientific community day by day, which makes the sensitivity and fragility of this ecosystem in response to global change widely recognized by scholars. On the basis of introducing the present research process on the degenerate mechanisim, measures of and approaches to recovery, carbon cycle and primary productivity toward global change, we put forward several propositions on studying the alpine grassland ecosystem in Northern Tibetan Plateau.

Early Cenozoic Mega Thrusting in the Qiangtang Block of the Northern Tibetan Plateau

Recent mapping and seismic survey reveal that intensive compression during the Early Cenozoic in the Qiangtang block of the central Tibetan Plateau formed an extensive complex of thrust sheets that moved relatively southward along several generally north-dipping great thrust systems. Those at the borders of the ~450 km wide block show it overrides the Lhasa block to the south and is overridden by the Hohxil-Bayanhar block to the north. The systems are mostly thin-skinned imbricate thrusts with associated folding. The thrust sheets are chiefly floored by Jurassic limestone that apparently slid over Triassic sandstone and shale, which is locally included, and ramped upward and over Paleocene-Eocene red-beds. Some central thrusts scooped deeper and carried up Paleozoic metamorphic rock, Permian carbonate and granite to form a central uplift that divides the Qiangtang block into two parts. These systems and their associated structures are unconformably overlain by little deformed Late Eocene-Oligocene volcanic rock or capped by Miocene lake beds. A thrust system in the northern part of the block, as well as one in the northern part of the adjacent Lhasa block, dip to the south and appear to be due to secondary adjustments within the thrust sheets. The relative southward displacement across this Early Cenozoic mega thrust system is in excess of 150 km in the Qiangtang block, and the average southward slip-rate of the southern Qiangtang thrusts ranged from 5.6 mm to 7.4 mm/a during the Late Eocene-Oligocene. This Early Cenozoic thrusting ended before the Early Miocene and was followed by Late Cenozoic crustal extension and strike-slip faulting within the Qiangtang block. The revelation and understanding of these thrust systems are very important for the evaluation of the petroleum resources of the region.

Carbon and oxygen isotopic constraints on paleoclimate and paleoelevation of the southwestern Qaidam basin, northern Tibetan Plateau

We investigate the growth of the northern Tibetan Plateau and associated climate change by applying oxygen and carbon isotopic compositions in Cenozoic strata in the southwestern Qaidam basin. The X-ray diffraction and isotopic studies reveal that the carbonate minerals are mainly authigenic and they do not preserve any evidence for detrital carbonate and diagenesis. The isotope data show large fluctuations in the δ^(18)O and δ^(13)C values in the middle-late Eocene, indicating relatively warm and seasonal dry climate.The positive correlation of the δ^(18)O and δ^(13)C values in the Oligocene and the positive shift of the δ^(13)C values from the Eocene to Oligocene suggest that the climate changed to arid in the Oligocene. However,the δ^(18) values show negative shift, which is closely related to the global cooling event. During the Miocene, the δ^(13)C values vary between-2‰ and-4‰, whereas the δ^(18)O values show continuous negative shift. The mean δ^(18) values decrease from-8.5‰, in the early Miocene to-10.0‰, in the late Miocene. The stable isotope-based paleoaltimetry results suggest that the elevation of the southwestern Qaidam basin was approximately 1500 m in the middle-late Eocene and Oligocene. Subsequently, during Miocene the crustal uplift process started and the elevation reached approximately 2000 m in the early Miocene and 2500 m in the late Miocene, which suggests large-scale growth of the northern Tibet Plateau during the Miocene.

Fluxes of CO_2,CH_4 and N_2O from alpine grassland in the Tibetan Plateau

Using stat ic chamber technique, fluxes of CO 2 , CH 4 and N 2 O were measured in the alpine grassland area from July 2000 to July 2001, d eterminations of mean fluxes showed that CO 2 and N 2 O were gene rally released from the soil, while the alpine grassland accounted for a weak CH 4 sink. Fluxes of CO 2 , CH 4 and N 2 O ranged widely. The highest CO 2 emission occurred in August, whereas a lmost 90% of the whole year emission occurred in the growing season. But the variations of CH 4 and N 2 O fluxes did not show any clear patterns over the one-year-experim ent. During a daily variation, the maximum CO 2 emission occurred at 16:00, and then decreased to the minimum emi ssion in the early morning. Daily pattern analyses indicated that the variation in CO 2 fluxes was positively related to air temperatures (R 2 =0.73) and soil temperatures at a depth of 5 cm (R 2 =0.86), whereas daily variations in CH 4 and N 2 O fluxes were poorly explained by soil temperatures and climatic va riables. CO 2 emissions in this area were much lower than other grasslands in plain areas .

Alpine grassland fPAR change over the Northern Tibetan Plateau from 2002 to 2011

In this study, two different methods including Digital Camera and Reference Panel （DCRP） and traditional in situ fPAR observation for measuring the in situ point fPAR of very short alpine grass vegetation were compared, and the Moderate Resolution Imaging Spectroradiometer （MODIS） fPAR products were evaluated and validated by in situ point data on the alpine grassland over the Northern Tibetan Plateau, which is sensitive to climate change and vulnerable to anthropogenic activities. Results showed that the MODIS alpine grassland fPAR product, examined by using DCRP, and traditional in situ fPAR observation had a significant relationship at the spatial and temporal scales. The decadal MODIS fPAR trend analysis showed that, average growing season fPAR increased by 1.2 × 10^-4 per year and in total increased 0.86% from 2002 to 2011 in alpine grassland, when most of the fPAR increments occurred in southeast and center of the Northern Tibetan Plateau, the alpine grassland tended to recover from degradation slightly. However, climatic factors have influenced the various alpine grassland vegetation fPAR over a period of 10 years; precipitation significantly affected the alpine meadow fPAR in the eastern region, whereas temperature considerably influenced the alpine desert steppe fPAR in the west region. These findings suggest that the regional heterogeneity in alpine grassland fPAR results from various environmental factors, except for vegetation characteristics, such as canopy structure and leaf area.

An off-line simulation of land surface processes over the northern Tibetan Plateau

In order to further understand the land surface processes over the northern Tibetan Plateau, this study produced an off-line simulated examination at the Bujiao site on the northern Tibetan Plateau from June 2002 to April 2004, using the Noah Land Surface Model （Noah LSM） and observed data from the CAMP/Tibet experiment. The observed data were neces- sarily corrected and the number of soil layers in the Noah LSM was changed from 4 to 10 to enable this off-line simulation and analysis. The main conclusions are as follows： the Noah LSM performed well on the northern Tibetan Plateau. The simulated net radiation, upward longwave radiation, and upward shortwave radiation demonstrated the same remarkable annual and seasonal variation as the observed data, especially the upward longwave radiation. The simulated soil temperatures were acceptably close to the observed temperatures, especially in the shallow soil layers. The simulated freezing and melting processes were shown to start from the surface soil layer and spread down to the deep soil layers, but they took longer than the observed processes. However, Noah LSM did not adequately simulate the soil moisture. Therefore, additional high-quality, long-term observations of land surface-atmosphere processes over the Tibetan Plateau will be a key factor in proper adiustments of the model parameters in the future.

Effects of freezing intensity on soil solution nitrogen and microbial biomass nitrogen in an alpine grassland ecosystem on the Tibetan Plateau,China

The change of freeze-thaw pattern of the Tibetan Plateau under climate warming is bound to have a profound impact on the soil process of alpine grassland ecosystem;however,the research on the impact of the freeze-thaw action on nitrogen processes of the alpine grassland ecosystem on the Tibetan Plateau has not yet attracted much attention.In this study,the impact of the freezing strength on the soil nitrogen components of alpine grassland on the Tibetan Plateau was studied through laboratory freeze-thaw simulation experiments.The 0–10 cm topsoil was collected from the alpine marsh meadow and alpine meadow in the permafrost region of Beilu River.In the experiment,the soil samples were cultivated at –10℃,–7℃,–5℃,–3℃ and –1℃,respectively for three days and then thawed at 2℃ for one day.The results showed that after the freeze-thaw process,the soil microbial biomass nitrogen significantly decreased while the dissolved organic nitrogen and inorganic nitrogen significantly increased.When the freezing temperature was below –7℃,there was no significant difference between the content of nitrogen components,which implied a change of each nitrogen component might have a response threshold toward the freezing temperature.As the freeze-thaw process can lead to the risk of nitrogen loss in the alpine grassland ecosystem,more attention should be paid to the response of the soil nitrogen cycle of alpine grasslands on the Tibetan Plateau to the freeze-thaw process.

Lithospheric Electrical Structure across the Eastern Segment of the Altyn Tagh Fault on the Northern Margin of the Tibetan Plateau

Project INDEPTH （InterNational DEep Profiling of Tibet and the Himalaya） is an interdisciplinary program designed to develop a better understanding of deep structures and mechanics of the Tibetan Plateau. As a component of magnetoteUuric （MT） work in the 4th phase of the project, MT data were collected along a profile that crosses the eastern segment of the Altyn Tagh fault on the northern margin of the plateau. Time series data processing used robust algorithms to give high quality responses. Dimensionality analysis showed that 2D approach is only valid for the northern section of the profile. Consequently, 2D inversions were only conducted for the northern section, and 3D inversions were conducted on MT data from the whole profile. From the 2D inversion model, the eastern segment of the Altyn Tagh fault only appears as a crustal structure, which suggests accommodation of strike slip motion along the Altyn Tagh fault by thrusting within the Qilian block. A large-scale off-proffie conductor within the mid-lower crust of the Qilian block was revealed from the 3D inversion model, which is probably correlated with the North Qaidam thrust belt. Furthermore, the unconnected conductors from the 3D inversion model indicate that deformations in the study area are generally localized.

Study on potential evapotranspiration and wet-dry condition in the seasonal frozen soil region of northern Tibetan Plateau

This study was based on the CEOP/CAMP-Tibet observed data at AWS （Automatic Weather Station） of MS3478 in the seasonal frozen soil region of northern Tibetan Plateau from March 2007 to February 2008. The variation characteristics of PE （potential evapotransph＇ation） were analyzed based on the Penman-Monteith method recommended by FAO （the Food and Agriculture Organization of the United Na- lions）. The contributions of dynamic, thermal and water factors to PE were discussed, and the wet-dry condition of the plateau region was further studied. The results indicated that daily PE was between 0.52 mm and 6.46 mm for the whole year. Monthly PE was over 107 mm from May to September, but decreased to less than 41 mm from November to February. Annual PE was 1,037.8mm. In the summer, thermal PE was significantly more than dynamic PE, but conversely in the winter. Annual variation of thermal PE was of sine wave pattern. In addition, drought and semi-drought climate lasted for a long time while semi-humid climate was short. The effect of water and dynamic factors on PE varied considerably with the seasons. Annual variation of thermal PE was of sine wave pattern.

Food-Chain Model of Grassland Degradation and Its Restoration Process in Northern Tibet Plateau:A Case Study in Nierong County

Based on the model of grassland climate ecological productivity, the process of grassland degradation and its restoration mechanism in northern Tibetan Plateau were discussed by the model of food-chain in which the environmental and human factors were corrected. The results of case study in Nierong County showed that： ① the climate trend of becoming warmer, more droughts and gales were conflicted with the restoration of grassland degradation, even under level of perfect management the climate ecological productivity was declined from 89. 3 kg/m^2 of 1983 to 71.8 kg/m^2 of 2003; ② from 1983 to 2003, the population increased fast, while the variation of livestock on hand was little, and the cost of its maintaining is rapid grassland degradation; ③ on the present condition of overgrazing, the livestock on hand can be maintained on the level of theoretical carrying capacity in 2033 by applying the mechanism of food-chain in grassland ecological system controlled with expected coefficients, so that to realize the policy of determining the quantity of livestock according to grass growth.

CO_2 processes in an alpine grassland ecosystem on the Tibetan Plateau

In this paper, the CO2 concentrations profile from 1.5 m depth in soil to 32 m height in atmosphere were measured from July 2000 to July 2001 in an alpine grassland ecosystem located in the permafrost area on the Tibetan Plateau, which revealed that CO2 concentrations varied greatly during this study period. Mean concentrations during the whole experiment in the atmosphere were absolutely lower than the CO2 concentrations in soil, which resulted in CO2 emissions from the alpine steppe soil to the atmosphere. The highest CO2 concentration was found at a depth of 1.5 m in soil while the lowest CO2 concentration occurred in the atmosphere. Mean CO2 concentrations in soil generally increased with depth. This was the compositive influence of the increasing soil moistures and decreasing soil pH, which induced the increasing biological activities with depth. Temporally, the CO2 concentrations at different layers in air remained a more steady state because of the atmospheric turbulent milking. During the seasonal variations, CO2 concentrations at surface soil interface showed symmetrical patterns, with the lowest accumulation of CO2 occurring in the late winter and the highest CO2 concentration in the growing seasons.

Air Temperature Estimation with MODIS Data over the Northern Tibetan Plateau

Time series of MODIS land surface temperature(Ts) and normalized difference vegetation index(NDVI) products,combined with digital elevation model(DEM) and meteorological data from 2001 to 2012,were used to map the spatial distribution of monthly mean air temperature over the Northern Tibetan Plateau(NTP). A time series analysis and a regression analysis of monthly mean land surface temperature(Ts) and air temperature(Ta) were conducted using ordinary linear regression(OLR) and geographical weighted regression(GWR). The analyses showed that GWR,which considers MODIS Ts,NDVI and elevation as independent variables,yielded much better results [RAdj2> 0.79; root-mean-square error(RMSE) =0.51℃–1.12℃] associated with estimating Tacompared to those from OLR(RAdj2= 0.40-0.78; RMSE = 1.60℃–4.38℃).In addition,some characteristics of the spatial distribution of monthly Taand the difference between the surface and air temperature(Td) are as follows. According to the analysis of the 0℃ and 10℃ isothermals,Tavalues over the NTP at elevations of 4000–5000 m were greater than 10℃ in the summer(from May to October),and Tavalues at an elevation of3200 m dropped below 0℃ in the winter(from November to April). Taexhibited an increasing trend from northwest to southeast. Except in the southeastern area of the NTP,T d values in other areas were all larger than 0℃ in the winter.

The Beila Ophiolite from the Bangong-Nujiang Suture Zone, Northern Tibetan Plateau

The Beila ophiolite is located in the middle part of the Bangong-Nujiang suture zone,northern Tibetan plateau.It is a complete ophiolite suite,and plays a key role in understanding the evolution of the Bangong-Nujiang suture zone,as well as the Meso-Tethys Ocean.The Beila ophiolite was composed of peridotite,serpentinite,gabbro,pillow basalt,and minor rodingite.Peridotites comprisemainlymedium–tocoarse–grained serpentinized harzburgites and minor plagioclase-bearing lherzolites and dunites.There are some felsic-ultramafic dykes within the peridotite and they are mainlypegmatoidal pyroxenites,coarse to fine-grained gabbros,and diabases.Gabbros included isotropic and cumulate gabbros,and they commonly contain minor pegmatoidal gabbros veins.Pillow basalts and basaltic andesites overlaid on the margin of the serpentinized peridotites.Rodingite occurs as lenses and/or dykes within the host serpentinized peridotites.Zircon SHRIMP U–Pb dating for two rodingite samples yielded the ages ranging from172 to 164 Ma.Whole-rock geochemical and zircon Hf isotopic data show that the Beila ophiolite shows SSZ-type ophiolite affinity.Finally,we suggest that the Beila ophiolite was generated in an initial subduction process at the middle Jurassic(164–172 Ma).

Sources of uncertainty in exploring rangeland phenology： A case study in an alpine meadow on the central Tibetan Plateau

Global climate change has been found to substantially influence the phenology of rangeland,especially on the Tibetan Plateau. However, there is considerable controversy about the trends and causes of rangeland phenology owing to different phenological exploration methods and lack of ground validation. Little is known about the uncertainty in the exploration accuracy of vegetation phenology.Therefore, in this study, we selected a typical alpine rangeland near Damxung national meteorological station as a case study on central Tibetan Plateau, and identified several important sources influencing phenology to better understand their effects on phenological exploration. We found man-made land use was not easily distinguished from natural rangelands, and therefore this may confound phenological response to climate change in the rangeland. Change trends of phenology explored by four methods were similar, but ratio threshold method(RTM) was more suitable for exploring vegetation phenology in terms of the beginning of growing season(BGS) and end of growing season(EGS). However, some adjustments are needed when RTM is used in extreme drought years. MODIS NDVI/EVI dataset was most suitable for exploring vegetation phenology of BGS and EGS. The discrimination capacities of vegetation phenology declined with decreasing resolution of remote sensing images from MODIS to GIMMS AVHRR datasets. Additionally, distinct trends of phenological change rates were indicated in different terrain conditions, with advance of growing season in high altitudes but delay of season in lower altitudes. Therefore, it was necessary to eliminate interference of complex terrain and man-made land use to ensure the representativeness of natural vegetation. Moreover, selecting the appropriate method to explore rangelands and fully considering the impact of topography are important to accurately analyze the effects of climate change on vegetation phenology.

Plant community dynamics during the growing season of typical ecosystems on the Tibetan Plateau

The relationships between vegetation and environmental factors have always been a core concern of ecologists.The dynamic characteristics of plant communities during the growing season can directly reflect these relation-ships,so we examined this issue for three typical ecosystems on the Tibetan Plateau.During the growing season,the dominant species remained stable while non-dominant species changed significantly in the alpine meadow and alpine steppe and a mono-dominant community was found in the temperate desert shrub.Due to the seasonal variations of temperature and soil water content,plant species diversity varied significantly during the growing season.Patrick richness,Pielou evenness and Simpson diversity indices differed significantly in the alpine meadow and alpine steppe.The total biomass of these three ecosystems was the largest during the middle growing season.Biomass was greater in the alpine meadow than the alpine steeps or temperature desert.The root-to-shoot ratio was the lowest during the middle growing season for the alpine meadow and alpine steppe and largest during the early growing season for temperate desert shrub.RDA showed the belowground and total biomass were greatly affected by soil physicochemical factors.Multiple linear stepwise regression showed the above ground biomass was greatly affected by relative atmospheric humidity and belowground and total biomass were greatly affected by soil organic carbon,total nitrogen at 0-20 cm soil depth and pH at 10-20 cm soil depth.These findings pro-vide insights into understanding the relationships between vegetation and environmental factors and promote the sustainable utilization of local grasslands on the Tibetan Plateau.

Mechanism of plant-soil feedback in a degraded alpine grassland on the Tibetan Plateau

Although biotic and abiotic factors have been confirmed to be critical factors that affect community dynamics,their interactive effects have yet to be fully considered in grassland degradation.Herein,we tested how soil nutrients and microbes regulated plant-soil feedback(PSF)in a degraded alpine grassland.Our results indicated that soil total carbon(STC;from 17.66 to 12.55 g/kg)and total nitrogen(STN;from 3.16 to 2.74 g/kg)exhibited significant(P<0.05)decrease from non-degraded(ND)to severely degraded(SD).Despite higher nutrients in ND soil generating significantly(P<0.05)positive PSF(0.52)on monocots growth when the soil was sterilized,a high proportion of pathogens(36%)in ND non-sterilized soil resulted in a strong negative PSF on monocots.In contrast,the higher phenotypic plasticity of dicots coupled with a higher abundance of mutualists and saprophytes(70%)strongly promoted their survival and growth in SD with infertile soil.Our findings identified a novel mechanism that there was a functional group shift from monocots with higher vulnerability to soil pathogens in the ND fertile soil to dicots with higher dependence on nutritional mutualists in the degraded infertile soil.The emerging irreversible eco-evolutionary in PSF after degradation might cause a predicament for the restoration of degraded grassland.

Deep structure of the Kumkol basin in the northern Tibetan Plateau and its resource environmental implications

The Kumkol basin is located in the northern Tibetan Plateau and is a closed plateau basin with an average altitude of>4000 m and an area of nearly 20000 km^(2). Its boundaries are limited by the Altyn Tagh fault, East Kunlun orogen and Qimantag orogen. Studying the deep structure of the Kumkol basin reveals 2 significant implications:(1) the basin has developed a large thickness of >7000 m Cenozoic continental sediments, recording the uplift history of the northern Tibetan Plateau, and(2)preliminary work indicates that the basin is likely to have oil and gas prospects. However, owing to the adverse natural conditions of the area and the strong tectonic activity in the Cenozoic, the latter of which was not conducive to hydrocarbon preservation,only regional geological mapping and petroleum exploration route surveys have been carried out, and there is no consensus on strata, structure and tectonic evolution. From 2021 to 2022, a deep seismic reflection profile implemented by the Second Tibetan Plateau Scientific Expedition and Research(STEP) project was the first high-resolution geophysical survey across the Kumkol basin. This study uses seismic reflection migration profiles, first-arrival wave tomographic imaging and previous research results to analyze the deep structure of the basin. The final merged model contains many features of tectonic and resource significance:(1) The Kumkol basin is ~90 km wide from north to south, with a basement depth of >9000 m. The main component is the Cenozoic continental deposits, which are divided into two major parts: the southern composite basin and the northern faulted basin. Owing to the later compression, the southern composite basin experienced significant deformation, but most parts still preserved their original sedimentary formations.(2) The structural deformation characteristics of the basin reveal a two-stage tectonic evolution process of the northern Tibetan Plateau in the Cenozoic: from the Oligocene to the Pliocene, the main mechanism was vertical differential uplift and subsidence, and after the Pliocene, it transformed to north-south compression and shortened deformation.(3) The strata, formation time, and source-reservoir-cap conditions of the Kumkol basin are similar to those of the Qaidam basin. If a breakthrough can be achieved, it is expected to expand the production capacity of the oil field in the Qaidam basin with a low-cost investment. Thus, further exploration is recommended.

Transformation and mechanisms of climate wet/dry change on the northern Tibetan Plateau under global warming:A perspective from paleoclimatology

The northern Tibetan Plateau is a climatically sensitive zone influenced by monsoon and westerly winds.In summer,water vapor transport can reach Qinghai Lake and the eastern section of the Qilian Mountains;in winter,westerly winds mainly control the climate.This article compares the wet/dry changes in the region during the mid-Holocene(MH)warm period,the medieval climate anomaly(MCA),the current warm period(CWP),and the future warm period from the perspective of paleoclimate.We found that the MH warm period was mainly affected by the orbit-controlled East Asian summer monsoon,and the region showed warm and humid climate characteristics.The MCA was mainly controlled by solar radiation,and there was a warm and dry phenomenon.The CWP and the future warm period are mainly controlled by the rise in temperature caused by the increase in greenhouse gases,and the climate is becoming more arid.The wet/dry patterns in the CWP and the future warm period in the next century on the northern Tibetan Plateau are similar to those in the MCA.Continued warming will lead to the expansion of the westerly belt and a gradually humid climate.The future wet/dry changes will be more similar to the MH warm period.

Effects of Temperature and Nitrogen Input on Nitrogen Mineralization in Alpine Soils on Tibetan Plateau

[Objective] The study aimed to investigate the effects of temperature and nitrogen input on nitrogen mineralization in alpine soils on the Tibetan Plateau.[Method]An incubation experiment in the laboratory was conducted using three alpine soils.These soils were collected from the top 10 cm depth in three grassland types：alpine meadow in Haibei,alpine steppe in Naqu and alpine wetland in Dangxiong.[Result] Temperature significantly affected nitrogen mineralization in alpine soils of three grassland types.The mineralization rate in alpine steppe soil rose with the rise of temperature,while the mineralization rate in the alpine meadow soil and alpine wetland soil decreased with the rise of temperature.Nitrogen input had no significant effect on nitrogen mineralization in the alpine meadow soil,but significantly increased nitrogen mineralization in the alpine steppe soil and the alpine wetland soil.Grassland types significantly affected nitrogen mineralization in alpine grasslands.[Conclusion] The effects of temperature and nitrogen input on nitrogen mineralization in alpine soils on the Tibetan Plateau were significant.And those different effects depended on different types of grassland.

Climatic and geographic factors affect ecosystem multifunctionality through biodiversity in the Tibetan alpine grasslands

Ecosystem multifunctionality(EMF), the simultaneous provision of multiple ecosystem functions, is often affected by biodiversity and environmental factors. We know little about how the interactions between biodiversity and environmental factors affect EMF. In this case study, a structural equation model was used to clarify climatic and geographic pathways that affect EMF by varying biodiversity in the Tibetan alpine grasslands. In addition to services related to carbon, nitrogen, and water cycling, forage supply, which is related to plantproductivity and palatability, was included in the EMF index. The results showed that 72% of the variation in EMF could be explained by biodiversity and other environmental factors. The ratio of palatable richness to all species richness explained 8.3% of the EMF variation. We found that air temperature, elevation, and latitude all affected EMF, but in different ways. Air temperature and elevation impacted the aboveground parts of the ecosystem, which included plant height, aboveground biomass, richness of palatable species, and ratio of palatable richness to all species richness. Latitude affected EMF by varying both aboveground and belowground parts of the ecosystem, which included palatable speciesrichness and belowground biomass. Our results indicated that there are still uncertainties in the biodiversity–EMF relationships related to the variable components of EMF, and climatic and geographic factors. Clarification of pathways that affect EMF using structural equation modeling techniques could elucidate the mechanisms by which environmental changes affect EMF.