Time variable Earth’s gravity field from SLR satellites

The time variable Earth’s gravity field contains information about the mass transport within the system Earth, i.e., the relationship between mass variations in the atmosphere, oceans, land hydrology, and ice sheets. For many years, satellite laser ranging (SLR) observations to geodetic satellites have provided valuable information of the low-degree coefficients of the Earth’s gravity field. Today, the Gravity Recovery and Climate Experiment (GRACE) mission is the major source of information for the time variable field of a high spatial resolution. We recover the low-degree coefficients of the time variable Earth’s gravity field using SLR observations up to nine geodetic satellites: LAGEOS-1, LAGEOS-2, Starlette, Stella, AJISAI, LARES, Larets, BLITS, and Beacon-C. We estimate monthly gravity field coefficients up to degree and order 10/10 for the time span 2003–2013 and we compare the results with the GRACE-derived gravity field coefficients. We show that not only degree-2 gravity field coefficients can be well determined from SLR, but also other coefficients up to degree 10 using the combination of short 1-day arcs for low orbiting satellites and 10-day arcs for LAGEOS-1/2. In this way, LAGEOS-1/2 allow recovering zonal terms, which are associated with long-term satellite orbit perturbations, whereas the tesseral and sectorial terms benefit most from low orbiting satellites, whose orbit modeling deficiencies are minimized due to short 1-day arcs. The amplitudes of the annual signal in the low-degree gravity field coefficients derived from SLR agree with GRACE K-band results at a level of 77 %. This implies that SLR has a great potential to fill the gap between the current GRACE and the future GRACE Follow-On mission for recovering of the seasonal variations and secular trends of the longest wavelengths in gravity field, which are associated with the large-scale mass transport in the system Earth.

A systematized workflow for assessing land resources using readily available spatial datasets

Land degradation as well as land conservation maps at a (sub-) national scale are critical for pro-ject planning for sustainable land management. It has long been recognized that online accessible and low-cost raster data sets (e.g. Landsat imagery, SRTM-DEM’s) provide a readily available basis for land resource assessments for developing countries. However, choice of spatial, tempo-ral and spectral resolution of such data is often limited. Furthermore, while local expert knowl-edge on land degradation processes is abundant, difficulties are often encountered when linking existing knowledge with modern approaches including GIS and RS. The aim of this study was to develop an easily applicable, standardized workflow for preliminary spatial assessments of land degradation and conservation, which also allows the integration of existing expert knowledge. The core of the developed method consists of a workflow for rule-based land resource assess-ment. In a systematic way, this workflow leads from predefined land degradation and conserva-tion classes to field indicators, to suitable spatial proxy data, and finally to a set of rules for clas-sification of spatial datasets. Pre-conditions are used to narrow the area of interest. Decision tree models are used for integrating the different rules. It can be concluded that the workflow presented assists experts from different disciplines in col-laboration GIS/RS specialists in establishing a preliminary model for assessing land degradation and conservation in a spatially explicit manner. The workflow provides support when linking field indicators and spatial datasets, and when determining field indicators for groundtruthing.

Using MDEFT MRI Sequences to Target the GPi in DBS Surgery.

OBJECTIVE Recent advances in different MRI sequences have enabled direct visualization and targeting of the Globus pallidus internus (GPi) for DBS surgery. Modified Driven Equilibrium Fourier Transform (MDEFT) MRI sequences provide high spatial resolution and an excellent contrast of the basal ganglia with low distortion. In this study, we investigate if MDEFT sequences yield accurate and reliable targeting of the GPi and compare direct targeting based on MDEFT sequences with atlas-based targeting. METHODS 13 consecutive patients considered for bilateral GPi-DBS for dystonia or PD were included in this study. Preoperative targeting of the GPi was performed visually based on MDEFT sequences as well as by using standard atlas coordinates. Postoperative CT imaging was performed to calculate the location of the implanted leads as well as the active electrode(s). The coordinates of both visual and atlas based targets were compared. The stereotactic coordinates of the lead and active electrode(s) were calculated and projected on the segmented GPi. RESULTS On MDEFT sequences the GPi was well demarcated in most patients. Compared to atlas-based planning the mean target coordinates were located significantly more posterior. Subgroup analysis showed a significant difference in the lateral coordinate between dystonia (LAT = 19.33 ± 0.90) and PD patients (LAT = 20.67 ± 1.69). Projected on the segmented preoperative GPi the active contacts of the DBS electrode in both dystonia and PD patients were located in the inferior and posterior part of the structure corresponding to the motor part of the GPi. CONCLUSIONS MDEFT MRI sequences provide high spatial resolution and an excellent contrast enabling precise identification and direct visual targeting of the GPi. Compared to atlas-based targeting, it resulted in a significantly different mean location of our target. Furthermore, we observed a significant variability of the target among the PD and dystonia subpopulation suggesting accurate targeting for each individual patient.

Understanding measured water rotational temperatures and column densities in the very innermost coma of Comet 73P/Schwassmann–Wachmann 3 B

Direct sublimation of a comet nucleus surface is usually considered to be the main source of gas in the coma of a comet. However, evidence from a number of comets including the recent spectacular images of Comet 103P/Hartley 2 by the EPOXI mission indicates that the nucleus alone may not be responsible for all, or possibly at times even most, of the total amount of gas seen in the coma. Indeed, the sublimation of icy grains, which have been injected into the coma, appears to constitute an important source. We use the fully-kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J., 685, 659−677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J., 732) to reproduce the measurements of column density and rotational temperature of water in Comet 73P-B/Schwassmann–Wachmann 3 obtained with a very high spatial resolution of ∼30 km using IRCS/Subaru in May 2006 (Bonev, B.P., Mumma, M.J., Kawakita, H., Kobayashi, H., Villanueva, G.L. [2008]. Icarus, 196, 241−248). For gas released solely from the cometary nucleus at a heliocentric distance of 1 AU, modeled rotational temperatures start at 110 K close to the surface and decrease to only several tens of degrees by 10–20 nucleus radii. However, the measured decay of both rotational temperature and column density with distance from the nucleus is much slower than predicted by this simple model. The addition of a substantial (distributed) source of gas from icy grains in the model slows the decay in rotational temperature and provides a more gradual drop in column density profiles. Together with a contribution of rotational heating of water molecules by electrons, the combined effects allow a much better match to the IRCS/Subaru observations. From the spatial distributions of water abundance and temperature measured in 73P/SW3-B, we have identified and quantified multiple mechanisms of release. The application of this tool to other comets may permit such studies over a range of heliocentric and geocentric distances.

Optical recording of calcium currents during impulse conduction in cardiac tissue

We explore the feasibility of obtaining a spatially resolved picture of Ca2+Ca2+ inward currents (ICaICa) in multicellular cardiac tissue by differentiating optically recorded Ca2+Ca2+ transients that accompany propagating action potentials. Patterned growth strands of neonatal rat ventricular cardiomyocytes were stained with the Ca2+Ca2+ indicators Fluo-4 or Fluo-4FF. Preparations were stimulated at 1 Hz, and Ca2+Ca2+ transients were recorded with high spatiotemporal resolution (50  μm50  μm, 2 kHz analog bandwidth) with a photodiode array. Signals were differentiated after appropriate digital filtering. Differentiation of Ca2+Ca2+ transients resulted in optically recorded calcium currents (ORCCs) that carried the temporal and pharmacological signatures of L-type Ca2+Ca2+ inward currents: the time to peak amounted to ∼2.1  ms∼2.1  ms (Fluo-4FF) and ∼2.4  ms∼2.4  ms (Fluo-4), full-width at half-maximum was ∼8  ms∼8  ms, and ORCCs were completely suppressed by 50  μmol/L50  μmol/LCdCl2CdCl2. Also, and as reported before from patch-clamp studies, caffeine reversibly depressed the amplitude of ORCCs. The results demonstrate that the differentiation of Ca2+Ca2+ transients can be used to obtain a spatially resolved picture of the initial phase of ICaICa in cardiac tissue and to assess relative changes of activation/fast inactivation of ICaICa following pharmacological interventions.

An automated method for mapping physical soil and water conservation structures on cultivated land using GIS and remote sensing techniques

An efficient and reliable automated model that can map physical Soil and Water Conservation (SWC) structures on cultivated land was developed using very high spatial resolution imagery obtained from Google Earth and ArcGIS, ERDAS IMAGINE, and SDC Morphology Toolbox for MATLAB and statistical techniques. The model was developed using the following procedures: (1) a high-pass spatial filter algorithm was applied to detect linear features, (2) morphological processing was used to remove unwanted linear features, (3) the raster format was vectorized, (4) the vectorized linear features were split per hectare (ha) and each line was then classified according to its compass direction, and (5) the sum of all vector lengths per class of direction per ha was calculated. Finally, the direction class with the greatest length was selected from each ha to predict the physical SWC structures. The model was calibrated and validated on the Ethiopian Highlands. The model correctly mapped 80% of the existing structures. The developed model was then tested at different sites with different topography. The results show that the developed model is feasible for automated mapping of physical SWC structures. Therefore, the model is useful for predicting and mapping physical SWC structures areas across diverse areas.

Palynostratigraphy of the last centuries in Switzerland based on 23 lake and mire deposits: chronostratigraphic pollen markers, regional patterns, and local histories

A total of 23 pollen diagrams [stored in the Alpine Palynological Data-Base (ALPADABA), Geobotanical Institute, Bern] cover the last 100 to over 1000 years. The sites include 15 lakes, seven mires, and one soil profile distributed in the Jura Mts (three sites), Swiss Plateau (two sites), northern Pre-Alps and Alps (six sites), central Alps (five sites), southern Alps (three sites), and southern Pre-Alps (four sites) in the western and southern part of Switzerland or just outside the national borders. The pollen diagrams have both a high taxonomic resolution and a high temporal resolution, with sampling distances of 0.5–3 cm, equivalent to 1 to 11 years for the last 100 years and 8 to 130 years for earlier periods. The chronology is based on absolute dating (14 sites: 210Pb 11 sites; 14C six sites; varve counting two sites) or on biostratigraphic correlation among pollen diagrams. The latter relies mainly on trends in Cannabis sativa, Ambrosia, Mercurialis annua, and Ostrya-type pollen. Individual pollen stratigraphies are discussed and sites are compared within each region. The principle of designating local, extra-local, and regional pollen signals and vegetation is exemplified by two pairs of sites lying close together. Trends in biostratigraphies shared by a major part of the pollen diagrams allow the following generalisations. Forest declined in phases since medieval times up to the late 19th century. Abies and Fagus declined consistently, whereas the behaviour of short-lived trees and trees of moist habitats differed among sites (Alnus glutinosa-type, Alnus viridis, Betula, Corylus avellana). In the present century, however, Picea and Pinus increased, followed by Fraxinus excelsior in the second half of this century. Grassland (traced by Gramineae and Plantago lanceolata-type pollen) increased, replacing much of the forest, and declined again in the second half of this century. Nitrate enrichment of the vegetation (traced by Urtica) took place in the first half of this century. These trends reflect the intensification of forest use and the expansion of grassland from medieval times up to the end of the last century, whereas subsequently parts of the grassland became used more intensively and the marginal parts were abandoned for forest regrowth. In most pollen diagrams human impact is the dominant factor in explaining inferred changes in vegetation, but climatic change plays a role at three sites.

Myelination of the right parahippocampal cingulum is associated with physical activity in young healthy adults.

Recent evidence suggests that individual differences in physical activity (PA) may be associated with individual differences in white matter microstructure and with grey matter volume of the hippocampus. Therefore, this study investigated the association between PA and white matter microstructure of pathways connecting to the hippocampus. A total of 33 young, healthy adults underwent magnetic resonance imaging (MRI). High angular resolution diffusion-weighted imaging and multi-component relaxometry MRI scans (multi-component driven equilibrium pulse observation of T1 and T2) were acquired for each participant. Activity levels (AL) of participants were calculated from 72-h actigraphy recordings. Tractography using the damped Richardson Lucy algorithm was used to reconstruct the fornix and bilateral parahippocampal cinguli (PHC). The mean fractional anisotropy (FA) and the myelin water fraction (MWF), a putative marker of myelination, were determined for each pathway. A positive correlation between both AL and FA and between AL and MWF were hypothesized for the three pathways. There was a selective positive correlation between AL and MWF in the right PHC (r = 0.482, p = 0.007). Thus, our results provide initial in vivo evidence for an association between myelination of the right PHC and PA in young healthy adults. Our results suggest that MWF may not only be more specific, but also more sensitive than FA to detect white matter microstructural alterations. If PA was to induce structural plasticity of the right PHC this may contribute to reverse structural alterations of the right PHC in neuropsychiatric disorder with hippocampal pathologies.

Geomorphology of the Imhotep region on comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Context. Since August 2014, the OSIRIS Narrow Angle Camera (NAC) onboard the Rosetta spacecraft has acquired high spatial resolution images of the nucleus of comet 67P/Churyumov-Gerasimenko, down to the decimeter scale. This paper focuses on the Imhotep region, located on the largest lobe of the nucleus, near the equator. Aims. We map, inventory, and describe the geomorphology of the Imhotep region. We propose and discuss some processes to explain the formation and ongoing evolution of this region. Methods. We used OSIRIS NAC images, gravitational heights and slopes, and digital terrain models to map and measure the morphologies of Imhotep. Results. The Imhotep region presents a wide variety of terrains and morphologies: smooth and rocky terrains, bright areas, linear features, roundish features, and boulders. Gravity processes such as mass wasting and collapse play a significant role in the geomorphological evolution of this region. Cometary processes initiate erosion and are responsible for the formation of degassing conduits that are revealed by elevated roundish features on the surface. We also propose a scenario for the formation and evolution of the Imhotep region; this implies the presence of large primordial voids inside the nucleus, resulting from its formation process.

NOVEL PHANTOMS AND POST-PROCESSING FOR DIFFUSION SPECTRUM IMAGING

High Angular Resolution Diffusion Imaging (HARDI) techniques, including Diffusion Spectrum Imaging (DSI), have been proposed to resolve crossing and other complex fiber architecture in the human brain white matter. In these methods, directional information of diffusion is inferred from the peaks in the orientation distribution function (ODF). Extensive studies using histology on macaque brain, cat cerebellum, rat hippocampus and optic tracts, and bovine tongue are qualitatively in agreement with the DSI-derived ODFs and tractography. However, there are only two studies in the literature which validated the DSI results using physical phantoms and both these studies were not performed on a clinical MRI scanner. Also, the limited studies which optimized DSI in a clinical setting, did not involve a comparison against physical phantoms. Finally, there is lack of consensus on the necessary pre- and post-processing steps in DSI; and ground truth diffusion fiber phantoms are not yet standardized. Therefore, the aims of this dissertation were to design and construct novel diffusion phantoms, employ post-processing techniques in order to systematically validate and optimize (DSI)-derived fiber ODFs in the crossing regions on a clinical 3T MR scanner, and develop user-friendly software for DSI data reconstruction and analysis. Phantoms with a fixed crossing fiber configuration of two crossing fibers at 90° and 45° respectively along with a phantom with three crossing fibers at 60°, using novel hollow plastic capillaries and novel placeholders, were constructed. T2-weighted MRI results on these phantoms demonstrated high SNR, homogeneous signal, and absence of air bubbles. Also, a technique to deconvolve the response function of an individual peak from the overall ODF was implemented, in addition to other DSI post-processing steps. This technique greatly improved the angular resolution of the otherwise unresolvable peaks in a crossing fiber ODF. The effects of DSI acquisition parameters and SNR on the resultant angular accuracy of DSI on the clinical scanner were studied and quantified using the developed phantoms. With a high angular direction sampling and reasonable levels of SNR, quantification of a crossing region in the 90°, 45° and 60° phantoms resulted in a successful detection of angular information with mean ± SD of 86.93°±2.65°, 44.61°±1.6° and 60.03°±2.21° respectively, while simultaneously enhancing the ODFs in regions containing single fibers. For the applicability of these validated methodologies in DSI, improvement in ODFs and fiber tracking from known crossing fiber regions in normal human subjects were demonstrated; and an in-house software package in MATLAB which streamlines the data reconstruction and post-processing for DSI, with easy to use graphical user interface was developed. In conclusion, the phantoms developed in this dissertation offer a means of providing ground truth for validation of reconstruction and tractography algorithms of various diffusion models (including DSI). Also, the deconvolution methodology (when applied as an additional DSI post-processing step) significantly improved the angular accuracy of the ODFs obtained from DSI, and should be applicable to ODFs obtained from the other high angular resolution diffusion imaging techniques.

External assessment of myocardial glucose metabolism with $\sp{18}$F-2 -deoxy-2 -fluoro-D-glucose

Despite the popularity of the positron emitting glucose analog, ($\sp{18}$F) -2-deoxy-2-fluoro-D-glucose (2FDG), for the noninvasive "metabolic imaging" of organs with positron emission tomography (PET), the physiological basis for the tracer has not been tested, and the potential of 2FDG for the rapid kinetic analysis of altered glucose metabolism in the intact heart has not been fully exploited. We, therefore, developed a quantitative method to characterize metabolic changes of myocardial glucose metabolism noninvasively and with high temporal resolution.^ The first objective of the work was to provide direct evidence that the initial steps in the metabolism of 2FDG are the same as for glucose and that 2FDG is retained by the tissue in proportion to the rate of glucose utilization. The second objective was to characterize the kinetic changes in myocardial glucose transport and phosphorylation in response to changes in work load, competing substrates, acute ischemia and reperfusion, and the addition of insulin. To assess changes in myocardial glucose metabolism isolated working rat hearts were perfused with glucose and 2FDG. Tissue uptake of 2FDG and the input function were measured on-line by external detection. The steady state rate of 2FDG phosphorylation was determined by graphical analysis of 2FDG time-activity curves.^ The rate of 2FDG uptake was linear with time and the tracer was retained in its phosphorylated form. Tissue accumulation of 2FDG decreased within seconds with a reduction in work load, in the presence of competing substrates, and during reperfusion after global ischemia. Thus, most interventions known to alter glucose metabolism induced rapid parallel changes in 2FDG uptake. By contrast, insulin caused a significant increase in 2FDG accumulation only in hearts from fasted animals when perfused at a sub-physiological work load. The mechanism for this phenomenon is not known but may be related to the existence of two different glucose transporter systems and/or glycogen metabolism in the myocardial cell.^ It is concluded that (1) 2FDG traces glucose uptake and phosphorylation in the isolated working rat heart; and (2) early and transient kinetic changes in glucose metabolism can be monitored with high temporal resolution with 2FDG and a simple positron coincidence counting system. The new method has revealed transients of myocardial glucose metabolism, which would have remained unnoticed with conventional methods. These transients are not only important for the interpretation of glucose metabolic PET scans, but also provide insights into mechanisms of glucose transport and phosphorylation in heart muscle. ^

Retreat of top cliff of Kurungnakh Island, Lena Delta, Siberia, Russia, 2010-2014, with links to shapefiles

Thawing-induced cliff top retreat in permafrost landscapes is mainly due to thermo-erosion. Ground-ice-rich permafrost landscapes are specifically vulnerable to thermo-erosion and may show high degradation rates. Within the HGF Alliance Remote Sensing and the FP7 PAGE21 permafrost programs we investigated how SAR and optical remote sensing can contribute to the monitoring of erosion rates of ice-rich cliffs in Arctic Siberia (Lena Delta, Russia). We produced two different vector products: i) Intra-annual cliff top retreat based on TerraSAR-X (TSX) satellite data (2012-2014): High-temporal resolution time series of TSX satellite data allow the inter-annual and intra-annual monitoring of the upper cliff-line retreat also under bad weather conditions and continuous cloud coverage. This published SAR product contains the retreating upper cliff lines of a 1.5 km long part of eroding ice-rich coast of Kurungnakh Island in the central Lena Delta. The upper cliff line was mapped using a thresholding approach for images acquired in the years 2012, 2013 and 2014 for the months June (2013, 2014), July (2013, 2014), August (2012, 2013, 2014) and September (2013, 2014). The cliff top retreat vector product is called 'upper_cliff_TerraSAR-X'. While the 2014 cliff lines show a clear retreat of 2 to 3 m/month, the cliff top lines for 2012 and 2013 are not chronologically ordered. However, lines from the end of the season of a year are always close to the lines from the beginning of the next summer season, indicating low cliff retreat in winter. ii) 4-year cliff top retreat based on optical satellite data (2010-2014): Long-term cliff top retreat could be assessed with two high-spatial resolution optical satellite images (GeoEye-1, 2010-08-05 and Worldview-1, 2014-08-19). The cliff top retreat vector product is called 'upper_cliff_optical'. Results: The long-term cliff top retreat derived from optical satellite data are 35 m cliff retreat within 4 years. The higher-temporal resolution SAR data equivalently show long-term rates of 18 m within 2 years and nearly now degradation activities in winter but maximum erosion rates in summer months.The Intra-seasonal cliff top retreat lines from 2014 show a rate of 2 to 3 m per month.

In-situ Iron isotope measurements and mineral compositions from ODP Hole 206-1256D and IODP Hole 335-U1256D

In-situ Fe isotope measurements have been carried out to estimate the impact of the hydrothermal metamorphic overprint on the Fe isotopic composition of Fe-Ti-oxides and Fe-sulfides of the different lithologies of the drilled rocks from IODP Hole 1256D (eastern equatorial Pacific; 15 Ma crust formed at the East Pacific Rise). Most igneous rocks normally have a very restricted range in their 56Fe/54Fe ratio. In contrast, Fe isotope compositions of hot fluids (> 300 °C) from mid-ocean-ridge spreading centers define a narrow range that is shifted to lower delta 56Fe values by 0.2 per mil - 0.5 per mil as compared to igneous rocks. Therefore, it is expected that mineral phases that contain large amounts of Fe are especially affected by the interaction with a fluid that fractionates Fe isotopes during exsolution/precipitation of those minerals. We have used a femtosecond UV-Laser ablation system to determine mineral 56Fe/54Fe ratios of selected samples with a precision of < 0.1 per mil (2 sigma level) at micrometer-scale. We have found significant variations of the delta 56Fe (IRMM-014) values in the minerals between different samples as well as within samples and mineral grains. The overall observed scale of delta 56Fe (magnetite) in 1256D rocks ranges from - 0.12 to + 0.64 per mil, and of delta 56Fe (ilmenite) from - 0.77 to + 0.01 per mil. Pyrite in the lowermost sheeted dike section is clearly distinguishable from the other investigated lithological units, having positive delta 56Fe values between + 0.29 and + 0.56 per mil, whereas pyrite in the other samples has generally negative delta 56Fe values from - 1.10 to - 0.59 permil. One key observation is that the temperature dependent inter-mineral fractionations of Fe isotopes between magnetite and ilmenite are systematically shifted towards higher values when compared to theoretically expected values, while synthesized, well equilibrated magnetite-ilmenite pairs are compatible with the theoretical predictions. Theoretical considerations including beta-factors of different aqueous Fe-chlorides and Rayleigh-type fractionations in the presence of a hydrous, chlorine-bearing fluid can explain this observation. The disagreement between observed and theoretical equilibrium fractionation, the fact that magnetite, in contrast to ilmenite shows a slight downhole trend in the delta 56Fe values, and the observation of small scale heterogeneities within single mineral grains imply that a general re-equilibration of the magnetite-ilmenite pairs is overprinted by kinetic fractionation effects, caused by the interaction of magnetite/ilmenite with hydrothermal fluids penetrating the upper oceanic crust during cooling, or incomplete re-equilibration at low temperatures. Furthermore, the observation of significant small-scale variations in the 56Fe/54Fe ratios of single minerals in this study highlights the importance of high spatial-resolution-analyses of stable isotope ratios for further investigations.

Physical properties of sediment cores from the Labrador Sea

Sediment core logs from six sediment cores in the Labrador Sea show millennial-scale climate variability during the last glacial by recording all Heinrich events and several major Dansgaard-Oeschger cycles. The same millennial-scale climate change is documented for surface-water d18O records of Neogloboquadrina pachyderma (left coiled); hence the surface-water d18O record can be derived from sediment core logging by means of multiple linear regression, providing a paleoclimate proxy record at very high temporal resolution (70 yrs). For the Labrador Sea, sediment core logs contain important information about deep-water current velocities and also reflect the variable input of IRD from different sources as inferred from grain-size analysis, benthic d18O, the relation of density and p-wave velocity, and magnetic susceptibility. For the last glacial, faster deep-water currents which correspond to highs in sediment physical properties, occurred during iceberg discharge and lasted for a several centuries to a few millennia. Those enhanced currents might have contributed to increased production of intermediate waters during times of reduced production of North Atlantic Deep Water. Hudson Strait might have acted as a major supplier of detrital carbonate only during lowered sea level (greater ice extent). During coldest atmospheric temperatures over Greenland, deep-water currents increased during iceberg discharge in the Labrador Sea, then surface water freshened shortly after, while the abrupt atmospheric temperature rise happened after a larger time lag of >=1 kyr. The correlation implies a strong link and common forcing for atmosphere, sea surface, and deep water during the last glacial at millennial time scales but decoupling at orbital time scales.

Osmium isotope record of Late Pleistocene sediments

We report high temporal resolution osmium isotopes records of bulk sediment and sediment leachates from DSDP Site 480 (Gulf of California) over the last 30 ka; from ODP Site 849 (Eastern equatorial Pacific) from the last 200 ka and from ODP Site 1002C (Cariaco Basin) across the 9-17 ka time interval in order to critically evaluate claims of a global 10% shift in the 187Os/188Os of seawater from glacial to interglacial intervals. We use organic-rich continental margin sites and carbonate-rich pelagic sites to isolate the temporal variations of the osmium seawater isotopic composition. Our results reveal that variations in 187Os/188Os fail to correlate with global changes in temperature across glacials/interglacials cycles as previously claimed. Instead, these results indicate differences of a few percent in the measured 187Os/188Os between each oceanic basin. We argue that these differences strongly suggest that seawater is not well homogenized with respect to its Os isotope composition.