Os and Re concentrations and Os isotope ratios for DSDP and ODP metalliferous sediments from the Pacific Ocean (Table 1)

We report new 187Os/186Os data and Re and Os concentrations in metalliferous sediments from the Pacific to construct a composite Os isotope seawater evolution curve over the past 80 m.y. Analyses of four samples of upper Cretaceous age yield 187Os/186Os values of between 3 and 6.5 and 187Re/186Os values below 55. Mass balance calculations indicate that the pronounced minimum of about 2 in the Os isotope ratio of seawater at the K-T boundary probably reflects the enormous input of cosmogenic material into the oceans by the K-T impactor(s). Following a rapid recovery to 187Os/186Os of 3.5 at 63 Ma, data for the early and middle part of the Cenozoic show an increase in 187Os/186Os to about 6 at 15 Ma. Variations in the isotopic composition of leachable Os from slowly accumulating metalliferous sediments show large fluctuations over short time spans. In contrast, analyses of rapidly accumulating metalliferous carbonates do not exhibit the large oscillations observed in the pelagic clay leach data. These results together with sediment leaching experiments indicate that dissolution of non-hydrogenous Os can occur during the hydrogen peroxide leach and demonstrate that Os data from pelagic clay leachates do not always reflect the Os isotopic composition of seawater. New data for the late Cenozoic further substantiate the rapid increase in the 187Os/186Os of seawater during the past 15 Ma. We interpret the correlation between the marine Sr and Os isotope records during this time period as evidence that weathering within the drainage basin of the Ganges-Brahmaputra river system is responsible for driving seawater Sr and Os toward more radiogenic isotopic compositions. The positive correlation between 87Sr/86Sr and U concentration, the covariation of U and Re concentrations, and the high dissolved Re, U and Sr concentrations found in the Ganges-Brahmaputra river waters supports this interpretation. Accelerating uplift of many orogens worldwide over the past 15 Ma, especially during the last 5 Ma, could have contributed to the rapid increase in 187Os/186Os from 6 to 8.5 over the past 15 Ma. Prior to 15 Ma the marine Sr and Os record are not tightly coupled. The heterogeneous distribution of different lithologies within eroding terrains may play an important role in decoupling the supplies of radiogenic Os and Sr to the oceans and account for the periods of decoupling of the marine Sr and Os isotope records.

Mg/Ca and stable oxygen measurement of the past 27 Ma of ODP Hole 120-747A

We explore the applicability of paired Mg/Ca and 18O/16O measurements on benthic foraminifera from Southern Ocean site 747 to paleoceanographic reconstructions on pre-Pleistocene timescales. We focus on the late Oligocene through Pleistocene (27-0 Ma) history of paleotemperatures and the evolution of the d18O values of seawater (d18Osw) at a temporal resolution of ~100-200 kyr. Absolute paleotemperature estimates depend on assumptions of how Mg/Ca ratios of seawater have changed over the past 27 Myr, but relative changes that occur on geologically brief timescales are robust. Results indicate that at the Oligocene to Miocene boundary (23.8 Ma), temperatures lag the increase in global ice-volume deduced from benthic foraminiferal d18O values, but the smaller-scale Miocene glaciations are accompanied by ocean cooling of -1°C. During the mid-Miocene phase of Antarctic ice sheet growth (~15-13 Ma), water temperatures cool by ~3°C. Unlike the benthic foraminiferal d18O values, which remain relatively constant thereafter, temperatures vary (by 3°C) and reach maxima at ~12 and ~8.5 Ma. The onset of significant Northern Hemisphere glaciation during the late Pliocene is synchronous with an ~4°C cooling at site 747. A comparison of our d18Osw curve to the Haq et al. (1987, doi:10.1126/science.235.4793.1156 ) sea level curve yields excellent agreement between sequence boundaries and times of increasing seawater 18O/16O ratios. At ~12-11 Ma in particular, when benthic foraminiferal d18O values do not support a further increase in ice volume, the d18Osw curve comes to a maximum that corresponds to a major mid-Miocene sea level regression. The agreement between the character of our Mg/Ca-based d18Osw curve and sequence stratigraphy demonstrates that benthic foramaniferal Mg/Ca ratios can be used to trace the d18Osw on pre-Pleistocene timescales despite a number of uncertainties related to poorly constrained temperature calibrations and paleoseawater Mg/Ca ratios. The Mg/Ca record also highlights that deep ocean temperatures can vary independently and unexpectedly from ice volume changes, which can lead to misinterpretations of the d18O record.

Organic geochemistry of water and bottom sediments during a seasonal flood in 2006 in the North Dvina River estuary

Heavy contaminant load released into the Northern Dvina River during flooding increased the concentrations of aliphatic (AHC) and polcyclic aromatic (PAH) hydrocarbons in water and bottom sediments. The composition of hydrocarbons was different from that of the summer low flow season. The concentrations of dissolved and particulate AHC ranged from 12 to 106 and from 192 to 599 µg/l, respectively, and bottom sediments contained from 26.2 to 329 µg/g AHC and 4 to 1785 ng/g PAH. As the transformation of AHC occurred at low spring temperatures, the alkane composition was shown to be dominated by terrigenous compounds, whereas more stable PAH showed elevated contents of petrogenic and pyrogenic compounds. It was also shown that the Northern Dvina-Dvina Bay geochemical barrier prevents contaminant input into the White Sea, i.e., acts as a marginal filter.

Concentrations and isotopic ratios of thorium and uranium in pelagic sediments

The predictable in situ production of 230Th from the decay of uranium in seawater, and its subsequent removal by scavenging onto falling particles, provides a valuable tool for normalizing fluxes to the seafloor. We describe a new application, determination of the 232Th that dissolves in the water column and is removed to the seafloor. 232Th is supplied to the ocean in continental minerals, dissolution of which leads to a measurable standing stock in the water column. Sedimentary adsorbed 232Th/230Th ratios have the potential to provide a proxy for estimating the amount of dissolved material that enters the ocean, both today and in the past. Ten core top samples were treated with up to eight different leaching techniques in order to determine the best method for the separating adsorbed from lattice bound thorium. In addition, separate components of the sediments were analyzed to test whether clay dissolution was an important contribution to the final measurement. There was no systematic correlation between the strength of acid used in the leach and the measured 232Th/230Th ratios. In all cases clean foraminifera produced the same ratio as leaches on bulk sediment. In three out of five samples leaches performed on non-carbonate detritus in the <63 µm size fraction were also identical. Without additional water column data it is not yet clear whether there is a simple one to one correlation between the expected deep-water 232Th/230Th and that produced by leaching, especially in carbonate-rich sediments. However, higher ratios, and associated high 232Th adsorbed fluxes, were observed in areas with high expected detrital inputs. The adsorbed fraction was ~35-50% of the total 232Th in seven out of ten samples. Our 230Th normalized 232Th fluxes are reasonable by comparison to global estimates of detrital inputs to the ocean. In nine cases out of ten, the total 230Th-normalized 232Th flux is greater than predicted from the annual dust fall at each specific location, but lower than the average global detrital input from all sources.

Uranium concentrations and isotope ratios Bahamas seawater and ODP holes pore-water

The geometry, timing, and rate of fluid-flow through carbonate margins and platforms is not well constrained. In this study, we use U concentrations and isotope ratios measured on small volumes of pore-water from Bahamas slope sediment, coupled with existing chlorinity data, to place constraints on the fluid-flow in this region and, by implication, other carbonate platforms. These data also allow an assessment of the behaviour of U isotopes in an unusually well constrained water-rock system. We report pore-water U concentrations which are controlled by dissolution of high-U organic material at shallow depths in the sediment and by reduction of U to its insoluble 4+ state at greater depths. The dominant process influencing pore-water (234U/238U) is alpha recoil. In Holocene sediments, the increase of pore-water (234U/238U) due to recoil provides an estimate of the horizontal flow rate of 11 cm/year, but with considerable uncertainty. At depths in the sediment where conditions are reducing, features in the U concentration and (234U/238U) profiles are offset from one another which constrains the effective diffusivity for U in these sediments to be c. 1-2 * 10**-8 cm**2/s. At depths between the Holocene and these reducing sediments, pore-water (234U/238U) values are unusually low due to a recent increase in the dissolution rate of grain surfaces. This suggests a strengthening of fluid flow, probably due to the flooding of the banks at the last deglaciation and the re-initiation of thermally-driven venting of fluid on the bank top and accompanying recharge on the slopes. Interpretation of existing chlorinity data, in the light of this change in flow rate, constrain the recent horizontal flow rate to be 10.6 ( 3.4) cm/year. Estimates of flow rate from (234U/238U) and Cl[-] are therefore in agreement and suggest flow rates close to those predicted by thermally-driven models of fluid flow. This agreement supports the idea that flow within the Bahamas Banks is mostly thermally driven and suggests that flow rates on the order of 10 cm/year are typical for carbonate platforms where such flow occurs.

Compositional and isotopic data on core MD84641

The organic carbon isotopic record of the sapropels(S1 and S3-S10) and intercalated marl oozes has been determined in a 12-m piston core from the eastern Mediterranean. The d13C_organic values are systematically lighter (mean=-21.0±0.82 per mil) in all sapropels and heavier (mean=-18.8±1.07 per mil) in the marl oozes. These differences are not due to variable marine and terrestrial organic matter mixtures because all values are heavier than modern plankton in the Mediterranean, there is no relationship between the C_organic/N ratios and the isotopic values, and published information on the abundance and distribution of organic biomarkers shows that terrestrial material constitutes a minor fraction of the total organic matter. Temperature effects on isotope fractionation are also discounted because the change in d13C_organic values between glacial and interglacial horizons is in the opposite sense. Diagenesis, which can produce relatively small changes in the carbon isotopic composition of sedimentary organic matter under certain circumstances, is unlikely to have caused the observed differences because this mechanism would cause an enrichmet in 12C, implying that all values were even heavier originally, and there is no secular trend in the d13C_organic record. The observed differences in d13C_organic between the two lithologies are probably produced by changes in the isotopic composition and the concentration of dissolved CO2. First, freshwater flooding during the formation of the sapropels caused the isotopic composition of the dissolved inorganic carbon in the surface waters of the Mediterranean to become lighter because of the 13C deficiency in fresh waters. Hence photosynthesis would have produced isotopically lighter organic material. Second, changes in atmospheric pCO2 between glacial and interglacial periods, as shown by the Vostok ice core, caused marked changes in the concentration of free dissolved CO2 in the mixed layer; lower values during glacial maxima caused a smaller fractionation of the carbon isotopes by phytoplankton, whereas levels were less limiting during the interglacials. Concentrations of dissolved CO2 could also have been much higher during the deposition of the sapropels because of the supply of regenerated CO2 to the mixed layer by upwelling, and this could have further lightened the d13C_organic values in the sapropels themselves. Carbon isotope records may provide an alternative method for estimating atmospheric pCO2 levels over longer time periods than can be obtained from ice cores.

Geochemistry of ODP Site 202-1233 and river sediments

Bulk sediment chemistry from three Chilean continental margin Ocean Drilling Program sites constrains regional continental erosion over the past 30,000 years. Sediments from thirteen rivers that drain the (mostly igneous) Andes and the (mostly metamorphic) Coast Range, along with existing rock chemistry datasets, define terrestrial provenance for the continental margin sediments. Andean river sediments have high Mg/Al relative to Coast-Range river sediments. Near 36°S, marine sediments have high-Mg/Al (i.e. more Andean) sources during the last glacial period, and lower-Mg/Al (less Andean) sources during the Holocene. Near 41°S a Ti-rich source, likely from coast-range igneous intrusions, is prevalent during Holocene time, whereas high-Mg/Al Andean sources are more prevalent during the last glacial period. We infer that there is a dominant ice-sheet control of sediment sources. At 36°S, Andean-sourced sediment decreased as Andean mountain glaciers retreated after ~17.6 ka, coincident with local oceanic warming and southward retreat of the Patagonian Forest and, by inference, westerly winds. At 41°S Andean sediment dominance peaks and then rapidly declines at ~19 ka, coincident with local oceanic warming and the earliest deglacial sea-level rise. We hypothesize that this decreased flux of Andean material in the south is related to rapid retreat of the marine-based portion of the Patagonian Ice Sheet in response to global sea-level rise, as the resulting flooding of the southern portion of the Central Valley created a sink for Andean sediments in this region. Reversal of the decreasing deglacial Mg/Al trend at 41°S from 14.5 to 13.0 ka is consistent with a brief re-advance of the Patagonian ice sheet coincident with the Antarctic Cold Reversal.

Sedimentation and mineral analysis on sediment cores from the Lena Delta

Core and outcrop analysis from Lena mouth deposits have been used to reconstruct the Late Quaternary sedimentation history of the Lena Delta. Sediment properties (heavy mineral composition, grain size characteristics, organic carbon content) and age determinations (14C AMS and IR-OSL) are applied to discriminate the main sedimentary units of the three major geomorphic terraces, which form the delta. The development of the terraces is controlled by complex interactions among the following four factors: (1) Channel migration. According to the distribution of 14C and IR-OSL age determinations of Lena mouth sediments, the major river runoff direction shifted from the west during marine isotope stages 5-3 (third terrace deposits) towards the northwest during marine isotope stage 2 and transition to stage 1 (second terrace), to the northeast and east during the Holocene (first terrace deposits). (2) Eustasy. Sea level rise from Last Glacial lowstand to the modern sea level position, reached at 6-5 ka BP, resulted in back-filling and flooding of the palaeovalleys. (3) Neotectonics. The extension of the Arctic Mid-Ocean Ridge into the Laptev Sea shelf acted as a halfgraben, showing dilatation movements with different subsidence rates. From the continent side, differential neotectonics with uplift and transpression in the Siberian coast ridges are active. Both likely have influenced river behavior by providing sites for preservation, with uplift, in particular, allowing accumulation of deposits in the second terrace in the western sector. The actual delta setting comprises only the eastern sector of the Lena Delta. (4) Peat formation. Polygenetic formation of ice-rich peaty sand (''Ice Complex'') was most extensive (7-11 m in thickness) in the southern part of the delta area between 43 and 14 ka BP (third terrace deposits). In recent times, alluvial peat (5-6 m in thickness) is accumulated on top of the deltaic sequences in the eastern sector (first terrace).

Distribution of Florisphaera profunda and Helicosphaera carteri in down core and surface sediments

A detailed assessment of the respective roles of production, export, and subsequent preservation of organic carbon (Corg) in the eastern Mediterranean (EMED) sediments during the formation of sapropels remains elusive. Here we present new micropaleontological results for both surface samples taken at several locations in the EMED and last interglacial sapropel S5 from core LC21 in the southeastern Aegean Sea. A strong exponential anticorrelation between relative abundances of the lower photic zone coccolithophore Florisphaera profundain the surface sediments and modern concentrations of chlorophyll a (Chl-a) at the sea surface suggests thatF. profunda percentages can be used to track past productivity changes in the EMED. Prior to S5 deposition, an abrupt and large increase of F. profunda percentages in LC21 coincided (within the multidecadal resolution of the records) with the marked freshening of EMED surface waters. This suggests a strong coupling between freshwater-bound surface to intermediate water (density) stratification and enhanced upward advection of nutrients to the base of the photic zone, fuelling a productive deep chlorophyll maximum (DCM) underneath a nutrient-starved surface layer. Our findings imply that (at least) at the onset of sapropel formation physical and biogeochemical processes likely operated in tandem, enabling high Corg accumulation at the seafloor.

Lead and strontium isotope composition, mass accumulation rates and rare earth elements of selected samples from DSDP Sites 92-597 to 92-601

Most of the Pb isotope data for the Leg 92 metalliferous sediments (carbonate-free fraction) form approximately linear arrays in the conventional isotopic plots, extending from the middle of the field for mid-ocean ridge basalts (MORB) toward the field for Mn nodules. These arrays are directed closely to the average values of Mn nodules, the composition of which reflects the Pb isotope composition of seawater (Reynolds and Dasch, 1971). Since the Leg 92 samples are almost devoid of continentally derived detritus, it can be inferred that the more radiogenic end-member is seawater. The less radiogenic end-member lies in the very middle of the MORB field, and hence can be considered to reflect the Pb isotope composition of typical ocean-ridge basalt. The array of data lying between these two end-members is most readily interpreted in terms of simple linear mixing of Pb from the two different end-member sources. According to this model, eight samples from Sites 599 to 601 contain 50 to 100% basaltic Pb. Five of these samples have compositions that are identical within the uncertainty of the analyses. We use the average of these five values to define our unradiogenic end-member in the linear mixing model. The ratios used for this average are 206Pb/204Pb = 18.425 ± 0.010; 207Pb/204Pb = 15.495 ± 0.018; 208Pb/204Pb = 37.879 ± 0.068. These values should approximate the average Pb isotope composition of discharging hydrothermal solutions, and therefore also that of the basaltic crust, over the period of time represented by these samples ( 4 m.y., from 4 to 8 Ma). Sr isotope ratios show a significant range of values, from 0.7082 to 0.7091. The lower ratios are well outside the value of 0.70910 ± 6 for modern-day seawater (Burke et al., 1982). However, most values correspond very closely to the curve of 87Sr/86Sr versus age for seawater, with older samples having progressively lower 87Sr/86Sr ratios. The simplest explanation for this progressive reduction is that recrystallization of the abundant biogenic carbonate in the sediments released older seawater Sr which was incorporated into ferromanganiferous phases during diagenesis. Leg 92 metalliferous sediments have total rare earth element (REE) contents that range on a carbonate-free basis from 131 to 301 ppm, with a clustering between 167 and 222 ppm. The patterns have strong negative Ce anomalies. Samples from Sites 599 to 601 display a slight but distinct enrichment in the heavy REE relative to the light REE, whereas those from Sites 597 to 598 show almost no heavy REE enrichment. The former patterns (those for Sites 599 to 601) are interpreted as indicating moderate diagenetic alteration of metalliferous sediments originating at the EPR axis; the latter reflect more complete diagenetic modification.

Chemical and isotope compositions of main phases of metalliferous sediments from holes of DSDP Leg 92

The carbonate-free metalliferous fraction of thirty-nine sediment samples from four DSDP Leg 92 sites has been analyzed for 12 elements, and a subset of 16 samples analyzed for Pb isotopic composition. The main geochemical features of this component are as follows: i) very high concentrations of Fe and Mn, typically 25-39% and 5-14%, respectively; ii) Al and Ca contents generally less than 2% and 5%, respectively; iii) high Cu (1000-2000 ppm), and Zn and Ni (500-1000 ppm) values; and iv) Co and Pb concentrations of 100-250 ppm. In terms of element partitioning within the metalliferous fraction, amorphous to poorly crystallized oxide-oxyhydroxides removed by the second leach carry virtually all of the Mn, and about 90% of the Ca, Sr and Ni. The well-crystallized goethite-rich material removed by the third leach carries the majority of Fe, Cu, and Pb. These relations hold for sediments as young as ~1-2 Ma, indicating early partitioning of hydrothermal Fe and Mn into separate phases. Calculated mass accumulation rates (MAR) for Fe, Mn, Cu, Pb, Zn and Ni in the bulk sediment show the same overall trends at three of the sites, with greatest MAR values near the basement, and a general decrease in MAR values towards the tops of the holes (for sediments deposited above the lysocline). These relations strongly support the concept of a declining hydrothermal contribution of these elements away from a ridge axis. Nevertheless, MAR values for these metals up to ~200 km from the ridge axis are orders of magnitude higher than on abyssal seafloor plains where there is no hydrothermal influence. Mn/Fe ratios throughout the sediment column at two sites indicate that the composition of the hydrothermal precipitates changed during transport through seawater, becoming significantly depleted in Mn beyond ~200-300 km from the axis, but maintaining roughly the same proportion of Fe. Most of the Pb isotope data for the Leg 92 metalliferous sediments form approximately linear arrays in the conventional isotopic plots, extending from the middle of the field for mid-ocean ridge basalts toward the field for Mn nodules. The array of data lying between these two end-members is most readily interpreted in terms of simple linear mixing of Pb derived from basaltic and seawater end-member sources. The least radiogenic sediments reflect the average Pb isotope composition of discharging hydrothermal solutions and ocean-ridge basalt at the EPR over the ~4-8 Ma B.P. interval. Pb in sediments deposited up to 250 km from the axis can be almost entirely of basaltic-hydrothermal origin. Lateral transport of some basaltic Pb by ocean currents appears to extend to distances of at least 1000 km west of the East Pacific Rise.

Geochemistry, isotopic ratios, granulometry and Uk'37 derived SST of cores obtained during Poseidon cruise POS287, off Portugal

Continental and marine conditions during the last millennium off Porto, Portugal (the southern pole of the North Atlantic Oscillation, NAO), are reconstructed from a sediment archive through a high-resolution multiproxy study and instrumental evidence. Results show multidecadal variability and sea surface temperatures (SSTs) that correlate well with previously published land and sea-based Northern Hemisphere temperature records, and appear to be responding to long-term solar insolation variability. Precipitation was negatively correlated with the NAO, whereas strong flooding events occurred at times of marked climate cooling (AD 1100-1150 and 1400-1470) and transitions in solar activity. AD 1850 marks a major shift in the phytoplankton community associated with a decoupling of d18O records of 3 planktonic foraminifera species. These changes are interpreted as a response to a reduction in the summer and/or annual upwelling and more frequent fall-winter upwelling-like events. This shift's coincidence with a decrease in SST and the increase in coherence between our data and the Atlantic Multidecadal Oscillation (AMO) confirms the connection of the upwelling variability to the North Atlantic Ocean's surface and thermohaline circulation on a decadal scale. The disappearance of this agreement between the AMO and our records beyond AD 1850 and its coincidence with the beginning of the recent rise in atmospheric CO2 supports the hypothesis of a strong anthropogenic effect on the last ~150 yr of the climate record. Furthermore, it raises an important question of the use of instrumental records as the sole calibration data set for climate reconstructions, as these may not provide the best analogue for climate beyond AD 1730.

Benthic foraminifera, stable isotope record and sedimentology of Holocene sediments in the Skagerrak

A high-resolution multi-proxy study of core MD99-2286 reveals a highly variable hydrographic environment in the Skagerrak from 9300 cal. yr BP to the present. The study includes foraminiferal faunas, stable isotopes and sedimentary parameters, as well as temperature and salinity reconstructions of a ca. 29 m long radiocarbon-dated core record. The multivariate technique fuzzy c-means was applied to the foraminiferal counts, and it was extremely valuable in defining subtle heterogeneities in the foraminiferal fauna data corresponding to hydrographic changes. The major mid-Holocene (Littorina) transgression, led to flooding of large former land areas in the North Sea, the opening of the English Channel and Danish straits and initiation of the modern circulation system. This is reflected by fluctuating C/N values and an explosive bloom of Hyalinea balthica. A slight indication of ameliorated conditions between 8000-5750 cal. yr BP is related to the Holocene Thermal Maximum. A subsequent increase in fresh water/Baltic water influence between 5750-4350 cal. yr BP is reflected by dominance of Bulimina marginata and depleted d18O-values. The Neoglacial cooling (after 4350 cal. yr BP) is seen in the Skagerrak as enhanced turbidity, increasing TOC-values and short-term changes in an overall Cassidulina laevigata dominated fauna suggesting a prevailing influence of Atlantic waters. This is in agreement with increased strength of westerly winds, as recorded for this period. The last 2000 years were also dominated by Atlantic Water conditions with generally abundant nutrient supply. However, during warm periods, particularly the Medieval Warm Period and the modern warming, the area was subject to a restriction in the supply of nutrients and/or the nutrient supply had a more refractory character.

GPS raw data (control points and ground control points) from the Liguria Region, Borghetto Santo Spirito, Italy

Monitoring the impact of sea storms on coastal areas is fundamental to study beach evolution and the vulnerability of low-lying coasts to erosion and flooding. Modelling wave runup on a beach is possible, but it requires accurate topographic data and model tuning, that can be done comparing observed and modeled runup. In this study we collected aerial photos using an Unmanned Aerial Vehicle after two different swells on the same study area. We merged the point cloud obtained with photogrammetry with multibeam data, in order to obtain a complete beach topography. Then, on each set of rectified and georeferenced UAV orthophotos, we identified the maximum wave runup for both events recognizing the wet area left by the waves. We then used our topography and numerical models to simulate the wave runup and compare the model results to observed values during the two events. Our results highlight the potential of the methodology presented, which integrates UAV platforms, photogrammetry and Geographic Information Systems to provide faster and cheaper information on beach topography and geomorphology compared with traditional techniques without losing in accuracy. We use the results obtained from this technique as a topographic base for a model that calculates runup for the two swells. The observed and modeled runups are consistent, and open new directions for future research.

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.