Seafloor seismic acquisition using autonomous underwater vehicles

Although narrow-azimuth towed-streamer data provide good image quality for structural interpretation, it is generally accepted that for wide-azimuth marine surveys seabed receivers deliver superior seismic reflection measurements and seismically derived reservoir attributes. However, seabed surveys are not widely used due to the higher acquisition costs when compared to streamer acquisition. In recent years, there have been significant engineering efforts to automate receiver deployment and retrieval in order to minimize the cost differential and conduct cost-efficient seabed receiver seismic surveys. These engineering efforts include industrially engineered nodes, nodes on a rope deployment schemes and even robotic nodes, which swim to and from the deployment location. This move to automation is inevitable, leading to robotization of seismic data acquisition for exploration and development activities in the oil and gas industry. We are developing a robotic-based technology, which utilizes autonomous underwater vehicles as seismic sensors without the need of using a remotely operated vehicle for deployment and retrieval. In this paper, we describe the autonomous underwater vehicle evolution throughout the project years from initial heavy and bulky nodes to fully autonomous light and flexible underwater receivers. Results obtained from two field pilot tests using different generations of autonomous underwater vehicles indicate that the seismic coupling, and navigation based on underwater acoustics are very reliable and robust.     

Mitigating the effects of sand dunes on seismic data from the Rub al Khali basin,Saudi Arabia

The Rub Al Khali region in Saudi Arabia is characterized by the presence of sand dunes separated by salt flats, also called Sabkhas. In general, the elevation of dunes in this region varies between 90 and 250 m above sea level. The presence of these sand dunes, along with the rapidly changing surface topography poses challenges for seismic data acquisition and processing. The high contrast in acoustic impedance between the dune base and the underlying formation often results in amplification of seismic waves that are recorded at stations located on the surface of sand dunes. Attempts to address the issue using conventional surface-consistent amplitude scaling methods without reducing these amplification effects generally fail, thus compromising the suitability of the processed data for amplitude versus offset analysis. In this study, we propose a new reference site technique to reduce the effects of sand dune amplification, enabling the production of data sets that are suitable for amplitude versus offset processing. The proposed technique uses a deterministic approach to derive surface-consistent, frequency-dependent de-amplification functions for shots and stations located on the dunes only. Two-dimensional synthetic and field data examples show that the technique significantly reduces the effects of sand dune amplification.     

Complex origin for the south-western Zamboanga metamorphic basement complex, Western Mindanao, Philippines

Abstract Amphibolites unconformably overlain by a metasedimentary sequence of quartz-muscovite-feldspar-kyanite schists, metagraywackes and epidote-bearing amphibolites occur in the northern portion of the south-western Zamboanga metamorphic basement complex, western Mindanao. These amphibolites (here identified as the Mount Dansalan amphibolites) display relict magmatic textures inherited from cumulate gabbro protoliths. Bulk-rock major and trace-element data are consistent with this hypothesis. Together with the chemistry of relict igneous clinopyroxenes, they indicate a magmatic arc-related signature for the gabbro protoliths. Geochemical data allow us to identify various sources for the associated metasediments: the gabbro themselves for the metagraywackes and a continental basement for the quartz-muscovite-feldspar-kyanite schists. Both sources contributed to the genesis of the epidote-amphibolite metasediments. The compositions of the metamorphic mineral assemblages suggest that the rocks have undergone metamorphism at temperatures ranging from 550°C to 700°C and pressures probably in the range of 5–9 kbar. ⁴⁰K–⁴⁰Ar isotopic study of amphibole separates from the Mount Dansalan samples document a metamorphic event dated at 24.6 ± 1.4, 22.2 ± 1.4 and 21.2 ± 1.2 Ma. Our results are in agreement with plate tectonic models which describe the south-western Zamboanga metamorphic basement as a continental terrane. However, its evolution was not as simple as it was usually considered. In particular the basement incorporated slivers of magmatic arc crust, which cannot be unambiguously related to any of the Tertiary arcs documented in the area.     

Uptake and tissue distribution of C4-C7 alkylphenols in Atlantic cod (Gadus morhua): relevance for biomonitoring of produced water discharges from oil production

The sensitivity of different tissues for assessment of chronic low-dose environmental exposure of fish to alkylphenols (APs) was investigated. We exposed Atlantic cod (Gadus morhua) in the laboratory to tritium labelled 4-tert-butylphenol, 4n-pentylphenol, 4n-hexylphenol, and 4n-heptylphenol via seawater (8 ng/l) and via contaminated feed (5 μg/kg fish per day). Measurements of different fish tissues during eight days of exposure and eight subsequent days of recovery revealed that APs administered via spiked seawater were readily taken up whereas the uptake was far less efficient when APs were administered in spiked feed. AP residues were mainly located in the bile fluid whereas the concentrations in liver were very low, indicating a rapid excretion and the liver-bile axis to be the major route of elimination. The biological half-life of APs in the exposed cod was short, between 10 and 20 h. Our study shows that in connection with biomonitoring of AP exposure in fish, assessment of AP metabolites in bile fluid is a more sensitive tool than detection of parent AP levels in liver or other internal tissues.     

Alteration geochemistry of the Nopal I uranium deposit (Sierra Peña Blanca,Mexico), a natural analogue for a radioactive waste repository in volcanic tuffs

Natural analogues provide an approach to characterize and test the long‐term modelling of a repository performance. This article presents geochemical information about the alteration conditions of the Nopal I uranium deposit, Mexico, an analogue for the proposed Yucca Mountain radioactive waste repository. Mineralization and hydrothermal alteration of volcanic tuffs are contemporaneous, according to petrographic observations. Trace element geochemistry (U, Th, REE) provides evidence for local mobilization of uranium under oxidizing conditions and further precipitation under reducing conditions. O‐ and H‐isotope geochemistry of kaolinite, smectite, opal and calcite suggests that argillic alteration proceeded at shallow depth with meteoric water at 25–75 °C, a low‐temperature context, unusual for volcanic‐hosted uranium deposits. This temperature range is compatible with some post‐closure evolution models of the proposed Yucca Mountain repository.     

Stratigraphy of the Cenomanian–Turonian Oceanic Anoxic Event OAE2 in shallow shelf sequences of NE Egypt

Two shallow water late Cenomanian to early Turonian sequences of NE Egypt have been investigated to evaluate the response to OAE2. Age control based on calcareous nannoplankton, planktic foraminifera and ammonite biostratigraphies integrated with δ¹³C stratigraphy is relatively good despite low diversity and sporadic occurrences. Planktic and benthic foraminiferal faunas are characterized by dysoxic, brackish and mesotrophic conditions, as indicated by low species diversity, low oxygen and low salinity tolerant planktic and benthic species, along with oyster-rich limestone layers. In these subtidal to inner neritic environments the OAE2 δ¹³C excursion appears comparable and coeval to that of open marine environments. However, in contrast to open marine environments where anoxic conditions begin after the first δ¹³C peak and end at or near the Cenomanian–Turonian boundary, in shallow coastal environments anoxic conditions do not appear until the early Turonian. This delay in anoxia appears to be related to the sea-level transgression that reached its maximum in the early Turonian, as observed in shallow water sections from Egypt to Morocco.     

How does micropetrography help us to understand the permeability and poromechanical behaviour of a rock?

Many poromechanical studies of rocks give little regard to the micropetrographic observation of their constituents. The aim of this study is to show how such a detailed study can help understand and provide more precision on the poromechanical behaviour of the material. The rock studied here is an oolitic limestone. Micropetrography reveals the nature of the constituents (calcite, mostly in the form of fine‐grained micrite), their abundance, their organization (micrite organized in ooids linked by a calcitic cement, small‐sized highly connected isotropic porosity, scarce unconnected porosity) and the structure of the rock (global isotropy, small heterogeneities due to sedimentological phenomena, no microcracks). Micropetrography completes the information obtained by poromechanics in that it allows the observation of elements deduced from macroscopical tests. The aim of this paper is to show the effect of microporosity on permeability and flow in a rock where microporous grains are in contact with each other.     

Contractivity-preserving explicit multistep Hermite–Obrechkoff series differential equation solvers

New optimal, contractivity-preserving (CP), explicit, d-derivative, k-step Hermite–Obrechkoff series methods of order p up to \(p=20\), denoted by CP HO(d, k, p), with nonnegative coefficients are constructed. These methods are used to solve nonstiff first-order initial value problems \(y'=f(t,y)\), \(y(t_0)=y_0\). The upper bound \(p_u\) of order p of HO(d, k, p) can reach, approximately, as high as 2.4 times the number of derivatives d. The stability regions of HO(d, k, p) have generally a good shape and grow with decreasing \(p-d\). We, first, note that three selected CP HO methods: 4-derivative 7-step HO of order 13, denoted by HO(4, 7, 13), 5-derivative 6-step HO of order 13, denoted by HO(5, 6, 13), and 9-derivative 2-step HO of order 13, denoted by CMDAHO(13) compare favorably with Adams–Cowell of order 13, denoted by AC(13), in solving standard N-body problems over an interval of 1000 periods on the basis of the relative error of energy as a function of the CPU time. Next, the three HO methods compare positively with AC(13) in solving standard N-body problems on the basis of the growth of relative positional error and relative energy error over 10, 000 periods of integration. Finally, these three methods compare also well with P-stable methods of Cash and Franco et al. on some quasi periodic, second-order linear and nonlinear problems. The coefficients of selected HO methods are listed in the appendix.     

Regional patterns of observed sea level change: insights from a 1/4° global ocean/sea-ice hindcast

A global eddy-admitting ocean/sea-ice simulation driven over 1958–2004 by daily atmospheric forcing is used to evaluate spatial patterns of sea level change between 1993 and 2001. In the present study, no data assimilation is performed. The model is based on the Nucleus for European Models of the Ocean code at the 1/4° resolution, and the simulation was performed without data assimilation by the DRAKKAR project. We show that this simulation correctly reproduces the observed regional sea level trend patterns computed using satellite altimetry data over 1993–2001. Generally, we find that regional sea level change is best simulated in the tropical band and northern oceans, whereas the Southern Ocean is poorly simulated. We examine the respective contributions of steric and bottom pressure changes to the total regional sea level changes. For the steric component, we analyze separately the contributions of temperature and salinity changes as well as upper and lower ocean contributions. Generally, the model results show that most regional sea level changes arise from temperature changes in the upper 750 m of the ocean. However, contributions of salinity changes and deep steric changes can be locally important. We also propose a map of ocean bottom pressure changes. Finally, we assess the robustness of such a model by comparing this simulation with a second simulation performed by MERCATOR-Ocean based on the same core model, but differing by its short length of integration (1992–2001) and its surface forcing data set. The long simulation presents better performance over 1993–2001 than the short simulation, especially in the Southern Ocean where a long adjustment time seems to be needed. In memory of my little brother Jean-Eudes, whose thirst for science filled out the rich discussions we had about my investigations and his job as user-service provider for MERCATOR-Ocean.     

Lagged rejuvenation of groundwater indicates internal flow structures and hydrological connectivity

Large proportions of rainwater and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time-varying processes. Chlorofluorocarbon (CFC)-based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km² forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from nine different wells with depths of 2–18 m close to the stream network. Immediately below the water table, CFC-based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age-depth-relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity.