Local changes in neocortical circuit dynamics coincide with the spread of seizures to thalamus in a model of epilepsy

During the generalization of epileptic seizures, pathological activity in one brain area recruits distant brain structures into joint synchronous discharges. However, it remains unknown whether specific changes in local circuit activity are related to the aberrant recruitment of anatomically distant structures into epileptiform discharges. Further, it is not known whether aberrant areas recruit or entrain healthy ones into pathological activity. Here we study the dynamics of local circuit activity during the spread of epileptiform discharges in the zero-magnesium in vitro model of epilepsy. We employ high-speed multi-photon imaging in combination with dual whole-cell recordings in acute thalamocortical (TC) slices of the juvenile mouse to characterize the generalization of epileptic activity between neocortex and thalamus. We find that, although both structures are exposed to zero-magnesium, the initial onset of focal epileptiform discharge occurs in cortex. This suggests that local recurrent connectivity that is particularly prevalent in cortex is important for the initiation of seizure activity. Subsequent recruitment of thalamus into joint, generalized discharges is coincident with an increase in the coherence of local cortical circuit activity that itself does not depend on thalamus. Finally, the intensity of population discharges is positively correlated between both brain areas. This suggests that during and after seizure generalization not only the timing but also the amplitude of epileptiform discharges in thalamus is entrained by cortex. Together these results suggest a central role of neocortical activity for the onset and the structure of pathological recruitment of thalamus into joint synchronous epileptiform discharges.

Computed Ultrasound Tomography in Echo mode (CUTE) of speed of sound for diagnosis and for aberration correction in pulse-echo sonography

Sound speed as a diagnostic marker for various diseases of human tissue has been of interest for a while. Up to now, mostly transmission ultrasound computed tomography (UCT) was able to detect spatially resolved sound speed, and its promise as a diagnostic tool has been demonstrated. However, UCT is limited to acoustically transparent samples such as the breast. We present a novel technique where spatially resolved detection of sound speed can be achieved using conventional pulse-echo equipment in reflection mode. For this purpose, pulse-echo images are acquired under various transmit beam directions and a two-dimensional map of the sound speed is reconstructed from the changing phase of local echoes using a direct reconstruction method. Phantom results demonstrate that a high spatial resolution (1 mm) and contrast (0.5 % of average sound speed) can be achieved suitable for diagnostic purposes. In comparison to previous reflection-mode based methods, CUTE works also in a situation with only diffuse echoes, and its direct reconstruction algorithm enables real-time application. This makes it suitable as an addition to conventional clinical ultrasound where it has the potential to benefit diagnosis in a multimodal approach. In addition, knowledge of the spatial distribution of sound speed allows full aberration correction and thus improved spatial resolution and contrast of conventional B-mode ultrasound. © (2014) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.

Reducing CT radiation dose with iterative reconstruction algorithms: the influence of scan and reconstruction parameters on image quality and CTDIvol

OBJECTIVES In this phantom CT study, we investigated whether images reconstructed using filtered back projection (FBP) and iterative reconstruction (IR) with reduced tube voltage and current have equivalent quality. We evaluated the effects of different acquisition and reconstruction parameter settings on image quality and radiation doses. Additionally, patient CT studies were evaluated to confirm our phantom results. METHODS Helical and axial 256 multi-slice computed tomography scans of the phantom (Catphan(®)) were performed with varying tube voltages (80-140kV) and currents (30-200mAs). 198 phantom data sets were reconstructed applying FBP and IR with increasing iterations, and soft and sharp kernels. Further, 25 chest and abdomen CT scans, performed with high and low exposure per patient, were reconstructed with IR and FBP. Two independent observers evaluated image quality and radiation doses of both phantom and patient scans. RESULTS In phantom scans, noise reduction was significantly improved using IR with increasing iterations, independent from tissue, scan-mode, tube-voltage, current, and kernel. IR did not affect high-contrast resolution. Low-contrast resolution was also not negatively affected, but improved in scans with doses <5mGy, although object detectability generally decreased with the lowering of exposure. At comparable image quality levels, CTDIvol was reduced by 26-50% using IR. In patients, applying IR vs. FBP resulted in good to excellent image quality, while tube voltage and current settings could be significantly decreased. CONCLUSIONS Our phantom experiments demonstrate that image quality levels of FBP reconstructions can also be achieved at lower tube voltages and tube currents when applying IR. Our findings could be confirmed in patients revealing the potential of IR to significantly reduce CT radiation doses.

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.

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.

Photometric monitoring of non-resolved space debris and databases of optical light curves.

The population of space debris increased drastically during the last years. Collisions involving massive objects may produce large number of fragments leading to significantly growth of the space debris population. An effective remediation measure in order to stabilize the population in LEO, is therefore the removal of large, massive space debris. To remove these objects, not only precise orbits, but also more detailed information about their attitude states will be required. One important property of an object targeted for removal is its spin period and spin axis orientation. If we observe a rotating object, the observer sees different surface areas of the object which leads to changes in the measured intensity. Rotating objects will produce periodic brightness vari ations with frequencies which are related to the spin periods. Photometric monitoring is the real tool for remote diagnostics of the satellite rotation around its center of mass. This information is also useful, for example, in case of contingency. Moreover, it is also important to take into account the orientation of non-spherical body (e.g. space debris) in the numerical integration of its motion when a close approach with the another spacecr aft is predicted. We introduce the two databases of light curves: the AIUB data base, which contains about a thousand light curves of LEO, MEO and high-altitude debris objects (including a few functional objects) obtained over more than seven years, and the data base of the Astronomical Observatory of Odessa University (Ukraine), which contains the results of more than 10 years of photometric monitoring of functioning satellites and large space debris objects in low Earth orbit. AIUB used its 1m ZIMLAT telescope for all light curves. For tracking low-orbit satellites, the Astronomical Observatory of Odessa used the KT-50 telescope, which has an alt-azimuth mount and allows tracking objects moving at a high angular velocity. The diameter of the KT-50 main mirror is 0.5 m, and the focal length is 3 m. The Odessa's Atlas of light curves includes almost 5,5 thousand light curves for ~500 correlated objects from a time period of 2005-2014. The processing of light curves and the determination of the rotation period in the inertial frame is challenging. Extracted frequencies and reconstructed phases for some interesting targets, e.g. GLONASS satellites, for which also SLR data were available for confirmation, will be presented. The rotation of the Envisat satellite after its sudden failure will be analyzed. The deceleration of its rotation rate within 3 years is studied together with the attempt to determine the orientation of the rotation axis.

GROMOS-C, a novel ground-based microwave radiometer for ozone measurement campaigns

Stratospheric ozone is of major interest as it absorbs most harmful UV radiation from the sun, allowing life on Earth. Ground-based microwave remote sensing is the only method that allows for the measurement of ozone profiles up to the mesopause, over 24 hours and under different weather conditions with high time resolution. In this paper a novel ground-based microwave radiometer is presented. It is called GROMOS-C (GRound based Ozone MOnitoring System for Campaigns), and it has been designed to measure the vertical profile of ozone distribution in the middle atmosphere by observing ozone emission spectra at a frequency of 110.836 GHz. The instrument is designed in a compact way which makes it transportable and suitable for outdoor use in campaigns, an advantageous feature that is lacking in present day ozone radiometers. It is operated through remote control. GROMOS-C is a total power radiometer which uses a pre-amplified heterodyne receiver, and a digital fast Fourier transform spectrometer for the spectral analysis. Among its main new features, the incorporation of different calibration loads stands out; this includes a noise diode and a new type of blackbody target specifically designed for this instrument, based on Peltier elements. The calibration scheme does not depend on the use of liquid nitrogen; therefore GROMOS-C can be operated at remote places with no maintenance requirements. In addition, the instrument can be switched in frequency to observe the CO line at 115 GHz. A description of the main characteristics of GROMOS-C is included in this paper, as well as the results of a first campaign at the High Altitude Research Station at Jungfraujoch (HFSJ), Switzerland. The validation is performed by comparison of the retrieved profiles against equivalent profiles from MLS (Microwave Limb Sounding) satellite data, ECMWF (European Centre for Medium-Range Weather Forecast) model data, as well as our nearby NDACC (Network for the Detection of Atmospheric Composition Change) ozone radiometer measuring at Bern.

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.

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.

Detection efficiency of microchannel plates for e(-) and pi(-) in the momentum range from 17.5 to 345 MeV/c

High-energy e(-) and pi(-) were measured by the multichannel plate (MCP) detector at the PiM1 beam line of the High Intensity Proton Accelerator Facilities located at the Paul Scherrer Institute, Villigen, Switzerland. The measurements provide the absolute detection efficiencies for these particles: 5.8% +/- 0.5% for electrons in the beam momenta range 17.5-300 MeV/c and 6.0% +/- 1.3% for pions in the beam momenta range 172-345 MeV/c. The pulse height distribution determined from the measurements is close to an exponential function with negative exponent, indicating that the particles penetrated the MCP material before producing the signal somewhere inside the channel. Low charge extraction and nominal gains of the MCP detector observed in this study are consistent with the proposed mechanism of the signal formation by penetrating radiation. A very similar MCP ion detector will be used in the Neutral Ion Mass (NIM) spectrometer designed for the JUICE mission of European Space Agency (ESA) to the Jupiter system, to perform measurements of the chemical composition of the Galilean moon exospheres. The detection efficiency for penetrating radiation determined in the present studies is important for the optimisation of the radiation shielding of the NIM detector against the high-rate and high-energy electrons trapped in Jupiter's magnetic field. Furthermore, the current studies indicate that MCP detectors can be useful to measure high-energy particle beams at high temporal resolution. (C) 2015 AIP Publishing LLC.

Mapping spectroscopic uncertainties into prospective methane retrieval errors from Sentinel-5 and its precursor

Sentinel-5 (S5) and its precursor (S5P) are future European satellite missions aiming at global monitoring of methane (CH4) column-average dry air mole fractions (XCH4). The spectrometers to be deployed onboard the satellites record spectra of sunlight backscattered from the Earth's surface and atmosphere. In particular, they exploit CH4 absorption in the shortwave infrared spectral range around 1.65 mu m (S5 only) and 2.35 mu m (both S5 and S5P) wavelength. Given an accuracy goal of better than 2% for XCH4 to be delivered on regional scales, assessment and reduction of potential sources of systematic error such as spectroscopic uncertainties is crucial. Here, we investigate how spectroscopic errors propagate into retrieval errors on the global scale. To this end, absorption spectra of a ground-based Fourier transform spectrometer (FTS) operating at very high spectral resolution serve as estimate for the quality of the spectroscopic parameters. Feeding the FTS fitting residuals as a perturbation into a global ensemble of simulated S5- and S5P-like spectra at relatively low spectral resolution, XCH4 retrieval errors exceed 0.6% in large parts of the world and show systematic correlations on regional scales, calling for improved spectroscopic parameters.

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.