Multi‐scale stratigraphic forward modelling of the Surat Basin for geological storage of CO2

Underground geological storage of CO₂ (GSC) requires a high level of subsurface understanding that is often hindered by a lack of data. This study demonstrates the use of stratigraphic forward modelling (SFM) in generating and characterising a static reservoir model using limited well data, with multiple potential applications within the GSC workflow. Sedsim SFM software was used to create a static model of the Surat Basin, including a high‐resolution nested model of the EPQ‐7 GSC tenement within the basin. Deposition and burial of the Jurassic Precipice Sandstone, Evergreen Formation and Hutton Sandstone were simulated. Modelling results show a close match with gamma‐ray well logs in the tenement area, and the model can be considered a credible model of the subsurface. The Sedsim‐predicted formation thicknesses and porosity and permeability distributions meet criteria set for GSC, suggesting that the EPQ‐7 tenement may be a prospective GSC location.     

The evolution of the Arabian lower crust and lithospheric mantle — Geochemical constraints from southern Syrian mafic and ultramafic xenoliths

Geochemical compositions of lower crustal and lithospheric mantle xenoliths found in alkali basaltic lavas from the Harrat Ash Shamah volcanic field in southern Syria place constraints on the formation of the Arabian–Nubian Shield in northern Arabia. Compositions of lower crustal granulites are compatible with a cumulate formation from mafic melts and indicate that they are not genetically related to their host rocks. Instead, their depletion in Nb relative to other incompatible elements points to an origin in a Neoproterozoic subduction zone as recorded by an average depleted mantle Sm–Nd model age of 630 Ma.Lithospheric spinel peridotites typically represent relatively low degree (< 10%) partial melting residues of spinel lherzolite with primitive mantle compositions as indicated by major and trace element modelling of clinopyroxene and spinel. The primary compositions of the xenoliths were subsequently altered by metasomatic reactions with low degree silicate melts and possibly carbonatites. Because host lavas lack these signatures any recent reaction of the lherzolites with their host magma can be ruled out. Sm–Nd data of clinopyroxene from Arabian lithospheric mantle lherzolites yield an average age of 640 Ma suggesting that the lithosphere was not replaced since its formation and supporting a common origin of the Arabian lower crustal and lithospheric mantle sections.The new data along with published Arabian mantle xenolith compositions are consistent with a model in which the lithospheric precursor was depleted oceanic lithosphere that was overprinted by metasomatic processes related to subduction and arc accretion during the generation of the Arabian–Nubian Shield. The less refractory nature of the northern Arabian lithosphere as indicated by higher Al, Na and lower Si and Mg contents of clinopyroxenes compared to the more depleted nature of the south Arabian lithospheric mantle, and the comparable low extent of melt extraction suggest that the northern Arabian lithosphere formed in a continental arc system, whereas the lithosphere in the southern part of Arabia appears to be of oceanic arc origin.     

Microfabrics of calcite ultramylonites as records of coaxial and non-coaxial deformation kinematics: Examples from the Rocher de l'Yret shear zone (Western Alps)

Microfabrics were analysed in calcite mylonites from the rim of the Pelvoux massif (Western Alps, France). WNW-directed emplacement of the internal Penninic units onto the Dauphinois domain produced intense deformation of an Eocene-age nummulitic limestone under lower anchizone metamorphic conditions (slightly below 300 °C). Two types of microfabrics developed primarily by dislocation creep accompanied by dynamic recrystallisation in the absence of twinning. Coaxial kinematics are inferred for samples exhibiting grain shape fabrics and textures with orthorhombic symmetry. Their texture (crystallographic preferred orientation, CPO) is characterised by two c-axis maxima, symmetrically oriented at 15° from the normal to the macroscopic foliation. Non-coaxial deformation is evident in samples with monoclinic shape fabrics and textures characterised by a single oblique c-axis maximum tilted with the sense of shear by about 15°. From the analysis of suitably oriented slip systems for the main texture components under given kinematics it is inferred that the orthorhombic textures, which developed in coaxial kinematics, favour activity of <10–11> and <02–21> slip along the f and r planes, respectively, with minor contributions of basal-<a> slip. In contrast, the monoclinic textures, which developed during simple shear, are most suited for duplex <a> slip along the basal plane. The transition between the dominating slip systems for the orthorhombic and monoclinic microfabrics is interpreted to be due to the effects of dynamic recrystallisation upon texture development. Since oblique c-axis maxima documented in the literature are most often rotated not with but against the shear sense, calcite textures alone should not be used as unequivocal shear sense indicators, but need to be complemented by microstructural criteria such as shape preferred orientations, grain size estimates and amount of twinning.     

On the variability of the Charnock constant and the functional dependence of the drag coefficient on wind speed: Part II-Observations

An analytical expression for the 10 m drag law in terms of the 10 m wind speed at the maximum in the 10 m drag coefficient, and the Charnock constant is presented, which is based on the results obtained from a model of the air-sea interface derived in Bye et al. (2010). This drag law is almost independent of wave age and over the mid-range of wind speeds (5?17 ms^?1) is very similar to the drag law based on observed data presented in Foreman and Emeis (2010). The linear fit of the observed data which incorporates a constant into the traditional definition of the drag coefficient is shown to arise to first-order as a consequence of the momentum exchange across the air-sea boundary layer brought about by wave generation and spray production which are explicitly represented in the theoretical model.     

Modelling the effects of temperature and food on individual growth and reproduction of Daphnia and their consequences on the population level

Individual based simulations of population dynamics require the availability of growth models with adequate complexity. For this purpose a simple-to-use model (non-linear multiple regression approach) is presented describing somatic growth and reproduction of Daphnia as a function of time, temperature and food quantity. The model showed a good agreement with published observations of somatic growth (r² = 0.954, n = 88) and egg production (r² = 0.898, n = 35). Temperature is the main determinant of initial somatic growth and food concentration is the main determinant of maximal body length and clutch size. An individual based simulation was used to demonstrate the simultaneous effects of food and temperature on the population level. Evidently, both temperature and food supply affected the population growth rate but at food concentrations above approximately 0.4 mg Cl⁻¹ Scenedesmus acutus temperature appeared as the main determinant of population growth.

Four simulation examples are given to show the wide applicability of the model: (1) analysis of the correlation between population birth rate and somatic growth rate, (2) contribution of egg development time and delayed somatic growth to temperature-effects on population growth, (3) comparison of population birth rate in simulations with constant vs. decreasing size at maturity with declining food concentrations and (4) costs of diel vertical migration. Due to its plausible behaviour over a broad range of temperature (2–20 °C) and food conditions (0.1–4 mg Cl⁻¹) the model can be used as a module for more detailed simulations of Daphnia population dynamics under realistic environmental conditions.     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.     

On CO<Subscript>2</Subscript> fluid flow and heat transfer behavior in the subsurface,following leakage from a geologic storage reservoir

Geologic storage of CO₂ is expected to produce plumes of large areal extent, and some leakage may occur along fractures, fault zones, or improperly plugged pre-existing wellbores. A review of physical and chemical processes accompanying leakage suggests a potential for self-enhancement. The numerical simulations presented here confirm this expectation, but reveal self-limiting features as well. It seems unlikely that CO₂ leakage could trigger a high-energy run-away discharge, a so-called “pneumatic eruption,” but present understanding is insufficient to rule out this possibility. The most promising avenue for increasing understanding of CO₂ leakage behavior is the study of natural analogues.     

Accuracy analysis, dem generation and validation using russian tk-350 stereo-images

TK-350 stereo-scenes of the Zonguldak testfield in the north-west of Turkey have been analysed. The imagery had a base-to-height ratio of 0·52 and covered an area of 200 km × 300 km, with each pixel representing 10 m on the ground. Control points digitised from 1:25 000 scale topographic maps were used in the test. A bundle orientation was executed using the University of Hanover program BLUH and PCI Geomatica OrthoEngine AE software packages. Tests revealed that TK-350 stereo-images can yield 3D geopositioning to an accuracy of about 10 m in planimetry and 17 m in height. A 40 m resolution digital elevation model (DEM) was generated by the PCI system and compared against a reference DEM, which was derived from digitised contour lines provided by 1:25 000 scale topographic maps. This comparison showed that accuracy depends mainly on the surface structure and the slope of the local terrain. Root mean square errors in height were found to be about 27 and 39 m outside and inside forested areas, respectively. The matched DEM demonstrated a systematic shift against the reference DEM visible as an asymmetric shift in the frequency distribution. This is perhaps caused by the presence of vegetation and buildings.