Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysics Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3Linked to original sources

Modeling and measuring the nocturnal drainage flow in a high-elevation, subalpine forest with complex terrain

The nocturnal drainage flow of air causes significant uncertainty in ecosystem CO 2 , H 2 O, and energy budgets determined with the eddy covariance measurement approach. In this study, we examined the magnitude, nature, and dynamics of the nocturnal drainage flow in a subalpine forest ecosystem with complex terrain. We used an experimental approach involving four towers, each with vertical profiling of wind speed to measure the magnitude of drainage flows and dynamics in their occurrence. We developed an analytical drainage flow model, constrained with measurements of canopy structure and SF 6 diffusion, to help us interpret the tower profile results. Model predictions were in good agreement with observed profiles of wind speed, leaf area density, and wind drag coefficient. Using theory, we showed that this one‐dimensional model is reduced to the widely used exponential wind profile model under conditions where vertical leaf area density and drag coefficient are uniformly distributed. We used the model for stability analysis, which predicted the presence of a very stable layer near the height of maximum leaf area density. This stable layer acts as a flow impediment, minimizing vertical dispersion between the subcanopy air space and the atmosphere above the canopy. The prediction is consistent with the results of SF 6 diffusion observations that showed minimal vertical dispersion of nighttime, subcanopy drainage flows. The stable within‐canopy air layer coincided with the height of maximum wake‐to‐shear production ratio. We concluded that nighttime drainage flows are restricted to a relatively shallow layer of air beneath the canopy, with little vertical mixing across a relatively long horizontal fetch. Insight into the horizontal and vertical structure of the drainage flow is crucial for understanding the magnitude and dynamics of the mean advective CO 2 flux that becomes significant during stable nighttime conditions and are typically missed during measurement of the turbulent CO 2 flux. The model and interpretation provided in this study should lead to research strategies for the measurement of these advective fluxes and their inclusion in the overall mass balance for CO 2 at this site with complex terrain.

Journal of Geophysical Research D: Atmospheres↗

Enhanced visualization for the interpretation of Magellan radar data: Supplement to the Magellan special issue

The differences of radar data from more familiar photographic images, coupled with the unique geologic characteristics of Venus's surface, can make interpretation of the standard Magellan data products difficult for both planetary scientists and nonspecialists. We describe a set of digital processing techniques for transforming individual remote‐sensing datasets in order to make the information they contain more apparent, and for combining multiple datasets of diverse resolution and information content into synthetic images that allow the viewer to explore correlations among the data. The creation of derived data products by techniques like those we describe is becoming increasingly widespread and important, both for scientific research and for use in educational and media presentations. We apply our processing techniques to Magellan synthetic aperture radar images, altimetry, and microwave emissivity measurements of four areas that exemplify some of the most important styles of geologic activity on Venus.

Journal of Geophysical Research E: Planets↗

Standardizing the nomenclature of Martian impact crater ejecta morphologies

The Mars Crater Morphology Consortium recommends the use of a standardized nomenclature system when discussing Martian impact crater ejecta morphologies. The system utilizes nongenetic descriptors to identify the various ejecta morphologies seen on Mars. This system is designed to facilitate communication and collaboration between researchers. Crater morphology databases will be archived through the U.S. Geological Survey in Flagstaff, where a comprehensive catalog of Martian crater morphologic information will be maintained.

Journal of Geophysical Research E: Planets↗

Modern thermokarst lake dynamics in the continuous permafrost zone, northern Seward Peninsula, Alaska

Quantifying changes in thermokarst lake extent is of importance for understanding the permafrost-related carbon budget, including the potential release of carbon via lake expansion or sequestration as peat in drained lake basins. We used high spatial resolution remotely sensed imagery from 1950/51, 1978, and 2006/07 to quantify changes in thermokarst lakes for a 700 km 2 area on the northern Seward Peninsula, Alaska. The number of water bodies larger than 0.1 ha increased over the entire observation period (666 to 737 or +10.7%); however, total surface area decreased (5,066 ha to 4,312 ha or -14.9%). This pattern can largely be explained by the formation of remnant ponds following partial drainage of larger water bodies. Thus, analysis of large lakes (>40 ha) shows a decrease of 24% and 26% in number and area, respectively, differing from lake changes reported from other continuous permafrost regions. Thermokarst lake expansion rates did not change substantially between 1950/51 and 1978 (0.35 m/yr) and 1978 and 2006/07 (0.39 m/yr). However, most lakes that drained did expand as a result of surface permafrost degradation before lateral drainage. Drainage rates over the observation period were stable (2.2 to 2.3 per year). Thus, analysis of decadal-scale, high spatial resolution imagery has shown that lake drainage in this region is triggered by lateral breaching and not subterranean infiltration. Future research should be directed toward better understanding thermokarst lake dynamics at high spatial and temporal resolution as these systems have implications for landscape-scale hydrology and carbon budgets in thermokarst lake-rich regions in the circum-Arctic.

Alaska↗

Ice core evidence of rapid air temperature increases since 1960 in alpine areas of the Wind River Range, Wyoming, United States

Site-specific transfer functions relating delta oxygen 18 (δ 18 O) values in snow to the average air temperature ( T A ) during storms on Upper Fremont Glacier (UFG) were used in conjunction with δ 18 O records from UFG ice cores to reconstruct long-term trends in air temperature from alpine areas in the Wind River Range, Wyoming. Transfer functions were determined by using data collected from four seasonal snowpacks (1989–1990, 1997–1998, 1998–1999, and 1999–2000). The timing and amount of each storm was determined from an automated snowpack telemetry (SNOTEL) site, 22 km northeast of UFG, and ∼1060 m in elevation below UFG. Statistically significant and positive correlations between δ 18 O values in the snow and T A were consistently found in three of the four seasonal snowpacks. The snowpack with the poor correlation was deposited in 1997–1998 during the 1997–1998 El Niño Southern Oscillation (ENSO). An ultrasonic snow-depth sensor installed on UFG provided valuable insights into site-specific storms and postdepositional processes that occur on UFG. The timing of storms recorded at the UFG and Cold Springs SNOTEL sites were similar; however, selected storms did not correlate. Snow from storms occurring after mid-October and followed by high winds was most susceptible to redeposition of snow. This removal of lower temperature snowfall could potentially bias the δ 18 O values preserved in ice core records to environmental conditions reflecting higher air temperatures and lower wind speeds. Transfer functions derived from seasonal snow cover on UFG were used to reconstruct T A values from δ 18 O values determined from two ice cores collected from UFG. Reconstructed air temperatures from the ice core data indicate an increase in T A of ∼3.5°C from the mid-1960s to the early 1990s in the alpine areas of northwestern Wyoming. Reconstructed T A from the ice core records between the end of the Little Ice Age (LIA), mid-1800s, and the early 1990s indicate a T A increase of ∼5°C. The historically reconstructed T A values from the UFG were significantly higher than the global average observed during the 20th Century but were in agreement with T A increases observed at selected, high-altitude and high-latitude sites in other parts of the world. Additional research is required to determine if part of the observed trend toward heavier δ 18 O values in ice from the UFG since the LIA (and increased T A ) is due to an increased proportion of snowfall from southerly storm tracks and moisture sources, as seems to have been the situation in 1997–1998.

Wyoming↗

Sediment and nutrient delivery from thermokarst features in the foothills of the North Slope, Alaska: Potential impacts on headwater stream ecosystems

Permafrost is a defining characteristic of the Arctic environment. However, climate warming is thawing permafrost in many areas leading to failures in soil structure called thermokarst. An extensive survey of a 600 km2 area in and around the Toolik Lake Natural Research Area (TLNRA) revealed at least 34 thermokarst features, two thirds of which were new since ???1980 when a high resolution aerial survey of the area was done. Most of these thermokarst features were associated with headwater streams or lakes. We have measured significantly increased sediment and nutrient loading from thermokarst features to streams in two well-studied locations near the TLNRA. One small thermokarst gully that formed in 2003 on the Toolik River in a 0.9 km2 subcatchment delivered more sediment to the river than is normally delivered in 18 years from 132 km2 in the adjacent upper Kuparuk River basin (a long-term monitoring reference site). Ammonium, nitrate, and phosphate concentrations downstream from a thermokarst feature on Imnavait Creek increased significantly compared to upstream reference concentrations and the increased concentrations persisted over the period of sampling (1999-2005). The downstream concentrations were similar to those we have used in a long-term experimental manipulation of the Kuparuk River and that have significantly altered the structure and function of that river. A subsampling of other thermokarst features from the extensive regional survey showed that concentrations of ammonium, nitrate, and phosphate were always higher downstream of the thermokarst features. Our previous research has shown that even minor increases in nutrient loading stimulate primary and secondary production. However, increased sediment loading could interfere with benthic communities and change the responses to increased nutrient delivery. Although the terrestrial area impacted by thermokarsts is limited, the aquatic habitat altered by these failures can be extensive. If warming in the Arctic foothills accelerates thermokarst formation, there may be substantial and wide-spread impacts on arctic stream ecosystems that are currently poorly understood. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗

Understanding the amplitudes of noise correlation measurements

Cross correlation of ambient seismic noise is known to result in time series from which station-station travel-time measurements can be made. Part of the reason that these cross-correlation travel-time measurements are reliable is that there exists a theoretical framework that quantifies how these travel times depend on the features of the ambient noise. However, corresponding theoretical results do not currently exist to describe how the amplitudes of the cross correlation depend on such features. For example, currently it is not possible to take a given distribution of noise sources and calculate the cross correlation amplitudes one would expect from such a distribution. Here, we provide a ray-theoretical framework for calculating cross correlations. This framework differs from previous work in that it explicitly accounts for attenuation as well as the spatial distribution of sources and therefore can address the issue of quantifying amplitudes in noise correlation measurements. After introducing the general framework, we apply it to two specific problems. First, we show that we can quantify the amplitudes of coherency measurements, and find that the decay of coherency with station-station spacing depends crucially on the distribution of noise sources. We suggest that researchers interested in performing attenuation measurements from noise coherency should first determine how the dominant sources of noise are distributed. Second, we show that we can quantify the signal-to-noise ratio of noise correlations more precisely than previous work, and that these signal-to-noise ratios can be estimated for given situations prior to the deployment of seismometers. It is expected that there are applications of the theoretical framework beyond the two specific cases considered, but these applications await future work.

Journal of Geophysical Research↗

Surface temperature patterns in complex terrain: Daily variations and long-term change in the central Sierra Nevada, California

A realistic description of how temperatures vary with elevation is crucial for ecosystem studies and for models of basin-scale snowmelt and spring streamflow. This paper explores surface temperature variability using temperature data from an array of 37 sensors, called the Yosemite network, which traverses both slopes of the Sierra Nevada in the vicinity of Yosemite National Park, California. These data indicate that a simple lapse rate is often a poor description of the spatial temperature structure. Rather, the spatial pattern of temperature over the Yosemite network varies considerably with synoptic conditions. Empirical orthogonal functions (EOFs) were used to identify the dominant spatial temperature patterns and how they vary in time. Temporal variations of these surface temperature patterns were correlated with large-scale weather conditions, as described by National Centers for Environmental Prediction-National Center for Atmospheric Research Reanalysis data. Regression equations were used to downscale larger-scale weather parameters, such as Reanalysis winds and pressure, to the surface temperature structure over the Yosemite network. These relationships demonstrate that strong westerly winds are associated with relatively warmer temperatures on the east slope and cooler temperatures on the west slope of the Sierra, and weaker westerly winds are associated with the opposite pattern. Reanalysis data from 1948 to 2005 indicate weakening westerlies over this time period, a trend leading to relatively cooler temperatures on the east slope over decadal timescale's. This trend also appears in long-term observations and demonstrates the need to consider topographic effects when examining long-term changes in mountain regions. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Introduction to special section on the California-Arizona crustal transect: Part II

In May 1988, a multidisciplinary conference was held in Flagstaff, Arizona, to synthesize recent research on the crustal structural and evolution of a transect across the southern Cordillera. This transect extends from the southern Colorado Plateau southwestward across the Arizona Transition Zone, through the Basin and Range province in the Mojave and Sonoran deserts, to the San Andreas fault system, and thence westward across the Imperial Valley and Peninsular Ranges and into the southern California Borderland. Most discussion concerned the eastern part of the transect, from the Colorado Plateau to the San Andreas fault [ Sass et al ., 1988; Howard et al ., 1990]. For this meeting, the name CACTIS (California‐Arizona Crustal Transect‐Interim Synthesis) was coined by J. Sass.

Journal of Geophysical Research↗

Integrating perspectives to understand lake ice dynamics in a changing world

Ice cover plays a critical role in physical, biogeochemical, and ecological processes in lakes. Despite its importance, winter limnology remains relatively understudied. Here, we provide a primer on the predominant drivers of freshwater lake ice cover and the current methodologies used to study lake ice, including in situ and remote sensing observations, physical based models, and experiments. We highlight opportunities for future research by integrating these four disciplines to address key knowledge gaps in our understanding of lake ice dynamics in changing winters. Advances in technology, data integration, and interdisciplinary collaboration will allow the field to move toward developing global forecasts of lake ice cover for small to large lakes across broad spatial and temporal scales, quantifying ice quality and ice thickness, moving from binary to continuous ice records, and determining how winter ice conditions and quality impact ecosystem processes in lakes over winter. Ultimately, integrating disciplines will improve our ability to understand the impacts of changing winters on lake ice.

Journal of Geophysical Research: Biogeosciences↗

Coseismic deformation of the Wrights tunnel during the 1906 San Francisco earthquake: A key to understanding 1906 fault slip and 1989 surface ruptures in the southern Santa Cruz Mountains, California

The Wrights tunnel is an abandoned railroad tunnel that crosses the San Andreas fault in the southern Santa Cruz Mountains in the vicinity of the 1989 Loma Prieta earthquake. The tunnel was damaged and deformed during the 1906 San Francisco earthquake and a plot showing postearthquake measurements made in the tunnel is given by Lawson [1908]. The amount of offset shown on this plot (1.5 m) has been used in several studies as being representative of the amount of fault offset along this segment of the San Andreas fault in 1906. Our historical research shows that different observers reported different amounts of fault offset in the tunnel and that the 1.5 m given on the plot is not a surveyed measurement. In addition, the plot of the tunnel has been interpreted in several previous studies as evidence of a broad (1.5 km) zone of faulting beneath Summit Ridge. Our analysis shows that this plot need not indicate a broad zone of deformation. Our historical research and modeling of the tunnel measurements indicate that faulting was confined to a zone less than 400 m wide and that 60–85% of the coseismic slip occurred across a single fault plane. There is no evidence for offset across a second shear zone beneath Summit Ridge in 1906. This implies that surface fractures reported on Summit Ridge in 1906 were not associated with significant deformation of the tunnel, implying that they were shallow, surficial features. By analogy, the very similar fractures that occurred on Summit Ridge in 1989 were also probably the result of shallow gravitational, rather than deep-seated tectonic, processes. Our modeling also indicates that total coseismic, near-surface slip across the San Andreas fault zone in the Wrights tunnel in 1906 was at least 1.7–1.8 m.

Journal of Geophysical Research B: Solid Earth↗

Arctic transitions in the Land - Atmosphere System (ATLAS): Background, objectives, results, and future directions

This paper briefly reviews the background, objectives, and results of the Arctic Transitions in the Land-Atmosphere System (ATLAS) Project to date and provides thoughts on future directions. The key goal of the ATLAS Project is to improve understanding of controls over spatial and temporal variability of terrestrial processes in the Arctic that have potential consequences for the climate system, i.e., processes that affect the exchange of water and energy with the atmosphere, the exchange of radiatively active gases with the atmosphere, and the delivery of freshwater to the Arctic Ocean. Three important conclusions have emerged from research associated with the ATLAS Project. First, associated with the observation that the Alaskan Arctic has warmed significantly in the last 30 years, permafrost is warming, shrubs are expanding, and there has been a temporary release of carbon dioxide from tundra soils. Second, the winter is a more important period of biological activity than previously appreciated. Biotic processes, including shrub expansion and decomposition, affect snow structure and accumulation and affect the annual carbon budget of tundra ecosystems. Third, observed vegetation changes can have a significant positive feedback to regional warming. These vegetation effects are, however, less strong than those exerted by land-ocean heating contrasts and the topographic constraints on air mass movements. The papers of this special section provide additional insights related to these conclusions and to the overall goal of ATLAS.

Journal of Geophysical Research D: Atmospheres↗

A detached eddy simulation model for the study of lateral separation zones along a large canyon-bound river

Lateral flow separation occurs in rivers where banks exhibit strong curvature. In canyon-boundrivers, lateral recirculation zones are the principal storage of fine-sediment deposits. A parallelized,three-dimensional, turbulence-resolving model was developed to study the flow structures along lateralseparation zones located in two pools along the Colorado River in Marble Canyon. The model employs thedetached eddy simulation (DES) technique, which resolves turbulence structures larger than the grid spacingin the interior of the flow. The DES-3D model is validated using Acoustic Doppler Current Profiler flowmeasurements taken during the 2008 controlled flood release from Glen Canyon Dam. A point-to-pointvalidation using a number of skill metrics, often employed in hydrological research, is proposed here forfluvial modeling. The validation results show predictive capabilities of the DES model. The model reproducesthe pattern and magnitude of the velocity in the lateral recirculation zone, including the size and position ofthe primary and secondary eddy cells, and return current. The lateral recirculation zone is open, havingcontinuous import of fluid upstream of the point of reattachment and export by the recirculation returncurrent downstream of the point of separation. Differences in magnitude and direction of near-bed andnear-surface velocity vectors are found, resulting in an inward vertical spiral. Interaction between therecirculation return current and the main flow is dynamic, with large temporal changes in flow direction andmagnitude. Turbulence structures with a predominately vertical axis of vorticity are observed in the shearlayer becoming three-dimensional without preferred orientation downstream.

Arizona↗

Chronological refinement of an ice core record at Upper Fremont Glacier in south central North America

The potential to use ice cores from alpine glaciers in the midlatitudes to reconstruct paleoclimatic records has not been widely recognized. Although excellent paleoclimatic records exist for the polar regions, paleoclimatic ice core records are not common from midlatitude locations. An ice core removed from the Upper Fremont Glacier in Wyoming provides evidence for abrupt climate change during the mid-1800s. Volcanic events (Krakatau and Tambora) identified from electrical conductivity measurements (ECM) and isotopic and chemical data from the Upper Fremont Glacier were reexamined to confirm and refine previous chronological estimates of the ice core. At a depth of 152 m the refined age-depth profile shows good agreement (1736 ± 10 A.D.) with the 14C age date (1729 ± 95 A.D.). The δ 18 O profile of the Upper Fremont Glacier (UFG) ice core indicates a change in climate known as the Little Ice Age (LIA). However, the sampling interval for δ 18 O is sufficiently large (20 cm) such that it is difficult to pinpoint the LIA termination on the basis of δ 18 O data alone. Other research has shown that changes in the δ18O variance are generally coincident with changes in ECM variance. The ECM data set contains over 125,000 data points at a resolution of 1 data point per millimeter of ice core. A 999-point running average of the ECM data set and results from f tests indicates that the variance of the ECM data decreases significantly at about 108 m. At this depth, the age-depth profile predicts an age of 1845 A.D. Results indicate the termination of the LIA was abrupt with a major climatic shift to warmer temperatures around 1845 A.D. and continuing to present day. Prediction limits (error bars) calculated for the profile ages are ±10 years (90% confidence level). Thus a conservative estimate for the time taken to complete the LIA climatic shift to present-day climate is about 10 years, suggesting the LIA termination in alpine regions of central North America may have occurred on a relatively short (decadal) timescale. Copyright 2000 by the American Geophysical Union.

Wyoming↗

Volcano monitoring using the Global Positioning System: Filtering strategies

Permanent Global Positioning System (GPS) networks are routinely used for producing improved orbits and monitoring secular tectonic deformation. For these applications, data are transferred to an analysis center each day and routinely processed in 24-hour segments. To use GPS for monitoring volcanic events, which may last only a few hours, real-time or near real-time data processing and subdaily position estimates are valuable. Strategies have been researched for obtaining station coordinates every 15 min using a Kalman filter; these strategies have been tested on data collected by a GPS network on Kilauea Volcano. Data from this network are tracked continuously, recorded every 30 s, and telemetered hourly to the Hawaiian Volcano Observatory. A white noise model is heavily impacted by data outages and poor satellite geometry, but a properly constrained random walk model fits the data well. Using a borehole tiltmeter at Kilauea's summit as ground-truth, solutions using different random walk constraints were compared. This study indicates that signals on the order of 5 mm/h are resolvable using a random walk standard deviation of 0.45 cm/√h. Values lower than this suppress small signals, and values greater than this have significantly higher noise at periods of 1–6 hours.

Hawai'i↗

Archie’s saturation exponent for natural gas hydrate in coarse-grained reservoirs

Accurately quantifying the amount of naturally occurring gas hydrate in marine and permafrost environments is important for assessing its resource potential and understanding the role of gas hydrate in the global carbon cycle. Electrical resistivity well logs are often used to calculate gas hydrate saturations, S h , using Archie's equation. Archie's equation, in turn, relies on an empirical saturation parameter, n . Though n = 1.9 has been measured for ice‐bearing sands and is widely used within the hydrate community, it is highly questionable if this n value is appropriate for hydrate‐bearing sands. In this work, we calibrate n for hydrate‐bearing sands from the Canadian permafrost gas hydrate research well, Mallik 5L‐38, by establishing an independent downhole S h profile based on compressional‐wave velocity log data. Using the independently determined S h profile and colocated electrical resistivity and bulk density logs, Archie's saturation equation is solved for n, and uncertainty is tracked throughout the iterative process. In addition to the Mallik 5L‐38 well, we also apply this method to two marine, coarse‐grained reservoirs from the northern Gulf of Mexico Gas Hydrate Joint Industry Project: Walker Ridge 313‐H and Green Canyon 955‐H. All locations yield similar results, each suggesting n ≈ 2.5 ± 0.5. Thus, for the coarse‐grained hydrate bearing ( S h > 0.4) of greatest interest as potential energy resources, we suggest that n = 2.5 ± 0.5 should be applied in Archie's equation for either marine or permafrost gas hydrate settings if independent estimates of n are not available.

Journal of Geophysical Research B: Solid Earth↗

Testing the accuracy of a 1-D volcanic plume model in estimating mass eruption rate

During volcanic eruptions, empirical relationships are used to estimate mass eruption rate from plume height. Although simple, such relationships can be inaccurate and can underestimate rates in windy conditions. One-dimensional plume models can incorporate atmospheric conditions and give potentially more accurate estimates. Here I present a 1-D model for plumes in crosswind and simulate 25 historical eruptions where plume height H obs was well observed and mass eruption rate M obs could be calculated from mapped deposit mass and observed duration. The simulations considered wind, temperature, and phase changes of water. Atmospheric conditions were obtained from the National Center for Atmospheric Research Reanalysis 2.5° model. Simulations calculate the minimum, maximum, and average values ( M min , M max , and M avg ) that fit the plume height. Eruption rates were also estimated from the empirical formula M empir = 140 H obs 4.14 ( M empir is in kilogram per second, H obs is in kilometer). For these eruptions, the standard error of the residual in log space is about 0.53 for M avg and 0.50 for M empir . Thus, for this data set, the model is slightly less accurate at predicting M obs than the empirical curve. The inability of this model to improve eruption rate estimates may lie in the limited accuracy of even well-observed plume heights, inaccurate model formulation, or the fact that most eruptions examined were not highly influenced by wind. For the low, wind-blown plume of 14–18 April 2010 at Eyjafjallajökull, where an accurate plume height time series is available, modeled rates do agree better with M obs than M empir .

Journal of Geophysical Research D: Atmospheres↗

Association between atmospheric circulation patterns and firn-ice core records from the Inilchek glacierized area, central Tien Shan, Asia

Glacioclimatological research in the central Tien Shan was performed in the summers of 1998 and 1999 on the South Inilchek Glacier at 5100-5460 m. A 14.36 m firn-ice core and snow samples were collected and used for stratigraphic, isotopic, and chemical analyses. The firn-ice core and snow records were related to snow pit measurements at an event scale and to meteorological data and synoptic indices of atmospheric circulation at annual and seasonal scales. Linear relationships between the seasonal air temperature and seasonal isotopic composition in accumulated precipitation were established. Changes in the ??18O air temperature relationship, in major ion concentration and in the ratios between chemical species, were used to identify different sources of moisture and investigate changes in atmospheric circulation patterns. Precipitation over the central Tien Shan is characterized by the lowest ionic content among the Tien Shan glaciers and indicates its mainly marine origin. In seasons of minimum precipitation, autumn and winter, water vapor was derived from the and and semiarid regions in central Eurasia and contributed annual maximal solute content to snow accumulation in Tien Shan. The lowest content of major ions was observed in spring and summer layers, which represent maximum seasonal accumulation when moisture originates over the Atlantic Ocean and Mediterranean and Black Seas. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗