Vibrations, Posture, and the Stabilization of Gaze: An Experimental Study on Impedance Control

Vibrations, Posture, and the Stabilization of Gaze: An Experimental Study on Impedance Control R. KREDEL, A. GRIMM & E.-J. HOSSNER University of Bern, Switzerland Introduction Franklin and Wolpert (2011) identify impedance control, i.e., the competence to resist changes in position, velocity or acceleration caused by environmental disturbances, as one of five computational mechanisms which allow for skilled and fluent sen-sorimotor behavior. Accordingly, impedance control is of particular interest in situa-tions in which the motor task exhibits unpredictable components as it is the case in downhill biking or downhill skiing. In an experimental study, the question is asked whether impedance control, beyond its benefits for motor control, also helps to stabi-lize gaze what, in turn, may be essential for maintaining other control mechanisms (e.g., the internal modeling of future states) in an optimal range. Method In a 3x2x4 within-subject ANOVA design, 72 participants conducted three tests on visual acuity and contrast (Landolt / Grating and Vernier) in two different postures (standing vs. squat) on a platform vibrating at four different frequencies (ZEPTOR; 0 Hz, 4 Hz, 8 Hz, 12 Hz; no random noise; constant amplitude) in a counterbalanced or-der with 1-minute breaks in-between. In addition, perceived exertion (Borg) was rated by participants after each condition. Results For Landolt and Grating, significant main effects for posture and frequency are re-vealed, representing lower acuity/contrast thresholds for standing and for higher fre-quencies in general, as well as a significant interaction (p < .05), standing for in-creasing posture differences with increasing frequencies. Overall, performance could be maintained at the 0 Hz/standing level up to a frequency of 8 Hz, if bending of the knees was allowed. The fact that this result is not only due to exertion is proved by the Borg ratings showing significant main effects only, i.e., higher exertion scores for standing and for higher frequencies, but no significant interaction (p > .40). The same pattern, although not significant, is revealed for the Vernier test. Discussion Apparently, postures improving impedance control not only turn out to help to resist disturbances but also assist in stabilizing gaze in spite of these perturbations. Con-sequently, studying the interaction of these control mechanisms in complex unpre-dictable environments seems to be a fruitful field of research for the future. References Franklin, D. W., & Wolpert, D. M. (2011). Computational mechanisms of sensorimotor control. Neuron, 72, 425-442.

Investigations on an egyptian coffin and mummified head from Burgdorf Museum, Switzerland - A case of a historical puzzle

In the Burgdorf Museum of Ethnology, a mummy rests in a coffin. According to the inventory book, it was purchased from the Cairo Egyptian museum in 1926. The coffin was now examined by Egyptologists and the mummy was radiocarbon dated and examined by Anthropologists. The aim of the study was to compare the results and to check whether mummy and coffin actually belong together. The skull was examined morphological-anthropologically and by CT as a “blank sample”. Coffin and skull imply that the individual was female. The coffin dates to the Ptolemaic period. Only skull bones are preserved, the ethmoid is damaged. CT images Show resinous substances, bone fragments and brain remnants inside the skull. The ethmoid bone was probably foraminated during the mummification process and thus ended up inside the skull. The individual was mummified between the New Kingdom and the Ptolemaic period. Due to its style, it is most probable that the coffin comes from the Gamhud necropolis. The Burgdorf museum of ethnology inventory book chronicles were largely falsified by the examinations. There is a time gap between coffin and the mummy, there are two possible interpretations: the body was mummified with older linen, or the mummy and the coffin do not belong together. The authors strongly advise further investigations.

Bone mineral density (BMD) and computer tomographic measurements of the equine proximal phalanx in correlation with breaking strength

Despite the fact that bone mineral density (BMD) is an important fracture risk predictor in human medicine, studies in equine orthopedic research are still lacking. We hypothesized that BMD correlates with bone failure and fatigue fractures of this bone. Thus, the objectives of this study were to measure the structural and mechanical properties of the proximal phalanx with dual energy X-ray absorptiometry (DXA), to correlate the data obtained from DXA and computer tomography (CT) measurements to those obtained by loading pressure examination and to establish representative region of interest (ROI) for in vitro BMD measurements of the equine proximal phalanx for predicting bone failure force. DXA was used to measure the whole bone BMD and additional three ROI sites in 14 equine proximal phalanges. Following evaluation of the bone density, whole bone, cortical width and area in the mid-diaphyseal plane were measured on CT images. Bones were broken using a manually controlled universal bone crusher to measure bone failure force and reevaluated for the site of fractures on follow-up CT images. Compressive load was applied at a constant displacement rate of 2 mm/min until failure, defined as the first clear drop in the load measurement. The lowest BMD was measured at the trabecular region (mean +/- SD: 1.52 +/- 0.12 g/cm2; median: 1.48 g/cm2; range: 1.38-1.83 g/cm2). There was a significant positive linear correlation between trabelcular BMD and the breaking strength (P = 0.023, r = 0.62). The trabecular region of the proximal phalanx appears to be the only significant indicator of failure of strength in vitro. This finding should be reassessed to further reveal the prognostic value of trabecular BMD in an in vivo fracture risk model.

Thoracic target volume delineation using various maximum-intensity projection computed tomography image sets for stereotactic body radiation therapy

The motion of lung tumors during respiration makes the accurate delivery of radiation therapy to the thorax difficult because it increases the uncertainty of target position. The adoption of four-dimensional computed tomography (4D-CT) has allowed us to determine how a tumor moves with respiration for each individual patient. Using information acquired during a 4D-CT scan, we can define the target, visualize motion, and calculate dose during the planning phase of the radiotherapy process. One image data set that can be created from the 4D-CT acquisition is the maximum-intensity projection (MIP). The MIP can be used as a starting point to define the volume that encompasses the motion envelope of the moving gross target volume (GTV). Because of the close relationship that exists between the MIP and the final target volume, we investigated four MIP data sets created with different methodologies (3 using various 4D-CT sorting implementations, and one using all available cine CT images) to compare target delineation. It has been observed that changing the 4D-CT sorting method will lead to the selection of a different collection of images; however, the clinical implications of changing the constituent images on the resultant MIP data set are not clear. There has not been a comprehensive study that compares target delineation based on different 4D-CT sorting methodologies in a patient population. We selected a collection of patients who had previously undergone thoracic 4D-CT scans at our institution, and who had lung tumors that moved at least 1 cm. We then generated the four MIP data sets and automatically contoured the target volumes. In doing so, we identified cases in which the MIP generated from a 4D-CT sorting process under-represented the motion envelope of the target volume by more than 10% than when measured on the MIP generated from all of the cine CT images. The 4D-CT methods suffered from duplicate image selection and might not choose maximum extent images. Based on our results, we suggest utilization of a MIP generated from the full cine CT data set to ensure a representative inclusive tumor extent, and to avoid geometric miss.

DUAL ENERGY COMPUTED TOMOGRAPHY FOR PROTON THERAPY TREATMENT PLANNING

Proton radiation therapy is gaining popularity because of the unique characteristics of its dose distribution, e.g., high dose-gradient at the distal end of the percentage-depth-dose curve (known as the Bragg peak). The high dose-gradient offers the possibility of delivering high dose to the target while still sparing critical organs distal to the target. However, the high dose-gradient is a double-edged sword: a small shift of the highly conformal high-dose area can cause the target to be substantially under-dosed or the critical organs to be substantially over-dosed. Because of that, large margins are required in treatment planning to ensure adequate dose coverage of the target, which prevents us from realizing the full potential of proton beams. Therefore, it is critical to reduce uncertainties in the proton radiation therapy. One major uncertainty in a proton treatment is the range uncertainty related to the estimation of proton stopping power ratio (SPR) distribution inside a patient. The SPR distribution inside a patient is required to account for tissue heterogeneities when calculating dose distribution inside the patient. In current clinical practice, the SPR distribution inside a patient is estimated from the patient’s treatment planning computed tomography (CT) images based on the CT number-to-SPR calibration curve. The SPR derived from a single CT number carries large uncertainties in the presence of human tissue composition variations, which is the major drawback of the current SPR estimation method. We propose to solve this problem by using dual energy CT (DECT) and hypothesize that the range uncertainty can be reduced by a factor of two from currently used value of 3.5%. A MATLAB program was developed to calculate the electron density ratio (EDR) and effective atomic number (EAN) from two CT measurements of the same object. An empirical relationship was discovered between mean excitation energies and EANs existing in human body tissues. With the MATLAB program and the empirical relationship, a DECT-based method was successfully developed to derive SPRs for human body tissues (the DECT method). The DECT method is more robust against the uncertainties in human tissues compositions than the current single-CT-based method, because the DECT method incorporated both density and elemental composition information in the SPR estimation. Furthermore, we studied practical limitations of the DECT method. We found that the accuracy of the DECT method using conventional kV-kV x-ray pair is susceptible to CT number variations, which compromises the theoretical advantage of the DECT method. Our solution to this problem is to use a different x-ray pair for the DECT. The accuracy of the DECT method using different combinations of x-ray energies, i.e., the kV-kV, kV-MV and MV-MV pair, was compared using the measured imaging uncertainties for each case. The kV-MV DECT was found to be the most robust against CT number variations. In addition, we studied how uncertainties propagate through the DECT calculation, and found general principles of selecting x-ray pairs for the DECT method to minimize its sensitivity to CT number variations. The uncertainties in SPRs estimated using the kV-MV DECT were analyzed further and compared to those using the stoichiometric method. The uncertainties in SPR estimation can be divided into five categories according to their origins: the inherent uncertainty, the DECT modeling uncertainty, the CT imaging uncertainty, the uncertainty in the mean excitation energy, and SPR variation with proton energy. Additionally, human body tissues were divided into three tissue groups – low density (lung) tissues, soft tissues and bone tissues. The uncertainties were estimated separately because their uncertainties were different under each condition. An estimate of the composite range uncertainty (2s) was determined for three tumor sites – prostate, lung, and head-and-neck, by combining the uncertainty estimates of all three tissue groups, weighted by their proportions along typical beam path for each treatment site. In conclusion, the DECT method holds theoretical advantages in estimating SPRs for human tissues over the current single-CT-based method. Using existing imaging techniques, the kV-MV DECT approach was capable of reducing the range uncertainty from the currently used value of 3.5% to 1.9%-2.3%, but it is short to reach our original goal of reducing the range uncertainty by a factor of two. The dominant source of uncertainties in the kV-MV DECT was the uncertainties in CT imaging, especially in MV CT imaging. Further reduction in beam hardening effect, the impact of scatter, out-of-field object etc. would reduce the Hounsfeld Unit variations in CT imaging. The kV-MV DECT still has the potential to reduce the range uncertainty further.

Scaling relations between trabecular bone volume fraction and microstructure at different skeletal sites

In this study, we investigated the scaling relations between trabecular bone volume fraction (BV/TV) and parameters of the trabecular microstructure at different skeletal sites. Cylindrical bone samples with a diameter of 8mm were harvested from different skeletal sites of 154 human donors in vitro: 87 from the distal radius, 59/69 from the thoracic/lumbar spine, 51 from the femoral neck, and 83 from the greater trochanter. μCT images were obtained with an isotropic spatial resolution of 26μm. BV/TV and trabecular microstructure parameters (TbN, TbTh, TbSp, scaling indices (< > and σ of α and αz), and Minkowski Functionals (Surface, Curvature, Euler)) were computed for each sample. The regression coefficient β was determined for each skeletal site as the slope of a linear fit in the double-logarithmic representations of the correlations of BV/TV versus the respective microstructure parameter. Statistically significant correlation coefficients ranging from r=0.36 to r=0.97 were observed for BV/TV versus microstructure parameters, except for Curvature and Euler. The regression coefficients β were 0.19 to 0.23 (TbN), 0.21 to 0.30 (TbTh), −0.28 to −0.24 (TbSp), 0.58 to 0.71 (Surface) and 0.12 to 0.16 (<α>), 0.07 to 0.11 (<αz>), −0.44 to −0.30 (σ(α)), and −0.39 to −0.14 (σ(αz)) at the different skeletal sites. The 95% confidence intervals of β overlapped for almost all microstructure parameters at the different skeletal sites. The scaling relations were independent of vertebral fracture status and similar for subjects aged 60–69, 70–79, and >79years. In conclusion, the bone volume fraction–microstructure scaling relations showed a rather universal character.

Monoenergetic computed tomography reconstructions reduce beam hardening artifacts from dental restorations

The aim of this study was to assess the potential of monoenergetic computed tomography (CT) images to reduce beam hardening artifacts in comparison to standard CT images of dental restoration on dental post-mortem CT (PMCT). Thirty human decedents (15 male, 58 ± 22 years) with dental restorations were examined using standard single-energy CT (SECT) and dual-energy CT (DECT). DECT data were used to generate monoenergetic CT images, reflecting the X-ray attenuation at energy levels of 64, 69, 88 keV, and at an individually adjusted optimal energy level called OPTkeV. Artifact reduction and image quality of SECT and monoenergetic CT were assessed objectively and subjectively by two blinded readers. Subjectively, beam artifacts decreased visibly in 28/30 cases after monoenergetic CT reconstruction. Inter- and intra-reader agreement was good (k = 0.72, and k = 0.73 respectively). Beam hardening artifacts decreased significantly with increasing monoenergies (repeated-measures ANOVA p < 0.001). Artifact reduction was greatest on monoenergetic CT images at OPTkeV. Mean OPTkeV was 108 ± 17 keV. OPTkeV yielded the lowest difference between CT numbers of streak artifacts and reference tissues (-163 HU). Monoenergetic CT reconstructions significantly reduce beam hardening artifacts from dental restorations and improve image quality of post-mortem dental CT.

Retrocrural space involvement on computed tomography as a predictor of mortality and disease severity in acute pancreatitis

BACKGROUND Because computed tomography (CT) has advantages for visualizing the manifestation of necrosis and local complications, a series of scoring systems based on CT manifestations have been developed for assessing the clinical outcomes of acute pancreatitis (AP), including the CT severity index (CTSI), modified CTSI, etc. Despite the internationally accepted CTSI having been successfully used to predict the overall mortality and disease severity of AP, recent literature has revealed the limitations of the CTSI. Using the Delphi method, we establish a new scoring system based on retrocrural space involvement (RCSI), and compared its effectiveness at evaluating the mortality and severity of AP with that of the CTSI. METHODS We reviewed CT images of 257 patients with AP taken within 3-5 days of admission in 2012. The RCSI scoring system, which includes assessment of infectious conditions involving the retrocrural space and the adjacent pleural cavity, was established using the Delphi method. Two radiologists independently assessed the RCSI and CTSI scores. The predictive points of the RCSI and CTSI scoring systems in evaluating the mortality and severity of AP were estimated using receiver operating characteristic (ROC) curves. PRINCIPAL FINDINGS The RCSI score can accurately predict the mortality and disease severity. The area under the ROC curve for the RCSI versus CTSI score was 0.962±0.011 versus 0.900±0.021 for predicting the mortality, and 0.888±0.025 versus 0.904±0.020 for predicting the severity of AP. Applying ROC analysis to our data showed that a RCSI score of 4 was the best cutoff value, above which mortality could be identified. CONCLUSION The Delphi method was innovatively adopted to establish a scoring system to predict the clinical outcome of AP. The RCSI scoring system can predict the mortality of AP better than the CTSI system, and the severity of AP equally as well.

Facial Nerve Image Enhancement from CBCT using Supervised Learning Technique

Facial nerve segmentation plays an important role in surgical planning of cochlear implantation. Clinically available CBCT images are used for surgical planning. However, its relatively low resolution renders the identification of the facial nerve difficult. In this work, we present a supervised learning approach to enhance facial nerve image information from CBCT. A supervised learning approach based on multi-output random forest was employed to learn the mapping between CBCT and micro-CT images. Evaluation was performed qualitatively and quantitatively by using the predicted image as input for a previously published dedicated facial nerve segmentation, and cochlear implantation surgical planning software, OtoPlan. Results show the potential of the proposed approach to improve facial nerve image quality as imaged by CBCT and to leverage its segmentation using OtoPlan.

Body Height Estimation from Femur Measurements in Postmortem Computed Tomography

After attending this presentation, attendees will: (1) understand how body height from computed tomography data can be estimated; and, (2) gain knowledge about the accuracy of estimated body height and limitations. The presentation will impact the forensic science community by providing knowledge and competence which will enable attendees to develop formulas for single bones to reconstruct body height using postmortem Computer Tomography (p-CT) data. The estimation of Body Height (BH) is an important component of the identification of corpses and skeletal remains. Stature can be estimated with relative accuracy via the measurement of long bones, such as the femora. Compared to time-consuming maceration procedures, p-CT allows fast and simple measurements of bones. This study undertook four objectives concerning the accuracy of BH estimation via p-CT: (1) accuracy between measurements on native bone and p-CT imaged bone (F1 according to Martin 1914); (2) intra-observer p-CT measurement precision; (3) accuracy between formula-based estimation of the BH and conventional body length measurement during autopsy; and, (4) accuracy of different estimation formulas available.1 In the first step, the accuracy of measurements in the CT compared to those obtained using an osteometric board was evaluated on the basis of eight defleshed femora. Then the femora of 83 female and 144 male corpses of a Swiss population for which p-CTs had been performed, were measured at the Institute of Forensic Medicine in Bern. After two months, 20 individuals were measured again in order to assess the intraobserver error. The mean age of the men was 53±17 years and that of the women was 61±20 years. Additionally, the body length of the corpses was measured conventionally. The mean body length was 176.6±7.2cm for men and 163.6±7.8cm for women. The images that were obtained using a six-slice CT were reconstructed with a slice thickness of 1.25mm. Analysis and measurements of CT images were performed on a multipurpose workstation. As a forensic standard procedure, stature was estimated by means of the regression equations by Penning & Riepert developed on a Southern German population and for comparison, also those referenced by Trotter & Gleser “American White.”2,3 All statistical tests were performed with a statistical software. No significant differences were found between the CT and osteometric board measurements. The double p-CT measurement of 20 individuals resulted in an absolute intra-observer difference of 0.4±0.3mm. For both sexes, the correlation between the body length and the estimated BH using the F1 measurements was highly significant. The correlation coefficient was slightly higher for women. The differences in accuracy of the different formulas were small. While the errors of BH estimation were generally ±4.5–5.0cm, the consideration of age led to an increase in accuracy of a few millimetres to about 1cm. BH estimations according to Penning & Riepert and Trotter & Gleser were slightly more accurate when age-at-death was taken into account.2,3 That way, stature estimations in the group of individuals older than 60 years were improved by about 2.4cm and 3.1cm.2,3 The error of estimation is therefore about a third of the common ±4.7cm error range. Femur measurements in p-CT allow very accurate BH estimations. Estimations according to Penning led to good results that (barely) come closer to the true value than the frequently used formulas by Trotter & Gleser “American White.”2,3 Therefore, the formulas by Penning & Riepert are also validated for this substantial recent Swiss population.

Interobserver variations of target volume delineation in multicatheter partial breast brachytherapy after open cavity surgery

PURPOSE To investigate interobserver variations of target volume delineations in accelerated partial breast irradiation with multicatheter brachytherapy (BT) and to assess the impact of guidelines on consistency of contouring. METHODS AND MATERIALS A contouring study with two phases in interstitial accelerated partial breast irradiation after open cavity surgery was conducted by the Groupe Européen de Curiethérapie-European Society for Radiotherapy and Oncology Breast Cancer Working Group. Contours of cavity and planning target volume (PTV) on preimplant and postimplant CT images were delineated. In Phase 1, nine radiation oncologists defined the target volumes of 5 patients, whereas in Phase 2, four observers draw the contours of 4 patients applying guidelines. In Phase 1, experience in breast BT after open cavity surgery was assessed. The delineations were compared between Phase 1 and Phase 2, the impact of guidelines was assessed, and cavity visualization score was related to consistency of delineations. RESULTS Significant interobserver variability in delineations of lumpectomy cavity and PTV was observed among the participants. Observers with BT experience after open cavity surgery outlined the cavity and PTV more consistently (conformity indexgen: 0.52 vs. 0.48 and 0.59 vs. 0.55 for preimplant and postimplant cavities). For all volumes, the mean Vmax/Vmin was 2.2 vs. 2.8. Having used guidelines all conformity indices increased significantly. For cavity, the increase was 14% and 11%, whereas for the PTV, 28% and 17% on the preimplant and postimplant CT images, respectively. A strong correlation was found between consistency of contours and cavity visualization score. CONCLUSIONS Simple guidelines on defining the lumpectomy cavity significantly increased the consistency of contouring. Reliable consistency of target volume definition can be expected only for good cavity visibility.

Computer-assisted virtual autopsy using surgical navigation techniques

UNLABELLED OBJECTIVE; Virtual autopsy methods, such as postmortem CT and MRI, are increasingly being used in forensic medicine. Forensic investigators with little to no training in diagnostic radiology and medical laypeople such as state's attorneys often find it difficult to understand the anatomic orientation of axial postmortem CT images. We present a computer-assisted system that permits postmortem CT datasets to be quickly and intuitively resliced in real time at the body to narrow the gap between radiologic imaging and autopsy. CONCLUSION Our system is a potentially valuable tool for planning autopsies, showing findings to medical laypeople, and teaching CT anatomy, thus further closing the gap between radiology and forensic pathology.

Assessment of cone beam computed tomography techniques for imaging lung damage in mice in vivo

Lung damage is a common side effect of chemotherapeutic drugs such as bleomycin. This study used a bleomycin mouse model which simulates the lung damage observed in humans. Noninvasive, in vivo cone-beam computed tomography (CBCT) was used to visualize and quantify fibrotic and inflammatory damage over the entire lung volume of mice. Bleomycin was used to induce pulmonary damage in vivo and the results from two CBCT systems, a micro-CT and flat panel CT (fpCT), were compared to histologic measurements, the standard method of murine lung damage quantification. Twenty C57BL/6 mice were given either 3 U/kg of bleomycin or saline intratracheally. The mice were scanned at baseline, before the administration of bleomycin, and then 10, 14, and 21 days afterward. At each time point, a subset of mice was sacrificed for histologic analysis. The resulting CT images were used to assess lung volume. Percent lung damage (PLD) was calculated for each mouse on both the fpCT (PLDfpcT) and the micro-CT (PLDμCT). Histologic PLD (PLDH) was calculated for each histologic section at each time point (day 10, n = 4; day 14, n = 4; day 21, n = 5; control group, n = 5). A linear regression was applied to the PLDfpCT vs. PLDH, PLDμCT vs. PLDH and PLDfpCT vs. PLDμCT distributions. This study did not demonstrate strong correlations between PLDCT and PLDH. The coefficient of determination, R, was 0.68 for PLDμCT vs. PLDH and 0.75 for the PLD fpCT vs. PLDH. The experimental issues identified from this study were: (1) inconsistent inflation of the lungs from scan to scan, (2) variable distribution of damage (one histologic section not representative of overall lung damage), (3) control mice not scanned with each group of bleomycin mice, (4) two CT systems caused long anesthesia time for the mice, and (5) respiratory gating did not hold the volume of lung constant throughout the scan. Addressing these issues might allow for further improvement of the correlation between PLDCT and PLDH. ^

Determinants of CD4 cell count decline after treatment interruption

Patients who had started HAART (Highly Active Anti-Retroviral Treatment) under previous aggressive DHHS guidelines (1997) underwent a life-long continuous HAART that was associated with many short term as well as long term complications. Many interventions attempted to reduce those complications including intermittent treatment also called pulse therapy. Many studies were done to study the determinants of rate of fall in CD4 count after interruption as this data would help guide treatment interruptions. The data set used here was a part of a cohort study taking place at the Johns Hopkins AIDS service since January 1984, in which the data were collected both prospectively and retrospectively. The patients in this data set consisted of 47 patients receiving via pulse therapy with the aim of reducing the long-term complications. ^ The aim of this project was to study the impact of virologic and immunologic factors on the rate of CD4 loss after treatment interruption. The exposure variables under investigation included CD4 cell count and viral load at treatment initiation. The rates of change of CD4 cell count after treatment interruption was estimated from observed data using advanced longitudinal data analysis methods (i.e., linear mixed model). Using random effects accounted for repeated measures of CD4 per person after treatment interruption. The regression coefficient estimates from the model was then used to produce subject specific rates of CD4 change accounting for group trends in change. The exposure variables of interest were age, race, and gender, CD4 cell counts and HIV RNA levels at HAART initiation. ^ The rate of fall of CD4 count did not depend on CD4 cell count or viral load at initiation of treatment. Thus these factors may not be used to determine who can have a chance of successful treatment interruption. CD4 and viral load were again studied by t-tests and ANOVA test after grouping based on medians and quartiles to see any difference in means of rate of CD4 fall after interruption. There was no significant difference between the groups suggesting that there was no association between rate of fall of CD4 after treatment interruption and above mentioned exposure variables. ^

Mineralogical and stable isotope analyses of carbonates from the eastern Black Sea

Authigenic carbonates associated with cold seeps provide valuable archives of changes in the long-term seepage activity. To investigate the role of shallow-buried hydrates on the seepage strength and fluid composition we analysed methane-derived carbonate precipitates from a high-flux hydrocarbon seepage area ("Batumi seep area") located on the south-eastern Black Sea slope in ca. 850 m. In a novel approach, we combined computerized X-ray tomography (CT) with mineralogical and isotope geochemical methods to get additional insights into the three-dimensional internal structure of the carbonate build-ups. X-ray diffractometry revealed the presence of two different authigenic carbonate phases, i.e. pure aragonitic rims associated with vital microbial mats and high-Mg calcite cementing the hemipelagic sediment. As indicated by the CT images, the initial sediment has been strongly deformed, first plastic then brittle, leading to brecciation of the progressively cemented sediment. The aragonitic rims on the other hand, represent a presumably recent carbonate growth phase since they cover the already deformed sediment. The stable oxygen isotope signature indicates that the high-Mg calcite cement incorporated pore water mixed with substantial hydrate water amounts. This points at a dominant role of high gas/fluid flux from decomposing gas hydrates leading to the deformation and cementation of the overlying sediment. In contrast, the aragonitic rims do not show an influence of 18O-enriched hydrate water. The differences in d18O between the presumably recent aragonite precipitates and the older high-Mg cements suggest that periods of hydrate dissociation and vigorous fluid discharge alternated with times of hydrate stability and moderate fluid flow. These results indicate that shallow-buried gas hydrates are prone to episodic decomposition with associated vigorous fluid flow. This might have a profound impact on the seafloor morphology resulting e.g. in the formation of carbonate pavements and pockmark-like structures but might also affect the local carbon cycle.