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.     

Recent ventifact development on the central Oregon coast,western USA

The unusual location of ventifacts, on a boulder‐built jetty at the mouth of the Siuslaw River, Oregon coast, western USA, allows ventifact age and wind abrasion rates to be estimated with some precision. The jetty was built mainly between 1892–1901 and extended throughout the twentieth century. Consideration of historical shoreline position and the history of jetty construction and repair suggests the ventifacts have formed since about 1930. Morphologically the ventifacts are aligned south‐to‐north reflecting winter winds and sediment transport from the adjacent beach. Wind‐parallel grooves and ridges with sharp, sinuous crests are developed on inclined boulder surfaces on top of the jetty and reflect suspended sand transport in wind vortices. Deeply pitted surfaces on steep boulder surfaces nearest the beach reflect impact by saltating sand grains. Based on present wind regimes (1992–2000) from three regional weather stations, southerly winds above the sand transport threshold occur for 21·9–29·6 per cent of the time. Based on estimated depth of loss from boulder surfaces, wind abrasion rates are calculated to be on the order of 0·24–1·63 mm a^?1. This is the first well‐constrained field estimate of ventifact age and ventifaction rate from a modern coastal environment. Copyright © 2003 John Wiley & Sons, Ltd.     

An experimental study on the shallow-level migmatization of ferrogabbros from the Fuerteventura Basal Complex, Canary Islands

Spectacular shallow-level migmatization of ferrogabbroic rocks occurs in a metamorphic contact aureole of a gabbroic pluton of the Tierra Mala massif (TM) on Fuerteventura (Canary Islands). In order to improve our knowledge of the low pressure melting behavior of gabbroic rocks and to constrain the conditions of migmatization of the TM gabbros, we performed partial melting experiments on a natural ferrogabbro, which is assumed as protolith of the migmatites. The experiments were performed in an internally heated pressure vessel (IHPV) at 200 MPa, 930–1150 °C at relatively oxidizing conditions. Distinct amounts of water were added to the charge.

From 930 to 1000 °C, the observed experimental phases are plagioclase (An_60–70), clinopyroxene, amphibole (titanian magnesiohastingsites), two Fe–Ti oxides, and a basaltic, K-poor melt. Above 1000 °C, amphibole is no longer stable. The first melts are very rich in normative plagioclase (>70 wt.%). This indicates that at the beginning of partial melting plagioclase is the major phase which is consumed to produce melt. In the experiments, plagioclase is stable up to high temperatures (1060 °C) showing increasing An content with temperature. This is not compatible with the natural migmatites, in which An-rich plagioclase is absent in the melanosomes, while amphibole is stable. Our results show that the partial melting of the natural rocks cannot be regarded as an “in-situ” process that occurred in a closed system. Considerable amounts of alkalis probably transported by water-rich fluids, derived from the mafic pluton underplating the TM gabbro, were necessary to drive the melting reaction out of the stability range of plagioclase. A partial melting experiment with a migmatite gabbro showing typical “in-situ” textures as starting material supports this assumption.

Crystallization experiments performed at 1000 °C on a glass of the fused ferrogabbro with different water contents added to the charge show that generally high water activities could be achieved (crystallization of amphibole), independently of the bulk water content, even in a system with very low initial bulk water content (0.3 wt.%). Increasing water contents produce plagioclase richer in An, reduces the modal proportion of plagioclase in the crystallizing assemblage and extends the melt fraction. High melt fractions of >30 wt.% could only be observed in systems with high bulk water contents (>2 wt.%). This indicates that the migmatites were generated under water-rich conditions (probably water-saturated), since those migmatites, which are characterized as “in-situ” formations, show generally high amounts of leucosomes (>30 wt.%).     

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.     

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.     

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.     

Origin of complex zoning in olivine from diverse,diamondiferous kimberlites and tectonic settings: Ekati (Canada), Alto Paranaiba (Brazil) and Kaalvallei (South Africa)

Olivine in kimberlites can provide unique insights into magma petrogenesis, because it is the most abundant xenocrystic phase and a stable magmatic product over most of the liquid line of descent. In this study we examined the petrography and chemistry of olivine in kimberlites from different tectonic settings, including the Slave craton, Canada (Ekati: Grizzly, Koala), the Brasilia mobile belt (Limpeza-18, Tres Ranchos-04), and the Kaapvaal craton, South Africa (Kaalvallei: Samada, New Robinson). Olivine cores display a wide range of compositions (e.g., Mg# = 78–95). The similarity in olivine composition, resorption of core zones and inclusions of mantle-derived phases, indicates that most olivine cores originated from the disaggregation of mantle peridotites, including kimberlite-metasomatised lithologies (i.e. sheared lherzolites and megacrysts). Olivine rims typically show a restricted range of Mg#, with decreasing Ni and increasing Mn and Ca contents, a characteristic of kimberlitic olivine worldwide. The rims host inclusions of groundmass minerals, which implies crystallisation just before and/or during emplacement. There is a direct correlation between olivine rim composition and groundmass mineralogy, whereby high Mg/Fe rims are associated with carbonate-rich kimberlites, and lower Mg/Fe rims are correlated with increased phlogopite and Fe-bearing oxide mineral abundances. There are no differences in olivine composition between explosive (Grizzly) and hypabyssal (Koala) kimberlites. Olivine in kimberlites also displays transitional zones and less common internal zones, between cores and rims. The diffuse transitional zones exhibit intermediate compositions between cores and rims, attributed to partial re-equilibration of xenocrystic cores with the ascending kimberlite melt. In contrast, internal zones form discrete layers with resorbed margins and restricted Mg# values, but variable Ni, Mn and Ca concentrations, which indicates a discrete crystallization event from precursor kimberlite melts at mantle depths. Overall, olivine exhibits broadly analogous zoning in kimberlites worldwide. Variable compositions for individual zones relate to different parental melt compositions rather than variations in tectonic setting or emplacement mechanism.

     

The age of carbonatites,kimberlites and lamprophyres from southern West Greenland: recurrent alkaline magmatism during 2500 million years

Alkaline rocks were intruded into the Precambrian basement in southern West Greenland during at least five separate episodes. The Tupertalik carbonatite intrusion was emplaced around or before 2650 m.y.; it is now metamorphosed in granulite facies but is recognisable by its trace element content. Lamprophyre dykes were intruded at ca. 1800 m.y. and again at ca. 1200 m.y. At ca. 600 m.y. the Sarfartôq carbonatite intrusion and extensive kimberlite dykes were emplaced in the Holsteiborg-Strømfjord region. At 225-115 m.y. the Qaqarssuk carbonatite complex was emplaced in the Sukkertoppen region, and numerous lamprophyre and kimberlite dykes were emplaced between Fishenasset and Ivigtut. All these episodes are correlatable with contemporaneous alkaline activity in Canada and Scandinavia and substantiate the periodicity of carbonatite emplacement. The Mesozoic suite of intrusions is related to rifting prior to continental break-up and the formation of the North Atlantic Ocean. The Eocambrian suite of 600 m.y. is thought to have formed under similar conditions during continental break-up and the formation of the Iapetus ocean.