Infaunal burrow ventilation and pore-water transport in muddy sediments

The ventilation of burrows by tube-dwelling benthos is understood to be important in determining rates of exchange of solutes between the sediment and overlying water. However, few models have attempted to link the burrow ventilation behavior of tube-dwelling organisms with their geochemical consequences. The classic cylinder model of bioirrigation in muddy sediments (Aller, R.C., 1980. Quantifying solute distributions in the bioturbated zone of marine sediments by defining an average microenvironment. Geochimica et Cosmochimica Acta 44, 1955–1965) links pore-water processes and burrow sizes and distributions in the sediment by assuming that burrows are fully flushed. The equivalence between the cylinder model and the more commonly used one-dimensional non-local exchange model depends upon this assumption. However, this assumption has seldom been tested in the field. We have extended the cylinder model of bioirrigation to include burrow ventilation activities of organisms. Burrow ventilation is modeled as a simple non-local exchange of burrow water with overlying water. Model simulations indicate that burrow ventilation has a large effect on vertical profiles and fluxes of solute tracers. We collected data on burrow geometry in the field by CT-scanning freshly collected sediment cores. At the same study sites, we measured activity profiles of ²²²Rn, a naturally occurring radionuclide tracer of pore-water transport. With model geometry independently constrained, we tested the model by comparing our model-predicted profiles with measured profiles. Our results demonstrate that burrows in the field are not fully flushed. Our estimated burrow ventilation rates compare favorably with previous laboratory measurements. The inclusion of realistic burrow ventilation in this pore-water transport model strongly affects modeled solute profiles and fluxes. We demonstrate how model parameters can be determined from field samples and present a model that more realistically simulates pore-water transport processes in muddy sediments.     

A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.     

Nickel-copper mineralization at Talnotry,Newton Stewart,Scotland

Summary Nickel-copper mineralization occurs near the base of a diorite intrusion close to its contact with hornfelsed Ordovician and Silurian shales and greywackes. The principal ore minerals are nickeline, gersdorffite, pyrrhotine, pentlandite and chalcopyrite with minor amounts of molybdenite, tellurobismutite, gold, sphalerite and argentopentlandite. Pyrite, marcasite, violarite and goethite also occur but are interpreted as later alteration products. Much of the pyrrhotine-rich mineralization at the base of the intrusion is in the form of blebs and interstitial aggregates with amphiboles, plagioclase feldspar, biotite, chlorite and quartz. Chalcopyrite-rich and nickeline-gersdorffite-rich mineralization occurs above this and immediately below unmineralized diorite in the form of patches, lenticular masses and stringers along joints and fractures.Whole rock and ore analyses and electron microprobe data on the silicates, sulphides and sulpharsenides are presented.The unmineralized diorite has low SiO₂ and high MgO contents compared to typical diorites and relatively high Cr, Ni and Ti trace element values. In the mineralized diorite, platinum-group elements occur in very low concentrations in the pyrrhotineand chalcopyrite-rich assemblages but Pt, Pd and Au show significant enrichment in the nickeline-gersdorffite-rich mineralization.A magmatic origin for the mineralization is proposed rather than formation by hydrothermal solutions or metasomatism.

---

Die nickel-kupfer vererzung von Talnotry, Newton Stewart, Schottland

Zusammenfassung Die Nickel-Kupfer Vererzung befindet sich im Kontaktbereich einer Dioritintrusion mit hornfelsdurchdrungenen silurischen und ordovizischen Grauwacken. Die am häufigsten auftretenden Erzmineralien sind Nickelin, Gersdorffit, Pyrrhotin, Pentlandit, und Chalcopyrit mit kleineren Mengen von Molybdenit, Tellurwismuth, Gold, Sphalerit und Argentopentlandit. Weiters treten als spätere Umwandlungsprodukte Pyrit, Markasit, Violarit und Goethit auf. Ein großer Teil der Pyrrhotin-reichen Vererzung am Fuße des eingedrungenen Diorits bildet fleckige, lückenfüllende Aggregate zwischen den Silikatphasen (Amphibol, Plagioklas, Biotit, Chlorit, Quarz). Die Chalcopyrit- und Nickelin-Gersdorffit-reiche Vererzung tritt zwischen dem unvererzten Diorit im Hangenden und der Pyrrhotin-reichen Vererzung im Liegenden als linsige, lappenartige Konzentration und entlang von Fugen und Brücken auf.Vollgesteins-Analysen sowie ausgewählte Erzanalysen und Mikrosondenergebnisse der Silikate, Sulfide und Sulfarsenide liegen vor.Im Unterschied zu den typischen Dioriten zeigt der unvererzte Diorit einen niedrigen SiO₂ und hohen MgO Gehalt und relativ hohe Cr, Ni und Ti Spurenelementwerte. Die Platingruppenelemente sind in den Pyrrhotin- und Chalcopyrit-reichen Vererzungen nur schwach vertreten, jedoch sind Pt, Pd und An in der Gersdorffit-reichen Vererzung stark angereichert.Ein magmatischer Ursprung dieser Vererzung, im Gegensatz zu einer metasomatischen Entstehung oder Bildung aus hydrothermalen Lösungen wird vermutet.

---

---

With 3 Figures     

Thermodynamic characterization of boltwoodite and uranophane: Enthalpy of formation and aqueous solubility

Boltwoodite and uranophane are uranyl silicates common in oxidized zones of uranium ore deposits. An understanding of processes that impact uranium transport in the environment, especially pertaining to the distribution of uranium between solid phases and aqueous solutions, ultimately requires determination of thermodynamic parameters for such crystalline materials. We measured formation enthalpies of synthetic boltwoodites, K(UO₂)(HSiO₄)·H₂O and Na(UO₂)(HSiO₄)·H₂O, and uranophane, Ca(UO₂)₂(HSiO₄)₂·5H₂O, by high temperature oxide melt solution calorimetry. We also studied the aqueous solubility of these phases from both saturated and undersaturated conditions at a variety of pH. The combined data permit the determination of standard enthalpies, entropies and Gibbs free energies of formation for each phase and analysis of its potential geological impact from a thermodynamic point of view.     

Influence of Environmental Factors on Spawning of the American Horseshoe Crab (<Emphasis Type="Italic">Limulus polyphemus</Emphasis>) in the Great Bay Estuary,New Hampshire,USA

The Great Bay Estuary, New Hampshire, USA is near the northern distribution limit of the American horseshoe crab (Limulus polyphemus). This estuary has few ideal beaches for spawning, yet it supports a modest population of horseshoe crabs. There is no organized monitoring program in the Great Bay Estuary, so it is unclear when and where spawning occurs. In this 2-year study (May through June, 2012 and 2013), >5,000 adult horseshoe crabs were counted at four sites in the estuary. The greatest densities of horseshoe crabs were observed at Great Bay sites in the upper, warmer reaches of the estuary. Peaks of spawning activity were not strongly correlated with the times of the new or full moons, and similar numbers of horseshoe crabs were observed mating during daytime and nighttime high tides. While many environmental factors are likely to influence the temporal and spatial patterns of spawning in this estuary, temperature appears to have the most profound impact.     

Extended Finite Element Method for the Analysis of Discontinuities in Rock Masses

The strength and deformability of rock mass primarily depend on the condition of joints and their spacing and partially on the engineering properties of rock matrix. Till today, numerical analysis of discontinuities e.g. joint, fault, shear plane and others is conducted placing an interface element in between two adjacent rock matrix elements. However, the applicability of interface elements is limited in rock mechanics problems having multiple discontinuities due to its inherent numerical difficulties often leading to non-convergent solution. Recent developments in extended finite element method (XFEM) having strong discontinuity imbedded within a regular element provide an opportunity to analyze discrete discontinuities in rock masses without any numerical difficulties. This concept is based on partition of unity principle and can be used for cohesive rock joints. This paper summarizes the mathematical frameworks for the implementation of strong discontinuities in 3 and 6 nodded triangular elements and also provides numerical examples of the application of XFEM in one and two dimensional problems with single and multiple discontinuities.     

Oscillatory zoning in garnet from the Willsboro Wollastonite Skarn, Adirondack Mts, New York: a record of shallow hydrothermal processes preserved in a granulite facies terrane

Oscillatory zoning in low δ¹⁸O skarn garnet from the Willsboro wollastonite deposit, NE Adirondack Mts, NY, USA, preserves a record of the temporal evolution of mixing hydrothermal fluids from different sources. Garnet with oscillatory zoning are large (1–3 cm diameter) euhedral crystals that grew in formerly fluid filled cavities. They contain millimetre‐scale oscillatory zoning of varying grossular–andradite composition (X_Adr = 0.13–0.36). The δ¹⁸O values of the garnet zones vary from 0.80 to 6.26‰ VSMOW and correlate with X_Adr. The shape, pattern and number of garnet zones varies from crystal to crystal, as does the magnitude of the correlated chemistry changes, suggesting fluid system variability, temporal and/or spatial, over the time of garnet growth. The zones of correlated Fe content and δ¹⁸O indicate that a high Fe³⁺/Al, high δ¹⁸O fluid mixed with a lower Fe³⁺/Al and δ¹⁸O fluid. The high δ¹⁸O, Fe enriched fluids were likely magmatic fluids expelled from crystallizing anorthosite. The low δ¹⁸O fluids were meteoric in origin. These are the first skarn garnet with oscillatory zoning reported from granulite facies rocks. Geochronologic, stable isotope, petrologic and field evidence indicates that the Adirondacks are a polymetamorphic terrane, where localized contact metamorphism around shallowly intruded anorthosite was followed by a regional granulite facies overprint. The growth of these garnet in equilibrium with meteoric and magmatic fluids indicates an origin in the shallow contact aureole of the anorthosite prior to regional metamorphism. The zoning was preserved due to the slow diffusion of oxygen and cations in the large garnet and protection from deformation and recrystallization in zones of low strain in thick, rigid, garnetite layers. The garnet provide new information about the hydrothermal system adjacent to the shallowly intruded massif anorthosite that predates regional metamorphism in this geologically complex, polymetamorphic terrane.     

Tuning a physically-based model of the air–sea gas transfer velocity

Air–sea gas transfer velocities are estimated for one year using a 1-D upper-ocean model (GOTM) and a modified version of the NOAA–COARE transfer velocity parameterization. Tuning parameters are evaluated with the aim of bringing the physically based NOAA–COARE parameterization in line with current estimates, based on simple wind-speed dependent models derived from bomb-radiocarbon inventories and deliberate tracer release experiments. We suggest that A = 1.3 and B = 1.0, for the sub-layer scaling parameter and the bubble mediated exchange, respectively, are consistent with the global average CO₂ transfer velocity k. Using these parameters and a simple 2nd order polynomial approximation, with respect to wind speed, we estimate a global annual average k for CO₂ of 16.4 ± 5.6 cm h^?1 when using global mean winds of 6.89 m s^?1 from the NCEP/NCAR Reanalysis 1 1954–2000. The tuned model can be used to predict the transfer velocity of any gas, with appropriate treatment of the dependence on molecular properties including the strong solubility dependence of bubble-mediated transfer. For example, an initial estimate of the global average transfer velocity of DMS (a relatively soluble gas) is only 11.9 cm h^?1 whilst for less soluble methane the estimate is 18.0 cm h^?1.     

Metals in eggs and hatchlings of Dermochelys coriacea (testudines) in Espírito Santo,Brazil and their relation to hatching and emergence success

Certain concentration of metal may influence survival in embryonic stages, reducing the birthrate and hatching emergence. This study aimed to evaluate concentration of metal in eggs and hatchlings of Dermochelys coriacea and correlate metal concentrations in eggs with hatching and emergence success. The samples were collected at Biological Reserve Station, ES, Brazil. At the nests' opening, unhatched eggs and stillborn hatchlings were collected to survey the successful hatching and emergence of D. coriacea. The eggs shell, egg content, and dead hatchlings metal concentrations were compared to hatching and emergence success. A positive correlation was found between the concentration of copper (Cu) in the hatchlings and the success of the emergence (r² = .28, p < .05), and a moderate positive correlation between the concentrations of iron (Fe) and barium (Ba; r² = .44, p < .05) and success hatching. The concentrations of metals in eggs and hatchlings of D. coriacea in the Espírito Santo are below the levels that could be considered harmful to the species. In summary, there is a positive correlation between metals levels and hatching and emergence success.     

A ternary solid solution model for omphacite and its application to geothermobarometry of eclogites from the Middle Adula nappe (Central Alps, Switzerland)

A ternary solid solution model for omphacite with the end-members jadeite (NaAlSi₂O₆), diopside (CaMgSi₂O₆) and hedenbergite (CaFeSi₂O₆) was derived from experimental data from the literature. The subregular solution model, fitted by linear programming, is best suited to omphacites with very little aegirine component in common eclogites. Applying this solution model to the calculation of equilibrium phase diagrams of eclogites from the Adula nappe (Central Alps, Switzerland) results in large stability fields for common eclogite assemblages (garnet+omphacite+quartz+H₂O±kyanite). Within this field the compositions of garnet and omphacite show very little variation. A precise determination of the peak-pressure and temperature is not possible. The occurrence of amphibole, overgrowing the peak-pressure assemblage in fresh eclogite, suggests retrograde re-equilibration, still under eclogite facies conditions. The computation of isopleths for garnet and pyroxene end-members allows the estimation of the pressure and temperature conditions of this re-equilibration event (19–21  kbar, c .  700 °C).