A moving point approach to model shallow ice sheets: a study case with radially-symmetrical ice sheets

Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly and the waiting time behaviours to be modelled efficiently. A finite difference moving point scheme is derived and applied in a simplified context (continental radially-symmetrical shallow ice approximation). The scheme, which is inexpensive, is validated by comparing the results with moving-margin exact solutions and steady states. In both cases the scheme is able to track the position of the ice sheet margin with high precision.

A moving-point approach to model shallow ice sheets: a study case with radially symmetrical ice sheets

Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving-point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly. Our approach is also well suited to capture waiting-time behaviour efficiently. A finite-difference moving-point scheme is derived and applied in a simplified context (continental radially symmetrical shallow ice approximation). The scheme, which is inexpensive, is verified by comparing the results with steady states obtained from an analytic solution and with exact moving-margin transient solutions. In both cases the scheme is able to track the position of the ice sheet margin with high accuracy.

Impact of Laptev Sea flaw polynyas on the atmospheric boundary layer and ice production using idealized mesoscale simulations

The interaction between polynyas and the atmospheric boundary layer is examined in the Laptev Sea using the regional, non-hydrostatic Consortium for Small-scale Modelling (COSMO) atmosphere model. A thermodynamic sea-ice model is used to consider the response of sea-ice surface temperature to idealized atmospheric forcing. The idealized regimes represent atmospheric conditions that are typical for the Laptev Sea region. Cold wintertime conditions are investigated with sea-ice–ocean temperature differences of up to 40 K. The Laptev Sea flaw polynyas strongly modify the atmospheric boundary layer. Convectively mixed layers reach heights of up to 1200 m above the polynyas with temperature anomalies of more than 5 K. Horizontal transport of heat expands to areas more than 500 km downstream of the polynyas. Strong wind regimes lead to a more shallow mixed layer with strong near-surface modifications, while weaker wind regimes show a deeper, well-mixed convective boundary layer. Shallow mesoscale circulations occur in the vicinity of ice-free and thin-ice covered polynyas. They are forced by large turbulent and radiative heat fluxes from the surface of up to 789 W m−2, strong low-level thermally induced convergence and cold air flow from the orographic structure of the Taimyr Peninsula in the western Laptev Sea region. Based on the surface energy balance we derive potential sea-ice production rates between 8 and 25 cm d−1. These production rates are mainly determined by whether the polynyas are ice-free or covered by thin ice and by the wind strength.

Observations of supercooling and frazil ice formation in the Laptev Sea coastal polynya

This paper examines a hydrographic response to the wind‐driven coastal polynya activity over the southeastern Laptev Sea shelf for April–May 2008, using a combination of Environmental Satellite (Envisat) advanced synthetic aperture radar (ASAR) and TerraSAR‐X satellite imagery, aerial photography, meteorological data, and SBE‐37 salinity‐temperature‐depth and acoustic Doppler current profiler land‐fast ice edgemoored instruments. When ASAR observed the strongest end‐of‐April polynya event with frazil ice formation, the moored instruments showed maximal acoustical scattering within the surface mixed layer, and the seawater temperatures were either at or 0.02°C below freezing. We also find evidence of the persistent horizontal temperature and salinity gradients across the fast ice edge to have the signature of geostrophic flow adjustment as predicted by polynya models.

Experimental Study on Friction between Saline Ice and Steel

Understanding the interaction of sea ice with offshore structures is of primary importance for the development of technology in cold climate regions. The rheological properties of sea ice (strength, creep, viscosity) as well as the roughness of the contact surface are the main factors influencing the type of interaction with a structure. A device was developed and designed and small scale laboratory experiments were carried out to study sea ice frictional interaction with steel material by means of a uniaxial compression rig. Sea-ice was artificially grown between a stainless steel piston (of circular cross section) and a hollow cylinder of the same material, coaxial to the former and of the same surface roughness. Three different values for the roughness were tested: 1.2, 10 and 30 μm Ry (maximum asperities height), chosen as representative values for typical surface conditions, from smooth to normally corroded steel. Creep tests (0.2, 0.3, 0.4 and 0.6 kN) were conducted at T = -10 ºC. By pushing the piston head towards the cylinder base, three different types of relative movement were observed: 1) the piston slid through the ice, 2) the piston slid through the ice and the ice slid on the surface of the outer cylinder, 3) the ice slid only on the cylinder surface. A cyclic stick-slip motion of the piston was detected with a representative frequency of 0.1 Hz. The ratio of the mean rate of axial displacement to the frequency of the stick-slip oscillations was found to be comparable to the roughness length (Sm). The roughness is the most influential parameter affecting the amplitude of the oscillations, while the load has a relevant influence on the their frequency. Guidelines for further investigations were recommended. Marco Nanetti - seloselo@virgilio.it

A refined TALDICE-1a age scale from 55 to 112 ka before present for the Talos Dome ice core based on high-resolution methane measurements

A precise synchronization of different climate records is indispensable for a correct dynamical interpretation of paleoclimatic data. A chronology for the TALDICE ice core from the Ross Sea sector of East Antarctica has recently been presented based on methane synchronization with Greenland and the EDC ice cores and δ18Oice synchronization with EDC in the bottom part (TALDICE-1). Using new high-resolution methane data obtained with a continuous flow analysis technique, we present a refined age scale for the age interval from 55–112 thousand years (ka) before present, where TALDICE is synchronized with EDC. New and more precise tie points reduce the uncertainties of the age scale from up to 1900 yr in TALDICE-1 to below 1100 yr over most of the refined interval and shift the Talos Dome dating to significantly younger ages during the onset of Marine Isotope Stage 3. Thus, discussions of climate dynamics at sub-millennial time scales are now possible back to 110 ka, in particular during the inception of the last ice age. Calcium data of EDC and TALDICE are compared to show the impact of the refinement to the synchronization of the two ice cores not only for the gas but also for the ice age scale.

Water isotopic ratios from a continuously melted ice core sample

A new technique for on-line high resolution isotopic analysis of liquid water, tailored for ice core studies is presented. We built an interface between a Wavelength Scanned Cavity Ring Down Spectrometer (WS-CRDS) purchased from Picarro Inc. and a Continuous Flow Analysis (CFA) system. The system offers the possibility to perform simultaneuous water isotopic analysis of δ18O and δD on a continuous stream of liquid water as generated from a continuously melted ice rod. Injection of sub μl amounts of liquid water is achieved by pumping sample through a fused silica capillary and instantaneously vaporizing it with 100% efficiency in a~home made oven at a temperature of 170 °C. A calibration procedure allows for proper reporting of the data on the VSMOW–SLAP scale. We apply the necessary corrections based on the assessed performance of the system regarding instrumental drifts and dependance on the water concentration in the optical cavity. The melt rates are monitored in order to assign a depth scale to the measured isotopic profiles. Application of spectral methods yields the combined uncertainty of the system at below 0.1‰ and 0.5‰ for δ18O and δD, respectively. This performance is comparable to that achieved with mass spectrometry. Dispersion of the sample in the transfer lines limits the temporal resolution of the technique. In this work we investigate and assess these dispersion effects. By using an optimal filtering method we show how the measured profiles can be corrected for the smoothing effects resulting from the sample dispersion. Considering the significant advantages the technique offers, i.e. simultaneuous measurement of δ18O and δD, potentially in combination with chemical components that are traditionally measured on CFA systems, notable reduction on analysis time and power consumption, we consider it as an alternative to traditional isotope ratio mass spectrometry with the possibility to be deployed for field ice core studies. We present data acquired in the field during the 2010 season as part of the NEEM deep ice core drilling project in North Greenland.

A sublimation technique for high-precision measurements of d13CO2 and mixing ratios of CO2 and N2O from air trapped in ice cores

In order to provide high precision stable carbon isotope ratios (δ13CO2 or δ13C of CO2) from small bubbly, partially and fully clathrated ice core samples we developed a new method based on sublimation coupled to gas chromatography-isotope ratio mass spectrometry (GC-IRMS). In a first step the trapped air is quantitatively released from ~30 g of ice and CO2 together with N2O are separated from the bulk air components and stored in a miniature glass tube. In an off-line step, the extracted sample is introduced into a helium carrier flow using a minimised tube cracker device. Prior to measurement, N2O and organic sample contaminants are gas chromatographically separated from CO2. Pulses of a CO2/N2O mixture are admitted to the tube cracker and follow the path of the sample through the system. This allows an identical treatment and comparison of sample and standard peaks. The ability of the method to reproduce δ13C from bubble and clathrate ice is verified on different ice cores. We achieve reproducibilities for bubble ice between 0.05 ‰ and 0.07 ‰ and for clathrate ice between 0.05 ‰ and 0.09 ‰ (dependent on the ice core used). A comparison of our data with measurements on bubble ice from the same ice core but using a mechanical extraction device shows no significant systematic offset. In addition to δ13C, the CO2 and N2O mixing ratios can be volumetrically derived with a precision of 2 ppmv and 8 ppbv, respectively.

Continuous measurements of methane mixing ratios from ice cores

This work presents a new, field-deployable technique for continuous, high-resolution measurements of methane mixing ratios from ice cores. The technique is based on a continuous flow analysis system, where ice core samples cut along the long axis of an ice core are melted continuously. The past atmospheric air contained in the ice is separated from the melt water stream via a system for continuous gas extraction. The extracted gas is dehumidified and then analyzed by a Wavelength Scanned-Cavity Ring Down Spectrometer for methane mixing ratios. We assess the performance of the new measurement technique in terms of precision (±0.8 ppbv, 1σ), accuracy (±8 ppbv), temporal (ca. 100 s), and spatial resolution (ca. 5 cm). Using a firn air transport model, we compare the resolution of the measurement technique to the resolution of the atmospheric methane signal as preserved in ice cores in Greenland. We conclude that our measurement technique can resolve all climatically relevant variations as preserved in the ice down to an ice depth of at least 1980 m (66 000 yr before present) in the North Greenland Eemian Ice Drilling ice core. Furthermore, we describe the modifications, which are necessary to make a commercially available spectrometer suitable for continuous methane mixing ratio measurements from ice cores.

Sympathetic stimulation using the cold pressor test increases coronary collateral flow

BACKGROUND Little is known about the vasomotor function of human coronary collateral vessels. The purpose of this study was to examine collateral flow under a strong sympathetic stimulus (cold pressor test, CPT). METHODS In 30 patients (62 +/- 12 years) with coronary artery disease, two subsequent coronary artery occlusions were performed with random CPT during one of them. Two minutes before and during the 1 minute-occlusion, the patient's hand was immerged in ice water. For the calculation of a perfusion pressure-independent collateral flow index (CFI), the aortic (Pao), the central venous (CVP) and the coronary wedge pressure (Poccl) were measured: CFI = (Poccl - CVP)/(Pao - CVP). RESULTS CPT lead to an increase in Pao from 98 +/- 14 to 105 +/- 15 mm Hg (p = 0.002). Without and with CPT, CFI increased during occlusion from 14% +/- 10% to 16% +/- 10% (p = 0.03) and from 17% +/- 9% to 19% +/- 9% (p = 0.006), respectively, relative to normal flow. During CPT, CFI was significantly higher at the beginning as well as at the end of the occlusion compared to identical instants without CPT. CFI at the end of the control occlusion did not differ significantly from the CFI at the beginning of occlusion with CPT. CONCLUSIONS During balloon occlusion, collateral flow increased due to collateral recruitment independent of external sympathetic stimulation. Sympathetic stimulation using CPT additionally augmented collateral flow. The collateral-flow-increasing effect of CPT is comparable to the recruitment effect of the occlusion itself. This may reflect a coronary collateral vasodilation mediated by the sympathetic nervous system.

The onset of Neoglaciation 6000 years ago in western Mongolia revealed by an ice core from the Tsambagarav mountain range

Glacier highstands since the Last Glacial Maximum are well documented for many regions, but little is known about glacier fluctuations and lowstands during the Holocene. This is because the traces of minimum extents are difficult to identify and at many places are still ice covered, limiting the access to sample material. Here we report a new approach to assess minimal glacier extent, using a 72-m long surface-to-bedrock ice core drilled on Khukh Nuru Uul, a glacier in the Tsambagarav mountain range of the Mongolian Altai (4130 m asl, 48°39.338′N, 90°50.826′E). The small ice cap has low ice temperatures and flat bedrock topography at the drill site. This indicates minimal lateral glacier flow and thereby preserved climate signals. The upper two-thirds of the ice core contain 200 years of climate information with annual resolution, whereas the lower third is subject to strong thinning of the annual layers with a basal ice age of approximately 6000 years before present (BP). We interpret the basal ice age as indicative of ice-free conditions in the Tsambagarav mountain range at 4100 m asl prior to 6000 years BP. This age marks the onset of the Neoglaciation and the end of the Holocene Climate Optimum. The ice-free conditions allow for adjusting the Equilibrium Line Altitude (ELA) and derive the glacier extent in the Mongolian Altai during the Holocene Climate Optimum. Based on the ELA-shift, we conclude that most of the glaciers are not remnants of the Last Glacial Maximum but were formed during the second part of the Holocene. The ice core derived accumulation reconstruction suggests important changes in the precipitation pattern over the last 6000 years. During formation of the glacier, more humid conditions than presently prevailed followed by a long dry period from 5000 years BP until 250 years ago. Present conditions are more humid than during the past millennia. This is consistent with precipitation evolution derived from lake sediment studies in the Altai.

A technique for continuous detection of drill liquid in ice cores

When drilling ice cores deeper than ∼100 m, drill liquid is required to maintain ice-core quality and to limit borehole closure. Due to high-pressure air bubbles in the ice, the ice core can crack during drilling and core retrieval, typically at 600–1200 m depth in Greenland. Ice from this 'brittle zone' can be contaminated by drill liquid as it seeps through cracks into the core. Continuous flow analysis (CFA) systems are routinely used to analyse ice for chemical impurities, so the detection of drill liquid is important for validating accurate measurements and avoiding potential instrument damage. An optical detector was constructed to identify drill liquid in CFA tubing by ultraviolet absorption spectroscopy at a wavelength of 290 nm. The set-up was successfully field-tested in the frame of the NEEM ice-core drilling project in Greenland. A total of 27 cases of drill liquid contamination were identified during the analysis of 175 m of brittle zone ice. The analyses most strongly affected by drill liquid contamination include insoluble dust particles, electrolytic conductivity, ammonium, hydrogen peroxide and sulphate. This method may also be applied to other types of drill liquid used at other drill sites.

Correction to Variations of Ice Bed Coupling Beneath and Beyond Ice Streams: The Force Balance

A geometrical force balance that links stresses to ice bed coupling along a flow band of an ice sheet was developed in 1988 for longitudinal tension in ice streams and published 4 years later. It remains a work in progress. Now gravitational forces balanced by forces producing tensile, compressive, basal shear, and side shear stresses are all linked to ice bed coupling by the floating fraction phi of ice that produces the concave surface of ice streams. These lead inexorably to a simple formula showing how phi varies along these flow bands where surface and bed topography are known: phi = h(O)/h(I) with h(O) being ice thickness h(I) at x = 0 for x horizontal and positive upslope from grounded ice margins. This captures the basic fact in glaciology: the height of ice depends on how strongly ice couples to the bed. It shows how far a high convex ice sheet (phi = 0) has gone in collapsing into a low flat ice shelf (phi = 1). Here phi captures ice bed coupling under an ice stream and h(O) captures ice bed coupling beyond ice streams.

Variations of Ice Bed Coupling Beneath and Beyond Ice Streams: The Force Balance

A geometrical force balance that links stresses to ice bed coupling along a flow band of an ice sheet was developed in 1988 for longitudinal tension in ice streams and published 4 years later. It remains a work in progress. Now gravitational forces balanced by forces producing tensile, compressive, basal shear, and side shear stresses are all linked to ice bed coupling by the floating fraction phi of ice that produces the concave surface of ice streams. These lead inexorably to a simple formula showing how phi varies along these flow bands where surface and bed topography are known: phi = h(O)/h(I) with h(O) being ice thickness h(I) at x = 0 for x horizontal and positive upslope from grounded ice margins. This captures the basic fact in glaciology: the height of ice depends on how strongly ice couples to the bed. It shows how far a high convex ice sheet (phi = 0) has gone in collapsing into a low flat ice shelf (phi = 1). Here phi captures ice bed coupling under an ice stream and h(O) captures ice bed coupling beyond ice streams.

Fracture and Back Stress Along the Byrd Glacier Flowband on the Ross Ice Shelf

East Antarctic ice discharged by Byrd Glacier continues as a flowband to the calving front of the Ross Ice Shelf. Flow across the grounding line changes from compressive to extensive as it leaves the fjord through the Transantarctic Mountains occupied by Byrd Glacier. Magnitudes of the longitudinal compressive stress that suppress opening of transverse tensile cracks are calculated for the flowband. As compressive back stresses diminish, initial depths and subsequent growth of these cracks, and their spacing, are calculated using theories of elastic and ductile fracture mechanics. Cracks are initially about one millimeter wide, with approximately 30 in depths and 20 in spacings for a back stress of 83 kPa at a distance of 50 kin beyond the fjord, where floating ice is 600 in thick. When these crevasses penetrate the whole ice thickness, they release tabular icebergs 20 kin to 100 kin wide, spaced parallel to the calving front of the Ross Ice Shelf