Method to measure tree-ring width,density,elemental composition,and stable carbon and oxygen isotopes using one sample

Tree-ring width(RW),density,elemental com-position,and stable carbon and oxygen isotope(δ^(13)C,δ^(18)O)are widely used as proxies to assess climate change,ecology,and environmental pollution;however,a specific pretreat-ment has been needed for each proxy.Here,we developed a method by which each proxy can be measured in the same sample.First,the sample is polished for ring width meas-urement.After obtaining the ring width data,the sample is cut to form a 1-mm-thick wood plate.The sample is then mounted in a vertical sample holder,and gradually scanned by an X-ray beam.Simultaneously,the count rates of the fluorescent photons of elements(for chemical characteriza-tion)and a radiographic grayscale image(for wood density)are obtained,i.e.the density and the element content are obtained.Then,cellulose is isolated from the 1-mm wood plate by removal of lignin,and hemicellulose.After producing this cellulose plate,cellulose subsamples are separated by knife under the microscope for inter-annual and intra-annual stable carbon and oxygen isotope(δ^(13)C,δ^(18)O)analysis.Based on this method,RW,density,elemental composition,δ^(13)C,and δ^(18)O can be measured from the same sample,which reduces sample amount and treatment time,and is helpful for multi-proxy comparison and combination research.

Enhancing forest insect outbreak detection by integrating tree-ring and climate variables

Annual tree rings are widely recognized as valuable tools for quantifying and reconstructing historical forest disturbances.However,the influence of climate can complicate the detection of disturbance signals,leading to limited accuracy in existing methods.In this study,we propose a random under-sampling boosting(RUB)classifier that integrates both tree-ring and climate variables to enhance the detection of forest insect outbreaks.The study focused on 32 sites in Alberta,Canada,which documented insect outbreaks from 1939 to 2010.Through thorough feature engineering,model development,and tenfold cross-validation,multiple machine learning(ML)models were constructed.These models used ring width indices(RWIs)and climate variables within an 11-year window as input features,with outbreak and non-outbreak occurrences as the corresponding output variables.Our results reveal that the RUB model consistently demonstrated superior overall performance and stability,with an accuracy of 88.1%,which surpassed that of the other ML models.In addition,the relative importance of the feature variables followed the order RWIs>mean maximum temperature(Tmax)from May to July>mean total precipita-tion(Pmean)in July>mean minimum temperature(Tmin)in October.More importantly,the dfoliatR(an R package for detecting insect defoliation)and curve intervention detec-tion methods were inferior to the RUB model.Our findings underscore that integrating tree-ring width and climate vari-ables as predictors in machine learning offers a promising avenue for enhancing the accuracy of detecting forest insect outbreaks.

A Study of the Distribution and Accumulation of Major Nutrient Elements in Different Densities of Eucalyptus grandis

This study was conducted to investigate the distribution and accumulation of major elements of 5-year-old Eucalyptus grandis forests with the densities of 556 trees/hm2 and 1 667 trees/hm2. The results showed that： （1） The distribution of ma- jor nutrient element contents in various organs of the E. grandis with the 2 densi- ties was nearly the same, and the ranked order was leaf〉tree trunk〉branch or branchlet〉fruit〉bark. （2） From the distribution of nutrient elements and organic mat- ter in the trunk, P, K, Mg and Ca contents presented basically a trend of increasing with the trunk height increasing. And the organic matter content increased with the trunk height increasing in the low density stand, while the result was opposite in the high density stand. Furthermore, N content increased with the trunk height increas- ing in the low density stand, but the content in the high density stand exhibited a trend of decreasing at first, increasing then and increasing at the last. （3） From the distribution of biomass, the ranked order was tree trunk〉bark or branch〉branchlet〉 leaf〉fruit. And （4） The accumulated amounts of major nutrient elements in various organs under the 2 stand densities was nearly same. Ca content was the highest in the 2 kinds of plantations, followed by N, K, Mg and P contents. Ca content was 601.78 kg/hm2 in the low density stand and 1 204. 43 kg/hm2 in the high density stand.

Effect of Different Sowing Dates and Densities on Individual Morphological Development of Super Short-season Insect-resistant Cotton

[Objective] The paper was to explore the effect of different sowing dates and densities on individual morphological development of super short-season insect-resistant cotton,confirm their effects on vegetative and reproductive growth of cotton,so as to provide theoretical and practical guidance for sowing date and density management of cotton planting in Jidong cotton growing region in Yellow River Basin.[Method] With super short-season insect-resistant cotton＂546＂as materials,the effects of different sowing dates（sowing dateⅠ：May 20;sowing dateⅡ：June 2;sowing date Ⅲ：June 14）and densities（low density：120 000 plants/hm2;middle density：150 000 plants/hm2;high density：180 000 plants/hm2）on individual morphological development of super short-season insect-resistant cotton were explored.[Result] Different sowing dates and density treatments significantly affected the individual morphological development of super short-season insect-resistant cotton＂546＂.The effectiveness of sowing date was higher than the effectiveness of density,and the effectiveness of sowing date on development of number of individual fruit branches was higher than that on plant height and stem diameter.[Conclusion] The regulation of sowing date and density during the cultivation process of super short-season insect-resistant cotton ＂546＂ in Jidong cotton growing region in Yellow River Basin could effectively promote vegetative and reproductive growth of cotton,strengthening its production base.

Effect of Chemical Regulation on Maize in Different Cultivation Densities

[Objective] This study aimed at investigating the effect of chemical regulation on maize cultivated in different densities.[Method] By taking maize hybrids Xundan 20 and Zhongnongda 4 as materials,the effect of chemical regulation on the yield,yield component factors and plant morphology of the two hybrids planting in Northern and Northeastern China were studied,respectively.[Result] Chemical Yuhuangjin could control maize stalks,lower the plant height and the ear height,make internodes below the spike more coarse and enhance its lodging-resistance.In addition,chemical regulation could increase the optimum cultivation density by 7 500 plants/hm2 more than traditional cultivation.[Conclusion] This paper provided theoretical reference for cultivating high-yield maize.

Effect of Fertilizing Level and Planting Densities on Yield and Nitrogen Utilization in Maize

Elite maize hybrid Guidan0810 was selected as material, and the effects of fertilizing level and planting densities on yield and nitrogen utilization were dis- cussed in the study. In field experiments as per double-cropping system, 4 main plots （fertilization levels） and 6 subplots （planting densities） were set in a split plot design. The results suggested that yield had close relationship with fertilization levels and planting densities. Different fertilization levels and planting densities significantly affected yield. With the increase of nitrogen fertilization, nitrogen use efficiency, nitrogen agronomic efficiency and nitrogen physiological efficiency declined. Under the same fertilization level, nitrogen use efficiency, nitrogen agronomic efficiency and nitrogen physiological efficiency grew a little with the increase of planting density, so nitrogen efficiency could be improved by regulating planting density. The results also showed that A2 （including N 225.0 kg/hm2, P205 75.0 kg/hm^2, K20 187.5 kg/hm^2） matching to B3 （52 500 plants/hm^2） or B4（60 000 plants/hm^2） was a better design, which could obtain a higher yield in the range of 7 913.2-8 207.8 kg/hm2, respectively.

Effects of Different Transplanting Densities on Yield and Yield Components of Japonica Longke 16

[Objective] This study was conducted to improve unit yield of rice in Longyang District and find the relationship between transplanting density and rice yield. [Method] Under the condition of field plot test, effects of cultivation density on yield and yield components were investigated with japonica rice Longke 16 as a test material. [Result] Under a planting density of 2.7-3.3×10^5 clusters/hm^2, equivalent to a density of seedlings with tillers of 5.4-6.6×10^5/hm^2, the yield reached 11 250-12 000 kg/hm^2. [Conclusion] Longke 16 is a japonica rice variety with good comprehensive characters, and rational dense planting and accurate quantitative cultivation could give full play to its high yield potential.

Analysis on Differences in Photosynthetic Characteristics of Korla Fragrant Pear among Different Planting Densities

[Objective] The aim of this study was to explore the differences in photosynthetic characteristics of Korla fragrant pear among different planting densities,providing a basis for the improvement of fruit yield and quality of Korla fragrant pear in production. [Method] The net photosynthetic rates, photoresponse curves and CO2 response curves of Korla fragrant pear under six different planting densities were determined using Li-6400 XT portable photosynthesis system. [Result] There were significant differences in photosynthetic parameters of Korla fragrant pear among six different planting densities. The net photosynthetic rate of Korla fragrant pear was highest under the planting density of 4.0 m × 6.0 m. At the planting density of4.0 m×6.0 m, when the light intensity reached 1 800 μmol/（m^2·s）, the net photosynthetic rate of Korla fragrant pear reached the peak [19.326 μmol/（m^2·s）], and the apparent quantum yield, carboxylation efficiency and dark respiration rate all reached the maximum values; at the planting density of 6.0 m ×7.0 m, Korla fragrant pear showed lower light compensation point and CO2 compensation point, but higher light saturation point. At the planting density of 3.0 m ×5.0 m, Korla fragrant pear had higher light saturation point. [Conclusion] Among the six different planting densities,Korla fragrant pear with planting density of 4.0 m ×6.0 m showed the highest net photosynthetic rate. It suggested that Korla fragrant pear, at the planting density of4.0 m ×6.0 m, had stronger utilization capacity for low light and low-concentration CO2. Therefore, the Korla fragrant pear with the planting density of 4.0 m ×6.0 m has the highest photosynthesis efficiency, and the planting density of 4.0 m ×6.0 m is the most ideal planting density for Korla fragrant pear.

Effects of Climatic Factors on Tree-ring Maximum Latewood Density of Picea schrenkiana in Xinjiang, China

Based on two tree-ring maximum latewood density （MXD） chronologies of Picea schrenkiana from the Manas River Basin, Xinjiang, the response characteristics of MXD to climate variation was discussed. Correlation analysis between MXD chronologies and instrumental records from Shihezi meteorological station showed that each chronology was significantly and positively correlated with the maximum monthly average temperature in July-August, and especially, the regional chronology （RC） was the most highly correlated variable （r=0.54, P〈0.001）. Afterwards, the maximum average temperature in July-August was reconstructed using RC. Comparison among reconstructed temperature, observed values, and the drought index （Is） confirmed that precipitation would affect MXD when the absolute value of Is is greater than 1.5σ （｜Is｜ 〉 2.5 during 1953-2008） or near to 1.5a over a 2-3 year period. The response characteristics are related to the semiarid climate of the study area. In dry years, lack of precipitation would limit the thickening of latewood cell walls and, as a result, impact MXD. Therefore, compared with relatively humid regions, the response of tree-ring MXD to air temperature similarly would be influenced by extreme moisture conditions in semiarid areas, and MXD, as a temperature proxy, should be used prudently on a limited scale.

QTL analysis of ear leaf traits in maize(Zea mays L.)under different planting densities

Modem maize varieties have become more productive than ever, owing largely to increased tolerance of high plant density. However, the genetics of ear leaf traits under different densities remains poorly understood. In this study, Zhongdan 909 recombinant inbred lines(RILs) derived from a cross between Z58 and HD568 were genotyped for 3072 single-nucleotide polymorphisms(SNPs), and phenotyped for leaf length(LL), leaf width(LW), and leaf angle(LA) of the uppermost ear leaf under three planting densities(52,500,67,500, and 82,500 plants ha-1, respectively). A genetic map was then constructed using1358 high-quality SNPs. The total length of the linkage map was 1985.2 cM and the average interval between adjacent markers 1.46 cM. With increasing density, LL and LW decreased from 63.68 to 63.02 cm and from 8.56 to 8.21 cm, respectively, while LA increased from19.42° to 19.66°. All three traits had high heritabilities, of 0.75, 0.78, and 0.84, respectively.Using inclusive composite interval mapping, 23, 25, and 17 quantitative trait loci(QTL) were detected for LL, LW, and LA, respectively. Of these, 35 were simultaneously detected under two or three plant densities, while 30 were detected under only one. Sixty-five individual QTL explained 2.41% to 16.53% of phenotypic variation, while eight accounted for >10%.These findings will help us understand the genetic basis of leaf traits in maize as well as the response of maize to increased plant density.

January to August temperature variability since 1776 inferred from tree-ring width of Pinus tabulaeformis in Helan Mountain

Tree-ring standardized chronologies are developed by 78 cores collected from the eastern and western Helan Mountain. Statistical analysis shows that both the STD and RES chronologies correlate negatively with the temperature of different periods of early half year, especially with January to August mean （JA） temperature, which means that JA temperature is one of the predominant limiting factors of tree growth in the Helan Mountain. Based on this analysis, we reconstructed JA temperature, and the explained variance is 43.3% （F=21.422, p〈0.001 ）. The comparatively high temperature periods in the reconstruction were： 1805-1818 1828-1857, 1899-1907, 1919-1931 and 1968-1995; and the comparatively low temperatu re periods happened in 1858-1872, 1883-1895 and 1935-1953. Ten-year moving average curve shows three slow uplifting trends： 1766-1853, 1862-1931 and 1944-1995. Each temperature increase was followed by a sudden temperature decrease about 10 years, that is to say, the JA temperature in the Helan Mountain is characterized by slow increase and sudden decrease. The 70- and 10.77-year periodicities detected in the temperature series correspond to the Gleissberg （80-year） and Schwabe （11-year） periodicities of solar activity respectively, the 2.11-2.62 years cycles are considered to be influenced by QBO （Quasi-Biennial-Oscillation） and the local environmental change.

Densities and viscosities for ionic liquids [BMIM][BF4] and [BMIM][Cl] and their binary mixtures at various temperatures and atmospheric pressure

The density and viscosity of 1-butyl-3-methylimidazolium tetrafluoroborate[BMIM][BF4]and 1-butyl-3-methylimidazolium chloride[BMIM][Cl]and their binary mixtures within the temperatures from 303.15 K to323.15 K and at ambient pressure were determined in this work.The temperature dependences of density and viscosity were satisfactorily described with the linear model and the Vogel-Tammann-Fulcher type equation,respectively.The molar volume and viscosity of binary IL mixtures were predicted through ideal mixing rules showing that almost null deviations for IL mixtures were observed and their mixing was remarkably close to linear ideal behavior in the molar volumes,while comparatively large errors in viscosity occurred.Additionally,the molar volume of the investigated pure ILs and their mixtures could well be predicted by a predictive model presented by Valderrama et al.(Fluid Phase Equilib.,275(2009)145).

Reconstruction of hydrological changes based on tree-ring data of the Haba River, northwestern China

Reconstructing the hydrological change based on dendrohydrological data has important implications for understanding the dynamic distribution and evolution pattern of a given river. The widespread, long-living coniferous forests on the Altay Mountains provide a good example for carrying out the dendrohydrological studies. In this study, a regional composite tree-ring width chronology developed by Lariat sibirica Ledeb. and Picea obovata Ledeb. was used to reconstruct a 301-year annual （from preceding July to succeeding June） streamflow for the Haba River, which originates in the southern Altay Mountains, Xinjiang, China. Results indicated that the reconstructed streamflow series and the observations were fitting well, and explained 47.5% of the variation in the observed streamflow of 1957-2008. Moreover, floods and droughts in 1949-2000 were precisely captured by the streamflow reconstruction. Based on the frequencies of the wettest/driest years and decades, we identified the 19th century as the century with the largest occurrence of hydrological fluctuations for the last 300 years. After applying a 21-year moving average, we found five wet （1724-1758, 1780-1810, 1822-1853, 1931-1967, and 1986-2004） and four dry （1759-1779, 1811-1821, 1854-1930, and 1968-1985） periods in the streamflow reconstruction. Furthermore, four periods （1770-1796, 1816-1836, 1884-1949, and 1973-1997） identified by the streamflow series had an obvious increasing trend. The increasing trend of streamflow since the 1970s was the biggest in the last 300 years and coincided with the recent warming-wetting trend in northwestern China. A significant correlation between streamflow and precipitation in the Altay Mountains indicated that the streamflow reconstruction contained not only local, but also broad-scale, hydro-climatic signals. The 24-year, 12-year, and 2.2-4.5-year cycles of the reconstruction revealed that the streamflow variability of the Haba River may be influenced by solar activity and the atmosphere-ocean system.

QTL Mapping for Stalk Related Traits in Maize (Zea mays L.) Under Different Densities

Stalk related traits, comprising plant height （PH）, ear height （EH）, internode number （IN）, average internode length （ALL）, stalk diameter （SD）, and ear height coefficient （EHC）, are significantly correlated with yield, density tolerance, and lodging resistance in maize. To investigate the genetic basis for stalk related traits, a doubled haploid （DH） population derived from a cross between NX531 and NX110 were evauluated under two densities over 2 yr. The additive quantitative trait loci （QTLs）, epistatic QTLs were detected using inclusive composite interval mapping and QTL-by-environment interaction were detected using mixed linear model. Differences between the two densities were significant for the six traits in the DH population. A linkage map that covered 1 721.19 cM with an average interval of 10.50 cM was constructed with 164 simple sequence repeat （SSR）. Two, two, seven, six, two, and eight additive QTLs for PH, IN, AIL, EH, SD, and EHC, respectively. The extend of their contribution to penotypic variation ranged from 10.10 to 31.93%. Seven QTLs were indentified simultaneously under both densities. One pair, two pairs and one pair of epistatic effects were detected for AIL, SD and EHC, respectively. No epistatic effects were detected for PH, EH, and IN. Nineteen QTLs with environment interactions were detected and their contribution to phenotypic variation ranged from 0.43 to 1.89%. Some QTLs were stably detected under different environments or genetic backgrounds comparing with previous studies. These QTLs could be useful for genetic improvement of stalk related traits in maize breeding.

Classification of densities and characteristics of curve of population centers in China by GIS

In this paper, with the spatial analysis functions in ArcGIS and the county-level census data of 2000 in China, the population density map was divided and shown by classes, meanwhile, the map system of population distribution and a curve of population centers were formed; in accordance with the geographical proximity principle, the classes of population densities were reclassified and a population density map was obtained which had the spatial clustering characteristic. The multi-layer superposition based on the population density classification shows that the population densities become denser from the Northwest to the Southeast; the multi-layer clustering phenomenon of the Chinese population distribution is obvious, the populations have a water-based characteristic gathering towards the rivers and coastlines. The curve of population centers shows the population densities transit from the high density region to the low one on the whole, while in low-density areas there are relatively dense areas, and in high-density areas there are relatively sparse areas. The reclassification research on the population density map based on the curve of population centers shows that the Chinese population densities can be divided into 9 classes, hereby, the geographical distribution of Chinese population can be divided into 9 type regions： the concentration core zone, high concentration zone, moderate concentration zone, low concentration zone, general transitional zone, relatively sparse area, absolute sparse area, extreme sparse area, and basic no-man＇s land. More than 3/4 of the population of China is concentrated in less than 1/5 of the land area, and more than half of the land area is inhabited by less than 2% of the population, the result reveals a better space law of China＇s population distribution.

Dynamically Exploring the QCD Matter at Finite Temperatures and Densities:A Short Review

We provide a concise review on recent theory advancements towards full-fledged(3+1)D dynamical descriptions of relativistic nuclear collisions at finite baryon density.Heavy-ion collisions at different collision energies produce strongly coupled matter and probe the QCD phase transition at the crossover,critical point,and first-order phase transition regions.Dynamical frameworks provide a quantitative tool to extract properties of hot QCD matter and map fireballs to the QCD phase diagram.Outstanding challenges are highlighted when confronting current theoretical frameworks with current and forthcoming experimental measurements from the RHIC beam energy scan programs.

Litter decomposition and the degradation of recalcitrant components in Pinus massoniana plantations with various canopy densities

To understand the decomposition characteristics of Pinus massoniana foliar litter and the degradation of its refractory compounds in plantations under five canopy densities,a litter bag experiment over a decomposition time of 392 days was carried out.The results show that canopy density significantly affected decomposition rates of litter and degradation rate of lignin and cellulose.Litter decomposition rates decreased significantly with decreasing canopy density.Both lignin and cellulose degradation rates were lower with canopy densities of 0.62 and 0.74 as compared with the three other densities.Lignin and cellulose losses were more rapid in the first 118 days.Soil fauna had significant impacts on litter decomposition and the degradation of refractory compounds.Canopy density had significant effects on factors such as soil properties and soil fauna community structure,which could be conducive to the decomposition of litter and the degradation of litter recalcitrant components.Canopy density between 0.6 and 0.7 might be a favorable management practice promoting litter decomposition and beneficial for the sustainable development of P.massoniana plantations.

Tree-ring response of Larix chinensis on regional climate and sea-surface temperature variations in alpine timberline in the Qinling Mountains

Tree-ring width chronologies of Larix chinensis were developed from the northern and southern slopes of the Qinling Mountains in Shaanxi Province,and climatic factors affecting the tree-ring widths of L.chinensis were examined.Correlation analysis showed that similar correlations between tree-ring width chronologies and climatic factors demonstrated that radial growth responded to climate change on both slopes.The radial growth of L.chinensis was mainly limited by temperature,especially the growing season.In contrast,both chronologies were negatively correlated with precipitation in May of the previous year and April of the current year.Spatial climate-correlation analyses with gridded land-surface climate data revealed that our tree-ring width chronologies contained a strong regional temperature signal over much of northcentral and eastern China.Spatial correlation with seasurface temperature fields highlights the influence of the Pacific Ocean,Indian Ocean,and North Atlantic Ocean.Wavelet coherence analysis indicated the existence of some decadal and interannual cycles in the two tree-ring width chronologies.This may suggest the influences of El Nin˜o-Southern Oscillation and solar activity on tree growth in the Qinling Mountains.These findings will help us understand the growth response of L.chinensis to climate change in the Qinling region,and they provide critical information for future climate reconstructions based on this species in semi-humid regions.

Performance of Different Cultivars in Direct Seeded Rice (<i>Oryza sativa</i>L.) with Various Seeding Densities

Optimum seeding density is necessary to achieve the production potential of a cultivar in all crops. To investigate the performance of different cultivars in direct seeded rice in response to various seeding densities, a field experiment was conducted at Agronomic Research Area, University of Agriculture, Faisalabad during the kharif 2012. Super basmati (fine rice) and KSK-133 (coarse rice) were sown using seeding densities of 30, 45, 60 and 75 kg&middotha-1. Results indicated that leaf area index and crop growth rate varied among rice cultivars as well as various seeding densities. The growth of rice at seeding density of 75 kg&middotha-1 was higher than rest of treatments. The response of rice yield and its components like number of productive tillers, number of kernels per panicle and 1000-kernal weight differed significantly under the influence of cultivars and various seeding densities. Seeding density of 75 kg&middotha-1 remained superior for both rice cultivars regarding final yield. This treatment furnished kernel yields of 3.83 and 5.49 t&middotha-1 in fine and coarse rice, respectively which were significantly higher than other treatments. However, Kernel quality was improved in terms of lower number of sterile and chalky kernels at lower seeding densities (30 kg&middotha-1). Based on these findings, it is concluded that using seeding density of 75 kg&middotha-1 for direct seeded fine and coarse rice cultivars is beneficial in terms of better growth and higher yield but kernel quality is reduced by increasing the seeding density.

Tree-ring-based reconstruction of temperature variability(1445–2011) for the upper reaches of the Heihe River Basin, Northwest China

Long-term temperature variability has significant effects on runoff into the upper reaches of inland rivers. This paper developed a tree-ring chronology of Qilian juniper (Sabina przewalskii Kom.) from the upper tree-line of the middle Qilian Mountains within the upper reaches of Heihe River Basin, Northwest China for a long-term reconstruction of temperature at the study site. In this paper, tree-ring chronology was used to examine climate-growth associations considering local climate data obtained from Qilian Meteorological Station. The results showed that temperatures correlated extremely well with standardized growth indices of trees (r=0.564, P<0.001). Tree-ring chronology was highest correlated with annual mean temperature (r=0.641, P<0.0001). Annual mean temperature which spans the period of 1445–2011 was reconstructed and explained 57.8% of the inter-annual to decadal temperature variance at the regional scale for the period 1961–2011. Spatial correlation patterns revealed that reconstructed temperature data and gridded temperature data had a significant correlation on a regional scale, indicating that the reconstruction represents climatic variations for an extended area surrounding the sampling sites. Analysis of the temperature reconstruction indicated that major cold periods occurred during the periods of 1450s–1480s, 1590s–1770s, 1810s–1890s, 1920s–1940s, and 1960s–1970s. Warm intervals occurred during 1490s–1580s, 1780s–1800s, 1900s–1910s, 1950s, and 1980s to present. The coldest 100-year and decadal periods occurred from 1490s–1580s and 1780s–1800s, respectively, while the warmest 100 years within the studied time period was the 20th century. Colder events and intervals coincided with wet or moist conditions in and near the study region. The reconstructed temperature agreed well with other temperature series reconstructed across the surrounding areas, demonstrating that this reconstructed temperature could be used to evaluate regional climate change. Compared to the tree-ring reconstructed temperature from nearby regions and records of glacier fluctuations from the surrounding high mountains, this reconstruction was reliable, and could aid in the evaluation of regional climate variability. Spectral analyses suggested that the reconstructed annual mean temperature variation may be related to large-scale atmospheric–oceanic variability such as the solar activity, Pacific Decadal Oscillation (PDO) and El Ni?o–Southern Oscillation (ENSO).