Lithospheric structure below the eastern Arabian Sea and adjoining West Coast of India based on integrated analysis of gravity and seismic data

Four uniformly spaced regional gravity traverses and the available seismic data across the western continental margin of India, starting from the western Indian shield extending into the deep oceanic areas of the eastern Arabian Sea, have been utilized to delineate the lithospheric structure. The seismically constrained gravity models along these four traverses suggest that the crustal structure below the northern part of the margin within the Deccan Volcanic Province (DVP) is significantly different from the margin outside the DVP. The lithosphere thickness, in general, varies from 110–120 km in the central and southern part of the margin to as much as 85–90 km below the Deccan Plateau and Cambay rift basin in the north. The Eastern basin is characterised by thinned rift stage continental crust which extends as far as Laxmi basin in the north and the Laccadive ridge in the south. At the ocean–continent transition (OCT), crustal density differences between the Laxmi ridge and the Laxmi basin are not sufficient to distinguish continental as against an oceanic crust through gravity modeling. However, 5-6 km thick oceanic crust below the Laxmi basin is a consistent gravity option. Significantly, the models indicate the presence of a high density layer of 3.0 g/cm³ in the lower crust in almost whole of the northern part of the region between the Laxmi ridge and the pericontinental northwest shield region in the DVP, and also below Laccadive ridge in the southern part. The Laxmi ridge is underlain by continental crust upto a depth of 11 km and a thick high density material (3.0 g/cm³) between 11–26 km. The Pratap ridge is indicated as a shallow basement high in the upper part of the crust formed during rifting. The 15 –17 km thick oceanic crust below Laccadive ridge is seen further thickened by high density underplated material down to Moho depths of 24–25 km which indicate formation of the ridge along Reunion hotspot trace.     

Tectonic core of a sedimentary drift: a potential ridge propagation feature beneath the Blake Outer Ridge

The Blake Outer Ridge is a 480–kilometer long linear sedimentary drift ridge striking perpendicular to the North American coastline. By modeling free-air gravity anomalies we tested for the presence of a crustal feature that may control the location and orientation of the Blake Outer Ridge. Most of our crustal density models that match observed gravity anomalies require an increase in oceanic crustal thickness of 1–3 km on the southwest side of the Blake Outer Ridge relative to the northeast side. Most of these models also require 1–4 km of crustal thinning in zone 20–30 km southwest of the crest of the Blake Outer Ridge. Although these features are consistent with the structure of oceanic fracture zones, the Blake Outer Ridge is not parallel to adjacent known fracture zones. Magnetic anomalies suggest that the ocean crust beneath this feature formed during a period of mid-ocean ridge reorganization, and that the Blake Outer Ridge may be built upon the bathymetric expression of an oblique extensional feature associated with ridge propagation. It is likely that the orientation of this trough acted as a catalyst for sediment deposition with the start of the Western Boundary Undercurrent in the mid-Oligocene.     

Evolution of errors in the altimetric bathymetry model used by Google Earth and GEBCO

We analyze errors in the global bathymetry models of Smith and Sandwell that combine satellite altimetry with acoustic soundings and shorelines to estimate depths. Versions of these models have been incorporated into Google Earth and the General Bathymetric Chart of the Oceans (GEBCO). We use Japan Agency for Marine-Earth Science and Technology (JAMSTEC) multibeam surveys not previously incorporated into the models as “ground truth” to compare against model versions 7.2 through 12.1, defining vertical differences as “errors.” Overall error statistics improve over time: 50th percentile errors declined from 57 to 55 to 49 m, and 90th percentile errors declined from 257 to 235 to 219 m, in versions 8.2, 11.1 and 12.1. This improvement is partly due to an increasing number of soundings incorporated into successive models, and partly to improvements in the satellite gravity model. Inspection of specific sites reveals that changes in the algorithms used to interpolate across survey gaps with altimetry have affected some errors. Versions 9.1 through 11.1 show a bias in the scaling from gravity in milliGals to topography in meters that affected the 15–160 km wavelength band. Regionally averaged (>160 km wavelength) depths have accumulated error over successive versions 9 through 11. These problems have been mitigated in version 12.1, which shows no systematic variation of errors with depth. Even so, version 12.1 is in some respects not as good as version 8.2, which employed a different algorithm.     

Gravity anomalies and deep structure of the Ninetyeast Ridge north of the equator,eastern Indian Ocean — a hot spot trace model

The Ninetyeast Ridge north of the equator in the eastern Indian Ocean is actively deforming as evidenced by seismicity and its eastward subduction below the Andaman Trench. Basement of the ridge is elevated nearly 2 km with respect to the Bengal Fan; seismic surveys demonstrate continuity of the ridge beneath sediment for 700 km north of 10° N where the ridge plunges below the Fan sediment. The ridge is characterised by a free-air gravity high of 50 mgal amplitude and 350 km wavelength, and along-strike continuity of 1500 km in a north-south direction, closely fringing (locally, even abutting) the Andaman arc-trench bipolar gravity field. Regression analysis between gravity and bathymetry indicates that the ridge gravity field cannot be explained solely by its elevation. The ridge gravity field becomes gradually subdued northwards where overlying Bengal Fan sediments have a smaller density contrast with the ridge material. Our gravity interpretation, partly constrained by seismic data, infers that the ridge overlies significant crustal mass anomalies consistent with the hot spot model for the ridge. The anomalous mass is less dense by about 0.27 g cm^–3 than the surrounding oceanic upper mantle, and acts as a cushion for isostatic compensation of the ridge at the base of the crust. This cushion is up to 8 km thick and 400–600 km wide. Additional complexities are created by partial subduction of the ridge below the Andaman Trench that locally modifies the arc-trench gravity field.     

Bathymetry enhancement by altimetry-derived gravity anomalies in the East Sea (Sea of Japan)

The gravity-geologic method (GGM) was used to enhance the bathymetry of the East Sea (Sea of Japan) with satellite altimetry-derived free-air gravity anomalies and shipborne depth measurements. By comparison with the bathymetry model of Smith and Sandwell’s (SAS) approach (1994), GGM was found to have an advantage with short wavelength (≤12 km) components, while SAS better predicts longer wavelength (≥25 km) components, despite its dependency on density contrast. To mitigate this limitation, a tuning density contrast of 10.25 g/cm³ between seawater and the seafloor was primarily estimated by the downward continuation method and then validated by the check points method with GGM. Similarly, SAS is limited by the “A” value in low-pass part of the Wiener filter, which defines the effective range of the wavelength components on bathymetry. As a final result, we present an enhanced GGM bathymetry model by integrating all available data.     

Gravity crustal models and heat flow measurements for the Eurasia Basin,Arctic Ocean

The Gakkel Ridge in the Arctic Ocean with its adjacent Nansen and Amundsen Basins is a key region for the study of mantle melting and crustal generation at ultraslow spreading rates. We use free-air gravity anomalies in combination with seismic reflection and wide-angle data to compute 2-D crustal models for the Nansen and Amundsen Basins in the Arctic Ocean. Despite the permanent pack-ice cover two geophysical transects cross both entire basins. This means that the complete basin geometry of the world’s slowest spreading system can be analysed in detail for the first time. Applying standard densities for the sediments and oceanic crystalline crust, the gravity models reveal an unexpected heterogeneous mantle with densities of 3.30 × 10³, 3.20 × 10³ and 3.10 × 10³ kg/m³ near the Gakkel Ridge. We interpret that the upper mantle heterogeneity mainly results from serpentinisation and thermal effects. The thickness of the oceanic crust is highly variable throughout both transects. Crustal thickness of less than 1 km dominates in the oldest parts of both basins, increasing to a maximum value of 6 km near the Gakkel Ridge. Along-axis heat flow is highly variable and heat flow amplitudes resemble those observed at fast or intermediate spreading ridges. Unexpectedly, high heat flow along the Amundsen transect exceeds predicted values from global cooling curves by more than 100%.     

Geoid height versus topography of the Northern Ninetyeast Ridge: implications on crustal compensation

In this paper we focused on understanding the isostatic compensation of the Ninetyeast Ridge in the overall context of the Bay of Bengal oceanic lithosphere and the interaction of the ridge system with the north Andaman subduction zone from north of 7–18°N. This region is characterized by the initial interaction of the Kerguelen hotspot with the Bay of Bengal oceanic lithosphere. We used satellite altimeter-derived marine geoid, as it should comprehensively reflect the compensations caused by large spatial wavelength dominated deeper anomaly sources in a hotspot affected lithospheric load like the Ninetyeast Ridge. Our analyses of the geoid-to-topography ratio (GTR), residual geoid, gravity-to-topographic kernel and upward continuation of anomalies show the existence of two different types of source compensation bodies beneath the northern (12–18°N) and southern (7–12°N) Ninetyeast Ridge. In the northern region, the geoid to topography ratio varies from 0.63 ± 0.05 to 0.44 ± 0.03, while in the southern region it ranges from 1.34 ± 0.09 to 1.31 ± 0.07 which resulted in a north to south increase in the apparent compensation depth from ~9 to 28 km. The presence of a shallow Moho, low GTR, broader gravity to topography kernel and the absence of a ridge anomaly from the mantle density dominated upward continued anomaly at z = 300 km indicates that at the northern segment the underplated low density crustal melt is the dominant isostatic compensating body. However, at the southern ridge segment the high GTR, strong gravity-to-topography kernel and the subsistence of the anomaly at long wavelengths, even at z = 300 km represents the existence of large volumes of hotspot related underplated dense material as the source of compensation. The proximity of the dense source compensating body of the southern Ninetyeast Ridge to the Andaman subduction zone affected the regional mantle driven density gradient flow, as observed from the z = 300 km continued gravity anomaly. The existence of a southern Ninetyeast Ridge in such a transpressional regime has caused the formation of a forearc sliver at its eastern flank, which is a major crustal deformational structure developed as a result of ridge-trench collision.     

Lava morphology mapping by expert system classification of high-resolution side-scan sonar imagery from the East Pacific Rise, 9°–10°N

Fine-scale lava morphology has been classified on the ridge crest of the East Pacific Rise between 9°15′N and 10°02′N using an expert system classification method. This method establishes the means to classify complicated seafloor environments by integrating textural and geometric feature attributes from a high-resolution side-scan sonar dataset where ground-reference data are available from near-bottom visual observations. The classification in this study focuses upon mapping the lava morphology distribution of sheet, lobate, and pillow flows along the East Pacific Rise. The reliability of the classification has been evaluated using an accuracy assessment. The study region yields a coverage area of 37,814 m² (44%) for lobate flows; 10,421 m² (12%) for pillow flows; 15,096 m² (18%) for sheet flows; 19,679 m² (23%) for fissured areas; and 2,967 m² (3%) for shadows or no data. The systematic distribution of lava morphology along the ridge found in this study supports the idea of using the regional distribution of surface morphology as an indicator of emplacement dynamics and supports an organization of the volcanic plumbing system at a third order segmentation scale beneath mid-ocean ridges.     

Utilization of high resolution satellite gravity over the Carlsberg Ridge

The Carlsberg Ridge lies between the equator and the Owen fracture zone. It is the most prominent mid-ocean ridge segment of the western Indian Ocean, which contains a number of earthquake epicenters. Satellite altimetry can be used to infer subsurface geological structures analogous to gravity anomaly maps generated through ship-borne survey. In this study, free-air gravity and its 3D image have been generated over the Carlsberg Ridge using a very high resolution data base, as obtained from Geosat GM, ERS-1, Seasat and TOPEX/POSEIDON altimeter data. As observed in this study, the Carlsberg Ridge shows a slow spreading characteristic with a deep and wide graben (average width ∼15 km). The transform fault spacing confirms variable slow to intermediate characteristics with first and second order discontinuities. The isostatically compensated region of the Carlsberg Ridge could be demarcated with near zero contour values in the free-air gravity anomaly images over and along the Carlsberg Ridge axes and over most of the fracture zone patterns. Few profiles have been generated across the Carlsberg Ridge and the characteristics of slow/intermediate spreading ridge of various orders of discontinuity could be identified. It has also been observed in zero contour image as well as in the characteristics of valley patterns along the ridge from NW to SE that different spreading rates, from slow to intermediate, are occurring in different parts of the Carlsberg ridge. It maintains the morphology of a slow spreading ridge in the NW, where the wide and deep axial valley (∼1.5–3 km) also implies the pattern of a slow spreading ridge. However, a change in the morphology/depth of the axial valley from NW to SE indicates the nature of the Carlsberg Ridge as a slow to intermediate spreading ridge. For the prevailing security restrictions, lat./lon. coordinates have been omitted in few images.     

The northern boundary of the Jan Mayen microcontinent, North Atlantic determined from ocean bottom seismic, multichannel seismic, and gravity data

The Jan Mayen microcontinent was as a result of two major North Atlantic evolutionary cornerstones—the separation of Greenland from Norway (~54 Ma), accompanied by voluminous volcanic activity, and the jump of spreading from the Aegir to the Kolbeinsey ridge (~33 Ma), which resulted in the separation of the microcontinent itself from Eastern Greenland (~24 Ma). The resulting eastern and western sides of the Jan Mayen microcontinent are respectively volcanic and non-volcanic rifted margins. Until now the northern boundary of the microcontinent was not precisely known. In order to locate this boundary, two combined refraction and reflection seismic profiles were acquired in 2006: one trending S–N and consisting of two separate segments south and north of the island of Jan Mayen respectively, and the second one trending SW–NE east of the island. Crustal P-wave velocity models were derived and constrained using gravity data collected during the same expedition. North of the West Jan Mayen Fracture Zone (WJMFZ) the models show oceanic crust that thickens from west to east. This thickening is explained by an increase in volcanic activity expressed as a bathymetric high and most likely related to the proximity of the Mohn ridge. East of the island and south of the WJMFZ, oceanic Layers 2 and 3 have normal seismic velocities but above normal average crustal thickness (~11 km). The similarity of the crustal thickness and seismic velocities to those observed on the conjugate M?re margin confirm the volcanic origin of the eastern side of the microcontinent. Thick continental crust is observed in the southern parts of both profiles. The northern boundary of the microcontinent is a continuation of the northern lineament of the East Jan Mayen Fracture Zone. It is thus located farther north than previously assumed. The crust in the middle parts of both models, around Jan Mayen island, is more enigmatic as the data suggest two possible interpretations—Icelandic type of oceanic crust or thinned and heavily intruded continental crust. We prefer the first interpretation but the latter cannot be completely ruled out. We infer that the volcanism on Jan Mayen is related to the Icelandic plume.     

Particulate organic carbon (POC) in surface sediments of the Baltic Sea

In this study, particulate organic carbon (POC) contents and their distribution pattern in surficial sediments of the Baltic Sea are presented for 1,471 sampling stations. POC contents range from approx. 0.1% in shallow sandy areas up to 16% in deep muddy basins (e.g. Gotland Basin). Some novel relationships were identified between sediment mass physical properties (dry bulk density (DBD), grain size) and POC levels. Notably, the highest POC concentrations (about 10–17 mg cm^–3) occur in sandy mud to mud (60–100% mud content) with intermediate POC contents of about 3–7% and DBDs of 0.1–0.4 g cm^–3. Areas with this range in values seem to represent the optimum conditions for POC accumulation in the Baltic Sea. The maximum POC contents (8–16%) are found in fluid mud of the central Baltic Sea characterized by extremely low DBDs (<0.1 g cm^–3) and moderate POC concentrations (4–7 mg cm^–3). Furthermore, sediment mass accumulation rates (MAR), based on ²¹⁰Pb and ¹³⁷Cs measurements and available for 303 sites of the Baltic Sea, were used for assessing the spatial distribution of POC burial rates. Overall, these vary between 14 and 35 g m^–2 year^–1 in the mud depositional areas and, in total, at least 3.5 (±2.9) Mt POC are buried annually. Distribution patterns of POC contents and burial rates are not identical for the central Baltic Sea because of the low MAR in this area. The presented data characterize Baltic Sea sediments as an important sink for organic carbon. Regional differences in organic carbon deposition can be explained by the origin and transport pathways of POC, as well as the environmental conditions prevailing at the seafloor (morphology, currents, redox conditions). These findings can serve to improve budget calculations and modelling of the carbon cycle in this large brackish-water marginal sea.     

Crustal structure of the ultra-slow spreading Knipovich Ridge,North Atlantic,along a presumed amagmatic portion of oceanic crustal formation

The ultra-slow, asymmetrically-spreading Knipovich Ridge is the northernmost part of the Mid Atlantic ridge system. In the autumn of 2002 a combined ocean-bottom seismometer multichannel seismic (OBS/MCS) and gravity survey along the spreading direction of the Knipovich Ridge was carried out. The main objective of the study was to gain an insight into the crustal structure and composition of what is assumed to be an amagmatic segment of oceanic crust. P-wave velocity and Vp/Vs models were built and complemented by a gravity model. The 190 km long transect reveals a much more complex crustal structure than anticipated. The magmatic crust is thinner than the global average of 7.1 ± 1.0 km. The young fractured portion of Oceanic Layer 2 has low seismic velocities while the older part has normal seismic velocities and is broken into several rotated fault blocks seen as thickness variations of Layer 2. The youngest part of Oceanic Layer 3 is also dominated by low velocities, indicative of fracturing, seawater circulation and thermal expansion. The remaining portion of Layer 3 exhibits inverse variations in thickness and seismic velocity. This is explained by a sequence of periods of faster spreading (estimated to be up to 8 mm/year from interpretation of magnetic anomalies) when more normal gabbroic crust was being generated and periods of slower spreading (5.5 mm/year) when amagmatic stretching and serpentinization of the upper mantle occurred, and crust composed of mixed gabbro and serpentinized mantle was generated. The volumetric changes and upward fluid migration, associated with the process of serpentinization in this part of the crust, caused disruption to the overlying sedimentary layers.     

Heat flow derived from BSR and its implications for gas hydrate stability zone in Shenhu Area of northern South China Sea

Bottom simulating reflectors (BSRs), known as the base of gas hydrate stability zone, have been recognized and mapped using good quality three-dimensional (3D) pre-stack migration seismic data in Shenhu Area of northern South China Sea. Additionally, seismic attribute technique has been applied to better constrain on the distribution of gas hydrate. The results demonstrate that gas hydrate is characterized by “blank” zone (low amplitude) in instantaneous amplitude attribute. The thickness of gas hydrate stability zone inferred from BSR ranges from 125 to 355 m with an average of 240 m at sea water depth from 950 to 1,600 m in this new gas hydrate province. The volume of gas in-place bound in hydrate is estimated from 1.7 × 10⁹ to 4.8 × 10⁹ m³, with the most likely value of around 3.3 × 10⁹ m³, using Monte Carlo simulation. Furthermore, geothermal gradient and heat flow are derived from the depths of BSRs using a conductive heat transfer model. The geothermal gradient varies from 35 to 95°C km⁻¹ with an average of 54°C km⁻¹. Corresponding heat flow values range from 43 to 105 mW m⁻² with an average of 64 mW m⁻². By comparison with geological characteristics, we suggest that the distribution of gas hydrate and heat flow are largely associated with gas chimneys and faults, which are extensively distributed in Shenhu Area, providing easy pathways for fluids migrating into the gas hydrate stability zone for the formation of gas hydrate. This study can place useful constraints for modeling gas hydrate stability zone from measured heat flow data and understanding the mechanism of gas hydrate formation in Shenhu Area.     

Flexure across a continent–ocean fracture zone: the northern Falkland/Malvinas Plateau, South Atlantic

 Bathymetry, satellite-derived gravity, and interpreted seismic reflection data across the northern Falkland/Malvinas Plateau fossil continent–ocean transform rim may record the degree of mechanical coupling across the boundary after ridge–transform intersection time. The rim comprises a broad microcontinental block in the east and a continental marginal fracture ridge 50–100 km wide elsewhere. Free-air gravity anomalies tentatively suggest that the fracture ridge is locked against oceanic elastic lithosphere both to the north (Argentine Basin) and south (Central Falkland Basin). Received: 18 January 1996 / Revision received: 25 March 1995     

Multi-scale pore structure of COx claystone: Towards the prediction of fluid transport

We aim to relate the morphology of the pore network of finely porous claystones to their fluid transport properties. By using Focused Ion Beam in combination with Scanning Electron Microscopy (FIB/SEM), we image the pore network of COx claystone from 2D image stacks and as 3D reconstructed volumes. Our FIB/SEM samples are representative of the mesoscopic matrix clay. Porosity resolvable by this technique is in the range 1.7–5.9% with peak pore sizes of 50–90 nm. 3D pore network skeletonization provides connected pore volumes between end surfaces, tortuosity, density, and shortest pore paths with their pore size distribution. At higher resolution, 2D transmission electron microscopy (TEM) reveals large amounts of smaller pores (2–20 nm) between clay aggregates, associated to a local porosity of 14–25%, and peak sizes of 4–6 nm. Liquid permeability predictions with Katz–Thompson model, at the FIB/SEM volume scale and at the TEM surface scale, are in good agreement with macroscopic measurements (on the order of 10⁻²⁰ m²), showing that both mesopore sizes (peaks at 50–90 nm and 4–6 nm), located within the clay matrix, contribute to liquid transport.     

Investigation of drag crisis phenomenon using CFD methods

When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 10⁴ and 7.57 × 10⁵. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem.     

Gravity anomalies of the ridge-transform system in the South Atlantic between 31 and 34.5° S: Upwelling centers and variations in crustal thickness

To decipher the distribution of mass anomalies near the earth's surface and their relation to the major tectonic elements of a spreading plate boundary, we have analyzed shipboard gravity data in the vicinity of the southern Mid-Atlantic Ridge at 31–34.5° S. The area of study covers six ridge segments, two major transforms, the Cox and Meteor, and three small offsets or discordant zones. One of these small offsets is an elongate, deep basin at 33.5° S that strikes at about 45° to the adjoining ridge axes.By subtracting from the free-air anomaly the three-dimensional (3-D) effects of the seafloor topography and Moho relief, assuming constant densities of the crust and mantle and constant crustal thickness, we generate the mantle Bouguer anomaly. The mantle Bouguer anomaly is caused by variations in crustal thickness and the temperature and density structure of the mantle. By subtracting from the mantle Bouguer anomaly the effects of the density variations due to the 3-D thermal structure predicted by a simple model of passive flow in the mantle, we calculate the residual gravity anomalies. We interpret residual gravity anomalies in terms of anomalous crustal thickness variations and/or mantle thermal structures that are not considered in the forward model. As inferred from the residual map, the deep, major fracture zone valleys and the median, rift valleys are not isostatically compensated by thin crust. Thin crust may be associated with the broad, inactive segment of the Meteor fracture zone but is not clearly detected in the narrow, active transform zone. On the other hand, the presence of high residual anomalies along the relict trace of the oblique offset at 33.5° S suggests that thin crust may have been generated at an oblique spreading center which has experienced a restricted magma supply. The two smaller offsets at 31.3° S and 32.5° S also show residual anomalies suggesting thin crust but the anomalies are less pronounced than that at the 33.5° S oblique offset. There is a distinct, circular-shaped mantle Bouguer low centered on the shallowest portion of the ridge segment at about 33° S, which may represent upwelling in the form of a mantle plume beneath this ridge, or the progressive, along-axis crustal thinning caused by a centered, localized magma supply zone. Both mantle Bouguer and residual anomalies show a distinct, local low to the west of the ridge south of the 33.5° S oblique offset and relatively high values at and to the east of this ridge segment. We interpret this pattern as an indication that the upwelling center in the mantle for this ridge is off-axis to the west of the ridge.     

Gravity and S-wave modelling across the Jan Mayen Ridge,North Atlantic; implications for crustal lithology

The horizontal components from fourteen Ocean Bottom Seismometers deployed along four profiles focused along the western margin of the Jan Mayen microcontinent, North Atlantic, have been modelled with regard to S-waves, based on P-wave models obtained earlier. The seismic models have furthermore been constrained by 2D gravity modelling. High V _p/V _s-ratios (2.3–7.9) within the Cenozoic sedimentary section are attributed to significant porosities, whereas V _p/V _s-ratios in the order of 1.9–2.2 for the Mesozoic and Paleozoic sedimentary rocks indicate shale-dominated lithology throughout the area. The eastern side of the Jan Mayen Ridge is interpreted as a passive, volcanic margin, based on relatively high crustal V _p/V _s-ratios (1.9), whereas lower V _p/V _s-ratios (1.75–1.8) suggest the presence of intermediate composition crust and non-volcanic margin on the western side of the ridge. In the westernmost part of the Jan Mayen Basin, slightly increased upper mantle V _p/V _s-ratios may indicate some degree of serpentization of upper mantle peridotites.     

Multidisciplinary investigations exploring indicators of gas hydrate occurrence in the Krishna–Godavari Basin offshore, east coast of India

We report some main results of multidisciplinary investigations carried out within the framework of the Indian National Gas Hydrate Program in 2002–2003 in the Krishna–Godavari Basin offshore sector, east coast of India, to explore indicators of likely gas hydrate occurrence suggested by preliminary multi-channel seismic reflection data and estimates of gas hydrate stability zone thickness. Swath bathymetry data reveal new evidence of three distinct geomorphic units representing (1) a delta front incised by several narrow valleys and mass flows, (2) a deep fan in the east and (3) a WNW–ESE-trending sedimentary ridge in the south. Deep-tow digital side-scan sonar, multi-frequency chirp sonar, and sub-bottom profiler records indicate several surface and subsurface gas-escape features with a highly resolved stratification within the upper 50 m sedimentary strata. Multi-channel seismic reflection data show the presence of bottom simulating reflections of continuous to discrete character. Textural analyses of 76 gravity cores indicate that the sediments are mostly silty clay. Geochemical analyses reveal decreasing downcore pore water sulphate (SO₄ ²⁻) concentrations (28.7 to <4 mM), increasing downcore methane (CH₄) concentrations (0–20 nM) and relatively high total organic carbon contents (1–2.5%), and microbial analyses a high abundance of microbes in top core sediments and a low abundance of sulphate-reducing bacteria in bottom core sediments. Methane-derived authigenic carbonates were identified in some cores. Combined with evidence of gas-escape features in association with bottom simulating reflections, the findings strongly suggest that the physicochemical conditions prevailing in the study area are highly conducive to methane generation and gas hydrate occurrence. Deep drilling from aboard the JOIDES Resolution during 2006 has indeed confirmed the presence of gas hydrate in the Krishna–Godavari Basin offshore.     

An integrated thermomechanical model for transform continental margin evolution

 A 2D numerical thermal model for transform continental margin evolution is presented that calculates thermally driven uplift and subsidence profiles across the margin, for any margin segment assuming both regional and local isostasy. Lateral variations in the magnitude of continental uplift along the transform are predicted. For a margin with a length of 900 km, with a spreading rate of 1 cm yr^-1, maximum continental uplift of 1300–1400 m is calculated, assuming local isostasy. Using a regional isostatic approximation, maximum uplift is reduced substantially to 335–470 m, and the exact magnitude, location, and timing of the maximum effect depends strongly on the assumption of a coupled or decoupled continent–ocean boundary. The length of time a margin point experiences continent–ocean shearing prior to ridge passing is also shown to be very significant. Received: 23 February 1995 / Revision received: 17 August 1995