The Geochemistry of Supraglacial Streams of Canada Glacier,Taylor Valley (Antarctica), and their Evolution into Proglacial Waters

We have investigated the geochemistry of supraglacial streams on the Canada Glacier, Taylor Valley, Antarctica during the 2001–2002 austral summer. Canada Glacier supraglacial streams represent the link between primary precipitation (i.e. glacier snow) and proglacial Lake Hoare. Canada Glacier supraglacial stream geochemistry is intermediate between glacier snow and proglacial stream geochemistry with average concentrations of 49.1 μeq L⁻¹ Ca2+, 19.9 μeq L⁻¹ SO₄²⁻, and 34.3 μeq L⁻¹ HCO₃⁻. Predominant west to east winds lead to a redistribution of readily soluble salts onto the glacier surface, which is reflected in the geochemistry of the supraglacial streams. Western Canada Glacier supraglacial streams have average SO₄²⁻:HCO₃⁻ equivalent ratios of 1.0, while eastern supraglacial streams average 0.5, suggesting more sulfate salts reach and dissolve in the western supraglacial streams. A graph of HCO₃⁻ versus Ca2+ for western and eastern supraglacial streams had slopes of 0.87 and 0.72, respectively with R2 values of 0.84 and 0.83. Low concentrations of reactive silicate (> 10 μmol L⁻¹) in the supraglacial streams suggested that little to no silicate weathering occurred on the glacier surface with the exception of cryoconite holes (1000 μmol L⁻¹). Therefore, the major geochemical weathering process occurring in the supraglacial streams is believed to be calcite dissolution. Proglacial stream, Anderson Creek, contains higher concentrations of major ions than supraglacial streams containing 5 times the Ca²⁺ and 10 times the SO₄²⁻. Canada Glacier proglacial streams also contain higher concentrations (16.6–30.6 μeq L⁻¹) of reactive silicate than supraglacial streams. This suggests that the controls on glacier meltwater geochemistry switch from calcite and gypsum dissolution to both salt dissolution and silicate mineral weathering as the glacier meltwater evolves. Our chemical mass balance calculations indicate that of the total discharge into Lake Hoare, the final recipient of Canada Glacier meltwater, 81.9% is from direct glacier runoff and 19.1% is from proglacial Andersen Creek. Although during a typical, low melt ablation season Andersen Creek contributes over 40% of the water added to Lake Hoare, its overall chemical importance is diluted by the direct inputs from Canada Glacier during high flow years. Decadal warming events, such as the 2001–2002 austral summer produce supraglacial streams that are a major source of water to Lake Hoare.     

Characteristics of winter mass balance of Glacier No.1 at the headwaters of the Urumqi River,Tianshan Mountains

This study analyzes the response of glacier to climate change during the past 49 years in Urumqi River source region, the Tianshan Mountains of China. The temporal and spatial variations of winter mass balance (bn-w) at different time scales were analyzed to identify their response to climate change during 1988–2006 (The observation of winter mass balance observation began in 1988) on the Glacier No.1 at the headwaters of the Urumqi River, Tianshan Mountains, China. The winter accumulation shows a significantly decreasing trend. The results show that the cumulative values on Glacier No.1 is 2,202 mm water equivalent during 1988–2006 and the mean values is 116 mm a⁻¹. Furthermore, the trend analysis of the winter mass balance indicates a rapid decrease since 1990, and the mean mass balance is only 79 mm a⁻¹ during 1997–2006. Winter mass balance correlates well negatively with the total evaporation from September to April (r = −0.68, α = 0.01), and positively with the total precipitation from September to April (r = 0.74, α = 0.01). However, winter mass balance shows a weak correlation with mean minimum air temperature during September to April (r = −0.35), and runoff on September (r = −0.13).     

冰川流域径流估算方法探索——以科其喀尔巴西冰川为例

利用科其喀尔巴西冰川2005年6月至2006年5月水文、气象资料,并结合15 m融合TM、1∶[KG-*2]50000地形图、FY 2C数值产品和NCEP/NCAR再分析资料等,构建了10个简单、具有一定自主创新意义的分布式冰川融水径流模型（空间分辨率60 m）,较好地模拟了研究冰川流域的日平均流量。结果表明,利用FY 2C总云量资料并结合辐射传输参数化方案能够较好地估算流域太阳入射短波辐射;单独利用总辐射有直接估算大型冰川流域某段时期融水径流的可能。气温与冰川末端流量呈指数关系,度日因子模型更适合于消融季节;提出的基于单元格气温和海拔的简单消融模型有望改进度日因子模型。在气温指数模型中加入太阳辐射调整系数,能够更好地估算冰川融水径流。简化分布式能量平衡模型能够反映大型冰川融水径流的变化;单层汇流方案在一定程度上能够概化托木尔型冰川的汇流过程。     

天山南坡科其卡尔巴契冰川度日因子变化特征研究

度日模型是估算冰川消融的一种简单而有效的方法.根据科其卡尔巴契冰川2003年的观测资料,分析了该冰川度日因子的空间变化规律及其影响因素.研究表明:各高度上的度日因子,介于2.0～9.7mm·℃^-1·d^-1之间变化,平均值为5.7mm·℃^-1·d^-1,与青藏高原各冰川及其它地区冰川相比较小;随着海拔的增高,度日因子随之递增;随平均气温的升高而随之递减.由于冰面状况复杂,度日因子变化幅度较大,裸冰区的度日因子明显大于表碛覆盖区.人为测量误差、反照率、地形等对度日因子的影响也不容忽视.     

The relationship between runoff and ground temperature in glacierized catchments in China

A single parameter index method, in which ground temperature and air temperature is the sole input variable, respectively, is used to evaluate and compare the glacial runoff in three typical glacier catchments, Dongkemadi glacier catchment in Tibetan plateau, Koxkar glacier catchment and the headwater catchment of Urumqi River catchment in Tianshan Mountains in West China. The method based on ground temperature is an attempt to evaluate glacier runoff in elevated terrains, as few studies have focused specifically on the association between glacier runoff and ground temperature. The results identify ground temperature versus a certain depth, which is a critical factor that affected glacier hydrological processes and showed that runoff data is much better correlated with ground temperature than air temperature. Especially, at the latter two catchments, the largest coefficients of exponential relationship R ² between glacier runoff and ground temperature are 0.9 and 0.83, respectively. The accuracy of the method makes it possible to estimate the glacier runoff with a certain depth ground temperature at a certain site, which may provide a new approach to evaluate the glacier runoff for areas where there is a lack of observation data.     

天山南坡科其喀尔冰川消融区运动特征分析

基于2009年5月至2011年10月科其喀尔冰川的花杆观测资料,对其消融区的表面运动特征进行分析. 结果表明：冰川消融区的年水平运动速度最大值为86.69 m·a^-1,年垂直运动速度最大值为15.34 m·a^-1,均出现在冰川海拔4 000~4 200 m的消融区上部;在靠近冰川末端的冰舌下部,受冰量补给减弱、厚层表碛覆盖等影响,冰川运动缓慢,年水平运动速度小于5 m·a^-1,而垂直运动速度值小于2 m·a^-1. 大多数横剖面的水平运动速度具有从中部向边缘逐渐减小的特征,而有的剖面却出现局部速度增大的区域. 整体而言,冰川水平及垂直运动速度随海拔降低而减小,符合冰川运动的一般规律,但主要受地形作用的影响,垂直运动速度随海拔的变化会出现波动. 消融期月水平运动速度与同期气温和降水的变化具有一定的相关性,可能反映出气候快速变化对冰川运动的影响.     

Glacier No. 4 of Sigong River over Mt. Bogda of eastern Tianshan, central Asia: thinning and retreat during the period 1962–2009

RTK-GPS data, aerial photographs and Aster images were used to quantify volume, surface elevation, terminus position and area changes of Glacier No. 4 of Sigong River over Mt. Bogda, Tianshan during the period from 1962 to 2009. Glacier surface elevation of the tongue area decreased by 15 ± 8 m (0.32 ± 0.17 m a⁻¹) and ice volume loss reached 0.014 ± 0.008 km³ (0.013 ± 0.007 km³ w. e.). The glacier terminus retreated at a rate of 8.0 m a⁻¹ and the area decreased by about 0.53 km², accounting for 15.8% of the glacier area in 1962 (3.33 km²). The changes can be primarily attributed to the significant increase in temperature in this region. A comparison with glacier changes by field measurements in other regions of eastern Tianshan showed obvious spatial differences in the magnitude of the changes, owing to a combination of regional climate change and topographical factors.     

Determining Spatial and Temporal Inputs of Freshwater,Including Submarine Groundwater Discharge,to a Subtropical Estuary Using Geochemical Tracers,Biscayne Bay,South Florida

Geochemical mixing models were used to decipher the dominant source of freshwater (rainfall, canal discharge, or groundwater discharge) to Biscayne Bay, an estuary in south Florida. Discrete samples of precipitation, canal water, groundwater, and bay surface water were collected monthly for 2 years and analyzed for salinity, stable isotopes of oxygen and hydrogen, and Sr²⁺/Ca²⁺ concentrations. These geochemical tracers were used in three separate mixing models and then combined to trace the magnitude and timing of the freshwater inputs to the estuary. Fresh groundwater had an isotopic signature (δ ¹⁸O = −2.66‰, δD −7.60‰) similar to rainfall (δ ¹⁸O = −2.86‰, δD = −4.78‰). Canal water had a heavy isotopic signature (δ ¹⁸O = −0.46‰, δD = −2.48‰) due to evaporation. This made it possible to use stable isotopes of oxygen and hydrogen to separate canal water from precipitation and groundwater as a source of freshwater into the bay. A second model using Sr²⁺/Ca²⁺ ratios was developed to discern fresh groundwater inputs from precipitation inputs. Groundwater had a Sr²⁺/Ca²⁺ ratio of 0.07, while precipitation had a dissimilar ratio of 0.89. When combined, these models showed a freshwater input ratio of canal/precipitation/groundwater of 37%:53%:10% in the wet season and 40%:55%:5% in the dry season with an error of ±25%. For a bay-wide water budget that includes saltwater and freshwater mixing, fresh groundwater accounts for 1–2% of the total fresh and saline water input.     

天山托木尔峰科其喀尔巴西冰川表面运动速度特征分析

天山托木尔峰科其喀尔巴西冰川是典型的树枝状山谷冰川,利用3组(6期)ASTER遥感影像通过COSI-corr软件反演了该冰川表面运动速度.与花杆测量数据进行对比,反演冰川表面运动速度平均绝对误差为3.1 m·a-1,相对误差为11.9%,二者在空间上的分布基本一致,表明其反演精度符合要求.在此基础上,分析冰川表面运动速...     

Solute dynamics of meltwater of Gangotri glacier,Garhwal Himalaya,India

The present study investigates solute dynamics of meltwater of Gangotri glacier system in terms of association of different chemical compounds with the geology of the area. In the meltwater, the presence of cations varied as c(Mg²⁺) > c(Ca²⁺) > c(Na⁺) > c(K⁺), while order of concentration of anions has been c(HCO₃ ⁻) > c(SO₄ ²⁻) > c(Cl⁻) > c(NO₃ ⁻) in years 2003 and 2004. The magnesium and calcium are found as the dominant cations along with bicarbonate and sulphate as dominant anions. The high ratios of c(Ca²⁺ + Mg²⁺)/total cations and c(Ca²⁺ + Mg²⁺)/c(Na⁺ + K⁺) indicate that the meltwater chemistry of the Gangotri glacier system catchment is mostly controlled by carbonate weathering. Attempts are made to develop rating curves for discharge and different cations. Sporadic rise in discharge without corresponding rise in concentration of most of cations is responsible for their loose correlation in a compound valley glacier like Gangotri glacier.     

Study on the change of Keqikaer Glacier during the last 30 years,Mt. Tuomuer,Western China

Using radio echo sounder, ice thickness of Keqikaer Glacier tongue was measured in 1981 and 2004. Data obtained by comparing topographical maps, aerial and satellite photographs at different times, illustrates changes of the thickness and advance/retreat of Keqikaer Glacier. Keqikaer Glacier has been in intensive retreat since the 1990s and become thinner since the 1980s. Measured thickness of the ice tongue indicates reducing with a speed of 0.5–1.5 m a⁻¹ since 1981. The shrinkage of the glacier terminus is less than 2% of the total length during the last 30 years; however, the retreat of terminus position and the thinning of the ice thickness provides significant information that these glaciers on the south slope of Mt. Tuomuer are in an intensively decreasing phase in recent time.     

西天山托木尔峰南麓大型山谷冰川冰舌区消融特征分析

基于对托木尔峰南麓托木尔型山谷冰川的野外考察和典型冰川的定位观测,对冰面被表碛广泛覆盖的所谓“托木尔型”冰川冰舌区表碛与冰面消融的关系进行了研究. 结果表明：表碛对冰面消融、冰川水文过程、冰川变化等均具有重要影响,当表碛厚度超过3 cm时,表碛对冰面消融就产生明显抑制作用,且随着厚度增加,冰面消融显明减弱. 科其喀尔冰川表面的观测表明,由末端向上,表碛厚度逐渐减薄. 受表碛影响,科其喀尔冰川区最大的消融量出现在海拔3 800~3 900 m之间、表碛物厚度小于10 cm的区域内;冰川消融强度由此向上随着海拔的升高而下降,向下随表碛厚度的增大而减弱. 冰面湖的发育是表碛覆盖冰川的又一主要特征,湖水对冰面的融蚀和快速排泄成为冰面产汇流的主要过程. 科其喀尔冰川研究表明,两三个冰面湖排泄形成的融蚀冰量就相当于冰川末端退缩造成的冰量损失. 因此,冰面湖等热喀斯特地形的形成、扩张融蚀、融穿排泄、形成湖区低地,这一周而复始的过程不仅是其主要消融方式之一,而且也强烈的影响着冰川水文及冰川变化. 托木尔峰南麓地区大型冰川变化主要以厚度减薄为主,而不是像大多数冰川显著的变化主要表现在末端和面积减少方面.     

Phytoplankton Size and Taxonomic Composition in a Temperate Estuary Influenced by Monsoon

Temporal and spatial variations in phytoplankton in Asan Bay, a temperate estuary under the influence of monsoon, were investigated over an annual cycle (2004). Phytoplankton blooms started in February (>20 μg chl l⁻¹) and continued until April (>13 μg chl l⁻¹) during the dry season, especially in upstream regions. The percentage contribution of large phytoplankton (micro-sized) was high (78–95%) during the blooms, and diatoms such as Skeletonema costatum and Thalassiosira spp. were dominant. The precipitation and freshwater discharge from embankments peaked and supplied nutrients into the bay during the monsoon event, especially in July. Species that favor freshwater, such as Oscillatoria spp. (cyanobacteria), dominated during the monsoon period. The phytoplankton biomass was minimal in this season despite nutrient concentrations that were relatively sufficient (enriched), and this pattern differed from that in tropical estuaries affected by monsoon and in temperate estuaries where phytoplankton respond to nutrient inputs during wet seasons. The flushing time estimated from the salinity was shorter than the doubling time in Asan Bay, which suggests that exports of phytoplankton maximized by high discharge directly from embankments differentiate this bay from other estuaries in temperate and tropical regions. This implies that the change in physical properties, especially in the freshwater discharge rates, has mainly been a regulator of phytoplankton dynamics since the construction of embankments in Asan Bay.     

Large-scale moraine deformation at the Athabasca Glacier,Jasper National Park,Alberta, Canada

In this paper the development of a large-scale gravitational deformation involving the eastern lateral moraine of the Athabasca Glacier in Jasper National Park, Alberta, Canada, is described. Interpretation and analysis of sequential aerial photographs indicates that a 540-m-wide segment of the eastern lateral moraine began to deform in the early 1950s; however, significant movement only began in the late 1960s. Since then, the moraine has undergone progressive gravitational deformation leading to a network of fractures, bulging, and the development of a large gap in the moraine crest. Geographic information system analysis of topographic changes between 1967 and 2006 indicates that the displaced volume of the moraine is approximately 9.0 × 10⁵ m³. In the last 39 years, the moraine crest has displaced 55 m (1.4 m yr⁻¹) down towards the glacier. The development of slope instability is linked to a combination of debuttressing from recent glacier recession, deformation of the moraine, as well as the movement of a large, mobile, debris-mantled slope impinging the upslope margin of the lateral moraine. This case study illustrates the importance of glacial conditioning and local geomorphological factors in creating conditions for large-scale moraine instability in recently deglacierized alpine basins.     

Insights into silicic melt generation using plagioclase, quartz and melt inclusions from the caldera-forming Rotoiti eruption, Taupo volcanic zone, New Zealand

The Rotoiti (~120 km³) and Earthquake Flat (~10 km³) eruptions occurred in close succession from the Okataina Volcanic Centre at ~50 ka. While accessory mineral geochronology points to long periods of crystallization prior to eruption (10⁴–10⁵ years) and separate thermal histories for the magmas, little was known about the rates and processes of the final melt production and eruption. Crystal zoning patterns in plagioclase and quartz reveal the thermal and compositional history of the magmatic system leading up to the eruption. The dominant modal phase, plagioclase, displays considerable within-crystal zonation: An_37–74, ~40–227 ppm MgO, 45–227 ppm TiO₂, 416–910 ppm Sr and 168–1164 ppm Ba. Resorption horizons in the crystals are marked by sharp increases (10–30%) in Sr, MgO and X_An that reflect changes in melt composition and are consistent with open system processes. Melt inclusions display further evidence for open system behaviour, some are depleted in Sr and Ba relative to accompanying matrix glass not consistent with crystallization of modal assemblage. MI also display a wide range in XH₂O that is consistent with volatile fluxing. Quartz CL images reveal zoning that is truncated by resorption, and accompanied by abrupt increases in Ti concentration (30–80 ppm) that reflect temperature increases ~50–110°C. Diffusion across these resorption horizons is restricted to zones of <20 μm, suggesting most crystallization within the magma occurred in <2000 years. These episodes are brief compared to the longevity (10⁴–10⁵ year) of the crystal mush zones. All textural and compositional features observed within the quartz and plagioclase crystals are best explained by periodic mafic intrusions repeatedly melting parts of a crystal-rich zone and recharging the system with silicic melt. These periodic influxes of silicic melt would have accumulated to form the large volume of magma that fed the caldera-forming Rotoiti eruption.     

Relationships among satellite chlorophylla, river inputs, and hypoxia on the Louisiana Continental shelf, Gulf of Mexico

SeaWiFS ocean color measurements were used to investigate interannual, monthly, and weekly variations in chlorophylla (chla) on the Louisiana shelf and to assess relationships with river discharge, nitrate load, and hypoxia. During the study period (2000–2003), interannual changes in shelf-wide chla concentrations averaged over January–July ranged from +57% to −33% of the 4-yr average, in accord with freshwater discharge changes of +20% to −29% and nitrate load changes of +20% to −35% from the Mississippi and Atchafalaya Rivers. Chla variations were largest on the shelf between the Mississippi and Atchafalaya Deltas. Within this region, which corresponds spatially to the area of most frequent hypoxia, lowest January–July mean chla concentrations (5.5 mg m⁻³ over 7,000 km²) occurred during 2000, the year of lowest freshwater discharge (16,136 m³ s⁻¹) and nitrate load (55,738 MT N d⁻¹) onto the shelf. Highest January–July mean chla concentrations (13 mg m⁻³ over 7,000 km²) were measured in 2002, when freshwater discharge (27,440 m³s⁻¹) and nitrate load (101,761 MT N d⁻¹) were highest and second highest, respectively. Positive correlations (R²=0.4–0.5) were found between chla and both fresh water and nitrate loads with 0 to 1 month lags, with the strongest relationships just west of the Mississippi Delta. In 2001, unusually clear skies allowed the identification of distinct spring and summer chla blooms west of the Mississippi Delta 4–5 wk after peaks in river discharge. East of the delta, the chla concentrations peaked in June and July, following the seasonal reversal in the coastal current. A clear linkage was not detected between satellite-measured chla and hypoxia during the 4-yr period, based on a time series of bottom oxygen concentrations at one station within the area of most frequent hypoxia. Clear relationships are confounded by the interaction of physical processes (wind stress effects) with the seasonal cycle of nutrient-enhanced productivity and are influenced by the prior year's nitrate load and carbon accumulation at the seabed.     

The significance of Cenozoic magmatism from the western margin of the eastern syntaxis, southeast Tibet

The major and trace-element geochemistry, Sr–Nd bulk-rock isotopes, U–Pb zircon chronology and Lu–Hf isotopic compositions are described for three granitic bodies which intrude the Nyingchi gneisses (Lhasa terrane) along the western margin of the eastern Himalayan syntaxis. The Bayi two-mica granite and Lunan granite–granodiorite were intruded at 22 ± 1 and 25.4 ± 0.3 Ma, respectively, whereas the Confluence biotite granite was emplaced at 49.1 ± 0.4 Ma. All share strong depletions in Y and HREE requiring a garnet-bearing source both during and following the Eocene collision of the Indian plate with the Lhasa terrane. The isotope geochemistry of these intrusives (ε _Nd(t) = −3 to −5, ⁸⁷Sr/⁸⁶Sr_(t) = 0.706–0.707) indicates a crustal source within the Lhasa terrane. Sr–Nd systematics of the garnet-bearing Nyingchi gneisses together with the U–Pb and Lu–Hf isotopic ratios of detrital zircons recovered from this unit identifies it as a potential melt source. The combined element and isotope geochemistry of the plutons indicate a mixed source; the gneisses provide the older component whereas the Gangdese batholith provides a younger, siliceous component. The involvement of garnet-bearing crustal material in melt sources from the Cretaceous (80 Ma) to the Miocene (20 Ma) is consistent with the presence of a thicker continental crust in the eastern Lhasa terrane, as is the presence of magmatic epidote in several plutons which indicates a regional deepening level of exposure eastwards. Post-collision crustal melting is synchronous with proposed slab break-off during the early Miocene, suggesting advective heating by rising asthenospheric melts.     

Trace element partitioning between orthopyroxene and anhydrous silicate melt on the lherzolite solidus from 1.1 to 3.2 GPa and 1,230 to 1,535°C in the model system Na<Subscript>2</Subscript>O–CaO–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>

We determined experimentally the Nernst distribution coefficient between orthopyroxene and anhydrous silicate melt for trace elements i in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂ (NCMAS) along the dry model lherzolite solidus from 1.1 GPa/1,230°C up to 3.2 GPa/1,535°C in a piston cylinder apparatus. Major and trace element composition of melt and orthopyroxene were determined with a combination of electron microprobe and ion probe analyses. We provide partitioning data for trace elements Li, Be, B, K, Sc, Ti, V, Cr, Co, Ni, Rb, Sr, Y, Zr, Nb, Cs, Ba, La, Ce, Sm, Nd, Yb, Lu, Hf, Ta, Pb, U, and Th. The melts were chosen to be boninitic at 1.1 and 2.0 GPa, picritic at 2.3 GPa and komatiitic at 2.7 and 3.2 GPa. Orthopyroxene is Tschermakitic with 8 mol% Mg-Tschermaks MgAl[AlSiO₆] at 1.1 GPa while at higher pressure it has 18–20 mol%. The rare earth elements show a continuous, significant increase in compatibility with decreasing ionic radius from D _La^opx−melt ∼ 0.0008 to D _Lu^opx−melt ∼ 0.15. For the high-field-strength elements compatibility increases from D _Th^opx−melt ∼ 0.001 through D _Nb^opx−melt ∼ 0.0015, D _U^opx−melt ∼ 0.002, D _Ta^opx−melt ∼ 0.005, D _Zr^opx−melt ∼ 0.02 and D _Hf^opx−melt ∼ 0.04 to D _Ti^opx−melt ∼ 0.14. From mathematical and graphical fits we determined best-fit values for D ₀^M1, D ₀^M2, r ₀^M1, r ₀^M2, E ₀^M1, and E ₀^M2 for the two different M sites in orthopyroxene according to the lattice strain model and calculated the intracrystalline distribution between M1 and M2. Our data indicate extreme intracrystalline fractionation for most elements in orthopyroxene; for the divalent cations D _i ^M2−M1 varies by three orders of magnitude between D _Co^M2−M1 = 0.00098–0.00919 and D _Ba^M2−M1 = 2.3–28. Trivalent cations Al and Cr almost exclusively substitute on M1 while the other trivalent cations substitute on M2; D _La^M2−M1 reaches extreme values between 6.5 × 10⁷ and 1.4 × 10¹⁶. Tetravalent cations Ti, Hf, and Zr almost exclusively substitute on M1 while U and Th exclusively substitute on M2. Our new comprehensive data set can be used for polybaric-polythermal melting models along the Earth’s mantle solidus. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

河西内陆河流域冰川融水近期变化

基于第一次中国冰川编目数据库,以国家气象台站的月降水与月气温资料为驱动数据,利用度日模型重建了河西内陆河流域各水系1961～2006年冰川物质平衡和融水径流序列,并与短期观测资料的对比验证表明,模型具有较高的可信度.河西内陆河流域1961～2006年平均冰川物质平衡为-49.5 mm/a,从东段石羊河流域到西段疏勒河、...