bluuebunny/crossref_metadata_2025 huggingface dataset:
Turbulent surface fluxes and air–ice coupling in the Baltic Air–Sea–Ice Study (BASIS)
Authors: Launiainen Jouko , Cheng Bin , Uotila Juha , Vihma Timo
Published: November 2007

Turbulent surface fluxes were studied using observations taken over sea ice in the Baltic Sea in March 1998. The fluxes of momentum and sensible heat were measured by a sonic anemometer and compared with fluxes derived from wind velocity and air-temperature profiles. The neutral 10 m drag coefficient showed no apparent dependence on wind speed (in the range 2–20 m s–1), resulting in a mean value of 1.0 × 10–3 for smooth snow-covered ice and 1.5 × 10−3 for deformed ice. The overall mean value was 1.28 × 10–3. The roughness length for temperature revealed a greater apparent dependence on wind speed and was slightly larger than the aerodynamic roughness for low wind speeds, and vice versa for moderate and high winds. We give an empirical expression that predicts how the scalar roughness depends on the aerodynamic roughness (drag coefficient) and wind speed. Agreement of the gradient-method results with the eddy-flux results supports the validity of the Monin-Obukhov similarity theory. Fluxes modelled by a coupled air-ice-sea model compared well with the eddy-flux and gradient methods. Surface temperature estimates by the three methods also agreed well. Tests and sensitivity analysis emphasize the need for especially accurate sensor calibration and strict information about the sensor heights for the gradient method.
Structure of sea ice in McMurdo Sound, Antarctica
Authors: Jones Stephen J. , Hill Brian T.
Published: November 2007

Sea-ice cores from 11 sites in McMurdo Sound, Antarctica, were collected in 1982 and their crystallography examined. All but one were first-year sea ice. The cores, approximately 2 m long, consisted typically of a thin layer of granular or snow ice (approximately 0.1 m) followed by columnar-grained ice in the top metre and platelet ice in the bottom metre. Salinity and temperature measurements are reported. The columnar-grained ice usually had a strong preferred c-axis orientation in the horizontal plane and also showed a change in this preferred direction with depth in the ice. The mean c-axis orientation, however, usually aligned well with measured or implied currents in the Sound. The platelets were usually aligned with c axis horizontal or close to horizontal, and did not exhibit as marked a preferred orientation as the columnar-grained ice.
Platelet ice and the land-fast sea ice of McMurdo Sound, Antarctica
Authors: Smith Inga J. , Langhorne Patricia J. , Haskell Timothy G. , Joe Trodahl H. , Frew Russell , Ross Vennell M.
Published: November 2007

Dendritic crystals of platelet ice appear beneath the columnar land-fast sea ice of McMurdo Sound, Antarctica. These leaf-like crystals are frozen into place by the advancing columnar growth. The platelets most probably begin to appear during July although in some parts of the Sound they may not appear at all. In addition, the amount and extent of platelet ice within the Sound varies from year to year. Previous authors have suggested that the formation of platelet ice is linked to the presence of the nearby ice shelf. It is a matter of debate whether these platelets form at depth and then float upwards or whether they grow in slightly supercooled water at the ice/water interface. The phenomenon is similar to that observed in the Weddell Sea region, but previous authors have suggested the two regions may experience different processes. This paper presents the results of field-work conducted in McMurdo Sound in 1999. Ice-structure analysis, isotopic analysis and salinity and temperature measurements near the ice/water interface are presented. Freezing points are calculated, and the possible existence of supercooling is discussed in relation to existing conjectures about the origin of platelets.
Characteristics of the surface heat budget during the ice-growth season in the southern Sea of Okhotsk
Authors: Toyota Takenobu , Wakatsuchi Masaaki
Published: November 2007

The heat budget over the ice-covered area of the southern Sea of Okhotsk is estimated from in situ meteorological and ice observation for four years, 1996−99. The data are from about 1 week in early February in each of four years. Ice-thickness distributions required for calculating the heat budget are quantitatively obtained from video analysis. A one-dimensional thermodynamical model is used to calculate the heat flux. The total heat flux is obtained by summing up the area-weighted heat flux of each ice-thickness category. In addition, to determine the characteristics of the heat budget in this region, we also calculated the heat budget in the northern Sea of Okhotsk using Special Sensor Microwave/Imager ice-extent data and European Centre for Medium-range Weather Forecasts meteorological data, and compared the results. Our investigations show the following characteristics in the southern Sea of Okhotsk: (1) Due to relatively thin ice thickness, the average turbulent heat flux is upward. (2) Thin ice and open water contribute significantly to the total turbulent heat flux. (3) Thermodynamic ice growth is limited to about 1 cm d−1 on average. (4) The heat budget is largely characterized by abundant solar radiation. The first, third and fourth results are characteristic of this region located at a relatively low latitude, while the second one is similar to that for polar regions.
Combined satellite- and ULS-derived sea-ice flux in the Weddell Sea, Antarctica
Authors: Drinkwater Mark R. , Liu Xiang , Harms Sabine
Published: November 2007

Several years of daily microwave satellite ice drift are combined with moored upward-looking sonar (ULS) ice drafts into an ice-volume flux record at points along a flux gate across the Weddell Sea, Antarctica. Monthly ice transport varies at the mooring locations from a maximum export of 0.4 m3 s−1 near joinville Island to −0.4 m3 s−1 imported along the Fimbul and Riiser-Larsen ice-shelf margins. Winter peaks are observed at each end of the flux gate, where high concentrations of deformed ice are advected in and out of the basin along the coastline. The central gyre, in contrast, exhibits negligible seasonality and much smaller volume transports. During the period of overlapping ULS operation, the mean monthly integrated ice export west of the gyre center is 59 × 103 m3 s−1, and the import in the East Wind Drift is −17 x103 m3s-1. ULS data are compared with ERS satellite observation of radar backscatter to obtain an empirical relationship between ice thickness and the rate of change of backscatter with incidence angle. Resulting proxy ice-thickness data are combined with Special Sensor Microwave/Imager-derived ice velocities to obtain seasonally varying estimates of net ice-volume flux for the period 1992−98. Significant interannual variability is observed in ice-volume flux expressed as fresh-water transport. A maximum annual mean of 0.054 Sv is observed in 1992; with a minimum of 0.015 Sv in 1996. A 6 year mean transport of 0.032 Sv is observed. Maximum seasonal ice export occurs in July 1992, with a minimum in November 1996. The 10 year mean area flux is 30 × 103 m2s–1 Interannual variations in net volume flux closely follow variations in area flux, with summer minima in 1990/91 and 1996/97. Maximum area transport occurs in 1991, and although this predates the ERS-1 scatterometer data, ice-thickness estimates by Harms and others confirm 1991 as a decadal peak in net integrated fresh-water transport.
Decadal variability in high northern latitudes as simulated by an intermediate-complexity climate model
Authors: Goosse H. , Selten F. M. , Haarsma R. J. , Opsteegh J. D.
Published: November 2007

A 2500 year integration has been performed with a global coupled atmospheric-sea-ice-ocean model of intermediate complexity with the main objective of studying the climate variability in polar regions on decadal time-scales and longer. The atmospheric component is the ECBILT model, a spectral T21 three-level quasi-geostrophic model that includes a representation of horizontal and vertical heat transfers as well as of the hydrological cycle. ECBILT is coupled to the CLIO model, which consists of a primitive-equation free-surface ocean general circulation model and a dynamic-thermodynamic sea-ice model. Comparison of model results with observations shows that the ECBILT-CLIO model is able to reproduce reasonably well the climate of the high northern latitudes. The dominant mode of coupled variability between the atmospheric circulation and sea-ice cover in the simulation consists of an annular mode for geopotential height at 800 hPa and of a dipole between the Barents and Labrador Seas for the sea-ice concentration which are similar to observed patterns of variability. In addition, the simulation displays strong decadal variability in the sea-ice volume, with a significant peak at about 18 years. Positive volume anomalies are caused by (1) a decrease in ice export through Fram Strait associated with more anticyclonic winds at high latitudes, (2) modifications in the freezing/melting rates in the Arctic due to lower air temperature and higher surface albedo, and (3) a weaker heat flux at the ice base in the Barents and Kara seas caused by a lower inflow of warm Atlantic water. Opposite anomalies occur during the volume-decrease phase of the oscillation.
Sea ice and the stability of north and northeast Greenland floating glaciers
Authors: Reeh Niels , Thomsen Henrik Højmark , Higgins Anthony K. , Weidick Anker
Published: November 2007

The interaction between sea ice and glaciers has been studied for the floating tongue of Nioghalvfjerdsfjorden glacier, northeast Greenland (79°30’N, 22° W). Information from glacial geological studies, expedition reports, aerial photographs and satellite imagery is used to document the glacier front position and fast-ice conditions on millennial to decadal time-scales. The studies indicate that the stability of the floating glacier margin is dependent on the presence of a protecting fast-ice cover in front of the glacier. In periods with a permanent fast-ice cover, no calving occurs, but after fast-ice break-up the glacier responds with a large calving activity, whereby several years of accumulated glacier-ice flux suddenly breaks away. Climate-induced changes of sea-ice conditions in the Arctic Ocean with seasonal break-up of the near-shore fast ice could lead to disintegration of the floating glaciers. The present dominant mass loss by bottom melting would then to a large extent be taken over by grounding-line calving of icebergs. The local influx of fresh water from the north Greenland glaciers to the sea would be reduced and the local iceberg production would increase.
Seasonal contrasts in snow-cover characteristics on Ross Sea ice floes
Authors: Morris Kim , Jeffries Martin O.
Published: November 2007

Snow-cover characteristics on ice floes in the Ross Sea, Antarctica, were examined during cruises in autumn 1998 and summer 1999. The autumn snow cover was shallower, colder and had higher and more variable salinity, and smaller single and composite grain-sizes than the summer snow cover. The autumn snow cover was dominated by rounded particles in chains of grains and clusters, while the summer snow cover was composed primarily of melt clusters. There was extensive flooding of the summer snow cover at the snow/ice interface. The summer snow cover was nearly isothermal and close to the melting point. It exhibited obvious signs of melting and refreezing in the form of ice lenses, pipes and superimposed ice, although no melt ponds were evident. Many of the ice lenses were located directly above the saline standing water found on most of the summer ice floes.
Arctic sea-ice oscillation: regional and seasonal perspectives
Authors: Wang Jia , Ikeda Moto
Published: November 2007

Variability of the sea-ice cover (extent) in the Northern Hemisphere (Arctic and subpolar regions) associated with the Arctic Oscillation (AO) is investigated using historical data from 1901 to 1997. A principal-component analysis (empirical orthogonal functions (EOFs)) was applied to sea-ice area (SIA) anomalies for the period 1953−95. The leading EOF mode for the SI A anomaly shows an in-phase fluctuation in response to the AO and is called the Arctic sea-ice oscillation (ASIO). Arctic sea ice experiences seasonal variations that differ in timing and magnitude. Four types of seasonal variation are identified in the Arctic sea ice, and are superimposed on long-term interannual to decadal variability. Consistent with the total Arctic SIA anomaly eight regional SIA anomalies have shown significant in-phase decrease (downward trend) since 1970, possibly part of a very long-term (century) cycle. Thus, it is recommended that SIA anomalies in the sensitive seasons be used to better capture interannual, interdecadal and longer (century) variability. Major decadal and interdecadal time-scales of SIA anomalies are found at 12−14 and 17−20 years. In the Sea of Okhotsk, a century time-scale is evident. The reduction rate (negative trend) of the total Arctic sea-ice cover in the last three decades is −4.5% per decade, with the summer rate being the highest (-10.2% per decade). The contribution to this total reduction varies from region to region, with sea-ice cover in the Greenland and Norwegian Seas experiencing the highest reduction rate of −20.2 % per decade.
Seasonal transformation of water masses, circulation and brine formation observed in Storfjorden, Svalbard
Authors: Haarpaintner Jörg , O’Dwyer Jane , Gascard Jean-Claude , Haugan Peter M. , Schauer Ursula , Østerhus Svein
Published: November 2007

The transformation of ocean water masses at high latitudes is closely related to the freezing and melting processes during the year. Downward salt fluxes from brine rejection during freezing increase the salinity and density of the water column underneath. Fresh-water input from river run-off and melting of sea ice reduces the density, mainly of the surface layer. Hydrographic profiles collected in Storfjorden, Svalbard, in spring and summer, show the strong seasonal and interannual variability of the water masses. Using, in addition, data from moorings, a ship-borne acoustic Doppler current profiler and drifting buoys, and results from models of ice drift, polynya evolution, ice formation and convection processes during freezing, we document the seasonal water-mass transformation and try to explain its interannual variability. The advection of ice and water through the two northern sounds and over the sill in the south of the fjord is examined. The release of brine-enriched bottom water over the sill must be balanced by advection from the Barents Sea in the upper layers. The interannual variability of the brine-enriched bottom layer is very high, and higher salinities are observed in a milder winter. The density anomaly resulting from freezing might depend more on the ice cover and geographical position of the polynya than on the total atmospheric forcing during winter.
The Antarctic Circumpolar Wave in a global, high-resolution, coupled ice-ocean model
Authors: Zhang Yuxia , Semtner Albert J.
Published: November 2007

The Antarctic Circumpolar Wave (ACW) is identified by White and Peterson (1996) as anomalies in sea-level pressure, meridional wind stress (MWS), sea-surface temperature (SST) and sea-ice extent (SIE) propagating eastward over the Southern Ocean. In this study, the ACW is examined using a global coupled ice-ocean model with an average horizontal grid size of 1/4°. The model is forced with 1979−93 daily average atmospheric data from the European Centre for Medium-range Weather Forecasts (ECMWF) re-analysis (ERA). The sea-ice model includes both dynamics and thermodynamics, and the ocean model is a primitive-equation, free-surface, z-coordinate model. Both standing and propagating oscillations are present in ERA surface net heat-flux (NHF) and MWS anomalies. The ocean and ice respond to such atmospheric forcing with similar standing and propagating oscillations. For the propagating mode, SIE, SST and sea-surface salinity anomalies propagate eastward with a period of about 4−5 years and take about 8−9 years to encircle the Antarctic continent. Thus, the simulated ACW is a wavenumber-2 phenomenon which agrees with the ACW identified by White and Peterson (1996). The correctly simulated strength of the Antarctic Circumpolar Current, which governs the phase speed of oceanic anomalies, in our high-resolution model is essential for obtaining the observed wavenumber-2 ACW mode in the ocean. The ACW signature is also present in ocean temperature and salinity anomalies down to about 1000 m depth with similar eastward-propagating speed. The anomalies in the interior ocean are more coherent and intense over the Pacific and Atlantic sectors than over the Indian sector. Northward (southward) MWS anomalies, northward (southward) SIE anomalies, cold (warm) SST anomalies and saltier (fresher) than normal salinity anomalies are in phase, while less (more) than normal NHF is 90° out of phase with them, indicating the ACW in sea ice and ocean is a response to that in the atmosphere.
Modeling the formation and evolution of oriented fractures in sea ice
Authors: Hibler W.D.
Published: November 2007

While sea-ice drift tends to be a relatively coherent field moving largely in response to spatially averaged wind, strain-rate observations show relative ice motion to be much more of a chaotic field with deformation occurring along narrow intersection strikes between effectively rigid plates. This relative motion also has considerable energy at high frequencies which are not coherent with the wind forcing. Plastic rheologies for a strongly interactive medium are inherently able to predict such fractures. To investigate this phenomenon, this paper presents a basic anisotropic theory for multiple-oriented weaknesses. When coupled with short-term evolution equations for ice strength, this model is found to yield strong zones of failure over the Arctic basin, with relatively rigid motion in between. These failure lines typically occur along geostrophic wind contours. These zones are also found to be created with appropriate "Coulombic" isotropic rheologies. When combined with a previously developed inertial embedding procedure for the oceanic boundary layer, the resulting deformation field tends to fluctuate with considerable power near the inertial period even though the wind field is constant. The mechanisms necessary for obtaining such oriented features are discussed.
Seasonal changes in Arctic sea-ice morphology
Authors: Perovich Donald K. , Richter-Menge Jacqueline A. , Tucker Walter B.
Published: November 2007

The morphology of the Arctic sea-ice cover undergoes large changes over an annual cycle. These changes have a significant impact on the heat budget of the ice cover, primarily by affecting the distribution of the solar radiation absorbed in the ice-ocean system. In spring, the ice is snow-covered and ridges are the prominent features. The pack consists of large angular floes, with a small amount of open water contained primarily in linear leads. By the end of summer the ice cover has undergone a major transformation. The snow cover is gone, many of the ridges have been reduced to hummocks and the ice surface is mottled with melt ponds. One surface characteristic that changes little during the summer is the appearance of the bare ice, which remains white despite significant melting. The large floes have broken into a mosaic of smaller, rounded floes surrounded by a lace of open water. Interestingly, this break-up occurs during summer when the dynamic forcing and the internal ice stress are small During the Surface Heat Budget of the Arctic Ocean (SHEBA) field experiment we had an opportunity to observe the break-up process both on a small scale from the ice surface, and on a larger scale via aerial photographs. Floe break-up resulted in large part from thermal deterioration of the ice. The large floes of spring are riddled with cracks and leads that formed and froze during fall, winter and spring. These features melt open during summer, weakening the ice so that modest dynamic forcing can break apart the large floes into many fragments. Associated with this break-up is an increase in the number of floes, a decrease in the size of floes, an increase in floe perimeter and an increase in the area of open water.
Investigating the anomalous sea-ice conditions in the Canadian High Arctic (Queen Elizabeth Islands) during summer 1998
Authors: Jeffers Sharon , Agnew Tom A. , Alt Bea Taylor , de Abreu Roger , McCourt Steve
Published: November 2007

The Queen Elizabeth Islands of the Canadian Arctic Archipelago exhibit one of the most complex sea-ice regimes in the Northern Hemisphere. Time series of minimum monthly passive-microwave sea-ice area (1979−98), minimum sea-ice extent, melting degree-days (1961−98) and minimum sea ice from the new Canadian Ice Service digital database (1968−98) are examined. The extreme nature of the amount of sea-ice melt in the summers of 1998 and 1962 is evident in these time series. The 38 year record of minimum ice, to date, shows no significant trend. Details of the sea-ice behavior during summer 1998 were then examined within 13 individual sea-ice regimes. The multi-year fast-ice plugs in both Sverdrup Channel and Nansen Sound broke up and became truly mobile in 1998. Discussion focuses on the areas surrounding the multi-year plugs, relating sea-ice conditions to weather. Results emphasize the importance of the timing of synoptic events in combination with strong thermal preconditioning in determining the sea-ice conditions in this area during summer 1998.
How a region of cracked sea ice affects ice-coupled wave propagation
Authors: Squire Vernon A. , Dixon Tony W.
Published: November 2007

By deriving the appropriate Green’s function, a model is developed that allows the interaction of normally-incident, ice-coupled waves with any number of cracks to be studied analytically. For a single crack a simple formula for the reflection and transmission coefficients, R and T, emerges that yields identical results to the computationally intensive work of Barrett and Squire (1996) but is much easier to apply. A crack is found to behave as a steep low-pass filter, allowing long waves through while inhibiting shorter waves, although there is also some fine structure to the response curve. The introduction of more cracks is straightforward. While in that case a formula for R and T is also possible in principle, it is easier to express the result as the solution of a simple matrix equation of order 2N, where N is the number of cracks. It is found that perfect transmission (|R| = 0) occurs at a set of discrete periods, hereinafter called a comb, for N > 1 and that the comb becomes finer as period decreases. For both periodically distributed cracks and ones that are randomly spaced, the gross shape of the response curve remains similar. The results suggest that it is improbable that waves travelling through the Arctic basin can be used as a remote-sensing agent to determine mean ice thickness. The Green’s function technique employed in this paper furnishes solutions to other problems of interest.
Arctic sea-ice morphological characteristics in summer 1996
Authors: Wadhams P. , Davis N.R.
Published: November 2007

A sea-ice thickness profile obtained in September 1996 from the Greenland Sea and the Eurasian Basin, extending as far as the North Pole, has shown an unusually open ice cover with low mean drafts, large amounts of open water and little deep pressure ridging. Comparisons with data obtained from the same region in October 1976 show that mean ice draft has declined by 43% and that the decline can largely be ascribed to a loss of the thickest ice.
Recent fluctuation of snow cover on mountainous areas in Japan
Authors: Shimizu Masujiro , Abe Osamu
Published: November 2007

To monitor the snow-cover distribution in relation to meteorological conditions on high mountainous areas in Japan, NIED constructed a snow-observation network which it has operated for approximately 10 years. The network consists of seven pairs of stations, each comprising a mountainous site and a low-lying flatland site, from Hokkaido district in the north to San-in in the southwest. Data obtained include snow depth, snow weight, air temperature and global solar radiation. This study presents recent fluctuation of snow cover on mountainous areas for several recent winters in Japan. Most winters were warmer than average, but winter 1995/96 was normal, and maximum snow depths were recorded in high-elevation areas. Seventy-seven avalanche accidents occurred in winter 1995/96. The relationship between meteorological and snow conditions in mountainous areas and flatland areas was analyzed.
Aeolian susceptibility maps: methodology and applications
Authors: Font D. , Furdada G. , Vilaplana J. M.
Published: November 2007

In this paper we present a new methodology for the study of drifting-snow phenomena at an Alpine site This method emphasizes the accurate determination of local wind directions and the distribution of snow denudation and accumulation zones through the observation, measurement and mapping of aeolian “snowforms” and drifts. The aim is to find out whether the local wind directions on extensive surfaces in mountain areas (of about a few km2) are independent of the main wind intensity and direction, and whether valid erosion and accumulation zones are spatially stable and can be established for dominant winds. The methodology proposed to study the spatial distribution of snow by the wind has been used in the redesign of an old ski slope at La Molina, a resort in the eastern Spanish Pyrenees. Results have been very useful for deciding the design of the slopes. Therefore, we believe that the prediction of spatial snow distribution by wind in mountain areas, without the need for expensive instruments, has a valuable part to play in the management of mountain and ski resorts and roads, and in avalanche forecasting.
A nearest-neighbor model for regional avalanche forecasting
Authors: Brabec Bernhard , Meister Roland
Published: November 2007

The application of nearest-neighbor algorithms to the task of regional avalanche forecasting in Switzerland is presented in this paper. The database used for the development of the model consists of snow and weather data from 60 manual weather stations and conventionally estimated avalanche-hazard levels. All these data are collected by the Swiss Federal Institute for Snow and Avalanche Research on a daily basis during winter. Data between 1987 and 1996 (10 winters) are used for our study. For the manual weather stations a nearest-neighbor model has been developed: NXD-VG calculates the 10 nearest-neighbor days by using a Euclidean weighted distance metric. A regional avalanche-hazard map is calculated by interpolating the results of NXD-VG between the stations. The avalanche forecasters can access the results of the model calculation directly because they are integrated into a program for bulletin construction. The model was validated using three complementary approaches. First, the database is cross-validated for all the winters available to estimate an unbiased prediction error of the models for two selected stations. Second, selected situations of the database are recalculated. Third, the model output is compared daily to the official forecasts published during winter 1999/2000.
The influence of the forest on night-time snow surface temperature
Authors: Höller Peter
Published: November 2007

Snow surface temperature (Ts) plays an important role in the formation of surface hoar or near-surface faceted crystals The goal of this study was to obtain detailed information on Ts in different forest stands nelr the timberline. The investigations were conducted during clear nights and showed that the snow surface temperature is influenced very strongly by the forest canopy. While the air temperature was very similar on the different experimental sites, Ts was higher in the forest than in the open field; on the south-facing slope the difference between the forest and the open field was 3–4.5°C, and on the north-facing slope approximately 3–7°C. Taking into account that εair is 0.7 and εtree is 0.94, the incoming radiation (I ↓) for the different experimental sites was calculated by the equation of Brunt (the canopy density was estimated using photographs taken with an 8 mm fish-eye). To calculate Ts, air temperature and averaged values of the net radiation (because the net radiation (I) has only a small range of variation during clear nights) were used. The results show that the calculated values were higher than the measured values (by approximately 2°C). However, a better correlation was found by using lower values of the emissivity (εair0.67 and εtree0.91).
Numerical modelling of sublimation on snow and comparison with field measurements
Authors: Hachikubo Akihiro
Published: November 2007

Various field observations of surface-hoar formation were carried out with measurements of vapor sublimation rate and meteorological conditions from 1994 to 1996 in a mountainous area of northern Hokkaido,Japan. The sublimation rate seemed to increase with wind speed when the relative humidity was high, whereas the sublimation rate changed from positive to negative as the wind speed increased when the relative humidity was low (60–70%RH). These results were explained with a numerical model (called the Simple model) for the heat balance at the snow surface; the results of the model showed that there is a specific wind speed which maximizes the sublimation rate when the relative humidity is low, and indicated that the sublimation rate becomes large under humid conditions (>90%RH) and if wind speeds are increasing in the range 0.5–3.5 m s–1. The heat balance at the snow surface reproduced by another snow model, Crocus, agreed well with the observation results. However, Crocus overestimated the sublimation rate by about 1.3 × 10–6 kg m–2 s–1 when surface hoar formed, and underestimated snow surface temperature by several degrees.
Application of statistical simulation for avalanche-risk evaluation
Authors: Chernouss P. A. , Fedorenko Yu.
Published: November 2007

Avalanche risk is considered as the probability of an avalanche event that could cause certain losses. A unified approach for avalanche-risk evaluation by statistical simulation is suggested. A chain of models, describing snow deposition, transformation, stability, avalanche dynamics and interaction with an obstacle, is used within such an approach. Each of the models evaluates a given situation in a deterministic manner, yielding a unique result value. Output data of each model can be the input for the next model in the chain. Uncertainty of input data is described in a probabilistic manner. Using the chain of the deterministic models and simulating the input data according to appropriate probability distributions with the Monte Carlo method, risk evaluations are obtained as the ratio of a number of certain types of outcome to the total number of tests. All kinds of information on weather, snow and avalanches can be used within this scheme. The simulation process can be started at any stage of the modelling. In this study it was started from the snow-cover stability simulation. Application of the statistical simulation in such a way gives an opportunity to reflect uncertainty of the initial data in the results obtained. The suggested scheme was used for producing a computer-assisted workplace for avalanche forecasting, “LAVINA”, which has been exploited at the Centre of Avalanche Safety of “Apatit”, Kirovsk, Russia, since the early 1990s. It is an integrated system that permits spatial and temporal estimations of snow stability and avalanche dynamics to be made. Assessments of the validity of the risk estimation made with LAVINA are presented.
Snow-cover depth, distribution and duration data from northeast Greenland obtained by continuous automatic digital photography
Authors: Christiansen Hanne H.
Published: November 2007

Detailed data on autumn, winter and spring snow-cover conditions from remote areas are often difficult or very expensive to obtain. Therefore, an inexpensive method of digital photography was tested in high-Arctic Greenland. Automatic digital photography has provided daily data on snow distribution and snow depth for > 1 year from the Zackenberg area (74°30’ N) in northeast Greenland. A standard digital hand-held camera (Kodak Digital Science DC50) was equipped to become automatic; it is supplied with an “automatic finger” and an external power supply and built into a protecting box with additional solar panels on top, in order to secure continuous operation throughout the year. The daily photograph covers a 100 m transect through a seasonal snowpatch, and thus on an annual basis also yields information on snow-cover duration in the different vegetation zones of the snowpatch. The camera was installed in mid-August 1998. Photographs from the period mid-August 1998 to early September 1999 were collected and analyzed. All photos are taken at astronomic noon in order to use the daylight as long as possible into the winter season before the 24 h winter darkness begins in mid-November.

The digital photographs yielded the following information for the year 1998/99: The first winter snowfall occurred on 18 October; small-scale snow redistribution by snow-drifting started when the winter snow cover was about 5 cm thick. The continuous winter snow cover lasted for a minimum of 170 days in the most snow-deprived areas downwind of the snowpatch, whereas snow in the centre of the snowpatch stayed on the ground for > 325 days, turning the snowpatch perennial in summer 1999.

Meteorological data obtained close to the photographed snowpatch site, in combination with the snow-cover depth and distribution data derived from the photographs, show that at wind speeds (at 2 m height) of up to around 6 m s−1, only small-scale snow-drifting took place, when the snow cover was thin. Intense snowdrifting, however, took place at wind speeds of 8–13 m s−1, particularly when enough snow was available in the upwind area. The automatic continuous photography technique demonstrated here could be particularly useful in remote areas at high risk of avalanches during winter. It is an alternative to traditional snow monitoring achieved mainly by sonic sensors, snow pillows or manual measurements of snow depths Likewise, it can provide better areal information than most standard methods, which give only point measurements.
Estimate of uncertainties in avalanche hazard mapping
Authors: Barbolini M. , Savi F.
Published: November 2007

The present work addresses the urgent demand for methods of quantifying the uncertainties inherent in the current procedures for avalanche hazard assessment. A Monte Carlo approach to hazard mapping is proposed for this purpose. This statistical sampling-analysis method allows us to evaluate the probability distributions of the relevant variables for avalanche hazard assessment –– essentially runout distance and impact pressure-once the release variables and the model parameters are expressed in terms of suitable probability distributions. In this way it is possible to explicitly account for uncertainties both in the input-data definition of the dynamic models and in the mapping results. The overall methodology is presented in detail and applied to a real-world avalanche mapping problem. The one-dimensional version of the VARA models is used for avalanche dynamics simulations.
On the characterization of seismic signals generated by snow avalanches for monitoring purposes
Authors: Suriñach E. , Furdada G. , Sabot F. , Biesca B. , Vilaplana J. M.
Published: November 2007

Seismic signals from artificially released avalanches were studied in an attempt to characterize them for avalanche-monitoring purposes. The seismic signals generated by different sizes and types of avalanches were recorded and analyzed in the time and frequency domains. Synchronized recordings of the corresponding seismic signals and the video images of the evolution of the avalanches were obtained together with a detailed cartography. Characteristic signatures in the frequency and time domains were found to depend on the characteristics of the avalanche path and measuring location, but to be mostly independent of avalanche size. The source of the different parts of the recorded seismic signals was determined. A relationship was observed between the avalanche size and the amplitude of the signals. Given the presence of local site effects, a prior seismic characterization of the avalanche path in relation to the recording sites is necessary for monitoring purposes. Moreover” it was found that sliding slabs in the early phase of acceleration produce little seismic energy, resulting in a time lapse between the observable start of the avalanche and the arrival of the detectable seismic waves at the receiving station.
Stable-isotope records from Dronning Maud Land, Antarctica
Authors: Graf Wolfgang , Oerter Hans , Reinwarth Oskar , Stichler Willibald , Wilhelms Frank , Miller Heinz , Mulvaney Robert
Published: November 2007

The European Project for Ice Coring in Antarctica (EPICA) includes a comprehensive pre-site survey on the inland ice plateau of Dronning Maud Land. This paper focuses on the investigation of the 18O content of shallow firn and ice cores. These cores were dated by profiles derived from dielectric-profiling and continuous flow analysis measurements. the individual records were stacked in order to obtain composite chronologies of 18O contents and accumulation rates with enhanced signal-to-noise variance ratios. These chronologies document variations in the last 200 and 1000 years. The18O contents and accumulation rates decreased in the 19th century and increased during the 20th century. Using the empirical relationships between stable isotopes, accumulation rates and the 10 m firn temperature, the variation of both parameters can be explained by the same temperature history. But other causes for these variations, such as the build-up of the snow cover, cannot be excluded. A marked feature of the 1000 year chronology occurs during the period AD 1180–1530 when the 18Ocontents remains below the long-term mean. Cross-correlation analyses between five cores from the Weddell Sea region and Dronning Maud Land show that 18O records can in some periods be positively correlated and in others negatively correlated, indicating a complex climatic history in time and space.
The irregular pattern of isotopic and ionic signals in the typical monsoon temperate-glacier area, Yulong mountain, China
Authors: Yuanqing He , Theakstone Wilfred H. , Tandong Yao , Tuo Chen , Zhang David D.
Published: November 2007

Sampling was carried out at Baishui glacier No. 1, the largest glacier on Yulong mountain, China, during the summers of 1999 and 2000, to investigate the spatial variations of oxygen isotopes in the atmosphere–glacier–river system. the results confirm that there is an apparent inverse relation between the oxygen isotopic composition of precipitation and air-temperature/precipitation amount in this region, with lower δ18O values when the amount of precipitation and air temperature in summer is higher, due to the influence of intense monsoon climate on the study area. There are marked differences in the δ18O values of winter-accumulated snow, glacial meltwater, summer precipitation and the glacier-fed river water. Spatial and temporal variations of isotopic composition are controlled by varied weather conditions at different altitudes. Isotopic depletion or fractionation occurs during phase changes, snow-to-ice and ice-to-meltwater transformations and runoff processes.Variations of stable isotopes in glacier runoff can indicate variations of sources of supply. Ionic changes occur as a result of meltwater contact with glacier bed materials.
Seasonal changes in snow surface roughness characteristics at Summit, Greenland: implications for snow and firn ventilation
Authors: Albert Mary R. , Hawley Robert L.
Published: November 2007

Air–snow transfer processes impact both ice-core interpretation and exchange affecting atmospheric chemistry. An understanding of seasonal differences in the character of the surface snow will facilitate evaluation of possible preferential seasonal exchange of reactive chemical species. Both diffusive processes and advective (ventilation) processes can serve to alter the physical, chemical and isotopic character of snow and firn. In this paper, we examine measurements of surface roughness over the course of a year at Summit, Greenland, and the implications for snow and firn ventilation. At Summit, during the winter-over experiment, summer and fall sastrugi amplitudes were approximately 5 cm and had smoothly curved profiles. the average amplitudes experienced mild increases in January, but by the end of February through March the amplitude increased to approximately 20 cm, and the profiles exhibited more abrupt geometries. Calculations are performed to show the potential impact of the changing roughness on interstitial ventilation rates in the snow, assuming that the permeability profile does not change in time. Under high winds, ventilation velocities in the near-surface snow can be up to 3 cms–1 in the winter, compared to 1 cms–1 in the summer. the frequency of 12 ms–1 winds in the summer, however, is less than in the winter. Under low-wind conditions, the summer roughness causes ventilation rates that are comparable to diffusion rates. However, in winter even 5 ms–1 wind conditions can cause the interstitial airflow due to ventilation to exceed the diffusion rates.
Mass balance and accumulation rate across Siple Dome, West Antarctica
Authors: Hamilton Gordon S.
Published: November 2007

Snow-accumulation rates and rates of ice-thickness change (mass balance) are studied at several sites on Siple Dome,West Antarctica. Accumulation rates are derived from analyses of gross beta radioactivity in shallow firn cores located along a 60 km transect spanning both flanks and the crest of the dome. There is a north–south gradient in snow-accumulation rate across the dome that is consistent with earlier radar mapping of internal stratigraphy. Orographic processes probably control this distribution. Mass balance is inferred from the difference between global positioning system (GPS)-derived vertical velocities and snow-accumulation rates for sites close to the firn-core locations. Results indicate that there is virtually no net thickness change at four of the five sites. the exception is at the northernmost site where a small amount of thinning is detected, that appears to be inconsistent with other studies. A possible cause of this anomalous thinning is recent retreat of the grounding line of Ice Stream D.
Potential for climate variability reconstruction from Andean glaciochemical records
Authors: Ginot Patrick , Schwikowski Margit , Schotterer Ulrich , Stichler Willibald , Gäggeler Heinz W. , Francou Bernard , Gallaire Robert , Pouyaud Bernard
Published: November 2007

In order to select glacier sites suitable for the reconstruction of short-term climate variability, chemical constituents of shallow cores from Chimborazo, Ecuador, Illimani, Bolivia, and Cerro Tapado, Chile, were analyzed to determine if they represent the climate and atmospheric conditions of the particular area. Pronounced variations were observed in the Chimborazo core, which were attributed to the seasonal occurrence of wet and dry periods. Using chemical tracers, the Pacific and the Amazon basin were identified as the two principal sources of atmospheric moisture. In the Illimani record, seasonality is obvious but is less regular than in the Chimborazo core. Chemical tracers point to the Amazon basin as a moisture source. Post-depositional sublimation masks a possible atmospheric signal in the glaciochemical records from Cerro Tapado. the arid conditions and a prevailing high condensation level also cause little variability in the stable-isotope content. Irregular, more negative δ18O values occur during the humid phases of El Niño events. Despite alterations due to sublimation, chemical tracers indicate the Pacific Ocean as a main moisture source for precipitation preserved on Glaciar CerroTapado.
Folding in strongly anisotropic layers near ice-sheet centers
Authors: Thorsteinsson Throstur , Waddington Edwin D.
Published: November 2007

Conditions for passive folding near ice-sheet centers are derived treating the ice as an anisotropic viscous medium. Vertical uniaxial compression at a dome, or pure shear stress near a ridge divide, both tend to stretch and flatten folds, while horizontal simple shear deformation tends to overturn folds. Overturned folding is likely in a given vicinity in a steady-state flow field, if the initial slope of a layer disturbance exceeds the ratio of compressive/extensive deformation to shear deformation. Analytical equations for particle tracks in steady state allow us to model the evolution of layers with initial slope disturbances (``wrinkles’’). the effects of anisotropy are explored using an analytical solution for the strain rate as a function of a vertically symmetric c-axis orientation distribution, called a cone fabric. Stronger anisotropy (small cone angle) makes the material softer in horizontal shear, and facilitates folding, i.e. elements with smaller slopes can be overturned for the same stress. the relation between anisotropy and folding is complicated by the fact that, for a range of cone angles, the material is also softer in compression, which opposes folding. Simulating a layer with spatially variable tilt of cone symmetry axes, and accounting for fabric development, demonstrates that variations in the fabric cause localized flow variations that could create the initial perturbations.
Changes in chemical and isotopic properties near infiltrated cracks in an ice core from Ürümqi glacier No. 1, Tien Shan, China
Authors: Xinqing Lee , Dahe Qin , Shugui Hou , Jiawen Ren , Keqin Duan , Hui Zhou
Published: November 2007

Ice cores recovered for paleoclimatic and/or paleoenvironmental reconstructions in the Tien Shan and Qinghai–Tibetan Plateau often encounter cracks. Although we expect that cracks opened to surface meltwater will inevitably change ice-core records, we do not know how and to what extent records are influenced. An ice core retrieved from glacier No. 1 at Ürümqi river head, Tien Shan, China, exhibits a crack nearly 2.5 m long that has admitted meltwater, forming secondary ice within the fracture. A small inclusion of the infiltrated ice in sampling is shown to reduce δ18O by an extent of Holocene vs Last Glacial Maximum while enhancing significantly the pH, conductivity and the following ionic species: CH3COO–, and CO(COO)22–. of the parameters increased, and HCOO– are the most affected, being enhanced nearly six-fold in the fractured section compared to the non-fractured sections, followed by CO(COO)22– and electrical conductivity measurement (ECM). Despite the alteration, primary fluctuations of some parameters are still recognizable. This suggests that if the infiltrated ice can be avoided in the sampling operation, ice cores with cracks may still provide authentic records. This shows the need to pay close attention to physical characteristics of ice cores in order to identify such secondary ice.