Mean and eddy motion in the Skagerrak/northern North Sea: insight from a numerical model

Insight regarding the mean and eddy motion in the Skagerrak/northern North Sea area is gained through an analysis of model-simulated currents, hydrography, kinetic energy and relative vorticity for the 2 years 2000 and 2001. In this a -coordinate ocean model is used. Since the tidal currents are generally strong in the area, care is exercised to distinguish the mesoscale (eddy) motion from higher-frequency motion such as tides, before computing the mean and eddy kinetic energy. The model-simulated response is first compared with available knowledge of the circulation in the area, and when available, also with sea-surface temperature obtained from satellite imagery. It is concluded that the model appears to faithfully reproduce most of what is known, in particularly the upper mixed layer circulation. An analysis of the mean and eddy kinetic energy reveals that many of the mesoscale structures found in the area are recurrent. This is particularly true for the structures off the southern tip of Norway. Also in general, areas of strong mean and eddy kinetic energy are co-located. The exception is the area off the southern tip of Norway, where the eddy kinetic energy is much larger than its mean counterpart. An analysis of the relative vorticity reveals that the variability found is due to the occurrence of recurrent anticyclonic eddies. It is hypothesized that these eddies are generated due to an offshore veering of the Norwegian coastal current (NCC) as it reaches the eastern end of the Norwegian Trench plateau. Here it becomes a free jet, which is then vulnerable to either barotropic instability caused by the horizontal shear in the jet-like structure of the NCC at this point, or a baroclinic (frontal) instability. The latter may come into play when the NCC veers offshore and its relatively fresh water meets the inflowing saline water of Atlantic origin, a frontogenesis that may become strong enough for cyclogenesis to take place. Due to the depth-independent nature of the model-generated eddies, the barotropic instability is the most likely candidate. It remains to resolve the reason for the offshore veering of the NCC. The most likely candidate mechanisms are vortex squeezing or simply that the coastline curvature is large enough for the NCC to separate from the coast in a hydraulic sense.Responsible Editor: Phil Dyke     

The age as a diagnostic of the dynamics of marine ecosystem models

The constituent-oriented age theory (CAT) worked out by Delhez et al. (1999) is a flexible tool that can be applied to diagnose complex models. It is shown here how this can be used to quantify the pace at which an ecosystem model works. At the cost of the introduction of one additional evolution equation for each compartment of the ecosystem model, the mean age of the biological material forming these compartments can be computed. The information obtained in this way complements the information provided by the concentration data; while the latter measures the standing stocks, the former provides an integrated assessment of the interaction rates and matter fluxes. The benefits of the method are demonstrated with a simple Lotka–Volterra system and a one-dimensional vertical model of the nitrogen cycle in the Ligurian Sea. The theory can be used to study the biological compartments individually or the ecosystem as a whole. In particular, the age is a valuable tool to quantify the overall cycling rate of nitrogen in the food web.Responsible Editor: Phil Dyke     

Spectra of poorly resolved ocean data

Standard temperature and pressure sensors on Aanderaa RCM8 current meters have a resolution of 0.024 °C and 0.6 bar, each equal to 1 digital number (value) over a range of 1024. It is shown that an 11-month deep-ocean temperature record using only four values can contain useful spectral information on internal wave motions. This is partially due to the modulation of high-frequency data by non-zero low-frequency (subinertial) variations. This result follows from the comparison of this record with artificial three- and four-value data constructed from temperature records observed in stronger stratified waters nearby. These artificial records show main features of the internal wave band similar to those observed in the original data spectra. Peaks at tidal harmonic frequencies and enhancements at sum-tidal-inertial interaction frequencies are preserved in the artificial data, but overall noise level (and thus the continuum spectral slope) is enhanced with respect to the properly resolved records (using 15 and 100 values). As a demonstration of the stable accuracy of the temperature sensors, the poorly resolved records provided an estimate of mean stratification to within 5% of the estimate using Seabird CTD data.Responsible Editor: Hans Burchard     

Seismic ground motion relationship in southern China based on stochastic finite-fault model

The characteristics of seismic ground motions in southern China are difficult to determine statistically due to a lack of strong ground motion data. In this study, a stochastic finite-fault ground motion model was adopted to simulate the seismic ground motions at bedrock for southern China, based on parameters derived from small and medium earthquakes that have occurred in the region. From these, the response spectra was estimated. A set of ground motion attenuation relations hipswas then developed based on simulated peak ground motions and response spectral parameters through regression, which would be applicable for use in engineering practice. Through comparisons, it was demonstrated that the proposed ground motion relationships are generally consistent with those obtained from other reported ground motion attenuation models for southern China.     

Complex Site Effects in Thessaloniki (Greece): I. Soil Structure and Comparison of Observations with 1D Analysis

This paper presents a 2D model of the geological structure of Thessaloniki city and results of empirical and theoretical approaches for the evaluation of site response due to complex site effects. The construction of the 2D model is based on the available geophysical and geotechnical data in terms of the most important parameters needed to model site response. The well-known subsoil structure, despite the existence of some local uncertainties, gave the possibility to correlate the geometry and the dynamic properties of the 2D model with the results of site response determined from the analysis of one event in frequency and time domains and 1D numerical modelling. The study of site response shows the effect of the lateral variations on ground motion and suggests that the contribution of locally generated surface waves to the resonant peak may be important. In this case history, the limitations of the 1D approximation to simulate ground motion under complex soil conditions in both frequency and time domains are also shown. This paper lays the ground for a companion article dealing with 2D site effects.     

Geochemical and isotopic characteristics of the Cretaceous Orikabe Plutonic Complex,Kitakami Mountains,Japan: magmatic evolution in a zoned pluton and significance of a subduction-related mafic parental magma

The Orikabe Plutonic Complex, northeast Japan, is a zoned pluton and one of the Cretaceous intrusions in the Circum-Pacific area. In the Main body, K-rich calc-alkaline rocks composed of marginal gabbro and a large amount of monzodiorite–quartz monzonite–monzogranite are intruded successively by innermost calc-alkaline rocks of granodiorite. The gabbro and monzodiorite–monzogranite have a continuous chemical variation, while the granodiorite has lower concentrations of K, Rb, Y, Zr, Nb and F at the same SiO₂ content. The gabbro and monzodiorite–quartz monzonite have a Rb-Sr whole-rock age of 119±12 Ma with an initial ⁸⁷Sr/⁸⁶Sr ratio of 0.70392±0.00007. The initial ⁸⁷Sr/⁸⁶Sr ratio of the innermost granodiorite is estimated to be about 0.7042. The ¹⁸O values of fresh rocks range from +6.7 to +8.3, indicating a positive correlation with SiO₂ contents. The K-rich calc-alkaline rocks were derived through fractional crystallization from a mafic parental magma with a slightly high ¹⁸O value, implying a major contribution of a sub-arc mantle at a continental margin. Trace element modeling indicates that the source could have been a fertile lherzolite enriched in LILE and depleted in HFSE. The innermost granodiorite was the differentiation product of a distinct parental magma, suggesting the involvement of a small amount of crustal component in the source and partial melting under a more hydrous condition.Editorial responsibility: J. Hoefs     

Carbonatite diversity in the Central Andes: the Ayopaya alkaline province,Bolivia

The Ayopaya province in the eastern Andes of Bolivia, 100 km NW of Cochabamba, hosts a Cretaceous alkaline rock series within a Palaeozoic sedimentary sequence. The alkaline rock association comprises nepheline-syenitic/foyaitic to ijolitic intrusions, carbonatite, kimberlite, melilititic, nephelinitic to basanitic dykes and diatremes, and a variety of alkaline dykes. The carbonatites display a wide petrographic and geochemical spectrum. The Cerro Sapo area hosts a small calciocarbonatite intrusion and a multitude of ferrocarbonatitic dykes and lenses in association with a nepheline-syenitic stock. The stock is crosscut by a spectacular REE-Sr-Th-rich sodalite-ankerite-baryte dyke system. The nearby Chiaracke complex represents a magnesiocarbonatite intrusion with no evidence for a relationship to igneous silicate rocks. The magnesiocarbonatite ( REE up to 1.3 wt%) shows strong HREE depletion, i.e. unusually high La/Yb ratios (520–1,500). Calciocarbonatites ( REE up to 0.5 wt%) have a flatter REE distribution pattern (La/Yb 95–160) and higher Nb and Zr contents. The sodalite-ankerite-baryte dyke system shows geochemical enrichment features, particularly in Na, Ba, Cl, Sr, REE, which are similar to the unusual natrocarbonatitic lavas of the recent volcano of Oldoinyo Lengai, Tanzania. The Cerro Sapo complex may be regarded as an intrusive equivalent of natrocarbonatitic volcanism, and provides an example for carbonatite genesis by late-stage crystal fractionation and liquid immiscibility. The magnesiocarbonatite intrusion of Chiaracke, on the other hand, appears to result from a primary carbonatitic mantle melt. Deep seated mantle magmatism/metasomatism is also expressed by the occurrence of a kimberlite dyke. Neodymium and strontium isotope data (_Nd 1.4–5.4, ⁸⁷Sr/⁸⁶ Sr<Bulk Earth) indicate a depleted mantle source for the alkaline magmatism. The magmatism of the Ayopaya region is attributed to failed rifting of western South America during the Mesozoic and represents the only occurrence of carbonatite and kimberlite rocks in the Andes.     

Distribution and petrogenetic behaviour of trace elements in granitic pegmatite quartz from South Norway

The present study documents that the trace-element distribution in granitic quartz is highly sensitive to CAFC processes in granitic melts. Igneous quartz efficiently records both the origin and the evolution of the granitic pegmatites. Aluminium, P, Li, Ti, Ge and Na in that order of abundance, comprises >95% of the trace elements. Most samples feature >1 ppm of any of these elements. The remnant 5% includes K, Fe, Be, B, Ba and Sr whereas the other elements are present at concentrations lower than the detection limit. Potassium, Fe, Be and Ti are relatively compatible hence obtain the highest concentrations in early formed quartz. Phosphorous, Ge, Li and Al are relatively incompatible and generally obtain the highest concentrations in quartz that formed at lower temperatures from more evolved granitic melts. The Ge/Ti, the Ge/Be, the P/Ge and the P/Be ratios of quartz are strongly sensitive to the origin and evolution of the granitic melts and similarly the Rb/Sr and the Rb/K ratios of K-feldspars may be utilised in petrogenetic interpretations. However, the quartz trace element ratios are better at distinguishing similarities and differences in the origin and evolution of granitic melts. After evaluating the different trace element ratios, the Ge/Ti ratio appears to be most robust during subsolidus processes in the igneous systems, hence probably should be the preferred ratio for analysing and understanding petrogenetic processes in granitic igneous rocks.Editorial responsibility: J. Hoefs     

A kinetic model of metamorphism: an application to siliceous dolomites

We use a kinetic model of a metamorphic system to study the effect of competing rates of reaction, fluid injection, and heating on the evolution of the reaction pathway in temperature/composition space at constant pressure. We show that for rocks in contact with mixed volatile (e.g., CO₂-H₂O) fluids the reaction path may be quite different from what is expected from equilibrium-based petrologic models. Equilibrium-based models, used to understand the development of rock systems undergoing mineral reactions during a metamorphic event, rely on the Gibbs phase rule and only consider stable phases. For constant pressure, the temperature-composition paths follow univariant curves and significant reactions may occur at invariant points. By contrast, the more general kinetic treatment is not constrained by equilibrium, although with the proper competing rates equilibrium is a possible endmember of the kinetic approach. The deviation from equilibrium depends on the competing rates of reaction, heating, and fluid injection. A key element required by the kinetic approach is the inclusion of metastable reactions in the formulation, whereas such reactions are irrelevant for equilibrium-based models. Metastable reactions are often involved in a complex interplay with common prograde stable metamorphic reactions. We present model results for the well-studied CaO-MgO-SiO₂-CO₂-H₂O (CMS) system to show how the system evolves under kinetic control. Our simulations and discussion focus on the behavior of the CMS system under a number of closed and open system conditions. Special attention is paid to closed system behavior in the vicinity of the (first) isobaric invariant point (with Dol, Qtz, Tlc, Cal, and Tr). Also, for open systems with massive fluid infiltration we consider heating rates varying from contact to regional metamorphic conditions. For some geologically reasonable rates of reactions, heating, and fluid injection, our results demonstrate that equilibrium conditions may be significantly overstepped in metamorphic systems. We used overall mineral reactions in this model with rates based on experimental results. Future models could rely on more fundamental dissolution and precipitation reactions. Such an extension would require additional kinetic rate data, as well as mineral solubilities in mixed volatile fluids.Editorial responsibility: J. Hoefs