Temperature, water level and bathymetry of thermokarst lakes in the continuous permafrost zone of northern Siberia - Lena River Delta, Siberia

Thermokarst lakes are typical features of the northern permafrost ecosystems, and play an important role in the thermal exchange between atmosphere and subsurface. The objective of this study is to describe the main thermal processes of the lakes and to quantify the heat exchange with the underlying sediments. The thermal regimes of five lakes located within the continuous permafrost zone of northern Siberia (Lena River Delta) were investigated using hourly water temperature and water level records covering a 3-year period (2009-2012), together with bathymetric survey data. The lakes included thermokarst lakes located on Holocene river terraces that may be connected to Lena River water during spring flooding, and a thermokarst lake located on deposits of the Pleistocene Ice Complex. Lakes were covered by ice up to 2 m thick that persisted for more than 7 months of the year, from October until about mid-June. Lake-bottom temperatures increased at the start of the ice-covered period due to upward-directed heat flux from the underlying thawed sediment. Prior to ice break-up, solar radiation effectively warmed the water beneath the ice cover and induced convective mixing. Ice break-up started at the beginning of June and lasted until the middle or end of June. Mixing occurred within the entire water column from the start of ice break-up and continued during the ice-free periods, as confirmed by the Wedderburn numbers, a quantitative measure of the balance between wind mixing and stratification that is important for describing the biogeochemical cycles of lakes. The lake thermal regime was modeled numerically using the FLake model. The model demonstrated good agreement with observations with regard to the mean lake temperature, with a good reproduction of the summer stratification during the ice-free period, but poor agreement during the ice-covered period. Modeled sensitivity to lake depth demonstrated that lakes in this climatic zone with mean depths > 5 m develop continuous stratification in summer for at least 1 month. The modeled vertical heat flux across the bottom sediment tends towards an annual mean of zero, with maximum downward fluxes of about 5 W/m**2 in summer and with heat released back into the water column at a rate of less than 1 W/m**2 during the ice-covered period. The lakes are shown to be efficient heat absorbers and effectively distribute the heat through mixing. Monthly bottom water temperatures during the ice-free period range up to 15 °C and are therefore higher than the associated monthly air or ground temperatures in the surrounding frozen permafrost landscape. The investigated lakes remain unfrozen at depth, with mean annual lake-bottom temperatures of between 2.7 and 4 °C.

Benthic and planktic foraminiferal stable isotope records combined with Nd isotope ratios of ODP Hole 198-1210B

The Pacific Ocean is the largest water body on Earth, and circulation in the Pacific contributed significantly to climate evolution in the latest Cretaceous, the culmination of a period of long-term cooling. Here, we present new high-resolution late Campanian to Maastrichtian benthic and planktic foraminiferal stable isotope data and a neodymium (Nd) isotope record obtained from sedimentary ferromanganese oxide coatings of Ocean Drilling Program Hole 1210B from the tropical Pacific Ocean (Shatsky Rise). These new records resolve 13 million years in the latest Cretaceous, providing insights into changes in surface and bottom water temperatures and source regions of deep to intermediate waters covering the carbon isotope excursions of the Campanian-Maastrichtian Boundary Event (CMBE) and the Mid-Maastrichtian event (MME). Our new benthic foraminiferal d18O and Nd isotope records together with published Nd isotope data show markedly parallel trends across the studied interval over a broad range of bathyal to abyssal water depths interpreted to reflect changes in the intensity of deep-ocean circulation in the tropical Pacific. In particular, we observe a three-million-year-long period of cooler conditions in the early Maastrichtian (72.5 to 69.5 Ma) when a concomitant change toward less radiogenic seawater Nd isotope signatures probably marks a period of enhanced admixture and northward flow of deep waters with Southern Ocean provenance. We suggest this change to have been triggered by intensified formation and convection of deep waters in the high southern latitudes, a process that weakened during the MME (69.5 to 68.5 Ma). The early Maastrichtian cold interval is closely related to the negative and positive carbon isotope trends of the CMBE and MME. The millions-of-years long duration of these carbon cycle perturbations suggests a tectonic forcing of climatic cooling, possibly related to changes in ocean basin geometry and bathymetry.

Flow separation and roughness lengths over large bedforms in a tidal environment

This study characterises the shape of the flow separation zone (FSZ) and wake region over large asymmetric bedforms under tidal flow conditions. High resolution bathymetry, flow velocity and turbulence data were measured along two parallel transects in a tidal channel covered with bedforms. The field data are used to verify the applicability of a numerical model for a systematic study using the Delft3D modelling system and test the model sensitivity to roughness length. Three experiments are then conducted to investigate how the FSZ size and wake extent vary depending on tidally-varying flow conditions, water levels and bathymetry. During the ebb, a large FSZ occurs over the steep lee side of each bedform. During the flood, no flow separation develops over the bedforms having a flat crest; however, a small FSZ is observed over the steepest part of the crest of some bedforms, where the slope is locally up to 15°. Over a given bedform morphology and constant water levels, no FSZ occurs for velocity magnitudes smaller than 0.1 m s**-1; as the flow accelerates, the FSZ reaches a stable size for velocity magnitudes greater than 0.4 m s**-1. The shape of the FSZ is not influenced by changes in water levels. On the other hand, variations in bed morphology, as recorded from the high-resolution bathymetry collected during the tidal cycle, influence the size and position of the FSZ: a FSZ develops only when the maximum lee side slope over a horizontal distance of 5 m is greater than 10°. The height and length of the wake region are related to the length of the FSZ. The total roughness along the transect lines is an order of magnitude larger during the ebb than during the flood due to flow direction in relation to bedform asymmetry: during the ebb, roughness is created by the large bedforms because a FSZ and wake develops over the steep lee side. The results add to the understanding of hydrodynamics of natural bedforms in a tidal environment and may be used to better parameterise small-scale processes in large-scale studies.

Eruption of the Håkon Mosby mud volcano recorded by the long-term observatory on mud-volcano eruptions (LOOME) between 2009 and 2010

Submarine mud volcanoes are considered an important source of methane to the water column. However, the temporal variability of their fluid transport including mud and methane emissions is largely unknown. Assuming that this transport was continuous and at steady state, methane emissions were previously proposed to result from a dynamic equilibrium between upward migration and consumption at the seabed by methane-consuming microbes. Here we have investigated non-steady state situations of vigorous mud movements and their reflection in fluid flow, seabed temperature and bathymetry. Time series of pressure, temperature, pH and seafloor photography were collected by a benthic observatory (LOOME) for 431 days at the active Håkon Mosby mud volcano. These new data document eruptions, which were accompanied by pulses of hot subsurface fluids and triggered rapid sediment uplift and lateral movement, as well as emissions of free gas.

Stable isotopic record of late Pliocene and Pleistocene foraminifera of ODP Site 114-704

We studied the stable isotopic and carbonate stratigraphy of ODP Hole 704A to reconstruct the paleoceanographic evolution of the eastern subantarctic sector of the South Atlantic Ocean. Site 704 is well positioned with respect to latitude (46°52.8'S, 7°25.3'E) and bathymetry (2532 m) to monitor past migrations in the position of Polar Front Zone (PFZ) and changes in deep-water circulation during the late Pliocene-Pleistocene. Several important changes occurred in proxy paleoceanographic indicators across the Gauss/Matuyama boundary at 2.47 Ma: (1) accumulation rates of biogenic sedimentary components increased by an order of magnitude (Froelich et al., this volume); (2) planktonic d1 8O values increased by an average of 0.5 per mil; (3) the amplitude of the benthic d18O signal increased; (4) the accumulation rate of ice-rafted detritus increased several fold (Warnke and Allen, this volume); and (5) carbon isotopic ratios of benthic foraminifers decreased by 0.5 per mil, as did the d13C of the fine-fraction carbonate by 1.5 per mil (Mead et al., 1991, doi:10.2973/odp.proc.sr.114.152.1991), but no change occurred in planktonic foraminiferal d13C values. Most of these changes are consistent with more frequent expansions and contractions of the PFZ over Site 704 after 2.47 Ma, bringing cold, nutrient-rich waters to 47°S that stimulated both carbonate and siliceous productivity. The synchronous increase in d18O values and ice-rafted detritus accumulation in Hole 704A indicates that the 2.4 Ma paleoceanographic event included ice volume growth on both Antarctica and Northern Hemisphere continents. The decrease in benthic d13C values indicates that the ventilation rate of Southern Ocean deep water decreased and the nutrient content increased during glacial events after 2.5 Ma. At the Gauss/Matuyama boundary, benthic d13C values of the Southern Ocean shifted toward those of the Pacific end member, indicating a decrease in the relative mixing ratio of Northern Component Water and Circumpolar Deep Water. During the early Matuyama (~2.3 to 1.7 Ma), the PFZ generally occupied a southerly position with respect to Site 704 and carbonate productivity prevailed. Exceptions to these general conditions occurred during strong glacial events of the early Matuyama (e.g., isotopic stages 82, 78, 74, and 70), when the PFZ migrated to the north and opal sedimentation predominated at Site 704. At 1.7 Ma, the PFZ migrated toward the equator and occupied a more northerly position for a prolonged interval between ~1.7 and 1.5 Ma. Beginning at ~1.5-1.4 Ma, surface and bottom water parameters (d18O, d13C, %CaCO3, and %opal) in the subantarctic South Atlantic became highly correlated such that glacial events (d18O maxima) corresponded to d13C and carbonate minima and opal maxima. This pattern is typical of the correlation found during the latest Pleistocene in the Southern Ocean (Charles and Fairbanks, in press). This event coincided with increased suppression of Northern Component Water during glacial events after 1.5 Ma (Raymo et al., 1990, doi:10.1016/0012-821X(90)90051-X), which may have influenced the climatology of the Southern Hemisphere by altering the flux of heat and salt to the Southern Ocean).

Distribution, density and abundance of Antarctic ice seals off Queen Maud Land and the eastern Weddell Sea

The Antarctic Pack Ice Seal (APIS) Program was initiated in 1994 to estimate the abundance of four species of Antarctic phocids: the crabeater seal Lobodon carcinophaga, Weddell seal Leptonychotes weddellii, Ross seal Ommatophoca rossii and leopard seal Hydrurga leptonyx and to identify ecological relationships and habitat use patterns. The Atlantic sector of the Southern Ocean (the eastern sector of the Weddell Sea) was surveyed by research teams from Germany, Norway and South Africa using a range of aerial methods over five austral summers between 1996-1997 and 2000-2001. We used these observations to model densities of seals in the area, taking into account haul-out probabilities, survey-specific sighting probabilities and covariates derived from satellite-based ice concentrations and bathymetry. These models predicted the total abundance over the area bounded by the surveys (30°W and 10°E). In this sector of the coast, we estimated seal abundances of: 514 (95 % CI 337-886) x 10**3 crabeater seals, 60.0 (43.2-94.4) x 10**3 Weddell seals and 13.2 (5.50-39.7) x 10**3 leopard seals. The crabeater seal densities, approximately 14,000 seals per degree longitude, are similar to estimates obtained by surveys in the Pacific and Indian sectors by other APIS researchers. Very few Ross seals were observed (24 total), leading to a conservative estimate of 830 (119-2894) individuals over the study area. These results provide an important baseline against which to compare future changes in seal distribution and abundance.

Seafloor spreading magnetic anomaly isochron map compilation for the Weddell Sea and Scotia Sea

The deep sea sedimentary record is an archive of the pre-glacial to glacial development of Antarctica and changes in climate, tectonics and ocean circulation. Identification of the pre-glacial, transitional and full glacial components in the sedimentary record is necessary for ice sheet reconstruction and to build circum-Antarctic sediment thickness grids for past topography and bathymetry reconstructions, which constrain paleoclimate models. A ~3300 km long Weddell Sea to Scotia Sea transect consisting of multichannel seismic reflection data from various organisations, were used to interpret new horizons to define the initial basin-wide seismostratigraphy and to identify the pre-glacial to glacial components. We mapped seven main units of which three are in the inferred Cretaceous-Paleocene pre-glacial regime, one in the Eocene-Oligocene transitional regime and three units in the Miocene-Pleistocene full glacial climate regime. Sparse borehole data from ODP leg 113 and SHALDRIL constrain the ages of the upper three units. Compiled seafloor spreading magnetic anomalies constrain the basement ages and the hypothetical age model. In many cases, the new horizons and stratigraphy contradict the interpretations in local studies. Each seismic sedimentary unit and its associated base horizon are continuous and traceable for the entire transect length, but reflect a lateral change in age whilst representing the same deposition process. The up to 1240 m thick pre-glacial seismic units form a mound in the central Weddell Sea basin and, in conjunction with the eroded flank geometry, support the interpretation of a Cretaceous proto-Weddell Gyre. The base reflector of the transitional seismic unit, which marks the initial ice sheet advances to the outer shelf, has a lateral model age of 26.6-15.5 Ma from southeast to northwest. The Pliocene-Pleistocene glacial deposits reveals lower sedimentations rates, indicating a reduced sediment supply. Sedimentation rates for the pre-glacial, transitional and full glacial components are highest around the Antarctic Peninsula, indicating higher erosion and sediment supply of a younger basement. We interpret an Eocene East Antarctic Ice Sheet expansion, Oligocene grounding of the West Antarctic Ice Sheet and Early Miocene grounding of the Antarctic Peninsula Ice Sheet.

Ocean temperature measured by elephant seal-CTDs in the Southern Ocean with links to datasets

The potential effects of ocean warming on marine predators are largely unknown, though the impact on the distribution of prey in vertical space may have far reaching impacts on diving predators such as southern elephant seals. We used data from satellite-tracked southern elephant seals from Marion Island to investigate the relationship between their dive characteristics (dive depths, dive durations and time-at-depth index values) and environmental variables (temperature at depth, depth of maximum temperature below 100 m, frontal zone and bathymetry) as well as other demographic and behavioural variables (migration stage, age-class, track day and vertical diel strategy). While other variables, such as bathymetry and vertical diel strategy also influenced dive depth, our results consistently indicated a significant influence of temperature at depth on dive depths. This relationship was positive for all groups of animals, indicating that seals dived to deeper depths when foraging in warmer waters. Female seals adjusted their dive depths proportionally more than males in warmer water. Dive durations were also influenced by temperature at depth, though to a lesser extent. Results from time-at-depth indices showed that both male and female seals spent less time at targeted dive depths in warmer water, and were presumably less successful foragers when diving in warmer water. Continued warming of the Southern Ocean may result in the distribution of prey for southern elephant seals shifting either poleward and/or to increasing depths. Marion Island elephant seals are expected to adapt their ranging and diving behaviour accordingly, though such changes may result in greater physiological costs associated with foraging.

Distribution of deep-sea foraminifera in surface sediments east of New Zealand

Paleobathymetric assessments of fossil foraminiferal faunas play a significant role in the analysis of the paleogeographic, sedimentary, and tectonic histories of New Zealand's Neogene marine sedimentary basins. At depths >100 m, these assessments often have large uncertainties. This study, aimed at improving the precision of paleodepth assessments, documents the present-day distribution of deep-sea foraminifera (>63 µm) in 66 samples of seafloor sediment at 90-700 m water depth (outer shelf to mid-abyssal), east of New Zealand. One hundred and thirty-nine of the 465 recorded species of benthic foraminifera are new records for the New Zealand region. Characters of the foraminiferal faunas which appear to provide the most useful information for estimating paleobathymetry are, in decreasing order of reliability: relative abundance of common benthic species; benthic species associations; upper depth limits of key benthic species; and relative abundance of planktic foraminifera. R mode cluster analysis on the quantitative census data of the 58 most abundant species of benthic foraminifera produced six species associations within three higher level clusters: (1) calcareous species most abundant at mid-bathyal to outer shelf depths (<1000 m); (2) calcareous species most abundant at mid-bathyal and greater depths (>600 m); (3) agglutinated species mostly occurring at deep abyssal depths (>3000 m). A detrended correspondence analysis ordination plot exhibits a strong relationship between these species associations and bathymetry. This is manifest in the bathymetric ranges of the relative abundance peaks of many of the common benthic species (e.g., Abditodentrix pseudothalmanni 500-2800 m, Bolivina robusta 200-650 m, Bulimina marginata f. marginata 20-600 m, B. marginata f. aculeata 400-3000 m, Cassidulina norvangi 1000-4500 m, Epistominella exigua 1000-4700 m, and Trifarina angulosa 10-650 m), which should prove useful in paleobathymetric estimates. The upper depth limits of 28 benthic foraminiferal species (e.g., Fursenkoina complanata 200 m, Bulimina truncana 450 m, Melonis affinis 550 m, Eggerella bradyi 750 m, and Cassidulina norvangi 1000 m) have potential to improve the precision of paleobathymetric estimates based initially on the total faunal composition. The planktic percentage of foraminiferal tests increases from outer shelf to upper abyssal depths followed by a rapid decline within the foraminiferal lysocline (below c. 3600 m). A planktic percentage <50% is suggestive of shelf depths, and >50% is suggestive of bathyal or abyssal depths above the CCD. In the abyssal zone there is dramatic taphonomic loss of most agglutinated tests (except some textulariids) at burial depths of 0.1-0.2 m, which negates the potential usefulness of these taxa in paleobathymetric assessments.