Elastic thickness structure of the Andaman subduction zone: Implications for convergence of the Ninetyeast Ridge

We use the Bouguer coherence (Morlet isostatic response function) technique to compute the spatial variation of effective elastic thickness (T-e) of the Andaman subduction zone. The recovered T-e map resolves regional-scale features that correlate well with known surface structures of the subducting Indian plate and the overriding Burma plate. The major structure on the India plate, the Ninetyeast Ridge (NER), exhibits a weak mechanical strength, which is consistent with the expected signature of an oceanic ridge of hotspot origin. However, a markedly low strength (0< T-e <3 km) in that region, where the NER is close to the Andaman trench (north of 10 N), receives our main attention in this study. The subduction geometry derived from the Bouguer gravity forward modeling suggests that the NER has indented beneath the Andaman arc. We infer that the bending stresses of the viscous plate, which were reinforced within the subducting oceanic plate as a result of the partial subduction of the NER buoyant load, have reduced the lithospheric strength. The correlation, T-e < T-s (seismogenic thickness) reveals that the upper crust is actively deforming beneath the frontal arc Andaman region. The occurrence of normal-fault earthquakes in the frontal arc, low Te zone, is indicative of structural heterogeneities within the subducting plate. The fact that the NER along with its buoyant root is subducting under the Andaman region is inhibiting the subduction processes, as suggested by the changes in trench line, interrupted back-arc volcanism, variation in seismicity mechanism, slow subduction, etc. The low T-e and thinned crustal structure of the Andaman back-arc basin are attributed to a thermomechanically weakened lithosphere. The present study reveals that the ongoing back-arc spreading and strike-slip motion along the West Andaman Fault coupled with the ridge subduction exerts an important control on the frequency and magnitude of seismicity in the Andaman region. (C) 2013 Elsevier Ltd. All rights reserved.

Micromorphology of Last Glacial Maximum Grounding line sediments around the Antarctic Peninsula

Sediments recovered from seven Last Glacial Maximum grounding lines sites, around the Antarctic Peninsula, were analyzed using micromorphology. This is the first evidence that grounding line sediments from around the Antarctic Peninsula have complex deformational histories and subglacial origins. It was determined that grounding zone wedge contain multiple units, or diamicton layers, with homogenized boundaries. The multiple diamicton units / layers are due to the accretionary formation of a grounding line wedge. All the sediments were deposited via deformation, and continual reincorporation, homogenization of lower diamicton layers by upper diamicton layers produced what macroscopically appeared to be a single massive diamicton unit. The morainal ridge that was sampled, alternatively, is composed of a single unit, or diamicton layer, that was subglacial in origin and believed to have been pushed out to form a ridge that was subsequently deformed via glacial push.

Minimum spanning trees for valley and ridge characterization in digital elevation maps

Texture synthesis employs neighbourhood matching to generate appropriate new content. Terrain synthesis has the added constraint that new content must be geographically plausible. The profile recognition and polygon breaking algorithm (PPA) [Chang et al. 1998] provides a robust mechanism for characterizing terrain as systems of valley and ridge lines in digital elevation maps. We exploit this to create a terrain characterization metric that is robust, efficient to compute and is sensitive to terrain properties.

Terrain regions are characterized as a minimum spanning tree derived from a graph created from the sample points of the elevation map which are encoded as weights in the edges of the graph. This formulation allows us to provide a single consistent feature definition that is sensitive to the pattern of ridges and valleys in the terrain Alternative formulations of these weights provide richer characteristicmeasures and we provide examples of alternate definitions based on curvature and contour measures.

We show that the measure is robust, with a significant portion derived directly from information local to the terrain sample. Global terrain characteristics introduce the issue of over- and underconnected valley/ridge lines when working with sub-regions. This is addressed by providing two graph construction strategies, which respectively provide an upper bound on connectivity as a single spanning tree, and a lower bound as a forest of trees.

Efficient minimum spanning tree algorithms are adapted to the context of terrain data and are shown to provide substantially better performance than previous PPA implementations. In particular, these are able to characterize valley and ridge behaviour at every point even in large elevation maps, providing a measure sensitive to terrain features at all scales.

The resulting graph based formulation provides an efficient and elegant algorithm for characterizing terrain features. The measure can be calculated efficiently, is robust under changes of neighbourhood position, size and resolution and the hybrid measure is sensitive to terrain features both locally and globally.

Mineral composition and major element oxides of gabbros from ODP Hole 179-1105A, southwest Indian Ridge

Hole 1105A penetrated 158 m of gabbros at a site offset 1.3 km east-northeast from Hole 735B on the Atlantis Bank near the Atlantis II Fracture Zone. A total of 118 m of dominantly medium- to coarse-grained intercalated Fe-Ti oxide gabbro and olivine gabbro was recovered from Hole 1105A that shows many petrographic features similar to those recovered from the upper part of Hole 735B. The main rock types are distinguished based on the constituent cumulus phases, with the most primitive gabbros consisting of olivine, plagioclase, and clinopyroxene. The inferred crystallization order is subsequently Fe-Ti oxides (ilmenite and titanomagnetite), followed by orthopyroxene, then apatite, and finally biotite. Orthopyroxene appears to replace olivine in a narrow middle interval. The magmatic evolution is likewise reflected in the mineral compositions. Plagioclase varies from An66 to An28. Olivine varies from Fo78 to Fo35. The gap in olivine crystallization occurs between Fo46 and Fo40 and coincides approximately with the appearance of orthopyroxene (~En50). The clinopyroxenes show large compositional variation in Mg/(Mg + Fe total) from 0.84 to 0.51. The nonquadrilateral cations of clinopyroxene similarly show large variations with Ti increasing for the olivine gabbros and decreasing for the Fe-Ti oxide gabbros with the decrease in Mg/(Mg + Fe total). The apatites are mainly flourapatites. The compositional variation in the gabbros is interpreted as a comagmatic suite resulting from fractional crystallization. Pyroxene geothermometry suggests equilibration temperatures from 1100°C and below. The coexisting Fe-Ti oxide minerals indicate subsolidus equilibration temperatures from 900°C for olivine gabbros to 700°C for the most evolved apatite-bearing gabbros. The cryptic variation in the olivine gabbros defines two or three lenses, 40 to 60 m thick, each characterized by a distinct convex zoning with a lower segment indicating upward reverse fractionation, a central maximum, and an upper segment showing normal fractionation. The Fe-Ti oxide gabbros show cryptic variations independent of the host olivine gabbros and reveal a systematic upward normal fractionation trend transgressing host olivine gabbro boundaries. Forward fractional crystallization modeling, using a likely parental magma composition from the Atlantis II Fracture Zone (MgO = 7.2 wt%; Mg/[Mg + Fe2+] = 0.62), closely matches the compositions of coexisting olivine, plagioclase, and clinopyroxene. This modeling suggests cosaturation of olivine, plagioclase, and clinopyroxene from 1155°C and the addition of Fe-Ti oxides from 1100°C. The liquid line of descent initially shows increasing FeO with moderately increasing SiO2. After saturation of Fe-Ti oxides, the liquid strongly decreases in FeO and TiO2 and increases in SiO2, reaching dacitic compositions at ~10% liquid remaining. The calculations indicate that formation of olivine gabbros can be accounted for by <65% fractionation and that only the residual 35% liquid was saturated in Fe-Ti oxides. The modeling of the solid fractionation products shows that both the olivine gabbro and the Fe-Ti oxide gabbros contain very small amounts of trapped liquid (<5%). The implications are that the gabbros represent crystal mush that originated in a recharging and tapping subaxial chamber. Compaction and upward melt migration in the crystal mush appear to have been terminated with relatively large amounts of interstitial liquid remaining in the upper parts of the cumulate mush. This termination may have been caused by tectonic disturbances, uplift, and associated withdrawal of magma into the subaxial dike and sill system. Prolonged compaction and cooling of the trapped melt in the mush formed small differentiated bodies and lenses by pressure release migration and crystallization along syntectonic channels. This resulted in differentiation products along lateral and vertical channelways in the host gabbro that vary from olivine gabbro, to Fe-Ti oxide gabbro, gabbronorite, and apatite gabbros and show large compositional variations independent of the host olivine gabbros.

Reconstructed ice-sheet limit in the Kara and Barents Sea area at about 50,000 years ago (Markhida Line)

Glacial landforms in northern Russia, from the Timan Ridge in the west to the east of the Urals, have been mapped by aerial photographs and satellite images supported by field observations. An east-west trending belt of fresh hummock-and-lake glaciokarst landscapes has been traced to the north of 67°N. The southern boundary of these landscapes is called the Markhida Line, which is interpreted as a nearly synchronous limit of the last ice sheet that affected this region. The hummocky landscapes are subdivided into three types according to the stage of postglacial modification: Markhida, Harbei and Halmer. The Halmer landscape on the Uralian piedmont in the east is the freshest, whereas the westernmost Markhida landscape is more eroded. The west- east gradient in morphology is considered to be a result of the time-transgressive melting of stagnant glacier ice and of the underlying permafrost. The pattern of ice-pushed ridges and other directional features reflects a dominant ice flow direction from the Kara Sea shelf. Traces of ice movement from the central Barents Sea are only discernible in the Pechora River left bank area west of 50°E. In the Polar Urals the horseshoe-shaped end moraines at altitudes of up to 560 m a.s.l. reflect ice movement up-valley from the Kara Ice Sheet, indicating the absence of a contemporaneous ice dome in the mountains. The Markhida moraines, superimposed onto the Eemian strata, represent the maximum ice sheet extent in the western part of the Pechora Basin during the Weichselian. The Markhida Line truncates the huge arcs of the Laya-Adzva and Rogovaya ice-pushed ridges protruding to the south. The latter moraines therefore reflect an older ice advance, probably also of Weichselian age. Still farther south, fluvially dissected morainic plateaus without lakes are of pre-Eemian age, because they plunge northwards under marine Eemian sediments. Shorelines of the large ice-dammed Lake Komi, identified between 90 and 110 m a.s.l. in the areas south of the Markhida Line, are radiocarbon dated to be older than 45 ka. The shorelines, incised into the Laya-Adzva moraines, morphologically interfinger with the Markhida moraines, indicating that the last ice advance onto the Russian mainland reached the Markhida Line during the Middle or Early Weichselian, before 45 ka ago.

Seismic Data from MSM21/4 to Knipovich Ridge on the western Svalbard margin

The seismic data were acquired north of the Knipovich Ridge on the western Svalbard margin during cruise MSM21/4. They were recorded using a Geometrics GeoEel streamer of either 120 channels (profiles p100-p208) or 88 channels (profiles p300-p805) with a group spacing of 1.56 m and a sampling rate of 2 kHz. A GI-Gun (2×1.7 l) with a main frequency of ~150 Hz was used as a source and operated at a shot interval of 6-8 s. Processing of profiles p100-p208 and p600-p805: Positions for each channel were calculated by backtracking along the profiles from the GI-Gun GPS positions. The shot gathers were analyzed for abnormal amplitudes below the seafloor reflection by comparing neighboring traces in different frequency bands within sliding time windows. To suppress surface-generated water noise, a tau-p filter was applied in the shot gather domain. Common mid-point (CMP) profiles were then generated through crooked-line binning with a CMP spacing of 1.5625 m. A zero-phase band-pass filter with corner frequencies of 60 Hz and 360 Hz was applied to the data. Based on regional velocity information from MCS data [Sarkar, 2012], an interpolated and extrapolated 3D interval velocity model was created below the digitized seafloor reflection of the high-resolution streamer data. This velocity model was used to apply a CMP stack and an amplitude-preserving Kirchhoff post-stack time migration. Processing of profiles p400-p500: Data were sampled at 0.5 ms and sorted into common midpoint (CMP) domain with a bin spacing of 5 m. Normal move out correction was carried out with a velocity of 1500 m s-1 and an Ormsby bandpass filter with corner frequencies at 40, 80, 600 and 1000 Hz was applied. The data were time migrated using the water velocity.

Population genetic and dispersal modeling data for Bathymodiolus mussels from the Mid-Atlantic Ridge

The zip folder comprises a text file and a gzipped tar archive. 1) The text file contains individual genotype data for 90 SNPs, 9 microsatellites and the mitochondrial ND4 gene that were determined in deep-sea hydrothermal vent mussels from the Mid-Atlantic Ridge (genus Bathymodiolus). Mussel specimens are grouped according to the population (pop)/location from which they have been sampled (first column). The remaining columns contain the respective allele/haplotype codes for the different genetic loci (names in the header line). The data file is in CONVERT format and can be directly transformed into different input files for population genetic statistics. 2) The tar archive contains NetCDF files with larval dispersal probabilities for simulated annual larval releases between 1998 and 2007. For each simulated vent location (Menez Gwen, Lucky Strike, Rainbow, Vent 1-10) two NetCDF files are given, one for an assumed pelagic larval duration of 1 year and the other one for an assumed pelagic larval duration of 6 months (6m).

(Table 1) Composition of volcanic glasses from the East Scotia Ridge

New U-series isotope data for lavas from the East Scotia backarc spreading center span a large range in (230Th/232Th) and (238U/232Th). Most of the backarc lavas have (238U/230Th) < 1, similar to the composition of mid-oceanic-ridge basalts (MORB). Lavas from two segments have (238U/230Th) > 1 and are enriched in fluid-mobile elements, implying a recent addition of a U-rich slab-derived component. The data from one segment suggest an influence from an aqueous fluid from altered MORB, whereas samples from the other slab-influenced segment show addition of a sediment melt. The slab-influenced samples form a distinct trend in the equiline diagram between aqueous fluid and sediment melt that is suggested to be a mixing line rather than an isochron.

(Table 1) Crystal sizes of measured samples from the Bush Hill region in the Gulf of Mexico and from ODP Leg 204 at the Hydrate Ridge offshore Oregon

The grain sizes of gas hydrate crystallites are largely unknown in natural samples. Single grains are hardly detectable with electron or optical microscopy. For the first time, we have used high-energy synchrotron diffraction to determine grain sizes of six natural gas hydrates retrieved from the Bush Hill region in the Gulf of Mexico and from ODP Leg 204 at the Hydrate Ridge offshore Oregon from varying depth between 1 and 101 metres below seafloor. High-energy synchrotron radiation provides high photon fluxes as well as high penetration depth and thus allows for investigation of bulk sediment samples. Gas hydrate grain sizes were measured at the Beam Line BW 5 at the HASYLAB/Hamburg. A 'moving area detector method', originally developed for material science applications, was used to obtain both spatial and orientation information about gas hydrate grains within the sample. The gas hydrate crystal sizes appeared to be (log-)normally distributed in the natural samples. All mean grain sizes lay in the range from 300 to 600 µm with a tendency for bigger grains to occur in greater depth. Laboratory-produced methane hydrate, aged for 3 weeks, showed half a log-normal curve with a mean grain size value of c. 40 µm. The grains appeared to be globular shaped.

Noble gas concentrations of oceanic crust and sediments from the Ninety East Ridge in the Indian Ocean

We have determined the concentrations and isotopic composition of noble gases in old oceanic crust and oceanic sediments and the isotopic composition of noble gases in emanations from subduction volcanoes. Comparison with the noble gas signature of the upper mantle and a simple model allow us to conclude that at least 98% of the noble gases and water in the subducted slab returns back into the atmosphere through subduction volcanism before they can be admixed into the earth's mantle. It seems that the upper mantle is inaccessible to atmospheric noble gases due to an efficient subduction barrier for volatiles.