Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.

Isotopic and chemical compositions of lower crust rocks and minerals from ODP Hole 176-735B

The composition of gabbroic rocks from the drill core of Hole 735B (ODP Leg 176) at the 11 Ma Atlantis II bank close to the slow spreading Southwest Indian Ridge (SWIR) has been analyzed for major and trace elements and Sr, Nd and Pb isotopic composition. The samples are thought to represent much of the mineralogical and geochemical variation in a vertical 1-km section (500-1500 m below the sea floor) of the lower ocean crust. Primitive troctolitic gabbros, olivine gabbros and gabbros that have Mg#=84-70, Ca#>61 and low Na# (Na/(Na+Al)) (8-17) are intruded by patches or veins of more evolved FeTi-oxide rich gabbroic and dioritic rocks with Mg# to 20, Ca# to 32, Na#=14-23, TiO2<7 wt.% and FeOtotal<18 wt.%. All rocks are acdcumulates, and incompatible element concentrations are low, e.g. Pb=0.1-0.7 ppm and Uand up to 2 ppm Pb and 0.2 ppm U in the evolved. The range of isotopic compositions of the unleached rocks is: 87Sr/86Sr=0.70280-0.70299, average 0.70287+/-0.00005 (1 S.D., N=30 samples) (except one felsic vein with 87Sr/86Sr=0.7045), 143Nd/144Nd=0.51304-0.51314, average 0.51310+/-0.00002 (1 S.D., N=28), 206Pb/204Pb=17.43-18.55, 207Pb/204Pb=15.40-15.61 and 208Pb/204Pb=37.19-38.28. The range of Sr and the almost constant Nd isotopic composition resemble that found in the upper 500 m of Hole 735B, while Pb ranges to more radiogenic compositions. In general, there is a decrease in isotopic variation of Sr and Pb as well as ? (238U/204Pb), U and Pb with depth, with a trend towards relatively unradiogenic compositions. This correlates with a decrease in alteration and frequency of evolved rock-types in the core. Leached samples generally have less radiogenic Pb with values trending towards 206Pb/204Pb=17.35, 207Pb/204Pb=15.35 and 208Pb/204Pb=37.0, while their 87Sr/86Sr ratios deviate less systematically from unleached rocks and reach both higher, 0.70307, and lower values, 0.70276. Separated clinopyroxene has elevated 87Sr/86Sr up to 0.7035, while plagioclase generally has close to whole rock Sr. Leaching reduced 87Sr/86Sr in clinopyroxene and in two (out of nine) cases leached separates and whole rock display isotopic equilibrium. Relatively minor hydrothermal seawater alteration is thought to have increased 87Sr/86Sr in the rocks, while a secondary high temperature percolation of a mantle-derived agent is thought to be the cause for the trend towards radiogenic Pb. This material had intermediate 87Sr/86Sr and may have originated from non-MORB off axis mantle. The main primary igneous isotopic variation of the gabbros is suggested to have been derived from the MORB-mantle and is defined mainly by leached samples from both ODP Leg 176 and Leg 118 and can be explained by two-component mixing of an end-member with composition like Central Indian Ridge basalts and an end-member with composition unlike any MORB. The latter is characterized by very unradiogenic Pb, in particular 207Pb/204Pb, and may have an origin with affinity to old depleted mantle (DM). The isotopic composition of the magmas parental to the FeTi-oxide rich rocks cannot be distinguished from the magmas parental to the primitive gabbros and an intimate relationship is indicated. The small-scale inhomogeneity indicated for the SWIR MORB-mantle at the Atlantis II Fracture Zone was probably inherited by the lower crustal rocks due to small-scale melting and monogenetic magma chambers at this slow spreading ridge.

Shear experiments from IODP Expedition 322 and 333 samples

The mechanical behavior of the plate boundary fault zone is of paramount importance in subduction zones, because it controls megathrust earthquake nucleation and propagation as well as the structural style of the forearc. In the Nankai area along the NanTroSEIZE (Kumano) drilling transect offshore SW Japan, a heterogeneous sedimentary sequence overlying the oceanic crust enters the subduction zone. In order to predict how variations in lithology, and thus mechanical properties, affect the formation and evolution of the plate boundary fault, we conducted laboratory tests measuring the shear strengths of sediments approaching the trench covering each major lithological sedimentary unit. We observe that shear strength increases nonlinearly with depth, such that the (apparent) coefficient of friction decreases. In combination with a critical taper analysis, the results imply that the plate boundary position is located on the main frontal thrust. Further landward, the plate boundary is expected to step down into progressively lower stratigraphic units, assisted by moderately elevated pore pressures. As seismogenic depths are approached, the décollement may further step down to lower volcaniclastic or pelagic strata but this requires specific overpressure conditions. High-taper angle and elevated strengths in the toe region may be local features restricted to the Kumano transect.

Core characteristics and events detected in ice cores from Svalbard

Ice coring and snow cover observations have been carried out at 3 sites in Nordaustlandet, Svalbard since 1995. The results of stratigraphic analyses, and chemical and d18O analyses from Vestfonna and Austfonna cores are presented here. The results from these sites show that most of the chemical constituents contained in the initial snow cover still remained in the ice cores, although re-distribution of them by melt water percolation had occurred. Anthropogenic increases in trace metals, sulfate and nitrate since about 1950 are detected. This suggests that ice-core chemistry records from Nordaustlandet, Svalbard, can be useful to reconstruct past atmospheric conditions. In addition to chemical records, records, that correlate well with the temperature records in Svalbard, can be used to reconstruct past temperature changes.

Wet tundra polygon centers and distance between polygon centroids on Samoylov Island, Lena Delta, Siberia, Russia, with links to shape files

Subgrid processes occur in various ecosystems and landscapes but, because of their small scale, they are not represented or poorly parameterized in climate models. These local heterogeneities are often important or even fundamental for energy and carbon balances. This is especially true for northern peatlands and in particular for the polygonal tundra, where methane emissions are strongly influenced by spatial soil heterogeneities. We present a stochastic model for the surface topography of polygonal tundra using Poisson-Voronoi diagrams and we compare the results with available recent field studies. We analyze seasonal dynamics of water table variations and the landscape response under different scenarios of precipitation income. We upscale methane fluxes by using a simple idealized model for methane emission. Hydraulic interconnectivities and large-scale drainage may also be investigated through percolation properties and thresholds in the Voronoi graph. The model captures the main statistical characteristics of the landscape topography, such as polygon area and surface properties as well as the water balance. This approach enables us to statistically relate large-scale properties of the system to the main small-scale processes within the single polygons.

137Cs and 90Sr depositions, and 137Cs/90Sr ratios of samples from the Kraton-3 underground nuclear explosion site, Russia

Strontium-90 activity concentrations in surface soils and areal deposition densities have been studied at a site contaminated by an accidental release to atmosphere from the underground nuclear explosion 'Kraton-3' conducted near the Polar Circle (65.9°N, 112.3°E) within the territory of the former USSR in 1978. In 2001-2002, the ground surface contamination at 14 plots studied ranged from 20 to 15000 kBq/m**2, which significantly exceeds the value of 0.44 kBq/m**2 deduced for three background plots. The zone with substantial radiostrontium contamination extends, at least, 2.5 km in a north-easterly direction from the borehole. The average 137Cs/90Sr ratio in the ground contamination originated from the 'Kraton-3' fallout was estimated to be 0.55, which is significantly different from the ratio of 2.05 evaluated for background plots contaminated mostly from global fallout. Although vertical migration of 90Sr in all undisturbed soil profiles studied is more rapid than that for 137Cs, the depth of percolation of both radionuclides into the ground is mostly limited to the top 10-20 cm, which may be explained, primarily, by permafrost conditions. The horizontal migration rate of radiostrontium in the aqueous phase exceeds the radiocaesium migration rate by many times. This phenomenon seems to be a reason for the significant enrichment of the soil surface layers by radiostrontium at some sites, with variations occurring in accordance with small-scale irregularities of landscape.

Physical properties and alpha-/gamma-HCH concentrations in sea-ice and snow samples obtained during a CCGS Amundsen cruise, eastern Beaufort Sea

The alpha- and gamma-hexachlorocyclohexanes (HCHs) are being scavenged from the atmosphere by falling snow, with the average total scavenging ratios (WT) of 3.8 x 10**4 and 9.6 x 10**3, respectively. After deposition, HCH snow concentrations can decrease by 40% because of snowpack ventilation and increase by 50% because of upward migration of brine from the ice. HCH vertical distribution in sufficiently cold winter sea ice, which maintains brine volume fractions <5%, reflects the ice growth history. Initially, the entrapment of brine (and HCHs) in ice depends on the rates of ice growth and desalination. However, after approximately the first week of ice formation, ice growth rate becomes dominant. Deviations of HCH concentrations from the values predicted by the ice bulk salinity (rate of brine entrapment) can be explained by spatial variability of HCHs in surface water. HCH burden in the majority of the ice column remains locked throughout most of the season until the early spring when snow meltwater percolates into the ice, delivering HCHs to the upper ocean via desalination by flushing. Percolation can lead to an increase in alpha- and gamma-HCH in the sea ice by up to 2%-18% and 4%-32%, respectively.

pH and calcium change in the microenvironment of a benthic foraminifer (Ammonia sp.) and its size during experiments

Calcareous foraminifera are well known for their CaCO3 shells. Yet, CaCO3 precipitation acidifies the calcifying fluid. Calcification without pH regulation would therefore rapidly create a negative feedback for CaCO3 precipitation. In unicellular organisms, like foraminifera, an effective mechanism to counteract this acidification could be the externalization of H+ from the site of calcification. In this study we show that a benthic symbiont-free foraminifer Ammonia sp. actively decreases pH within its extracellular microenvironment only while precipitating calcite. During chamber formation events the strongest pH decreases occurred in the vicinity of a newly forming chamber (range of gradient about 100 µm) with a recorded minimum of 6.31 (< 10 µm from the shell) and a maximum duration of 7 h. The acidification was actively regulated by the foraminifera and correlated with shell diameters, indicating that the amount of protons removed during calcification is directly related to the volume of calcite precipitated. The here presented findings imply that H+ expulsion as a result of calcification may be a wider strategy for maintaining pH homeostasis in unicellular calcifying organisms.

Grain-size distribution and frictional behaviour of sediment cores from the Nice Slope

The upper shelf of the landslide-prone Ligurian Margin (Western Mediterranean Sea) off Nice well-known for the 1979 Airport Landslide is a natural laboratory to study preconditioning factors and trigger mechanisms for submarine landslides. For this study low-stress ring shear experiments have been carried out on a variety of sediments from >50 gravity cores to characterise the velocity-dependent frictional behaviour. Mean values of the peak coefficient of friction vary from 0.46 for clay-dominated samples (53 % clay, 46 % silt, 1 %) sand up to 0.76 for coarse-grained sediments (26 % clay, 57 % silt, 17 % sand). The majority of the sediments tested show velocity strengthening regardless of the grain size distribution. For clayey sediments the peak and residual cohesive strength increases with increasing normal stress, with values from 1.3 to 10.6 kPa and up to 25 % of all strength supported by cohesive forces in the shallowmost samples. A pseudo-static slope stability analysis reveals that the different lithologies (even clay-rich material with clay content >=50 %) tested are stable up to slope angles <26° under quasi-drained conditions.