Fault system and thermal regime in the vicinity of site NGHP-01-10, Krishna–Godavari basin, Bay of Bengal

Drilling/coring activities onboard JOIDES Resolution for hydrate resource estimation have confirmed gas hydrate in the continental slope of Krishna-Godavari (KG) basin, Bay of Bengal and the expedition recovered fracture filled gas hydrate at the site NGHP-01-10. In this paper we analyze high resolution multi-channel seismic (MCS), high resolution sparker (HRS), bathymetry, and sub-bottom profiler data in the vicinity of site NGHP-01-10 to understand the fault system and thermal regime. We interpreted the large-scale fault system (>5 km) predominantly oriented in NNW-SSE direction near NGHP-01-10 site, which plays an important role in gas hydrate formation and its distribution. The increase in interval velocity from the baseline velocity of 1600 m/s to 1750–1800 m/s within the gas hydrate stability zone (GHSZ) is considered as a proxy for the gas hydrate occurrence, whereas the drop in interval velocity to 1400 m/s suggest the presence of free gas below the GHSZ. The analysis of interval velocity suggests that the high concentration of gas hydrate occurs close to the large-scale fault system. We conclude that the gas hydrate concentration near site NGHP-01-10, and likely in the entire KG Basin, is controlled primarily by the faults and therefore has high spatial variability.We also estimated the heat flow and geothermal gradient (GTG) in the vicinity of NGHP-01-10 site using depth and temperature of the seafloor and the BSR. We observed an abnormal GTG increase from 38 °C/km to 45 °C/km at the top of the mound, which remarkably agrees with the measured temperature gradient at the mound (NGHP-01-10) and away from the mound (NGHP-01-03). We analyze various geological scenarios such as topography, salinity, thermal non-equilibrium of BSR and fluid/gas advection along the fault system to explain the observed increase in GTG. The geophysical data along with the coring results suggest that the fluid advection along the fault system is the primary mechanism that explains the increase in GTG. The approximate advective fluid flux estimated based on the thermal measurement is of the order of few tenths of mm/yr (0.37–0.6 mm/yr).     

Modelling of thermal convection in sedimentary basins and its relevance to diagenetic reactions

Mathematical calculations of thermal convection have been carried out using a porous three-player model to simulate pore-water flow in a sedimentary basin with layers of different permeabilities. The calculated flow lines demonstrate that even very thin layers (< 1 m) with low permeability, like shales in a more permeable sandstone sequence, will split potentially larger convection cells into smaller units of clean sandstone beds which may then be too small to exceed the critical Rayleigh number. In a situation where we have horizontal isotherms, a sandstone bed with 1 Darcy permeability must have a thickness of more than 330 m without shale interbeds for convection to occur. This situation is rarely met in sedimentary basins. In the case of sloping isotherms, non-Rayleigh convection will always occur but the flow rate is proportional to the effective thickness of the bed and the slope of the isotherms. Calculations suggest that the flow rates are insignificant for diagnetic reactions except in situations where we have steeply sloping isotherms eg, around salt domes and igneous or hydrothermal intrusions.Analyses of formation waters from sedimentary basins like the North Sea often show evidence of a crude stratification of the pore water with respect to salinity and isotopic composition. This is evidence suggesting that large scale convection, or other types of mixing, does not take place and positive salinity depth gradient may help to physically stabilize the formation water. The fact that many of the shallow reservoir rocks from the North Sea have formation water with low salinity and negative δ ¹⁸O values suggest that this is modified meteoric water and that migration of petroleum from deeper parts of the basin occurred as a separate phase along restricted migration pathways and was not associated with a high flux of pore water. In the absence of thermal convection the total pore water flux through sandstones will be rather small except locally where we may have focused compactional flow through small cross-sections. Diagenetic reactions will, therefore, normally be relatively isochemical during deeper burial.     

Flow simulation with reactive transport applied to carbonate rock diagenesis

Significant volumes of the known hydrocarbon reserves are found in carbonate rocks, many of these dolomitized. The spatial distribution of diagenesis on these rocks is one of the main challenges in oil reservoir modeling. Reactive transport models can be a powerful tool to understand the active diagenetic processes and their effects on the quality of these reservoirs. In this study it was used, for the first time, the CMG-GEM simulator to model diagenetic evolution of a carbonate sequence, subjected to compaction-driven and geothermal flow in a simulated period of 200 thousand years. It was simulated carbonate cementation, dolomitization and dissolution, with and without presence of faults. Among the analyzed variables, the volume of circulating fluid was the most important factor. For both mechanisms, flow simulated velocities obtained had magnitudes smaller than 10⁻⁶ m/day. Diagenesis was insignificant for these low speeds at simulated time interval. Only dolomitized facies presented relevant diagenesis in form of calcite dissolution and dolomite precipitation. Simulations with flow rates of 1 m/day revealed a considerable increase in observed diagenesis, especially in carbonate cementation and in porosity enhancement. Diagenesis was more pronounced in more permeable sediments, highlighting the role of fluid flow in diagenetic reactions. Relative dissolution was greatly reduced during simulations performed with absence of dolomite and dolomitization reactions. The presence of faults strongly influences spatial distribution of diagenesis, especially relatively to dissolution. More permeable facies were more dissolved near fault, decreasing with increasing distance. Low permeability facies, as mudstones, are not dissolved, even near fault. Spatial distribution of diagenesis would then be dependent mainly on the quality of original pore structure, of fault presence and mineral composition of rock.     

Litter production and litter elemental composition in two rehabilitated Kandelia obovata mangrove forests in Jiulongjiang Estuary,China

Spatial and seasonal variations in litter production and C, N, and P concentrations were compared between the 24 and 48 year old Kandelia obovata mangrove forests in the Jiulongjiang estuary, China. The 24 yr forest had significantly higher production of total, leaf and branch litter, but lower flower and fruit litter than the 48 yr forest. Total, leaf and branch litter production were significantly positively correlated to monthly temperature and rainfall. Spatial patterns of litter production among the inner, mid and outer zones in the same forest were similar to those of tree heights. C, N and P concentrations of leaf litter showed significant seasonality but varied little among these three forest zones. C/N and N/P ratios of leaf litter were significantly lower for the 24 yr forest than those for the 48 yr forest. During the entire sampling year, total litter of the 24 and 48 yr forests contained 590.31 and 437.31 g C m⁻² yr⁻¹, 8.46 and 5.47 g N m⁻² yr⁻¹, 1.92 and 1.16 g P m⁻² yr⁻¹, respectively.     

A comparative study between present and palaeo-heat flow in the Qiangtang Basin,northern Tibet,China

The Qiangtang Basin is a significant prospective area for hydrocarbon and gas hydrate resources in the Tibetan Plateau, China. However, relatively little work has been performed to characterise heat flow in this basin, which has restricted petroleum and gas hydrate exploration. In this study, we compare present and palaeo-heat flow in the Qiangtang Basin to provide information on geothermal regime, hydrocarbon generation and permafrost that is necessary for further petroleum and gas hydrate exploration. We base our study on temperature data from a thermometer well, thermal conductivity tests, vitrinite reflectance data, homogenisation temperature data from fluid inclusions, stratigraphic information and a time-independent modelling approach. Our results indicate that in the central Qiangtang Basin, the present thermal gradient is approximately 15.5 °C/km, and heat flow is approximately 46.69 mW/m². Heat flow in the Qiangtang Basin is not relatively stable since the Early Jurassic, as previous research has suggested, and it is generally decreasing with time. Additionally, there is a clear difference between the hottest thermal regime of the southern and northern Qiangtang Depressions during Cretaceous to Pleistocene time. In the southern Qiangtang Depression, the palaeogeothermal gradient is approximately 32.0 °C/km, and palaeo-heat flow is approximately 70 mW/m². However, in the northern Qiangtang Depression, the palaeogeothermal gradient exceeds 81.8 °C/km, and palaeo-heat flow is greater than 172.09 mW/m². The high thermal regime in the northern Qiangtang Depression is driven mainly by hydrothermal convection. Gas reservoirs are possible targets for hydrocarbon exploration in this depression. Currently, the northwestern part of the northern Qiangtang Depression is the most favourable area for gas hydrate exploration in the Qiangtang Basin.     

Computations of the energy flux to mixing processes via baroclinic wave drag on barotropic tides

The geographical distribution of barotropic to baroclinic transfer of tidal energy by baroclinic wave drag in the abyssal ocean is estimated. Using tidal velocities from a state-of-the-art numerical tidal model, the total loss of barotropic tidal energy in the deep ocean (between 70°S and 70°N and at depths greater than 1000 m) is estimated to be about 0.7 TW (M₂) corresponding to a mean value of the energy flux (e) of 2.4×10⁻³ W/m². The distribution of e is however highly skewed with a median of about 10⁻⁶ W/m². Only 10% of the area is responsible for more than 97% of the total energy transfer.To assess the possible influence of the relatively coarse bathymetry representation upon the present estimate, complementary calculations using better resolved sea floor topography are carried out over a control area around the Hawaiian Ridge. There are no major differences between the results achieved using the two different bathymetry databases. Fluxes of about 16 GW or 6×10⁻³ W/m² are computed in both cases, and the main contributions to the total fluxes originate in the same range of e-values and cover equally large parts of the total area.It is not clear whether the present model is valid at flat or subcritical bottom slopes. However, for the Hawaiian region, only 2% of the total energy flux as calculated in the present study originates in areas of critical and subcritical slopes.     

Single- and two-phase fluid flow properties of cataclastic fault rocks in porous sandstone

Understanding the impact of faults on fluid flow in the subsurface is important for the extraction of oil, gas and groundwater as well as the geological storage of waste products. We address two problems present in current industry-standard workflows for fault seal analysis that may lead to fault rocks not being represented adequately in computational fluid flow models. Firstly, fluid flow properties of fault rocks are often measured only for small-scale faults with throws not exceeding a few centimetres. Large seismic-scale faults (throws >20 m) are likely to act as baffles or conduits to flow but they are seldom recovered from subsurface cores and consequently fault rock data for them is sparse. Secondly, experimental two-phase fluid flow data is lacking for fault rocks and, consequently, uncertainties exist when modelling flow across faults in the presence of two or more immiscible phases. We present a data set encompassing both single- and two-phase fluid flow properties of fault and host rocks from the 90-Fathom fault and its damage zone at Cullercoats Bay, NE England. Measurements were made on low-throw single and zones of deformation bands as well as on slip-surface cataclasites present along the ～120 m throw main fault. Samples were analysed using SEM and X-ray tomography prior to petrophysical measurements. We show that single deformation bands, deformation band zones and slip-surface cataclasites exhibit dissimilar single- and two-phase fluid flow properties. This is due to grain-size reduction being more pronounced in slip-surface cataclasites and changes in microstructure being fault-parallel for deformation bands but mostly fault-perpendicular for slip-surface cataclasites. A trend of fault rocks with low absolute permeabilities exhibiting lower relative permeabilities than more permeable rocks at the same capillary pressure is evident.     

Advective cooling within sedimentary rift basins—application to the Upper Rhinegraben (Germany)

Measurements of dispersed vitrinite along several exploration wells within the northern Rhinegraben are indicative of a thermal graben history that is influenced by a combination of basal conductive and groundwater-flow related convective heat transfer. To determine the conductive/convective components of heat transfer within the rift today, a series of 2D numerical groundwater flow and heat models are developed along a cross-sectional transect across the northern Rhinegraben. Fault zone permeability is varied in the simulations of these models to determine the possible fluid pathways and the effects of circulating groundwater on the graben temperature field. Depending on the fault permeability, negative thermal anomalies always develop in areas of cold recharging groundwater along the graben flanks regardless of fault permeability, whereas hot discharging groundwater near the topographic low of the graben only results in positive thermal anomalies under the assumption of high fault permeability. Simulation results suggest that the modern groundwater flow system has an overall net cooling effect on the temperature field of the rift.Without convective cooling by groundwater, vitrinite reflectance levels in wells would be expected to be much higher on average than are observed. Although relatively high heat flow densities (100 mW/m²) are documented in the Rhinegraben, an average of only 65 mW/m² would be sufficient to produce the observed vitrinite reflectance levels. Thus, a long-lived (>10 My) cooling convective fluid flow in combination with a high basal heat flow seems to be active.     

A geochemical traverse along the “Sperchios Basin – Evoikos Gulf” graben (Central Greece): Origin and evolution of the emitted fluids

The studied area is a 130 km long fast spreading graben in Central Greece. Its complex geodynamical setting includes both the presence of a subduction slab at depth responsible for the recent (Quaternary) volcanic activity in the area and the western termination of a tectonic lineament of regional importance (the North-Anatolian fault). A high geothermal gradient is made evident by the presence of many thermal springs with temperatures from 19 to 82 °C, that discharge along the normal faults bordering the graben.In the period 2004–2012, 58 gas and 69 water samples were collected and their chemical and isotopic analysis revealed a wide range of compositions.Two main groups of thermal waters can be distinguished on the basis of their chemical composition. The first, represented by dilute waters (E.C. <0.6 mS/cm) of the westernmost sites, is characterised by the presence of CH₄-rich and mixed N₂–CH₄ gases. The second displays higher salinities (E.C. from 12 to 56 mS/cm) due to mixing with a modified marine component. Reservoir temperatures of 150–160 °C were estimated with cationic geothermometers at the easternmost sites.Along the graben, from west to east, the gas composition changes from CH₄- to CO₂-dominated through mixed N₂–CH₄ and N₂–CO₂ compositions, while at the same time the He isotopic composition goes from typical crustal values (<0.1 R/R_A) up to 0.87 R/R_A, showing in the easternmost sites a small (3–11%) but significant mantle input. The δ¹³C values of the CO₂-rich samples suggest a mixed origin (mantle and marine carbonates).     

Seismic imaging of a fractured gas hydrate system in the Krishna–Godavari Basin offshore India

Gas hydrate was discovered in the Krishna–Godavari (KG) Basin during the India National Gas Hydrate Program (NGHP) Expedition 1 at Site NGHP-01-10 within a fractured clay-dominated sedimentary system. Logging-while-drilling (LWD), coring, and wire-line logging confirmed gas hydrate dominantly in fractures at four borehole sites spanning a 500 m transect. Three-dimensional (3D) seismic data were subsequently used to image the fractured system and explain the occurrence of gas hydrate associated with the fractures. A system of two fault-sets was identified, part of a typical passive margin tectonic setting. The LWD-derived fracture network at Hole NGHP-01-10A is to some extent seen in the seismic data and was mapped using seismic coherency attributes. The fractured system around Site NGHP-01-10 extends over a triangular-shaped area of ∼2.5 km² defined using seismic attributes of the seafloor reflection, as well as “seismic sweetness” at the base of the gas hydrate occurrence zone. The triangular shaped area is also showing a polygonal (nearly hexagonal) fault pattern, distinct from other more rectangular fault patterns observed in the study area. The occurrence of gas hydrate at Site NGHP-01-10 is the result of a specific combination of tectonic fault orientations and the abundance of free gas migration from a deeper gas source. The triangular-shaped area of enriched gas hydrate occurrence is bound by two faults acting as migration conduits. Additionally, the fault-associated sediment deformation provides a possible migration pathway for the free gas from the deeper gas source into the gas hydrate stability zone. It is proposed that there are additional locations in the KG Basin with possible gas hydrate accumulation of similar tectonic conditions, and one such location was identified from the 3D seismic data ˜6 km NW of Site NGHP-01-10.     

New heat flow measurements in the Ulleung Basin, East Sea (Sea of Japan): relationship to local BSR depth, and implications for regional heat flow distribution

In July 2007, new marine heat flow data were collected at ten sites (HF01–10) in the central and southwestern sectors of the Ulleung Basin (East Sea or Sea of Japan) as part of regional gas hydrate research. In addition, cores were collected at five of these sites for laboratory analysis. The results show that the geothermal gradient ranged from 103–137 mK/m, and the in-situ thermal conductivity from 0.82–0.95 W/m·K. Laboratory measurements of thermal conductivity were found to deviate by as much as 40% from the in-situ measurements, despite the precautions taken to preserve the cores. Based on the in-situ conductivity, the heat flow was found to increase with water depth toward the center of the basin, ranging from 84–130 mW/m². Using a simple model, we estimated the heat flow from the depths of the BSR, and compared this with the observed heat flow. In our study area, the two sets of values were quite consistent, the observed heat flows being slightly higher than the BSR-derived ones. The evaluation of regional pre-1994 data revealed that the heat flow varied widely from 51–157 mW/m² in and around the basin. Due to a large scatter in these older data, a clear relationship between heat flow and water depth was not evident, in contrast to what would be expected for a rifted sedimentary basin. This raises the question as to whether the pre-1994 data represent the true background heat flow from the underlying basin crust since the basin opening, and/or whether they contain large measurement errors. In fact, evidence in support of the latter explanation exists. BSRs are generally found in the deep parts of the basin, and vary by only ±15 m in depth below the seafloor. From the average BSR depth, we inferred the background heat flow using a simple model, which in the case of the Ulleung Basin is approximately 120 and 80 mW/m² for 2.5 and 1 km below sea level, respectively.     

The thermohaline structure and evolution of the deep waters in the Canada Basin,Arctic Ocean

Below the sill depth (at about 2400 m) of the Alpha-Mendeleyev ridge complex, the waters of the Canada Basin (CB) of the Arctic Ocean are isolated, with a ¹⁴C isolation age of about 500 yr. The potential temperature θ decreases with depth to a minimum θ_m≈−0.524°C near 2400 m, increases with depth through an approximately 300 m thick transition layer to θ_h≈−0.514°C, and then remains uniform from about 2700 m to the bottom at 3200–4000 m. The salinity increases monotonically with depth through the deep θ_m and transition layer from about 34.952 to about 34.956 and then remains uniform in the bottom layer. A striking staircase structure, suggestive of double-diffusive convection, is observed within the transition layer. The staircase structure is observed for about 1000 km across the basin and has been persistent for more than a decade. It is characterized by 2–3 mixed layers (10–60 m thick) separated by 2–16 m thick interfaces. Standard formulae, based on temperature and salinity jumps, suggest a double-diffusive heat flux through the staircase of about 40 mW m⁻², consistent with the measured geothermal heat flux of 40–60 mW m⁻². This is to be expected for a scenario with no deep-water renewal at present as we also show that changes in the bottom layer are too small to account for more than a small fraction of the geothermal heat flux. On the other hand, the observed interfaces between mixed layers in the staircase are too thick to support the required double-diffusive heat flux, either by molecular conduction or by turbulent mixing, as there is no evidence of sufficiently vigorous overturns within the interfaces. It therefore seems, that while the staircase structure may be maintained by a very weak heat flux, most of the geothermal heat flux is escaping through regions of the basin near lateral boundaries, where the staircase structure is not observed. The vertical eddy diffusivity required in these near-boundary regions is O(10⁻³) m² s⁻¹. This implies Thorpe scales of order 10 m. We observe what may be Thorpe scales of this magnitude in boundary-region potential temperature profiles, but cannot tell if they are compensated by salinity. The weak stratification of the transition layer means that the large vertical mixing rate implies a local dissipation rate of only O(10⁻¹⁰) W kg⁻¹, which is not ruled out by plausible energy budgets. In addition, we discuss an alternative scenario of slow, continuous renewal of the CB deep water. In this scenario, we find that some of the geothermal heat flux is required to heat the new water and vertical fluxes through the transition layer are reduced.     

Triangulation-to-GPS and GPS-to-GPS geodesy in Trinidad, West Indies: Neotectonics, seismic risk, and geologic implications

We estimated horizontal velocities at 25 sites first surveyed in a 1901-1903 British Ordnance Survey triangulation and then resurveyed with GPS in 1994-1995 to identify Trinidad's principal active on-land faults, quantify fault-slip-rates, and test for elastic locking. Our best-fit single-fault elastic dislocation model put 12 ± 3 mm/yr of dextral strike-slip on the Central Range Fault (1-2 km locking depth), an apparently aseismic active fault. The estimated motions also showed statistically insignificant horizontal motion (2.2 ± 1.8 mm/yr of right-slip; 2.7 ± 2.0 mm/yr of N-S shortening) on the eastward on-strike extension of the El Pilar Fault, known to be the active transform fault in Venezuela. Repeat GPS measurements made between 1994 and 2005 at two sites spanning the island north to south showed a 14 ± 3 mm/yr eastward (plate-motion-parallel) dextral velocity differential, consistent with our best-fit historic (1901-1995) fault-slip-rate. Paleoseismology trenching demonstrates that the Central Range Fault cuts <5000-year-old sediment and is capped by ∼550-year-old sediment, suggesting that it may be locked and may have ruptured at least once during this time interval. About ∼5 mm/yr of slip could be taken up on the Los Bajos Fault and additional faults in the offshore south of Trinidad. The existing 1901-1995 and 1994-2005 geodetic data alone cannot resolve whether the Central Range Fault is essentially creeping (≤1-2 km locking depth) or locked to a more standard depth of 10 km.     

The nested structure of Arctic thermohaline intrusions

Arctic thermohaline intrusions have a “nested” temperature/salinity structure characterized by the lining up of widely separated profiles along a series of well-defined lines in the T–S plane. The nesting pattern is coherent across much of the Arctic Basin (>2000 km), and roughly 90% of the water column between 150 and 350 m depth is found to lie along these nesting lines. We propose the nested structure results from a type of slanted convection occurring within the thick, salt-finger stratified layers in the intrusions. The convection cells are reminiscent of those investigated in laboratory experiments done by Thorpe, Hutt, and Soulsby, and it is estimated that the cells have horizontal dimensions ranging from 50 to 100 km at depths close to the Atlantic Water temperature maximum. Simple theoretical ideas suggest the convective cells may appear when the intrusions reach a critical amplitude, driving them toward a nested configuration. Our analysis provides a new estimate of the effective lateral diffusivity due to the intrusions (≃50 m² s⁻¹), as well as an estimate of the vertical diffusivity near the core of the Atlantic Water layer (≃0.8×10⁻⁵ m² s⁻¹).     

Habitat,growth and physiological ecology of a basaltic community of Ridgeia piscesae from the Juan de Fuca Ridge

The vestimentiferan tubeworm Ridgeia piscesae is an ecosystem-structuring organism in the hydrothermal vent environments of the Northeast Pacific. During this study, a single representative aggregation of the long-skinny morphotype of R. piscesae from the main endeavor segment was monitored for 3 yr before being collected in its entirety with a hydraulically actuated collection device manipulated in situ by a research vehicle. Vestimentiferan growth rates in this aggregation were determined by staining the exterior of the tubes and measuring newly deposited tube sections. The average growth rate of R. piscesae in this aggregation was very low in both years of the growth study (3.2 mm yr⁻¹). Although the incidence of plume damage from partial predation was very high (>95%), mortality was very low (<4% yr⁻¹). The distribution and the very tight clustering of recently recruited individuals indicated gregarious settlement behavior that is hypothesized to be partly due to biotic cues from settled larvae. Coupled measurements of vent fluid sulfide concentration and temperature were used to calculate the exposure of the vestimentiferans to sulfide from short- and long-term temperature monitoring. Plume-level temperature records indicate that most of the time individuals in this aggregation were exposed to extremely low levels of vent fluid, and therefore sulfide (<0.1 μM), while their posterior sections were consistently exposed to sulfide concentrations in the 100 μM range. A rootball-like structure formed the common base of the aggregation. In contrast to the anterior sections of the tubeworm tubes, the portions of the tubes within the “rootball” were freely permeable to sulfide. The results of this study show that R. piscesae, unlike vestimentiferans from the East Pacific Rise, can survive and grow in areas of low diffuse vent flow with very low plume-level exposure to sulfide. We propose that this morphotype of R. piscesae has the ability to acquire sulfide from sources near their posterior ends, similar to some species of cold seep vestimentiferans from the Gulf of Mexico. The ability of this single species of vestimentiferan to survive low exposure to vent flow with low mortality coupled with sulfide uptake across posterior tube sections may help explain the occurrence of a single vent vestimentiferan species in a wide variety of habitat conditions at hydrothermal vent sites in the Northeast Pacific.     

Reynolds stress and TKE production in an estuary with a tidal bore

We report new measurements of the turbulent properties of the flow in a tidally energetic estuarine channel of almost uniform cross-section. A high-frequency (1.2 MHz), bottom-mounted Acoustic Doppler Current Profiler (ADCP) has been used to observe the velocity field at a sampling rate of 10 Hz in parallel with measurements of the surface elevation by tide gauges. Our data have been analysed using the Variance Method to determine turbulent kinetic energy (TKE), shear stress and TKE production over the tidal cycle with a time resolution of 60 s. During the highly energetic but brief flood period, when the surface axial velocity reaches 2 m s⁻¹, we observed large values of stress (>2 Pa) and shear production (5 W m⁻³). TKE is also input through the release of energy in the bore itself which results in a brief but intense injection of energy at the bore front with large transient TKE levels (100 J m⁻³). Subsequent input by shear production maintains TKE levels which are generally lower (20 J m⁻³) than the strong peak associated with the bore for the rest of the flood. On the ebb, the flow is relatively tranquil with maximum speeds 0.5 m s⁻¹ and peak TKE production rates of 0.1 W m⁻³.The flow and elevation data have also been used to estimate the energy fluxes into and out of the estuary. Short (1 h), intense energy inputs (8 MW at springs) on the flood flow are largely balanced by longer, less intense seaward energy flow on the ebb. The net energy input is found to be 0.1 MW at springs which is consistent with estimates of upstream dissipation. Peak dissipation in the bore itself may exceed the mean energy input but it is active only for a small fraction of the tidal cycle and its average contribution does not exceed 12% of total dissipation.     

Seasonal oxygen and carbon isotope variability in euthecosomatous pteropods from the Sargasso Sea

We examine seasonal variations in the stable carbon and oxygen isotopic composition of individual shells of the pteropods Limacina inflata and Styliola subula, collected from Oceanic Flux Program sediment traps (at 500 m depth) near Bermuda in the western Sargasso Sea. Calcification depths estimated from L. inflata δ¹⁸O vary between 200 and 650 m in late winter and spring, and between 50 and 250 m in late summer and fall. S. subula shows similar seasonal variability with calcification depths between 250 and 600 m in late winter and spring and 50–400 m in late summer and fall. These results suggest that both species calcify across a greater range of depths than indicated by previous geochemical studies. Furthermore, the data indicate that these species change their calcification depth in conjunction with changes in thermal stratification of the water column. Pteropod shell δ¹³C values vary inversely with δ¹³C_DIC but show a positive correlation with seawater [CO₃²⁻] and temperature after depth differences in δ¹³C_DIC are accounted for. We hypothesize that either the influence of temperature on metabolic CO₂ incorporation during shell growth and/or the influence of ambient [CO₃²⁻] on shell geochemistry can explain these relationships. Taken together, the individual shell δ¹⁸O and δ¹³C data suggest that shell calcification, and by inference the life cycle, of these pteropods is several months or less. Individual pteropod shell analyses have potential for contributing to our understanding of the environmental parameters that play a role in seasonal calcification depth shifts, as well as to our knowledge of past upper ocean thermal structure.     

An estimate of Atlantic Ocean thermal energy conversion (OTEC) resources

An interhemispheric box model of the Atlantic thermohaline circulation (THC) is modified by replacing the tropical box with two vertically resolved sub-domains. Seawater flows from large-scale ocean thermal energy conversion (OTEC) are allowed in one of the tropical sub-domains. Under present conditions and standardized OTEC operations, steady-state net power production density would reach a maximum of about 80 kW/km² (corresponding to 1.8 TW) with a cold seawater withdrawal per unit area of about 14 m/yr. This maximum reflects the impact of large OTEC flows on the oceanic thermal structure, although the THC would not be significantly affected. It is larger than a recently suggested worldwide value of the order of 30 kW/km² because of the relative strength of the Atlantic THC. Under asymmetric high-latitude warming scenarios potentially representative of current climatic trends, a substantial weakening or a reversal of the THC are possible. In the former case, recoverable OTEC resources could practically vanish. In the latter case, the emergence of a stronger reverse THC eventually could boost OTEC resources. Such events are hypothetical and would unfold over centuries, but the mere possibility of their occurrence challenges the accepted notion that OTEC resources are forever renewable.     

Source rock characteristics and hydrocarbon generation modelling of Upper Cretaceous Mukalla Formation in the Jiza-Qamar Basin,Eastern Yemen

The Upper Cretaceous Mukalla coals and other organic-rich sediments which are widely exposed in the Jiza-Qamar Basin and believed to be a major source rocks, were analysed using organic geochemistry and petrology. The total organic carbon (TOC) contents of the Mukalla source rocks range from 0.72 to 79.90% with an average TOC value of 21.50%. The coals and coaly shale sediments are relatively higher in organic richness, consistent with source rocks generative potential. The samples analysed have vitrinite reflectance in the range of 0.84–1.10 %Ro and pyrolysis T_max in the range of 432–454 °C indicate that the Mukalla source rocks contain mature to late mature organic matter. Good oil-generating potential is anticipated from the coals and coaly shale sediments with high hydrogen indices (250–449 mg HC/g TOC). This is supported by their significant amounts of oil-liptinite macerals are present in these coals and coaly shale sediments and Py-GC (S₂) pyrograms with n-alkane/alkene doublets extending beyond nC₃₀. The shales are dominated by Type III kerogen (HI < 200 mg HC/g TOC), and are thus considered to be gas-prone.One-dimensional basin modelling was performed to analysis the hydrocarbon generation and expulsion history of the Mukalla source rocks in the Jiza-Qamar Basin based on the reconstruction of the burial/thermal maturity histories in order to improve our understanding of the of hydrocarbon generation potential of the Mukalla source rocks. Calibration of the model with measured vitrinite reflectance (Ro) and borehole temperature data indicates that the present-day heat flow in the Jiza-Qamar Basin varies from 45.0 mW/m² to 70.0 mW/m² and the paleo-heat flow increased from 80 Ma to 25 Ma, reached a peak heat-flow values of approximately 70.0 mW/m² at 25 Ma and then decreased exponentially from 25 Ma to present-day. The peak paleo-heat flow is explained by the Gulf of Aden and Red Sea Tertiary rifting during Oligocene-Middle Miocene, which has a considerable influence on the thermal maturity of the Mukalla source rocks. The source rocks of the Mukalla Formation are presently in a stage of oil and condensate generation with maturity from 0.50% to 1.10% Ro. Oil generation (0.5% Ro) in the Mukalla source rocks began from about 61 Ma to 54 Ma and the peak hydrocarbon generation (1.0% Ro) occurred approximately from 25 Ma to 20 Ma. The modelled hydrocarbon expulsion evolution suggested that the timing of hydrocarbon expulsion from the Mukalla source rocks began from 15 Ma to present-day.