Dispersion of surface plumes in the Southern North Sea

Summary Continuous releases of fluorescent dye were made in the southern North Sea and the power law dependence of the plume width, represented by _y , against the diffusion time determined. The spreading was non-Fickian and could be represented in the form _y ² =B ² t whereB was a diffusion velocity of magnitude 1.4×10^–2 m/s. Such spreading was reproduced in both particle tracking and finite difference plume models by allowing the horizontal diffusivity,K _H, to depend linearly on diffusion time. The weakness of this method is that it is not clear how the diffusion parameters can be extrapolated to weather conditions that are different to those prevailing at the time of the experiments. However, comparable spreading rates were obtained by combining a Fickian diffusion model with an advective field that represented the near-surface current shears due to wind and waves. The resulting shear diffusion effect produced realistic simulation of the observed dispersion rates. An advantage of this approach is that it enables predictions to be made over a range of weather conditions provided that the wind and wave shears can be accurately parameterized.

---

Ausbreitung von Tracerwolken an der Oberfläche in der südlichen Nordsee

Zusammenfassung Fluoreszierender Farbstoff wurde in der südlichen Nordsee kontinuierlich freigesetzt und die Breite der Wolke, repräsentiert durch _y , in Abhägigkeit von der Ausbreitungszeit bestimmt. Die Ausbreitung folgte nicht dem Fickschen Ansatz und konnte in der Form _y ² =B ² t dargestellt werden, wobeiB eine Ausbreitungsgeschwindigkeit von der Größenordnung 1.4×10^–2 m/s darstellt. Eine derartige Ausbreitung wurde sowohl mit Modellen des Partikel-Tracking als auch mit Modellen der finiten Differenzen reproduziert, wobei die horizontale AusbreitungsgrößeK _H linear von der Diffusionszeit abhängig war. Die Schwäche dieser Methode ist, daß es nicht klar ist, wie die Diffusions-Parameter für andere Wetterbedingungen als zur Zeit des Experiments zu extrapolieren wären. Vergleichbare Ausbreitungsverhältnisse werden durch die Kombination eines Fickschen Diffusionsmodells mit einem advektiven Feld erreicht, welches die oberflächennahen durch Wind und Wellen verursachten Stromscherungen berücksichtigt. Der so resultierende Scher-Diffusions-Effekt erbrachte ein realistische Simulation der beobachteten Ausbreitung. Ein Vorteil dieses Vorgehens ist, daß es Vorhersagen über einen Bereich von Wetterbedingungen ermöglicht, vorausgesetzt, daß die durch Wind und Wellen verurschte Scherung zutreffend parameterisiert werden kann.

---

La dispersion de la panache de surface dans la partie sud de la Mer du Nord

Résumé A partir de déversements continus de produits fluorescents effectués dans la partie sud de la Mer du Nord, on a pu déterminer l'expression de la largeur du panache _y , par son carré en fonction du temps de diffusion. La dispersion n'était pas du type Fickian et était représentée sous la forme _y ² =B ² t aB était une vitesse de diffusion de grandeur égale 1.4×10^–2 m/s. Une telle dispersion peut être reproduite par une modélisation du panache à la fois par suivi de particules et aux différences finies en admettant que le coefficient de diffusion horizontaleK _H dépendait linéairement du temps. L'inconvénient de cette méthode est qu'il n'est pas facile de savoir comment les paramètres de diffusion peuvent être extrapolés à des conditions météorologiques différentes de celles prévalant au moment des expériences. Quoi qu'il en soit, des taux de dispersion comparables ont été obtenus en combinant un modèle de diffusion Fickian avec un champ d'advection qui représentait les cisaillements de courant de surface dus au vent et aux vagues. L'effet de diffusion du cisaillement résultant produisit une simulation réaliste des taux de dispersion observés. Un avantage de cette approche est de rendre possible des prédictions réalisées à travers un éventail de conditions météorologiques à condition que les cisaillements liés au vent et aux vagues puissent être paramétrisés avec précision.

---

     

http://www.sciencedirect.com/science/article/pii/S1674987114001601

Data from a migmatised metapelite raft enclosed within charnockite provide quantitative constraints on the pressure-temperature-time[P-T-t) evolution of the Nagercoil Block at the southernmost tip of peninsular India.An inferred peak metamorphic assemblage of garnet,K-feldspar.sillimanite,plagioclase,magnetite,ilmenite,spinel and melt is consistent with peak metamorphic pressures of 6-8 kbar and temperatures in excess of 900℃.Subsequent growth of cordierite and biotite record high-temperature retrograde decompression to around 5 kbar and 800 C.SHRIMP U-Pb dating of magmatic zircon cores suggests that the sedimentary protoliths were in part derived from felsic igneous rocks with Palaeoproterozoic crystallisation ages.New growth of metamorphic zircon on the rims of detrital grains constrains the onset of melt crystallisation,and the minimum age of the metamorphic peak,to around560 Ma.The data suggest two stages of monazite growth.The first generation of REE-enriched monazite grew during partial melting along the prograde path at around 570 Ma via the incongruent breakdown of apatite.Relatively REE-depleted rims,which have a pronounced negative europium anomaly,grew during melt crystallisation along the retrograde path at around 535 Ma.Our data show the rocks remained at suprasolidus temperatures for at least 35 million years and probably much longer,supporting a long-lived high-grade metamorphic history.The metamorphic conditions,timing and duration of the implied clockwise P-T-t path are similar to that previously established for other regions in peninsular India during the Ediacaran to Cambrian assembly of that part of the Gondwanan supercontinent.     

Isochron‐burial dating of glaciofluvial deposits: First results from the Swiss Alps

In this study, we use isochron‐burial dating to date the Swiss Deckenschotter, the oldest Quaternary deposits of the northern Alpine Foreland. Concentrations of cosmogenic ¹⁰Be and ²⁶Al in individual clasts from a single stratigraphic horizon can be used to calculate an isochron‐burial age based on an assumed initial ratio and the measured ²⁶Al/¹⁰Be ratio. We suggest that, owing to deep and repeated glacial erosion, the initial isochron ratio of glacial landscapes at the time of burial varies between 6.75 and 8.4. Analysis of 22 clasts of different lithology, shape, and size from one 0.5 m thick gravel bed at Siglistorf (Canton Aargau) indicates low nuclide concentrations: <20 000 ¹⁰Be atoms/g and <150 000 ²⁶Al atoms/g. Using an ²⁶Al/¹⁰Be ratio of 7.6 (arithmetical mean of 6.75 and 8.4), we calculate a mean isochron‐burial age of 1.5 ± 0.2 Ma. This age points to an average bedrock incision rate between 0.13 and 0.17 mm/a. Age data from the Irchel, Stadlerberg, and Siglistorf sites show that the Higher Swiss Deckenschotter was deposited between 2.5 and 1.3 Ma. Our results indicate that isochron‐burial dating can be successfully applied to glaciofluvial sediments despite very low cosmogenic nuclide concentrations. Copyright © 2017 John Wiley & Sons, Ltd.     

Assessing the contribution of off‐fault deformation to slip‐rate estimates within the Taupo Rift,New Zealand,using 3‐D ground‐penetrating radar surveying and trenching

We demonstrate that conventional palaeoseismic trenching and mapping techniques that do not account for the effects of off‐fault deformation can significantly underestimate a fault’s slip rate. Using combined interpretations of 3‐D ground‐penetrating radar (GPR) and palaeoseismic trench data, we show that drag folding and hangingwall and footwall horizontal‐axis rotations have accommodated up to 41% of total extension across a normal fault within the Taupo Rift, New Zealand, over the past 24.6 ± 1.0 cal. ka BP. Our results may explain why geologically determined fault‐slip rates for the central and southern Taupo Rift are anomalously low when compared with geodetic estimates. We suggest that a combination of GPR surveying and palaeoseismic trenching may help resolve differences between geodetically and geologically determined strain rates observed across active extensional regimes worldwide.     

Simultaneous X-ray and IR variability in the quasar 3C 273

From a combination of high-quality X-ray observations from the NASA Rossi X-ray Timing Explorer ( RXTE ) and IR observations from the UK Infrared Telescope (UKIRT) we show that the medium-energy X-ray (3–20 keV) and near-IR fluxes in the quasar 3C 273 are highly correlated. It is widely believed that the X-ray emission in quasars like 3C 273 arises from Compton scattering of low-energy seed photons, and our observations provide the first reliable detection of correlated variations in 3C 273 between the X-ray band and any lower energy band. For a realistic electron distribution we demonstrate that it is probable that each decade of the seed-photon distribution from the mm to IR wavebands contributes roughly equally to the medium-energy X-ray flux. However, the expected mm variations are too small to be detected above the noise, probably explaining the lack of success of previous searches for a correlation between X-ray and mm variations. In addition, we show that the IR leads the X-rays by  0.75±0.25 d  . These observations rule out the 'External Compton' emission process for the production of the X-rays.     

A parallel domain decomposition algorithm for coastal ocean circulation models based on integer linear programming

This paper presents a new parallel domain decomposition algorithm based on integer linear programming (ILP), a mathematical optimization method. To minimize the computation time of coastal ocean circulation models, the ILP decomposition algorithm divides the global domain in local domains with balanced work load according to the number of processors and avoids computations over as many as land grid cells as possible. In addition, it maintains the use of logically rectangular local domains and achieves the exact same results as traditional domain decomposition algorithms (such as Cartesian decomposition). However, the ILP decomposition algorithm may not converge to an exact solution for relatively large domains. To overcome this problem, we developed two ILP decomposition formulations. The first one (complete formulation) has no additional restriction, although it is impractical for large global domains. The second one (feasible) imposes local domains with the same dimensions and looks for the feasibility of such decomposition, which allows much larger global domains. Parallel performance of both ILP formulations is compared to a base Cartesian decomposition by simulating two cases with the newly created parallel version of the Stevens Institute of Technology’s Estuarine and Coastal Ocean Model (sECOM). Simulations with the ILP formulations run always faster than the ones with the base decomposition, and the complete formulation is better than the feasible one when it is applicable. In addition, parallel efficiency with the ILP decomposition may be greater than one.     

Mineralogy and petrology of the Abee enstatite chondrite breccia and its dark inclusions

The Abee E4 enstatite chondrite breccia consists of clasts (many rimmed by metallic Fe, Ni), dark inclusions and matrix. The clasts and matrix were well equilibrated by thermal metamorphism, as evidenced by uniform mineral compositions, recrystallized chondrules, low MnO content of enstatite and high abundance of orthoenstatite. The clasts acquired their metal-rich rims prior to this metamorphic episode. The occurrence in Abee of relatively unmetamorphosed dark inclusions, clasts with nearly random magnetic orientations and a matrix with a uniform magnetic orientation [18,19] indicates that clast and matrix metamorphism occurred prior to the agglomeration of the breccia.The dark inclusions are an unusual kind of enstatite chondritic material, distinguished from the clasts and matrix by their relative enrichments in REE [21–23], low relative abundances of kamacite, total metallic Fe, Ni and silica, lower niningerite/(total sulfide) ratios, high relative abundances of oldhamite and martensite, smaller euhedral enstatite, more heterogeneous enstatite and metallic Fe, Ni, more calcic enstatite and more nickeliferous schreibersite.We propose the following model for the petrogenesis of the Abee breccia: The maximum metamorphic temperature of breccia parent material was?- 840°C (the minimum temperature of formation of Abee niningerite) and perhaps near 950–1000°C (the Fe-Ni-S eutectic temperature). Euhedral enstatite crystals in metallic Fe, Ni- and sulfide-rich areas grew at these metamorphic temperatures into pliable metal and sulfide. Breccia parent material was impact-excavated from depth, admixed with dark inclusions and rapidly cooled (700 to 200°C in about 2 hours) [15]. During this cooling, clast and matrix material acquired thermal remanent magnetization. Random conglomeration of clasts and unconsolidated matrix materials caused the clasts to have random magnetic orientations and the matrix areas to have net magnetic intensities of zero (due to the cancellation of numerous randomly oriented magnetic vectors of equal intensity in the matrix). A subsequent ambient magnetic field imparted a uniform net magnetic orientation to the matrix and caused the magnetic orientations of the clasts to be somewhat less random. The Abee breccia was later consolidated, possibly by shock or by shallow burial and very long-period/low-temperature (< 215°C) metamorphism.     

A parallel mesh-free contaminant transport model based on the Analytic Element and Streamline Methods

This paper introduces a new method for simulating large-scale subsurface contaminant transport that combines an Analytic Element Method (AEM) groundwater flow solution with a split-operator Streamline Method for modeling reactive transport. The key feature of the method is the manner in which the vertically integrated AEM flow solution is used to construct three-dimensional particle tracks that define the geometry of the Streamline Method. The inherently parallel nature of the algorithm supports the development of reactive transport models for spatial domains much larger than current grid-based methods. The applicability of the new approach is verified for cases with negligible transverse dispersion through comparisons to analytic solutions and existing numerical solutions, and parallel performance is demonstrated through a realistic test problem based on the regional-scale transport of agricultural contaminants from spatially distributed sources.     

Simulation of a Stably Stratified Atmospheric Boundary Layer Using One-Dimensional Turbulence

One-dimensional turbulence (ODT) is a single-column simulation in which vertical motions are represented by an unsteady advective process, rather than their customary representation by a diffusive process. No space or time averaging of mesh-resolved motions is invoked. Molecular-transport scales can be resolved in ODT simulations of laboratory-scale flows, but this resolution of these scales is prohibitively expensive in ODT simulations of the atmospheric boundary layer (ABL), except possibly in small subregions of a non-uniform mesh.Here, two methods for ODT simulation of the ABL on uniform meshes are described and applied to the GABLS (GEWEX Atmospheric Boundary Layer Study; GEWEX is the Global Energy and Water Cycle Experiment) stable boundary-layer intercomparison case. One method involves resolution of the roughness scale using a fixed eddy viscosity to represent subgrid motions. The other method, which is implemented at lower spatial resolution, involves a variable eddy viscosity determined by the local mesh-resolved flow, as in multi-dimensional large-eddy simulation (LES). When run at typical LES resolution, it reproduces some of the key high-resolution results, but its fidelity is lower in some important respects. It is concluded that a more elaborate empirically based representation of the subgrid physics, closely analogous to closures currently employed in LES of the ABL, might improve its performance substantially, yielding a cost-effective ABL simulation tool. Prospects for further application of ODT to the ABL, including possible use of ODT as a near-surface subgrid closure framework for general circulation modeling, are assessed.