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.     

The Agulhas Ridge, South Atlantic: The Peculiar Structure of a Fracture Zone

The Agulhas Ridge is a prominent topographic feature that parallels the Agulhas-Falkland Fracture Zone (AFFZ). Seismic reflection and wide angle/refraction data have led to the classification of this feature as a transverse ridge. Changes in spreading rate and direction associated with ridge jumps, combined with asymmetric spreading within the Agulhas Basin, modified the stress field across the fracture zone. Moreover, passing the Agulhas Ridge’s location between 80 and 69 Ma, the Bouvet and Shona Hotspots may have supplied excess material to this part of the AFFZ thus altering the ridge’s structure. The low crustal velocities and overthickened crust of the northern Agulhas Ridge segment indicate a possible continental affinity that suggests it may be formed by a small continental sliver, which was severed off the Maurice Ewing Bank during the opening of the South Atlantic. In early Oligocene times the Agulhas Ridge was tectono-magmatically reactivated, as documented by the presence of basement highs disturbing and disrupting the sedimentary column in the Cape Basin. We consider the Discovery Hotspot, which distributes plume material southwards across the AAFZ, as a source for the magmatic material.     

Nocdap results: On the relationship of extreme waves and currents for design of offshore structures

Simultaneous wind, wave, and current data during 21 storms spanning four winters at Tromsøflaket (230 m depth) were analysed to determine joint probabilities of occurrence. Waves were measured with a Waverider Buoy, winds with a recording anemometer onboard a vessel and currents at up to five depths with Aandera RCM-4 current meters. Measured currents were filtered to separate tidal currents from residual currents. In most of the statistical analyses, the actual current profiles were replaced with a simplified “equivalent” profile (constant above 50 m depth) whose magnitude was chosen to give the same drag load on a single pile as the actual profile when both profiles were combined with wave orbital velocities. The data suggest that the equivalent current is weakly correlated with the wave height. Within the range of the data analysed, an in-line equivalent residual current of 30 cm sec⁻¹ is adequate for design drag force calculations.The data presented here are used to illustrate a procedure for data analysis and are not suitable for use as design criteria.     

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.     

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.     

U–Pb ages of Neoarchean granitoids from the Black Hills, South Dakota, USA: implications for crustal evolution in the Archean Wyoming province

Zircons in basement rocks from the eastern Wyoming province (Black Hills, South Dakota, USA) have been analyzed by ion microprobe (SHRIMP) in order to determine precise ages of Archean tectonomagmatic events. In the northern Black Hills (NBH) near Nemo, Phanerozoic and Proterozoic (meta)sedimentary rocks are nonconformably underlain by Archean biotite–feldspar gneiss (BFG) and Little Elk gneissic granite (LEG), both of which intrude older schists. The Archean granitoid gneisses exhibit a pervasive NW–SE-trending fabric, whereas an earlier NE–SW-trending fabric occurs sporadically only in the BFG, which is intruded by the somewhat younger LEG. Zircon crystals obtained from the LEG and BFG exhibit double terminations, oscillatory zoning, and Th/U ratios of 0.6±0.3—thereby confirming a magmatic origin for both lithologies. In situ analysis of the most U–Pb concordant domains yields equivalent ²⁰⁷Pb/²⁰⁶Pb ages (upper intercept, U–Pb concordia) of 2559±6 and 2563±6 Ma (both ±2σ) for the LEG and BFG, respectively, which constrains a late Neoarchean age for sequential pulses of magmatism in the NBH. Unzoned (in BSE) patches of 2560 Ma zircon commonly truncate coeval zonation in the same crystals with no change in Th/U ratio, suggesting that deuteric, fluid-assisted recrystallization accompanied post-magmatic cooling. A xenocrystic core of magmatic zircon observed in one LEG zircon yields a concordant age of 2894±6 Ma (±2σ). This xenocryst represents the oldest crustal material reported thus far in the Black Hills. Whether this older zircon originated as unmelted residue of 2900 Ma crust that potentially underlies the Black Hills or as detritus derived from 2900 Ma crustal sources in the Wyoming province cannot be discerned. In the southern Black Hills (SBH), the peraluminous granite at Bear Mountain (BMG) of previously unknown age intrudes biotite–plagioclase schist. Zircon crystals from the BMG are highly metamict and altered, but locally preserve small domains suitable for in situ analysis. A U–Pb concordia upper intercept age of 2596±11 Ma (±2σ) obtained for zircon confirms both the late Neoarchean magmatic age of the BMG and a minimum age for the schist it intrudes. Taken together, these data indicate that the Neoarchean basement granitoids were emplaced at 2590–2600 Ma (SBH) and 2560 Ma (NBH), most likely in response to subduction associated with plate convergence (final assembly of supercontinent Kenorland?). In contrast, thin rims present on some LEG–BFG zircons exhibit strong U–Pb discordance, high common Pb, and low Th/U ratios—suggesting growth or modification under hydrothermal conditions, as previously suggested for similar zircons from SE Wyoming. The LEG–BFG zircon rims yield a nominal upper intercept date of 1940–2180 Ma, which may represent a composite of multiple rifting events known to have affected the Nemo area between 2480 and 1960 Ma. Together, these observations confirm the existence of a Paleoproterozoic rift margin along the easternmost Wyoming craton. Moreover, the 2480–1960 Ma time frame inferred for rifting in the Black Hills (Nemo area) corresponds closely to a 2450–2100 Ma time frame previously inferred for the fragmentation of supercontinent Kenorland.     

CO2-H2O mixtures in the geological sequestration of CO2. II. Partitioning in chloride brines at 12-100°C and up to 600 bar

Correlations presented by Spycher et al. (2003) to compute the mutual solubilities of CO₂ and H₂O are extended to include the effect of chloride salts in the aqueous phase. This is accomplished by including, in the original formulation, activity coefficients for aqueous CO₂ derived from several literature sources, primarily for NaCl solutions. Best results are obtained when combining the solubility correlations of Spycher et al. (2003) with the activity coefficient formulation of Rumpf et al. (1994) and Duan and Sun (2003), which can be extended to chloride solutions other than NaCl. This approach allows computing mutual solubilities in a noniterative manner with an accuracy typically within experimental uncertainty for solutions up to 6 molal NaCl and 4 molal CaCl₂.     

Deciphering late Devonian–early Carboniferous P–T–t path of mylonitized garnet-mica schists from Prins Karls Forland,Svalbard

Quartz-in-garnet inclusion barometry integrated with trace element thermometry and calculated phase relations is applied to mylonitized schists of the Pinkie unit cropping out on the island of Prins Karls Forland, western part of the Svalbard Archipelago. This approach combines conventional and novel techniques and allows deciphering of the pressure–temperature (P–T) evolution of mylonitic rocks, for which the P–T conditions could not have been easily deciphered using traditional methods. The results obtained suggest that rocks of the Pinkie unit were metamorphosed under amphibolite facies conditions at 8–10 kbar and 560–630°C and mylonitized at ~500 to 550°C and 9–11 kbar. The P–T results are coupled with in-situ Th–U-total Pb monazite dating, which records amphibolite facies metamorphism at c. 359–355 Ma. This is the very first evidence of late Devonian–early Carboniferous metamorphism in Svalbard and it implies that the Ellesmerian Orogeny on Svalbard was associated with metamorphism up to amphibolite facies conditions. Thus, it can be concluded that the Ellesmerian collision between the Franklinian margin of Laurentia and Pearya and Svalbard caused not only commonly accepted brittle deformation and weak greenschist facies metamorphism, but also a burial and deformation of rock complexes at much greater depths at elevated temperatures.