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.

Heat flow ratios of stations from the Iberian abyssal plain (Table 1)

New heat flow observations have been made in the Iberia abyssal plain off the Galicia margin along the transeat of Ocean Drilling Program Leg 149 drill sites. in order to investigate the nature of this unusually wide and deep continet-ocean transition region. Our results indicate the presence of three separate zones. Average values of 47.5 +/- 3 mW/m in the westernmost zone III agree with predictions of standard oceanic lithospheric models for its estimated age of 126 Ma. In contrast, the heat flow within zone II is 5-15 mW/m higher than predicted. assuming that the mantle heat flow remains constant across the basin. This region of high values is coincident with the location of a major intra-crustal "S"-type reflector east of ODP Site 900. and the anomaly is consistent with the presence of 2-3 km of primarily upper continental crust above the reflector, with concentrations of radiogenic components similar to those from granodiorite samplles dredged off Galicia Bank. It is not, however, consistent with the low values of heat production measured on gabbroic sanhples from its western end at ODP Site 900. In zone I, detailed measurements across the tilted fault block south of ODP Site 901 show consistent variations which closely match predictions due to the effects of basement structure and sediment deposition. There is no evidence for variations due to vertical convective transport along the dipping basement fault block. Once corrected for these variations. measurements in zone I yield average values that agree quite well with previous measurements across Calicia Bank. indicating no systematic landward increase in heat flow with decreasing amounts of continental, extension.

(Appendix) Petrographic information, trace element and Sr isotope ratios of Pito Deep samples

Fluid flow through the axial hydrothermal system at fast spreading ridges is investigated using the Sr-isotopic composition of upper crustal samples recovered from a tectonic window at Pito Deep (NE Easter microplate). Samples from the sheeted dike complex collected away from macroscopic evidence of channelized fluid flow, such as faults and centimeter-scale hydrothermal veins, show a range of 87Sr/86Sr from 0.7025 to 0.7030 averaging 0.70276 relative to a protolith with 87Sr/86Sr of ~0.7024. There is no systematic variation in 87Sr/86Sr with depth in the sheeted dike complex. Comparison of these new data with the two other localities that similar data sets exist for (ODP Hole 504B and the Hess Deep tectonic window) reveals that the extent of Sr-isotope exchange is similar in all of these locations. Models that assume that fluid-rock reaction occurs during one-dimensional (recharge) flow lead to significant decreases in the predicted extent of isotopic modification of the rock with depth in the crust. These model results show systematic misfits when compared with the data that can only be avoided if the fluid flow is assumed to be focused in isolated channels with very slow fluid-rock exchange. In this scenario the fluid at the base of the crust is little modified in 87Sr/86Sr from seawater and thus unlike vent fluids. Additionally, this model predicts that some rocks should show no change from the fresh-rock 87Sr/86Sr, but this is not observed. Alternatively, models in which fluid-rock reaction occurs during upflow (discharge) as well as downflow, or in which fluids are recirculated within the hydrothermal system, can reproduce the observed lack of variation in 87Sr/86Sr with depth in the crust. Minimum time-integrated fluid fluxes, calculated from mass balance, are between 1.5 and 2.6 * 10**6 kg/m**2 for all areas studied to date. However, new evidence from both the rocks and a compilation of vent fluid compositions demonstrates that some Sr is leached from the crust. Because this leaching lowers the fluid 87Sr/86Sr without changing the rock 87Sr/86Sr, these mass balance models must underestimate the time-integrated fluid flux. Additionally, these values do not account for fluid flow that is channelized within the crust.

Age models of northern North Atlantic sediment cores

Benthic foraminiferal stable isotope records from four high-resolution sediment cores, forming a depth transect between 1237 m and 2303 m on the South Iceland Rise, have been used to reconstruct intermediate and deep water paleoceanographic changes in the northern North Atlantic during the last 21 ka (spanning Termination I and the Holocene). Typically, a sampling resolution of ~100 years is attained. Deglacial core chronologies are accurately tied to North Greenland Ice Core Project (NGRIP) ice core records through the correlation of tephra layers and changes in the percent abundance of Neogloboquadrina pachyderma (sinistral) with transitions in NGRIP. The evolution from the glacial mode of circulation to the present regime is punctuated by two periods with low benthic d13C and d18O values, which do not lie on glacial or Holocene water mass mixing lines. These periods correlate with the late Younger Dryas/Early Holocene (11.5-12.2 ka) and Heinrich Stadial 1 (14.7-16.8 ka) during which time freshwater input and sea-ice formation led to brine rejection both locally and as an overflow exported from the Nordic seas into the northern North Atlantic, as earlier reported by Meland et al. (2008). The export of brine with low ?13C values from the Nordic seas complicates traditional interpretations of low d13C values during the deglaciation as incursions of southern sourced water, although the spatial extent of this brine is uncertain. The records also reveal that the onset of the Younger Dryas was accompanied by an abrupt and transient (~200-300 year duration) decrease in the ventilation of the northern North Atlantic. During the Holocene, Iceland-Scotland Overflow Water only reached its modern flow strength and/or depth over the South Iceland Rise by 7-8 ka, in parallel with surface ocean reorganizations and a cessation in deglacial meltwater input to the North Atlantic.

Age models, sedimentology and geochemistry from sediment cores in Drake Passage throughflow

The Drake Passage (DP) is the major geographic constriction for the Antarctic Circumpolar Current (ACC) and exerts a strong control on the exchange of physical, chemical, and biological properties between the Atlantic, Pacific, and Indian Ocean basins. Resolving changes in the flow of circumpolar water masses through this gateway is, therefore, crucial for advancing our understanding of the Southern Ocean's role in global ocean and climate variability. Here, we reconstruct changes in DP throughflow dynamics over the past 65,000 y based on grain size and geochemical properties of sediment records from the southernmost continental margin of South America. Combined with published sediment records from the Scotia Sea, we argue for a considerable total reduction of DP transport and reveal an up to ~40% decrease in flow speed along the northernmost ACC pathway entering the DP during glacial times. Superimposed on this long-term decrease are high-amplitude, millennial-scale variations, which parallel Southern Ocean and Antarctic temperature patterns. The glacial intervals of strong weakening of the ACC entering the DP imply an enhanced export of northern ACC surface and intermediate waters into the South Pacific Gyre and reduced Pacific-Atlantic exchange through the DP ("cold water route"). We conclude that changes in DP throughflow play a critical role for the global meridional overturning circulation and interbasin exchange in the Southern Ocean, most likely regulated by variations in the westerly wind field and changes in Antarctic sea ice extent.

Age models for ODPO Leg 130 sites

The Neogene carbonate stratigraphy of five sites drilled on Ontong Java Plateau during Leg 130 reveals a number of patterns which are unexpected, and which we refer to as loss paradox, equatorial insensitivity, and climate paradox. They denote the following unresolved questions. 1 The loss of carbonate at depth (as derived from differences in accumulation rates) is much greater than suggested by the change in carbonate percentages (calculated under the assumption that carbonate dissolution is the cause of loss). This indicates an important role for redeposition processes, such as winnowing (bottom currents), sifting (resuspension and catabatic flow) and episodic sloughing or solifluction (presumably stimulated by earthquakes). 2 Accumulation rates are not markedly increased at the time a site crosses the equator. There are several possible reasons. Equatorial upwelling may be unimportant in controlling sedimentation rates this far in the western Pacific, or its output may be spread over a considerable distance from the equator. Alternatively, increased supply below the equator is compensated for by increased removal (e.g. from resuspension by bioturbation, combined with catabatic flow). It is conceivable that errors in the timescale could also produce the effect seen. 3 There is an overall tendency for agreement between the stratigraphic patterns of carbonate content and of accumulation rates, but neither pattern is readily explained by reference to changes in climate (represented by benthic delta18O) or in sea-level (as derived from sequence stratigraphy).

Debris flow magnitudes: a global catalog

The episodic occurrence of debris flow events in response to stochastic precipitation and wildfire events makes hazard prediction challenging. Previous work has shown that frequency-magnitude distributions of non-fire-related debris flows follow a power law, but less is known about the distribution of post-fire debris flows. As a first step in parameterizing hazard models, we use frequency-magnitude distributions and cumulative distribution functions to compare volumes of post-fire debris flows to non-fire-related debris flows. Due to the large number of events required to parameterize frequency-magnitude distributions, and the relatively small number of post-fire event magnitudes recorded in the literature, we collected data on 73 recent post-fire events in the field. The resulting catalog of 988 debris flow events is presented as an appendix to this article. We found that the empirical cumulative distribution function of post-fire debris flow volumes is composed of smaller events than that of non-fire-related debris flows. In addition, the slope of the frequency-magnitude distribution of post-fire debris flows is steeper than that of non-fire-related debris flows, evidence that differences in the post-fire environment tend to produce a higher proportion of small events. We propose two possible explanations: 1) post-fire events occur on shorter return intervals than debris flows in similar basins that do not experience fire, causing their distribution to shift toward smaller events due to limitations in sediment supply, or 2) fire causes changes in resisting and driving forces on a package of sediment, such that a smaller perturbation of the system is required in order for a debris flow to occur, resulting in smaller event volumes.

Surface velocity fields, digital elevation models and ice front positions derived from multi-mission SAR remote sensing data at Sjögren Inlet glaciers, Antarctic Peninsula

Substantial retreat or disintegration of numerous ice shelves have been observed on the Antarctic Peninsula. The ice shelf in the Prince Gustav Channel retreated gradually since the late 1980's and broke-up in 1995. Tributary glaciers reacted with speed-up, surface lowering and increased ice discharge, consequently contributing to sea level rise. We present a detailed long-term study (1993-2014) on the dynamic response of Sjögren Inlet glaciers to the disintegration of Prince Gustav Ice Shelf. We analyzed various remote sensing datasets to observe the reactions of the glaciers to the loss of the buttressing ice shelf. A strong increase in ice surface velocities was observed with maximum flow speeds reaching 2.82±0.48 m/d in 2007 and 1.50±0.32 m/d in 2004 at Sjögren and Boydell glaciers respectively. Subsequently, the flow velocities decelerated, however in late 2014, we still measured about two times the values of our first measurements in 1996. The tributary glaciers retreated 61.7±3.1 km² behind the former grounding line of the ice shelf. In regions below 1000 m a.s.l., a mean surface lowering of -68±10 m (-3.1 m/a) was observed in the period 1993-2014. The lowering rate decreased to -2.2 m/a in recent years. Based on the surface lowering rates, geodetic mass balances of the glaciers were derived for different time steps. High mass loss rate of -1.21±0.36 Gt/a was found in the earliest period (1993-2001). Due to the dynamic adjustments of the glaciers to the new boundary conditions the ice mass loss reduced to -0.59±0.11 Gt/a in the period 2012-2014, resulting in an average mass loss rate of -0.89±0.16 Gt/a (1993-2014). Including the retreat of the ice front and grounding line, a total mass change of -38.5±7.7 Gt and a contribution to sea level rise of 0.061±0.013 mm were computed. Analysis of the ice flux revealed that available bedrock elevation estimates at Sjögren Inlet are too shallow and are the major uncertainty in ice flux computations. This temporally dense time series analysis of Sjögren Inlet glaciers shows that the adjustments of tributary glaciers to ice shelf disintegration are still going on and provides detailed information of the changes in glacier dynamics.