(Supplementary Table 1) Iceberg-rafted debris (IBRD) mass accumulation rates (MAR) of IODP site 318-U1361

The Pliocene and Early Pleistocene, between 5.3 and 0.8 million years ago, span a transition from a global climate state that was 2-3 °C warmer than present with limited ice sheets in the Northern Hemisphere to one that was characterized by continental-scale glaciations at both poles. Growth and decay of these ice sheets was paced by variations in the Earth's orbit around the Sun. However, the nature of the influence of orbital forcing on the ice sheets is unclear, particularly in light of the absence of a strong 20,000-year precession signal in geologic records of global ice volume and sea level. Here we present a record of the rate of accumulation of iceberg-rafted debris offshore from the East Antarctic ice sheet, adjacent to the Wilkes Subglacial Basin, between 4.3 and 2.2 million years ago. We infer that maximum iceberg debris accumulation is associated with the enhanced calving of icebergs during ice-sheet margin retreat. In the warmer part of the record, between 4.3 and 3.5 million years ago, spectral analyses show a dominant periodicity of about 40,000 years. Subsequently, the powers of the 100,000-year and 20,000-year signals strengthen. We suggest that, as the Southern Ocean cooled between 3.5 and 2.5 million years ago, the development of a perennial sea-ice field limited the oceanic forcing of the ice sheet. After this threshold was crossed, substantial retreat of the East Antarctic ice sheet occurred only during austral summer insolation maxima, as controlled by the precession cycle.

The use of direct boundary element method for gaining insight into complex seismic site response

The boundary element method is specially well suited for the analysis of the seismic response of valleys of complicated topography and stratigraphy. In this paper the method’s capabilities are illustrated using as an example an irregularity stratified (test site) sedimentary basin that has been modelled using 2D discretization and the Direct Boundary Element Method (DBEM). Site models displaying different levels of complexity are used in practice. The multi-layered model’s seismic response shows generally good agreement with observed data amplification levels, fundamental frequencies and the high spatial variability. Still important features such as the location of high frequencies peaks are missing. Even 2D simplified models reveal important characteristics of the wave field that 1D modelling does not show up.