Palaeo-carbonate seep structures above an oil reservoir, Gryphon Field, Tertiary, North Sea

Petrographic and geochemical analyses performed on a North Sea core from the Gryphon Field reveal the presence of palaeo-degassing features surrounded by injected sandstones in the Eocene interval. The injected sandstones are oil-stained and poorly cemented by carbonate and quartz. ¹⁸O isotope analyses indicate that carbonate cementation occurred during shallow burial (likely less than about 300 m). Depleted ¹³C (around –30 V-PDB) carbonate cement suggests that bicarbonate was derived from the microbial oxidation of oil and gas. Late quartz overgrowths enclose oil present in the injected units. The tubular degassing conduits are composed of zoned cements and have ¹⁸O and ¹³C isotope values similar to the injected sandstones, indicating that oil and gas seepage induced the precipitation of authigenic carbonate in the shallow subsurface. Oil inclusions in inter- and intra-crystal cement sites in both injected sandstones and degassing conduits indicate that oil seepage was an ongoing feature at shallow burial. A proposed model involves oil and gas seepage and the formation of the degassing conduits, followed by a sand injection phase. It seems likely that oil and gas continued to leak towards the seabed by exploiting the network of permeable injected sandstones and the horizons of porous degassing features.     

A review of eastern tropical Pacific oceanography: Summary

The collection of articles in this volume reviewing eastern tropical Pacific oceanography is briefly summarized, and updated references are given. The region is an unusual biological environment as a consequence of physical characteristics and patterns of forcing – including a strong and shallow thermocline, the ITCZ and coastal wind jets, equatorial upwelling, the Costa Rica Dome, eastern boundary and equatorial current systems, low iron input, inadequate ventilation of subthermocline waters, and dominance of ENSO-scale temporal variability. Remaining unanswered questions are presented.     

The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

Parameter sensitivity of three Kalman filter schemes for assimilation of water levels in shelf sea models

In applications of data assimilation algorithms, a number of poorly known assimilation parameters usually need to be specified. Hence, the documented success of data assimilation methodologies must rely on a moderate sensitivity to these parameters. This contribution presents a parameter sensitivity study of three well known Kalman filter approaches for the assimilation of water levels in a three dimensional hydrodynamic modelling system. The filters considered are the ensemble Kalman filter (EnKF), the reduced rank square root Kalman filter (RRSQRT) and the steady Kalman filter. A sensitivity analysis of key parameters in the schemes is undertaken for a setup in an idealised bay. The sensitivity of the resulting root mean square error (RMSE) is shown to be low to moderate. Hence the schemes are robust within an acceptable range and their application even with misspecified parameters is to be encouraged in this perspective. However, the predicted uncertainty of the assimilation results are sensitive to the parameters and hence must be applied with care. The sensitivity study further demonstrates the effectiveness of the steady Kalman filter in the given system as well as the great impact of assimilating even very few measurements.     

Trends in the solubility of iron,aluminium, manganese and phosphorus in aerosol collected over the Atlantic Ocean

The solubility of iron, aluminium, manganese and phosphorus has been determined in aerosol samples collected between 49°N and 52°S during three cruises conducted in the Atlantic Ocean as part of the European Union funded IRONAGES programme. Solubilities (defined at pH 4.7) determined for Fe and Al in samples of Saharan dust were significantly lower (medians 1.7% and 3.0%, respectively) than the solubilities of these metals in aerosols from other source regions (whole dataset medians 5.2% and 9.0%, respectively). Mn solubility also varied with aerosol source, but the median solubility of Mn in Saharan dust was very similar to the median for the dataset as a whole (55% and 56%, respectively). The observed solubility of aerosol P was ∼ 32%, with P solubility in Saharan aerosol perhaps as low as 10%. Laboratory studies have indicated that aerosol Fe solubility is enhanced by acid processing. No relationship could be found between Fe solubility and the concentrations of acid species (non-seasalt SO₄²⁻, NO₃⁻) nor the net acidity of the aerosol, so we are unable to confirm that this process is significant in the atmosphere. In terms of the supply of soluble Fe to oceanic ecosystems on a global scale, the observed higher solubility for Fe in non-Saharan aerosols is probably not significant because the Sahara is easily the dominant source of Fe to the Atlantic. On a smaller scale however, higher solubility for aerosol Fe may alter our understanding of Fe cycling in regions such as the remote Southern Ocean.     

Stress Sharing Behavior and Its Mechanism during Consolidation Process in Composite Ground Improved by Sand Compaction Piles with Low Replacement Area Ratio

In order to accurately design a sand compaction pile (SCP) with low replacement area ratio, it is important to understand the mechanical interaction between the sand pile and clay ground and its mechanism during consolidation process in composite ground. In this article, therefore, a series of numerical analyses on composite ground improved by SCP with low replacement area ratio were carried out. The applicability of numerical analyses, in which an elasto-viscoplastic consolidation finite element method was applied, were confirmed by comparing the results obtained from a series of laboratory model tests with the composite ground improved by SCP. Through the results of the numerical analyses, mechanical behavior of the sand pile and clay in composite ground during consolidation is elucidated, together with a stress sharing mechanism between sand pile and clay.     

Large aperture seismic imaging at a convergent margin: Techniques and results from the Costa Rica seismogenic zone

In March and April 1995 a cooperative German, Costa Rican, and United States research team recorded onshore-offshore seismic data sets along the Pacific margin of Costa Rica using the R/V Ewing. Off the Nicoya Peninsula we used a linear array of ocean bottom seismometers and hydrophones (OBS/H) with onshore seismometers extending across much of the isthmus. In the central area we deployed an OBS/H areal array consisting of 30 instruments over a 9 km by 35-km area and had land stations on the Nicoya Peninsula adjacent to this marine array and also extending northeast on the main Costa Rican landmass. Our goal in these experiments was to determine the crustal velocity structure along different portions of this convergent margin and to use the dense instrument deployments to create migrated reflection images of the plate boundary zone and the subducting Cocos Plate. Our specific goal in the central area was to determine whether a subducted seamount is present at the location of the 1990, M 7 earthquake off the Nicoya Peninsula and can thus be linked to its nucleation. Subsequently we have processed the data to improve reflection signals, used the data to calculate crustal velocity models, and developed several wide-aperture migration techniques, based on a Kirchhoff algorithm, to produce reflection images. Along the northern transect we used the ocean bottom data to construct a detailed crustal velocity model, but reflections from the plate boundary and top and bottom of the subducting Cocos plate are difficult to identify and have so far produced poor images. In contrast, the land stations along this same transect recorded clear reflections from the top of the subducting plate or plate boundary, within the seismogenic zone, and we have constructed a clear image from this reflector beneath the Nicoya shelf. Data from the 3-D seismic experiment suffer from high-amplitude, coherent noise (arrivals other than reflections), and we have tried many techniques to enhance the signal to noise ratio of reflected arrivals. Due to the noise, an apparent lack of strong reflections from the plate boundary zone, and probable structural complexity, the resulting 3-D images only poorly resolve the top of the subducting Cocos Plate. The images are not able to provide compelling evidence of whether there is a subducting seamount at the 1990 earthquake hypocenter. Our results do show that OBS surveys are capable of creating images of the plate boundary zone and the subducting plate well into the seismogenic zone if coherent reflections are recorded at 1.8 km instrument spacing (2-D) and 5 km inline by 1 km crossline spacing for 3-D acquisition. However, due to typical high amplitude coherent noise, imaging results may be poorer than expected, especially in unfavorable geologic settings such as our 3-D survey area. More effective noise reduction in acquisition, possibly with the use of vertical hydrophone arrays, and in processing, with advanced multiple removal and possibly depth filtering, is required to achieve the desired detailed images of the seismogenic plate boundary zone.