# Interpretation of the Chemical and Physical Time-Series Retrieved from Sentik Glacier, Ladakh Himalaya, India

Paul Mayewski, W. Berry Lyons, N. Ahmad, Gordon Smith, M. Pourchet

► To cite this version:

Paul Mayewski, W. Berry Lyons, N. Ahmad, Gordon Smith, M. Pourchet. Interpretation of the Chemical and Physical Time-Series Retrieved from Sentik Glacier, Ladakh Himalaya, India. Journal of Glaciology, 1984, 30 (104), pp.66-76. 10.3189/S0022143000008509. insu-03443630

HAL Id: insu-03443630

https://insu.hal.science/insu-03443630

Submitted on 23 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Fig. 1. Location of core site from a portion of NASA Landsat Image E-30531-04491-C displaying the Nun Kun massif in the lower center of the photograph. Peaks in massif center attain elevations as high as 7135 m (Nun) and 7077 m (Kun). Glacier with banded moraine flows directly north off the massif. Moraine bands as seen on the photograph are approximately 6 km in length. The Suru River meanders past the snout of this glacier. S marks the core site (4908 m) on Sentik Glacier.

to south-east near the ground surface and west to east at the 500 mbar level (Boucher, 1975). Summer air-mass circulation includes south-west to north-east to north-west flow near the ground surface and at the 500 mbar level air flow is eastward along the northern portion of the Himalaya and westward along the southern portion of the range (Boucher, 1975). An anticyclone rests for a portion of the summer season over the Himalaya and divides the eastward- and westward-flowing 500 mbar level air system (Boucher, 1975). The lower air mass is believed to be relatively thin (<2 km) and overlain by an arid continental air mass (Sequeira and Kelkar, 1978).

Primary maximum precipitation comes to the region of the study site during the period July to August and primary minimum during November (Boucher, 1975). Estimated total rainfall for valley bottoms in the Nun Kun area is 0.90 m a⁻¹ water equivalent (Qazi, 1973). Rainfall records are available from several stations in the region, however the data are scant and contain numerous inconsistencies. Discharge records for the Suru River (Fig. 1), which flows in the valley directly north of the core site, suggest that the warmest time of the year is June to August (Qazi, 1973).

The main core site was chosen at an elevation of 4908 m on a relatively flat, single-flow portion of a 3 km long alpine glacier, Sentik Glacier (S on Fig. 1), that flows north-east off the Nun Kun massif. Several other core sites were investigated in the massif area, at elevations as high as 5512 m, but these were either

affected by avalanches or as in the case of the Nun Kun plateau itself primarily composed of wind-blown snow.

Interpretation of down-hole glaciochemistry recovered from the Sentik Glacier core is based on the assumption that neither melting nor percolation, to depths significantly less than seasonal accumulation levels, is present at the core site. Although, as already noted, few meteorological data are available from the region, radiosonde measurements from Srinagar (1700 m a.s.l. 120 km east of Nun Kun) reveal that even during the warmest months of the year, June-August, at only 500-600 m above Srinagar, mean monthly temperatures barely rise to 0°C. Furthermore, a set of down-hole temperatures made at the study site in July suggest that although the upper few centimeters of the glacier are close to 0°C, immediately beneath the surface, temperatures drop until finally at a depth of 9 m the core attains and maintains a temperature of -3.0°C to a depth of at least 16.6m.

As an additional test for the investigation of melt-water potential, plexiglass trays designed to capture melt water were inserted in snow-pits close to the sample site. Although this monitoring was conducted in late July, close to the peak of the summer season and at several depths (5, 10, 20, 30, 40 and 50 cm) no melt water was found in the trays. The presence of several thin ice bands, comprising <5% of the total core, suggests that melt events, although not common, do occur at certain times. Hastenrath (1978) has demonstrated that high-altitude areas,
such as the summit of Quelccaya ice cap, Peru, may absorb too little net radiation to allow melting. Although we do not have radiation data to compare to Hastenrath's, the same situation may have existed at the Sentik Glacier core site for much of the time period represented by the core.

## METHODOLOGY

Samples collected from a 16.6 m core on Sentik Glacier were analyzed for density, microparticles, and chemical properties (chloride, sodium, reactive iron, reactive silicate, reactive phosphate, ammonium, δD, δ(180), total β-activity, and pH). Although the collection interval was continuous at 8 cm, most analyses were performed at 16 cm intervals due to sample volume constraints. Extreme care was taken at all times during sample collection, handling, and analyses to prevent sample contamination as required by the nature of the chemical studies. Cores were taken using a manually operated pre-cleaned teflon-coated SIPRE-type auger. Plastic suits, plastic gloves, and particle masks were worn by all members of the party who came in direct contact with the exposed core.

Each sample had its outer 1-2 cm scraped immediately after sectioning from the cores using pre-cleaned plexiglass scrapers, and the scrapings were discarded to reduce possible contamination by the auger. Once scraped, samples were placed in pre-cleaned plastic bags. Pre-cleaning of bags, scrapers, and the auger involved a soap and water wash with rinsing several times in distilled-deionized water. After melting in the plastic bags each sample was transferred to two 30 ml and one 125 ml (LPE) container, pre-cleaned as described above, and one 125 ml (LPE) container that was pre-cleaned by acid-

washing in concentrated hydrochloric acid followed by rinsing several times in distilled-deionized water.

The portion of the sample in the soap-and-water cleaned 125 ml container was used for analysis of chloride (Zall and others, 1956) and of reactive phosphate employing AutoAnalyzer techniques (Glibert and Loder, 1977). The portion of each sample in the acid-cleaned 125 ml container was used for analysis of reactive silicate and reactive iron by AutoAnalyzer using the techniques of Glibert and Loder (1977) and Stookey (1970), respectively. These samples were acidified with ultrapure HNO3 and allowed to sit at room temperature for at least two weeks prior to analysis. Sodium was analyzed via atomic absorption spectrophotometry. All of the aforementioned were analyzed at the University of New Hampshire. Duplicate analyses of standards were used to determine the precision of analysis for chloride, reactive iron, reactive silicate, reactive phosphate and sodium. Analytical precision of the replicate analyses of standards was ±1% at 7.0 µM, ±1% at 4.5 µM, ±1% at 4 µM, ±4% at 0.5 µM, and ±3% at 13.0 µM, respectively.

The two 30 ml (LPE) containers mentioned previously were used for collecting portions of the melted samples that were analyzed for δD, δ(180), and total β-activity measurements were made in the CNRS Laboratoire de Glaciologie, France, on one set of 30 ml containers that were sealed with wax to prevent diffusion. Volumes required for total β-activity measurement required lumping of samples. pH was analyzed in the field on the other set of 30 ml samples using an Orion Ion-analyzer T.M. model 399A portable pH meter with an Orion 91-06 combination pH electrode standardized using NBS buffers of pH = 4.01 and pH = 7.00. pH analyses were made without stirring. Microparticles were analyzed at The Ohio State University using the 30 ml sample set from pH. The technique employed a Coulter

<table><caption>TABLE I. MEAN VALUES FOR TIME-SERIES SETS IN APPENDIX A</caption><thead><tr><th>Time-series set</th><th>Mean</th><th>Sample number</th></tr></thead><tbody><tr><td>Chloride</td><td>2.94 µM</td><td>114</td></tr><tr><td>Sodium</td><td>3.14 µM</td><td>114</td></tr><tr><td>Reactive iron</td><td>13.91 µM</td><td>228</td></tr><tr><td>Reactive silicate</td><td>8.35 µM</td><td>228</td></tr><tr><td>Reactive phosphate</td><td>0.39 µM</td><td>228</td></tr><tr><td>Ammonium</td><td>1.07 µM</td><td>114</td></tr><tr><td>δD</td><td>-97.90 °/oo</td><td>114</td></tr><tr><td>δ(180)</td><td>-15.07 °/oo</td><td>19</td></tr><tr><td>pH</td><td>4.12-5.85 pH units*</td><td>114</td></tr><tr><td>Total β-activity</td><td>224.76 dph kg<sup>-1</sup></td><td>21</td></tr><tr><td>Density</td><td>570 kg m<sup>-3</sup></td><td>209</td></tr><tr><td colspan="3">Microparticles (examples of 3 out of the 15 size grades available appear below)</td></tr><tr><td>0.50-0.63 µm</td><td>34 x 10<sup>3</sup> µm size particles/500 µl of sample</td><td>114</td></tr><tr><td>1.00-1.25 µm</td><td>20 x 10<sup>2</sup> µm size particles/500 µl of sample</td><td>114</td></tr><tr><td>8.00-10.00 µm</td><td>3 µm size particles/500 µl of sample</td><td>114</td></tr></tbody></table>

* range rather than mean is specified due to the log scale used to measure pH.
Counter (L.G. Thompson, 1977) and yielded 15 size grades in the range 0.50 to 12.70 µM.

## DATA ARRAY AND INTERPRETATION

Time-series sets of down-hole chemical properties (chloride, sodium, reactive iron, reactive silicate, reactive phosphate, δD, δ(¹⁸O), pH, and total β-activity) and physical properties (debris and ice band locations, density, and three selected microparticle size grades out of the 15 size grades analyzed) appear in Appendix A. Table I contains a list of mean values for all of the data sets appearing in Appendix A. Results stemming from the analysis of this data array follow.

## Core chronology

Total β-activity measurements although integrated over more than one sample are available for the entire length of the core. The resultant total β-activity record has a built-in ±0.3 error in down-hole position due to the sampling interval required for this analysis. In the depth range 15.40-15.96 m a maximum total β-activity value for the core was encountered, 639 dph kg⁻¹, which is believed to coincide with the 1963 thermonuclear test level. The time period represented by the core is therefore ≈ 17 ± 0.3 years which is

equivalent to an average net mass balance of ≈ 0.92 m a⁻¹. Using an average core density of 570 kg m⁻³ this yields an average net mass balance of 0.62 m a⁻¹ water equivalent.

## Periodic events

With the exception of total β-activity, δ(¹⁸O), and debris and ice band locations all of the data sets in Appendix A are nearly regularly spaced. The latter yield an array of 23 time series sets each with a value every 0.16 m down core and some (reactive iron, reactive silicate, and reactive phosphate) with values every 0.08 m down core. The volume of data and potential combinations presented by these data require synthesis using a statistical technique capable of revealing trends such as periodicity. Spectral analysis involving the computation of spectral density (power) at chosen frequencies provides just such a statistical tool. Power at a given frequency is the measure of how much of the variance of a time-series, or in this case a depth-series exists at a particular frequency (period).

Power, frequency, period, and confidence interval, at specified frequencies, appear in Figure 2a-c. For most of these depth-series several power peaks are recognizable. Reliability of peaks was tested using simplified filtering and by varying lag number in sev-

<figure><img src="image_3.png" /><figcaption>2a - includes time-series sets displaying annual and subannual power peaks at 10 lags (n = 114).</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>2b - includes time-series sets displaying annual and subannual power peaks at 20 lags (n = 228).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>2c - includes time-series sets at 10 lags displaying just subannual power.</figcaption></figure>

Fig. 2. Power (unit²/cy m⁻¹), period (m) and frequency (cy m⁻¹) for time-series sets. Units equal to µM for chemical species, °/oo for δD, pH units for pH, kg m⁻³ for density and particles/500 µl of sample for specified µm microparticle size grades. Vertical bar represents 80% confidence interval for power peaks at 12 degrees of freedom for n = 114 and 24 degrees of freedom for n = 228 (after Koopmans, 1974). The lower limit of the vertical bar defines whether the peak departs significantly enough from baseline to the classified as a peak. Note that confidence interval is valid only at frequency (period) specified. Confidence interval for period = ± 0.16 m using a 2σ sampling interval.

eral runs. These peaks occur at periodicities of: ≥3.22 m, <3.22 to >0.32 m and ≤0.32 m. Using the chronology provided by the total β-activity measurements and the resultant average net mass balance, 0.92 m a⁻¹, a depth to time conversion yields power peaks at ≥3.5 a, <3.5 to <0.35 a, and ≤0.35 a.

Power peaks with periodicity ≥3.5 a (≥3.22 m) are displayed by all time-(depth-)series. These events cannot, however, be precisely constrained with the sample number and length of record available in this study. Power at ≥3.5 a is both a product of the length of the time-series used in this study and of "real" trends in the data. The latter was tested by

increasing the number of lags examined in the spectral study. The presence of power in the multi-annual period suggests that collection of longer time-series, i.e. deeper cores, may yield valuable information concerning events with periodicities ≥3.5 a.

Power peaks with periodicity ≤0.35 a (≤0.35 m) are displayed by several time-series. However, events represented by power in this periodicity are almost indistinguishable from noise and thus at the sampling interval used in this study further resolution of these events is not statistically justifiable.

Assuming an average net mass balance for the core of 0.92 m a⁻¹ and an error in sample position of
± 0.16 m events with periodicities of 0.92 ± 0.15 m and 0.46 ± 0.16 m are annual and subannual, respectively. All time-series sets display power in the subannual period while only chloride, reactive phosphate, microparticles in the size range <0.50 to 0.63 µm, and reactive silicate display significant power peaks in the annual range.

### Non-periodic events

Since spectral analysis provides a view of the periodic elements in a time-series assemblage this technique yields a smoothed display of the data. Non-periodic components of the time-series sets, although comprising a small portion of the record in this study, need to be examined on an individual basis. The most apparent non-periodic anomaly in concentration is found at the depth 3.0 to 3.5 m. It is displayed by density, debris bands, and the time-series for chloride, sodium, reactive iron, reactive silicate, reactive phosphate, ammonium, and microparticles in the range 0.50 to 0.80 µm and 3.17 to 10.00 µm as relatively high values. The presence of a major density transition at this depth suggests that densification and/or ablation may exert some control on the distribution of concentrations down core. Alternatively, relatively high concentrations of chemical species may imply deposition during a period of intensified atmospheric circulation. The lack of spectral power in the annual period for the density time-series and the general non-synchroneity of the spectral signature (Fig. 2c) of this time-series with others in this study suggests that trends in the density time-series are not associated with trends in the chemical species and/or microparticle time-series.

Notably although several of the chemical time-series record annual and subannual events, the number of debris bands encountered down core does not match the total β-activity of the core. The presence of debris bands is therefore not usable as a stratigraphic marker but rather as an indication of events such as unusual ablation and/or intensified transport of debris to the core site.

### Sources

Chemical and physical constituents in the core are introduced from several source types, size source areas, and distances. Source input timing and atmospheric circulation interpretations, where appropriate, for the chemical and physical constituents analyzed are discussed in the following order: chloride, sodium, reactive iron, reactive silicate, reactive phosphate, ammonium, pH, microparticles, δD, and δ(18O). Volcanism although a potential source for several of the preceding has been deleted from the discussion that follows since it is unlikely as a periodically occurring source at the study site.

The primary source for the chloride in this core is marine and the closest marine source to Nun Kun during the period of the summer monsoon is the Arabian Sea, a distance of ≈ 1250 km from the core site. Values of chloride from fresh surface snows analyzed at several sites in the Nun Kun area (Mayewski and others, 1983) compared to those of fresh ran samples collected 75 km down-wind from Nun Kun at Gulmarg (Sequeira and Keilkar, 1978), suggest progressive depletion of chloride as air travels off the Arabian Sea to the core site (Mayewski and others, 1983). The annual power peak derived from the spectral analysis of the chloride time-series is larger than the sub-annual peak which is consistent with the view that the major input of chloride comes once each year. Since concentrations of chloride in fresh summer snow (Mayewski and others, 1983) are very similar to those representing annual pulses in the time-series record, the high core values are assumed to be summer events. Subannual input of chloride probably results from the inflow of oceanic air which to a varying degree enters this site throughout the year from sources more distant than the Arabian Sea.

Sodium is brought to the core site either in com-

bination with chloride from a sea-salt source or independently from crustal weathering. Since sodium has a higher mean concentration than chloride it probably comes from an admixture of these sources. Sodium has a relatively strong subannual spectral power peak and a weak annual peak. Crustal weathering would be a source from any wind direction and, therefore, any season thus explaining the stronger peak in the sub-annual period.

The primary source for periodically occurring reactive iron and reactive silicate is crustal weathering. Reactive iron has the higher mean concentration and a dominant subannual power peak suggesting perhaps a more local and/or larger source. Annual power in the reactive silicate spectral signature is less prominent than subannual power suggesting that annual inputs of this species come from a source that is limited or long-travelled by comparison with the sub-annual input or inputs of reactive iron. Subannual input of the reactive silicate may come from the same local weathering source as the reactive iron while the annual reactive silicate pulse may be separate. Measurements of reactive iron and reactive silicate on fresh summer snow from Nun Kun (Mayewski and others, 1983) suggest that these species are derived, in general, from a relatively low-elevation air mass directed off either the Indian subcontinent or the Tibetan Plateau.

Reactive phosphate has potentially both biogenic and anthropogenic sources, and for this study displays spectral power in both the annual and subannual ranges. Examination of temporal and spatial distributions measured in this study of reactive phosphate reveal: (1) that several of the high reactive phosphate values in the time-series correspond with debris band positions while the mean reactive phosphate value of samples having no visible debris bands is lower than that of rain from remote northern hemispheric areas <0.06 µM (Graham and Duce, 1979) and (2) similarities in distribution of reactive phosphate and reactive iron concentrations in fresh summer snow collected from several elevations on Nun Kun (Mayewski and others, 1983). The foregoing are suggestive of phosphate influx due to adsorption onto wind-blown materials such as iron oxide coatings as might be expected from an agricultural site where exposed terrain has been subjected to fertilizer usage (Mayewski and others, 1983). Aerosol loading from eroded soil material is a notable component of the atmosphere in north-west India (Bryson and Swain, 1981). Since summer reactive phosphate concentrations, determined from the fresh snow study (Mayewski and others, 1983) are similar to the maximum values displayed in the reactive phosphate time-series the annual power peak probably represents summer input. While transport of wind-blown material may be a likely source for summer inputs of reactive phosphate, subannual inputs may also be derived from this source as well as biogenic emissions and fuel-burning (Davidson and others, 1981).

Ammonium present in the core samples has as its most probable sources in this area anthropogenic and biogenic activity. The ammonium time-series display the same relative power in both the annual and sub-annual period as might be expected considering potential multiple sources and multiple timing of inputs. No unique sources can be defined and thus this species may not be useful as a stratigraphic marker. However, analysis of ammonium from freshly fallen summer snow collected in the study area reveals that the relatively high values, >0.5 µM, are found at elevations below 5300 m (Mayewski and others, 1983). Examination of Nepalese aerosols by Ikegami and others (1980) reveals a similar distribution in ammonium concentration with ammonium sulfate particles at elevations <3000 m and sulfate particles at >5000 m.

pH has similar spectral power in both annual and subannual period which is probably indicative of the fact that it is a composite measure of strong acids (i.e. nitric and sulfuric) and weak acids (i.e. carbonic) from multiple sources. While pH measurements
on fresh summer snow from Nun Kun in conjunction with other measurements have been used to help differentiate air masses (Mayewski and others, 1983) their use as a climate reconstruction tool in areas such as the Himalaya is suspect not only because of potential multiple sources but also because a number of workers (Scholander and others, 1961; Berner and others, 1978; Stauffer and Berner, 1978; Oeschger and others, 1982) have shown that even polar ice and snow are easily "contaminated" with CO₂.

The fifteen size grades of microparticles measured as part of this study display a regular decrease in concentration per size grade as expected from Junge (1963, p. 111-201). This relationship indicates that sorting in transport and source area exert control on the microparticle distributions measured in this study. Microparticle sources are not only varied but they can be additive and include, as summarized by E.M. Thompson (1980), chemical reactions (natural and anthropogenic) and incorporation of volcanic, extraterrestrial, terrestrial, and marine components.

Junge (1963, p. 111-201) divides microparticles up into: Aitken (≤0.1 µm), large (0.1 to 1.0 µm), and giant (>1.0 µm) sizes. Preservation of this grouping was investigated for the microparticle size grades in this study as a first approximation prior to interpretation of the data. The microparticle size grades used in this study, <0.50 to 12.70 µm, do not provide sufficient detail to differentiate Junge's (1963) finer grade, Aitken size particles, but they are divisible by comparison of their spectral properties into three size-grade groups (Fig. 2a and b): <0.50 to 0.80 µm, 0.80 to 2.00 µm, and 2.00 to 12.70 µm. This division roughly spans Junge's (1963) large and giant range with an intermediate "gray" group at 0.80 to 2.00 µm. The <0.50 to 0.80 µm microparticle size group displays power in both annual and subannual periods (Fig. 2a) while the 0.80 to 2.00 µm group has a spectral signature in which power in the subannual period dominates over power in the annual period (Fig. 2a). The 2.00 to 12.70 µm group displays power in the subannual period but none in the annual period (Fig. 2c).

As noted by Lamb (1970) and Shaw (1979), microparticles <1.00 µm are believed to have long residence times in the stratosphere and/or high troposphere. Residence time, transport distance, and potential multiple sources for these microparticles prevents determination of unique source types and/or source areas. Although no elemental analysis of the microparticles in our study was undertaken, such studies at South Pole (E.M. Thompson, 1980) reveal that a large portion of these microparticles contain iron and silicate. This is substantiated in our study by the similarity in spectral signatures of the time-series for reactive silicate and the <0.50 to 0.80 µm size grade of microparticles.

The 2.00 to 12.70 µm microparticle size group is probably derived from a more local source than the 0.50 to 0.80 µm group based purely on size. The site-specific local derivation of the 2.00 to 12.70 µm microparticles is further substantiated by their high concentration in the Nun Kun area relative to concentrations of microparticles in other high-altitude alpine areas such as Mt Kenya (L.G. Thompson, [1981] and Quelccaya, Peru (Thompson and others, 1979).

Values of δD were measured regularly every 0.16 m down core whereas δ(18O) was only measured at selected depths. To test δD to δ(18O) relationships down-core values and a suite of fresh summer snowfall samples collected in the elevation range 5130 to 5512 m (Mayewski and others, 1983) were examined. Best-fit linear regressions for δD and δ(18O) down-hole and surface samples are δD = 7.1 δ(18O) + 5.2 and δD = 8.1 × δ(18O) + 13.4, respectively. By comparison the best fit for northern hemisphere continental stations is δD = 8.0 δ(18O) + 10 (Dansgaard, 1964) and for suites of samples relatively close to Nun Kun in the central Hindu Kush and Garhwals δD = 8.1 δ(18O) + 12 and δD = 8.0 δ(18O) + 10, respectively (Niewodniczański and

others, 1981). Although it should be noted that the surface and down-hole samples in our study come from different elevations, comparison of these with fits from Dansgaard (1964) and Niewodniczański and others (1981) suggest that our down-hole samples may have experienced some degree of post-depositional alteration perhaps due to slight diffusion and/or melting. Although post-depositional alteration may have had some effect on the distribution down core of δD values it has not affected the distribution of other chemical species and the microparticles. Glaciochemical analysis of a core collected from a low elevation Norwegian glacier (Davies and others, 1982) reveals that ions migrate down core by order of solubility. The coincidence down core of reactive iron, reactive silicate, and microparticles all of low solubility, with higher solubility chloride, reactive phosphate, and ammonium suggests that post-depositional alteration by melting and/or diffusion if present is minimal at the core site.

Spectral power for the δD values is apparent for the subannual period but not annually even though monsoonal circulation in the area results in major shifts in the temperature of the air masses entering this site. The study of the fresh surface snow reveals why a major annual δD power peak may not exist since it can be used to demonstrate that during any single event two adjacent air masses may have quite different temperatures and thus individual storm events may be as effective as seasonal air mass shifts in creating δD distributions.

### Synchronicity of selected time-series sets

Time-dependent relationships between two time-series sets can be expressed using cross-spectral analysis in terms of phase shift and the strength of their dependence at various frequencies (Edwards and Thornes, 1973). Coherence between base and cross sets as close to unity as possible are sought although knowledge of the zero level of significance for coherence values (Koopmans, 1974) provides information concerning the reliability of coherence values less than unity. Phase shift (lead or lag) of selected time-series sets that are useful as annual and/or subannual indicators in this study, and that display power at statistically similar periods are examined (Table II). Since air mass circulation in the Himalaya is complicated, it cannot be assumed that all time-series sites with annual or subannual spectral power necessarily enter the core site simultaneously. However, several examples of time-series sets which do enter simultaneously are dealt with below.

Phase shift for the chemical-species time series chloride, sodium, and reactive phosphate (Table II) suggest that all three species enter the core site in-phase annually. Since the annual pulse of high chloride concentration enters the core site in summer, the cross-spectral analysis demonstrates that the annual power peaks of sodium and reactive phosphate correspond to summer events.

Phase shift for examples of the spectrally defined groups equivalent to Junge's (1963) large and giant microparticles are also presents (Table II). Based on their cross-spectral relationships the large and giant particles enter the core site in phase subannually.

## SUMMARY AND DISCUSSION

Analysis of selected chemical properties (chloride, sodium, reactive iron, reactive silicate, reactive phosphate, δD, δ(18O), total β-activity, and pH) and physical properties (density, debris band and ice lens locations, and microparticles in 15 size grades from 0.50 to 12.70 µm) reveal details concerning the net mass balance, timing of precipitation events and source of precipitation for a core collected at 4908 m on Sentik Glacier. The core is believed to have experienced minimal post-depositional alter-
<table><caption>TABLE II. CROSS-SPECTRAL RELATIONSHIPS FOR SELECTED TIME-SERIES SETS</caption><thead><tr><th>Base set</th><th>Cross set</th><th>Period<sup>1</sup> m</th><th>Amplitude of power peak<sup>2</sup></th><th>Coherence</th><th>Zero significance level for coherence<sup>3</sup></th><th>Phase<sup>4</sup> m</th></tr></thead><tbody><tr><td>Chloride</td><td>Sodium</td><td>0.81</td><td>5.03</td><td>0.85</td><td>0.66</td><td>0.05 (0.01)</td></tr><tr><td>Chloride</td><td>Sodium</td><td>1.08</td><td>5.92</td><td>0.86</td><td>0.66</td><td>1.00 (1.08)</td></tr><tr><td>Chloride</td><td>Reactive Phosphate</td><td>0.81</td><td>1.73</td><td>0.70</td><td>0.44</td><td>0.91 (0.74)</td></tr><tr><td>Chloride</td><td>Reactive Phosphate</td><td>1.08</td><td>1.70</td><td>0.70</td><td>0.44</td><td>0.91 (0.99)</td></tr><tr><td>0.50 to 0.63 µm</td><td>0.63 to 0.80 µm</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Micro-particles</td><td>Micro-particles</td><td>0.45</td><td>40 x 10<sup>6</sup></td><td>0.87</td><td>0.66</td><td>0.99 (0.45)</td></tr><tr><td>0.63 to 8.00 µm</td><td>8.00 to 10.00 µm</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Micro-particles</td><td>Micro-particles</td><td>0.45</td><td>23.90</td><td>0.98</td><td>0.66</td><td>1.00 (0.45)</td></tr></tbody></table>

1. 0.81 and 1.08 m are used to span the annual period.

2. chemical species in (µM)²/cy m⁻¹; microparticles in (µm)²/cy m⁻¹.

3. from Koopmans (1974) at 90% confidence level.

4. in fractions of a circle.

ation and thus to be a true monitor of the combined wet and dry deposition reaching this site.

Based on a total β-activity core chronology the 16.6 m core represents an ≈ 17 ± 0.3 year record. The average net mass balance computed from this chronology is 0.92 m a⁻¹ which using an average core density of 570 kg m⁻³ is 0.62 m a⁻¹ water equivalent. The latter value is somewhat less than the precipitation rate of 0.90 m a⁻¹ water equivalent estimated for the valley bottom adjacent to the core site (Qazi, 1973) suggesting a decrease in precipitation with elevation in this area.

Spectral analysis of chemical and physical time-series sets developed from the core reveal recognizable periodic events at ≲0.32 m (noise), ≲3.22 to ≳0.32 m, and ≳3.22 m. The periodicity ≲3.22 to ≳0.32 m includes annual and subannual events. The existence of spectrally defined periodic events is of particular significance in this study since they help to prove that selected high-altitude temperate glaciers in the Himalaya can provide stratigraphic records for use in the delineation of former climate. Discussion of these periodic events and the source and atmospheric circulation implications of the chemical and physical time-series are summarized below.

Chloride concentrations are highest in the summer and are derived at this time from a marine source, the Arabian Sea. Subannual chloride inputs probably come from more distant marine sources. Sodium coveries with chloride since much of it is marine-derived. However, some sodium is probably brought into the core site from crustal weathering during both the summer and winter monsoon events. Reactive phosphate is probably carried to the core site during the summer adsorbed onto wind-blown materials derived from exposed agricultural sites. Subannual inputs of this species may come from combined biologic and anthropogenic sources. Ammonium has multiple sources which are not differentiable in this study. Reactive iron and reactive silicate are both derived from crustal weathering and are deposited at the core site during at least the summer and winter monsoons and silicate comes to the core site from either a distinct limited source or over a long distance annually. pH distribution with time may be a composite monitor of strong and weak acid contributions to the site and interpre-

tation of its distribution down core may be too complicated to allow development of a useful record. Microparticles have multiple sources and are spectrally divisible into two major groups with a "gray" intermediate group. The smallest particles, ≲0.50 to 0.80 µm, have a strong annual input and a weaker sub-annual input, are believed to be long-travelled and derived from a relatively high-altitude air mass and composed of at least silicate. Microparticle size grades 2.00 to 12.70 µm are probably of more local derivation than the ≲0.50 to 0.80 µm group. δD displays a weak subannual periodicity either due to dampening by melting and/or diffusion or because temperature differences season to season are no greater than those found during subseasonal events.

Multiannual periodicity is displayed by all time-series. Events with periodicity ≳3.5 years require far more investigation because of their potential value to monsoon studies. This information could be gained by deeper drilling.

ACKNOWLEDGEMENTS

Members of the 1980 Nun Kun Expedition included P.A. Mayewski (leader), W.B. Lyons, P. Axelson, P.A. Jeschke, and J. Sevigny of the University of New Hampshire, N. Ahmad, M.N. Kaul, and K. Kahn of Aligarh Muslim University, India, and P.T. Davis of Mount Holyoke College. We would like to thank Ghulam Ali and the people of Tongul, Ladakh, for all of their help in this project. Ted Loder (University of New Hampshire) gave invaluable advice concerning our chemical studies. Isotope analyses were graciously provided by J. Jouzel and L. Merlivat (Saclay, France) and C. Lorius (Grenoble, France) and we are deeply indebted. L.G. Thompson (Ohio State University) very kindly provided all of the microparticle analyses. J. Borkland and P. Vornberger assisted in the preparation of the statistical analyses and N. Pettigrew and M. Foster provided helpful insight into statistical interpretations. Critical reviews of the manuscript were provided by C. Lorius, R.M. Koerner, and L.G. Thompson. This study is part of a joint U.S. - Indian study of Himalayan glaciers funded by U.S. National Science Foundation grant INT80-03175. A permit for this study was provided by the Government of
India and their cooperation and help in this program
is very much appreciated.

REFERENCES

Berner, W.B., and others. 1978. Past atmospheric com-
position and climate: gas parameters measured on
ice cores, by W.B. Berner, B. Stauffer, and
H. Oeschger. *Nature*, Vol. 276, No. 5683, p. 53-55.

Boucher, K. 1975. *Global climate*. New York, John Wiley and Sons, Inc.

Bryson, R.A., and Swain, M. 1981. Holocene variations of monsoonal rainfall in Rajasthan. *Quaternary Research*, Vol. 16, No. 2, p. 135-45.

Dansgaard, W. 1964. Stable isotopes in precipitation.
*Tellus*, Vol. 16, No. 4, p. 436-68.

Davidson, C.I., and others. 1981. Airborne lead and other elements derived from local fires in the Himalayas, by C.I. Davidson, T.C. Grimm, and M.A. Nasta. *Science*, Vol. 214, No. 4527, p. 1344-46.

Davies, T.D., and others. 1982. Preferential elution of strong acid from a Norwegian ice cap, by T.D. Davies, C.E. Vincent, and P. Brimblecombe. *Nature*, Vol. 300, No. 5888, p. 161-63.

Edwards, A.M.C., and Thornes, J.B. 1973. Annual cycle in river water quality: a time-series approach. *Water Resources Research*, Vol. 19, No. 5, p. 1286-95.

Glibert, P.M., and Loder, T. 1977. Automated analysis of nutrients in sea water: a manual of techniques. *Woods Hole Technical Report* 77-47.

Graham, W.F., and Duce, R.A. 1979. Atmospheric path-
way, of the phosphorus cycle. *Geochimica et Cos-
mochimica Acta*, Vol. 43, No. 8, p. 1185-1208.

Hahn, D.G., and Shukla, J. 1976. An apparent relation-
ship between Eurasian snow cover and Indian monsoon
rainfall. *Journal of the Atmospheric Sciences*,
Vol. 33, No. 12, p. 2461-62.

Hastenrath, S. 1978. Heat-budget measurements on the Quelccaya ice cap, Peruvian Andes. *Journal of Glaciology*, Vol. 20, No. 82, p. 85-97.

Ikegami, K., and others. 1980. Preliminary report on the vertical distribution of aerosol particles over the Nepal Himalayas, by K. Ikegami, J. Inoue, K. Higuchi, and A. Ono. *Seppyo*, Vol. 42, Special Issue, p. 86-89.

Junge, C.E. 1963. *Air chemistry and radioactivity.*
New York, Academic Press.

Koopmans, L.H. 1974. *The spectral analysis of time-
series.* New York, Academic Press. (Probability and
Mathematical Statistics, 22.)

Lamb, H.H. 1970. Volcanic dust in the atmosphere, with a chronology and assessment of its meteorological significance. *Philosophical Transactions of the Royal Society of London, Ser. A*, Vol. 266, No. 1178, p. 425-533.

Mayewski, P.A., and Jeschke, P.A. 1979. Himalayan and
trans-Himalayan glacier fluctuations since AD 1812.
*Arctic and Alpine Research*, Vol. 11, No. 3, p. 267-
87.

Mayewski, P.A., and others. 1981. Reconnaissance
glacio-chemical studies in the Indian Himalayas, by
P.A. Mayewski, W.B. Lyons, and N. Ahmad. *Proceed-
ings of the Eastern Snow Conference*, 38th annual
meeting, p. 45-48.

Mayewski, P.A., and others. 1983. Chemical composition of a high altitude fresh snowfall in the Ladakh Himalayas, by P.A. Mayewski, W.B. Lyons, and N. Ahmad. *Geophysical Research Letters*, Vol. 10, No. 1, p. 105-08.

Mukherjee, B.K., and others. 1981. High-level tempera-
tures and winds over tropics and Indian summer mon-
soon, by B.K. Mukherjee, R.S. Reddy, and
Bh. V. Ramana Murty. *Journal of Geophysical Research*,
Vol. 86, No. C10, p. 9688-92.

Niewodniczański, J., and others. 1981. The altitude effect on the isotopic composition of snow in high mountains, by J. Niewodniczański, J. Grabczak, L. Barański, and J. Rzepka. *Journal of Glaciology*, Vol. 27, No. 95, p. 99-111.

Oeschger, H., and others. 1982. Atmospheric CO₂ con-
tent in the past deduced from ice-core analyses, by
H. Oeschger, B. Stauffer, A. Neftel, J. Schwander,
and R. Zumbrunn. *Annals of Glaciology*, Vol. 3,
p. 227-32.

Qazi, N.A. 1973. *Water resources development in Suru Basin (Ladakh)*. Srinagar, Vacant Press.

Rangarajan, C., and Eapen, C.D. 1981. Estimates of interhemispheric transport of radioactive debris by the East African low-level jet stream. *Journal of Geophysical Research*, Vol. 86, No. C12, p. 12153-54.

Reiter, E.R. 1981. The Tibet connection. *Natural History*, Vol. 90, No. 1, p. 64-71.

Scholander, P.F., and others. 1961. Composition of gas bubbles in Greenland icebergs, by P.F. Scholander, E.A. Hemmingsen, L.K. Coachman, and D.C. Nutt. *Journal of Glaciology*, Vol. 3, No. 29, p. 813-22.

Sequeira R., and Keilkar, D. 1978. Geochemical impli-
cations of summer monsoonal rainwater composition
over India. *Journal of Applied Meteorology*, Vol. 17,
No. 9, p. 1390-96.

Shaw, G.E. 1979. Considerations on the origin and
properties of the Antarctic aerosol. *Reviews of
Geophysics and Space Physics*, Vol. 17, No. 8,
p. 1983-98.

Stauffer, B., and Berner, W. 1978. CO₂ in natural ice.
*Journal of Glaciology*, Vol. 21, No. 85, p. 291-300.

Stookey, L.L. 1970. Ferrozine - a new spectrophoto-
metric reactant from iron. *Analytical Chemistry*,
Vol. 42, No. 7, p. 778-881.

Thompson, E.M. 1980. 911 years of microparticle de-
position at the South Pole: a climatic interpre-
tation. *Ohio State University. Institute of Polar
Studies. Report No. 73.*

Thompson, L.G. 1977. Microparticles, ice sheets and climates. *Ohio State University. Institute of Polar Studies. Report No. 64.*

Thompson, L.G. [1981.] Ice core studies from Mt Kenya, Africa, and their relationships to other tropical ice core studies. [*Union Géodésique et Géophysique Internationale. Association Internationale des Sciences Hydrologiques.*] Sea level, ice, and climatic change. *Proceedings of the symposium held 7-8 December 1979 during the 17th general assembly of the International Union of Geodesy and Geophysics, Canberra, p. 55-62.* (IAHS Publication No. 131.)

Thompson, L.G., and others. 1979. Climatic ice core records from the tropical Quelccaya ice cap, by L.G. Thompson, S. Hastenrath, and B. Morales Arnao. *Science*, Vol. 203, No. 4386, p. 1240-43.

Wissmann, H. von. 1960. Die heutige Vergletscherung und Schneegrenze in Hochasien, mit Hinweisen auf die Vergletscherung der letzten Eiszeit. *Akademie der Wissenschaften und der Literatur in Mainz. Abhandlungen der Mathematisch-naturwissenschaftlichen Klasse*. Jahrg. 1959, Nr. 14, p. 1103-1407.

Zall, D.M., and others. 1956. Photometric determina-
tion of chloride in water, by D.M. Zall, D. Fisher,
and M.D. Garner. *Analytical Chemistry*, Vol. 28,
No. 11, p. 1665-68.
APPENDIX A. DOWN-HOLE DATA SETS

Symbol

Data set (units)

BETA = Total β-activity (dph kg⁻¹)

DEBRIS AND ICE BANDS = Debris bands (2 = diffuse; 4 = semi-concentrated) and ice bands (= 1)

DENSITY = Density (x 10³kg m⁻³)

CL = Chloride (µM); note scale change

NA = Sodium (µM); note scale change

FE = Reactive iron (µM); note two scales

SI02 = Reactive silicate (µM); note two scales

Po4 = Reactive phosphate (µM); note two scales

NH4 = Ammonium (µM)

DELTA D = δD (°/oo)

DELTA 018 = δ(¹⁸O) (°/oo)

pH = pH (pH units)

MICRQ = 0.63 x 10³ = 0.50-0.63 x 10³ = Microparticles 0.63 µ x 10³ (particles/500 µl of sample)

MICRQ = 1.25 x 10² = 1.00-1.25 x 10² = Microparticles 1.25 µ x 10² (particles/500 µl of sample)

MICRQ = 10.0 = 8.00-10.00 = Microparticles 10.0 µm (particles/500 µl of sample)

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

MS. received 7 April 1983 and in revised form 4 August 1983