Robust trends of landscape dynamics in the Arctic Lena Delta with temporally dense Landsat time-series stacks, with links to GeoTIFFs

Arctic permafrost landscapes are among the most vulnerable and dynamic landscapes globally, but due to their extent and remoteness most of the landscape changes remain unnoticed. In order to detect disturbances in these areas we developed an automated processing chain for the calculation and analysis of robust trends of key land surface indicators based on the full record of available Landsat TM, ETM +, and OLI data. The methodology was applied to the ~ 29,000 km**2 Lena Delta in Northeast Siberia, where robust trend parameters (slope, confidence intervals of the slope, and intercept) were calculated for Tasseled Cap Greenness, Wetness and Brightness, NDVI, and NDWI, and NDMI based on 204 Landsat scenes for the observation period between 1999 and 2014. The resulting datasets revealed regional greening trends within the Lena Delta with several localized hot-spots of change, particularly in the vicinity of the main river channels. With a 30-m spatial resolution various permafrost-thaw related processes and disturbances, such as thermokarst lake expansion and drainage, fluvial erosion, and coastal changes were detected within the Lena Delta region, many of which have not been noticed or described before. Such hotspots of permafrost change exhibit significantly different trend parameters compared to non-disturbed areas. The processed dataset, which is made freely available through the data archive PANGAEA, will be a useful resource for further process specific analysis by researchers and land managers. With the high level of automation and the use of the freely available Landsat archive data, the workflow is scalable and transferrable to other regions, which should enable the comparison of land surface changes in different permafrost affected regions and help to understand and quantify permafrost landscape dynamics.

Animation of micro-computed tomography of bioerosion trace Entobia megastoma made by East Australian bioeroding sponge Cliona celata

Endolithic bioerosion is difficult to analyse and to describe, and it usually requires damaging of the sample material. Sponge erosion (Entobia) may be one of the most difficult to evaluate as it is simultaneously macroscopically inhomogeneous and microstructurally intricate. We studied the bioerosion traces of the two Australian sponges Cliona celata Grant, 1826 (sensu Schönberg 2000) and Cliona orientalis Thiele, 1900 with a newly available radiographic technology: high resolution X-ray micro-computed tomography (MCT). MCT allows non-destructive visualisation of live and dead structures in three dimensions and was compared to traditional microscopic methods. MCT and microscopy showed that C. celata bioerosion was more intense in the centre and branched out in the periphery. In contrast, C. orientalis produced a dense, even trace meshwork and caused an overall more intense erosion pattern than C. celata. Extended pioneering filaments were not usually found at the margins of the studied sponge erosion, but branches ended abruptly or tapered to points. Results obtained with MCT were similar in quality to observations from transparent optical spar under the dissecting microscope. Microstructures could not be resolved as well as with e.g. scanning electron microscopy (SEM). Even though sponge scars and sponge chips were easily recognisable on maximum magnification MCT images, they lacked the detail that is available from SEM. Other drawbacks of MCT involve high costs and presently limited access. Even though MCT cannot presently replace traditional techniques such as corrosion casts viewed by SEM, we obtained valuable information. Especially for the possibility to measure endolithic pore volumes, we regard MCT as a very promising tool that will continue to be optimised. A combination of different methods will produce the best results in the study of Entobia.

Microfossils detected in a block of ancient dense clay coated with a ferromanganese crust, Clarion-Clipperton Province, East Equatorial Pacific

The area in study is characterized by a regional stratigraphic hiatus from Early Miocene to Quaternary. Deposits from Late Eocene to Early Miocene occur on the bottom surface or under a thin sedimentary cover. Ferromanganese nodules, mostly of Oligocene age, formed on surface layers of Tertiary or Quaternary sediments. A detailed micropaleontological study of a block of dense ancient clay coated with a ferromanganese crust was carried out. Composition of found radiolarian and diatomaceous complexes proved that the crust formed in Quaternary on an eroded surface of Late Oligocene clay. In Quaternary Neogene sediments were eroded and washed away by bottom currents. It is likely that the erosion began 0.9-0.7 Ma at the beginning of the "Glacial Pleistocene". The erosion could be initiated by loosening and resuspension of surface sediments resulting from seismic activity generated by strong earthquakes in the Central America subduction zone. The same vibration maintained residual nodules at the seafloor surface. Thus, for the area in study a common reason and a common Quaternary interval for formation of the following features is supposed: a regional stratigraphic hiatus, formation of residual nodule fields, and position of ancient nodules on the surface of Quaternary sediments.

Gas hydrate structures and C1-C5 hydrocarbons at individual sites in the Gulf of Mexico

Gas hydrate samples from various locations in the Gulf of Mexico (GOM) differ considerably in their microstructure. Distinct microstructure characteristics coincide with discrete crystallographic structures, gas compositions and calculated thermodynamic stabilities. The crystallographic structures were established by X-ray diffraction, using both conventional X-ray sources and high-energy synchrotron radiation. The microstructures were examined by cryo-stage Field-Emission Scanning Electron Microscopy (FE-SEM). Good sample preservation was warranted by the low ice fractions shown from quantitative phase analyses. Gas hydrate structure II samples from the Green Canyon in the northern GOM had methane concentrations of 70-80% and up to 30% of C2-C5 of measured hydrocarbons. Hydrocarbons in the crystallographic structure I hydrate from the Chapopote asphalt volcano in the southern GOM was comprised of more than 98% methane. Fairly different microstructures were identified for those different hydrates: Pores measuring 200-400 nm in diameter were present in structure I gas hydrate samples; no such pores but dense crystal surfaces instead were discovered in structure II gas hydrate. The stability of the hydrate samples is discussed regarding gas composition, crystallographic structure and microstructure. Electron microscopic observations showed evidence of gas hydrate and liquid oil co-occurrence on a micrometer scale. That demonstrates that oil has direct contact to gas hydrates when it diffuses through a hydrate matrix.

Particel size distribution and microstructures of sediments from ODP Leg 205 sites (Table 1)

Ocean Drilling Program Legs 170 and 205 offshore Costa Rica provide structural observations which support a new model for the geometry and deformation response to the seismic cycle of the frontal sedimentary prism and decollement. The model is based on drillcore, thin section, and electron microscope observations. The decollement damage zone is a few tens of meters in width, it develops mainly within the frontal prism. A clear cm-thick fault core is observed 1.6 km from the trench. The lower boundary of the fault core is coincident with the lithological boundary between the frontal prism and the hemipelagic and pelagic sediment of the Cocos plate. Breccia clast distributions in the upper portion of the decollement damage zone were studied through fractal analysis. This analysis shows that the fractal dimension changes with brecciated fragment size, implying that deformation was not accommodated by self-similar fracturing. A higher fractal dimensionality correlates with smaller particle size, which indicates that different or additional grain-size reduction processes operated during shearing. The co-existence of two distinct fracturing processes is also confirmed by microscopic analysis in which extension fracturing in the upper part of the damage zone farthest from the fault core is frequent, while both extension and shear fracturing occur approaching the fault core. The coexistence of extensional and shear fracturing seems to be best explained by fluid pressure variations in response to variations of the compressional regime during the seismic cycle. During the co-seismic event, sub-horizontal compression and fluid pressure increase, triggering shear fracturing and fluid expulsion. Fractures migrate upward with fluids, contributing to the asymmetric shape of the decollement, while slip propagates. In the inter-seismic interval the frontal prismrelaxes and fluid pressure drops. The frontal prismgoes into diffuse extension during the intervalwhen plate convergence is accommodated by creep along the ductile fault core. The fault core is typically a barrier to deformation, which is explained by its weak, but impermeable, nature. The localized development of a damage zone beneath the fault core is characterized by shear fracturing that appears as the result of local strengthening of the detachment.

Sediment microstructures of ODP Leg 105 holes

Possible genetic relationships between syn- and post-depositional processes and sediment microstructure were investigated. Samples from cores at Sites 646 and 647 of Ocean Drilling Program (ODP) Leg 105 included examples of bottom current deposition (contourites), turbidity current deposition, consolidation, and diagenesis. Examination of nearly 200 micrographs of 14 samples from Site 646 and 13 samples from Site 647 leads to the conclusion that sedimentation processes do not appear to have an obvious influence on fabric. The effects of post-depositional processes, such as bioturbation, coring disturbance, and even remolding, appear to be less significant than one might expect as a result of the relatively coarse grain size of the sediments studied. Consolidation resulting from increased overburden stress results in increased particle alignment and compression of fabric elements with depth. The transition from open, random fabric in shallow samples to preferred orientation at depth represents the only change in these sediments that can be ascribed directly to a specific depositional or post-depositional process. Mineralogical variations, owing to changes in weathering processes and growth of authigenic/diagenetic minerals, also have a pronounced effect on sediment fabric.

Dissolution of a spatially variable DNAPL (dense non aqueous phase liquid) source

Dissolution of non-aqueous phase liquids (NAPLs) or gases into groundwater is a key process, both for contamination problems originating from organic liquid sources, and for dissolution trapping in geological storage of CO2. Dissolution in natural systems typically will involve both high and low NAPL saturations and a wide range of pore water flow velocities within the same source zone for dissolution to groundwater. To correctly predict dissolution in such complex systems and as the NAPL saturations change over time, models must be capable of predicting dissolution under a range of saturations and flow conditions. To provide data to test and validate such models, an experiment was conducted in a two-dimensional sand tank, where the dissolution of a spatially variable, 5x5 cm**2 DNAPL tetrachloroethene source was carefully measured using x-ray attenuation techniques at a resolution of 0.2x0.2 cm**2. By continuously measuring the NAPL saturations, the temporal evolution of DNAPL mass loss by dissolution to groundwater could be measured at each pixel. Next, a general dissolution and solute transport code was written and several published rate-limited (RL) dissolution models and a local equilibrium (LE) approach were tested against the experimental data. It was found that none of the models could adequately predict the observed dissolution pattern, particularly in the zones of higher NAPL saturation. Combining these models with a model for NAPL pool dissolution produced qualitatively better agreement with experimental data, but the total matching error was not significantly improved. A sensitivity study of commonly used fitting parameters further showed that several combinations of these parameters could produce equally good fits to the experimental observations. The results indicate that common empirical model formulations for RL dissolution may be inadequate in complex, variable saturation NAPL source zones, and that further model developments and testing is desirable.