Geochemistry and petrogenesis of a peralkaline granite complex from the Midian Mountains, Saudi Arabia

A zoned intrusion with a biotite granodiorite core and arfvedsonite granite rim represents the source magma for an albitised granite plug near its eastern margin and radioactive siliceous veins along its western margin. A study of selected REE and trace elements of samples from this complex reveals that the albitised granite plug has at least a tenfold enrichment in Zr, Hf, Nb, Ta, Y, Th, U and Sr, and a greatly enhanced heavy/light REE ratio compared with the peralkaline granite. The siliceous veins have even stronger enrichment of these trace elements, but a heavy/light REE ratio and negative eu anomaly similar to the peralkaline granite. It is suggested that the veins were formed from acidic volatile activity and the plug from a combination of highly fractionated magma and co-existing alkaline volatile phase. The granodiorite core intrudes the peralkaline granite and has similar trace element geochemistry. The peralkaline granite is probably derived from the partial melting of the lower crust in the presence of halide-rich volatiles, and the granodiorite from further partial melting under volatile-free conditions.     

Land-derived pollutant hydrocarbons

Evidence is given for the derivation of oil pollution from on-shore sources and for a possible relation between its concentration in estuarine sediment and in the lipids of the sediment biomass.     

Lie group variational integrators for the full body problem in orbital mechanics

Equations of motion, referred to as full body models, are developed to describe the dynamics of rigid bodies acting under their mutual gravitational potential. Continuous equations of motion and discrete equations of motion are derived using Hamilton’s principle. These equations are expressed in an inertial frame and in relative coordinates. The discrete equations of motion, referred to as a Lie group variational integrator, provide a geometrically exact and numerically efficient computational method for simulating full body dynamics in orbital mechanics; they are symplectic and momentum preserving, and they exhibit good energy behavior for exponentially long time periods. They are also efficient in only requiring a single evaluation of the gravity forces and moments per time step. The Lie group variational integrator also preserves the group structure without the use of local charts, reprojection, or constraints. Computational results are given for the dynamics of two rigid dumbbell bodies acting under their mutual gravity; these computational results demonstrate the superiority of the Lie group variational integrator compared with integrators that are not symplectic or do not preserve the Lie group structure.     

Effects of atmospheric dynamics and ocean resolution on bi-stability of the thermohaline circulation examined using the Grid ENabled Integrated Earth system modelling (GENIE) framework

We have used the Grid ENabled Integrated Earth system modelling (GENIE) framework to undertake a systematic search for bi-stability of the ocean thermohaline circulation (THC) for different surface grids and resolutions of 3-D ocean (GOLDSTEIN) under a 3-D dynamical atmosphere model (IGCM). A total of 407,000 years were simulated over a three month period using Grid computing. We find bi-stability of the THC despite significant, quasi-periodic variability in its strength driven by variability in the dynamical atmosphere. The position and width of the hysteresis loop depends on the choice of surface grid (longitude-latitude or equal area), but is less sensitive to changes in ocean resolution. For the same ocean resolution, the region of bi-stability is broader with the IGCM than with a simple energy-moisture balance atmosphere model (EMBM). Feedbacks involving both ocean and atmospheric dynamics are found to promote THC bi-stability. THC switch-off leads to increased import of freshwater at the southern boundary of the Atlantic associated with meridional overturning circulation. This is counteracted by decreased freshwater import associated with gyre and diffusive transports. However, these are localised such that the density gradient between North and South is reduced tending to maintain the THC off state. THC switch-off can also generate net atmospheric freshwater input to the Atlantic that tends to maintain the off state. The ocean feedbacks are present in all resolutions, across most of the bi-stable region, whereas the atmosphere feedback is strongest in the longitude–latitude grid and around the transition where the THC off state is disappearing. Here the net oceanic freshwater import due to the overturning mode weakens, promoting THC switch-on, but the atmosphere counteracts this by increasing net freshwater input. This increases the extent of THC bi-stability in this version of the model. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Quantifying strain birefringence halos around inclusions in diamond

The pressure and temperature conditions of formation of natural diamond can be estimated by measuring the residual stress that an inclusion remains under within a diamond. Raman spectroscopy has been the most commonly used technique for determining this stress by utilising pressure-sensitive peak shifts in the Raman spectrum of both the inclusion and the diamond host. Here, we present a new approach to measure the residual stress using quantitative analysis of the birefringence induced in the diamond. As the analysis of stress-induced birefringence is very different from that of normal birefringence, an analytical model is developed that relates the spherical inclusion size, R _i, host diamond thickness, L, and measured value of birefringence at the edge of the inclusion, \Updelta n(R_\texti )_\textav \Updelta n(R_{\text{i}} )_{\text{av}} , to the peak value of birefringence that has been encountered; to first order \Updelta n_\textpk = (3/4)(L/R_\texti )  \Updelta n(R_\texti )_\textav \Updelta n_{\text{pk}} = (3/4)(L/R_{\text{i}} ) \, \Updelta n(R_{\text{i}} )_{\text{av}} . From this birefringence, the remnant pressure (P _i) can be calculated using the photoelastic relationship \Updelta n_\textpk = - (3/4)n³ q_\textiso P_\texti \Updelta n_{\text{pk}} = - (3/4)n^{3} q_{\text{iso}} P_{\text{i}} , where q _iso is a piezo-optical coefficient, which can be assumed to be independent of crystallographic orientation, and n is the refractive index of the diamond. This model has been used in combination with quantitative birefringence analysis with a MetriPol system and compared to the results from both Raman point and 2D mapping analysis for a garnet inclusion in a diamond from the Udachnaya mine (Russia) and coesite inclusions in a diamond from the Finsch mine (South Africa). The birefringence model and analysis gave a remnant pressure of 0.53 ± 0.01 GPa for the garnet inclusion, from which a source pressure was calculated as 5.7 GPa at 1,175°C (temperature obtained from IR analysis of the diamond host). The Raman techniques could not be applied quantitatively to this sample to support the birefringence model; they were, however, applied to the largest coesite inclusion in the Finsch sample. The remnant pressure values obtained were 2.5 ± 0.1 GPa (birefringence), 2.5 ± 0.3 GPa (2D Raman map), and 2.5–2.6 GPa (Raman point analysis from all four inclusions). However, although the remnant pressures from the three methods were self-consistent, they led to anomalously low source pressure of 2.9 GPa at 1,150°C (temperature obtained from IR analysis) raising serious concerns about the use of the coesite-in-diamond geobarometer.     

Thick sections of layered ultramafic cumulates in the Oman ophiolite revealed by an airborne hyperspectral survey: Petrogenesis and relationship to mantle diapirism

Using the HyMap instrument, we have acquired visible and near infrared hyperspectral data over the Maqsad area of the Oman ophiolite (~ 15 × 60 km). This survey allowed us to identify and map the distribution of clinopyroxene-rich cumulates (inter-layered clinopyroxenites and wehrlites) whose occurrence was previously undocumented in this area. The cumulates reach several hundred meters in thickness and crop out at distances exceeding 15 km on both sides of the Maqsad former spreading centre. They occur either in mantle harzburgites, as km-sized layered intrusions surrounded by fields of pegmatitic dykes consisting of orthopyroxene-rich pyroxenite and gabbronorites, or at the base of the crustal section where they are conformably overlain by cumulate gabbros. These ultramafic cumulates crystallized from silica- and Mg-rich melts derived from a refractory mantle source (e.g. high Cr#, low [Al₂O₃], low [TiO₂]). These melts are close to high-Ca boninites, although, strictly speaking, not perfect equivalents of present-day, supra-subduction zone, boninites. Chemical stratigraphy reveals cycles of replenishment, mixing and fractional crystallization from primitive (high Mg#) melts, typical of open magma chambers and migration of inter-cumulus melts. The TiO₂ content of clinopyroxene is always low (≤ 0.2 wt.%) but quite variable compared to the associated pegmatites that are all derived from a source ultra-depleted in high field strength elements (HFSE). This variability is not caused by fractional crystallization alone, and is best explained by hybridization between the ultra-depleted melts (parent melts of the pegmatites) and the less depleted mid-ocean ridge basalts (MORB) parent of the dunitic–troctolitic–gabbroic cumulates making up the crustal section above the Maqsad diapir.We propose that, following a period of magma-starved spreading, the Maqsad mantle diapir, impregnated with tholeiitic melts of MORB affinity, reached shallow depths beneath the ocean ridge. This diapir induced melting of the formerly accreted and hydrothermally altered lithosphere. At this stage, these boninitic-like lithospheric melts crystallized as pegmatitic dykes. As the diapir continued to rise, the amount of MORB reaching shallow depths increased, together with the surrounding temperature, leading to the formation of magma chambers where the crystallization of layered cumulates became possible. These cumulates remained rich in pyroxene and devoid of plagioclase as long as the contribution of MORB-derived melts was moderate relative to the lithospheric-derived melts. As the contribution of MORB to the refilling of the magma chamber increased, gabbroic cumulates started to crystallize.     

The Use of Geographically Weighted Regression for Spatial Prediction: An Evaluation of Models Using Simulated Data Sets

Increasingly, the geographically weighted regression (GWR) model is being used for spatial prediction rather than for inference. Our study compares GWR as a predictor to (a) its global counterpart of multiple linear regression (MLR); (b) traditional geostatistical models such as ordinary kriging (OK) and universal kriging (UK), with MLR as a mean component; and (c) hybrids, where kriging models are specified with GWR as a mean component. For this purpose, we test the performance of each model on data simulated with differing levels of spatial heterogeneity (with respect to data relationships in the mean process) and spatial autocorrelation (in the residual process). Our results demonstrate that kriging (in a UK form) should be the preferred predictor, reflecting its optimal statistical properties. However the GWR-kriging hybrids perform with merit and, as such, a predictor of this form may provide a worthy alternative to UK for particular (non-stationary relationship) situations when UK models cannot be reliably calibrated. GWR predictors tend to perform more poorly than their more complex GWR-kriging counterparts, but both GWR-based models are useful in that they provide extra information on the spatial processes generating the data that are being predicted.     

Changes to the drivers of fire weather with a warming climate – a case study of southeast Tasmania

Projected changes to the global climate system have great implications for the incidence of large infrequent fires in many regions. Here we examine the synoptic-scale and local-scale influences on the incidence of extreme fire weather days and consider projections of the large-scale mean climate to explore future fire weather projections. We focus on a case study region with periodic extreme fire dangers; southeast Tasmania, Australia. We compare the performance of a dynamically downscaled regional climate model with Global Climate Model outputs as a tool for examining the local-scale influences while accounting for high regional variability. Many of the worst fires in Tasmania and the southeast Australian region are associated with deep cold fronts and strong prefrontal winds. The downscaled simulations reproduce this synoptic type with greater fidelity than a typical global climate model. The incidence of systems in this category is projected to increase through the century under a high emission scenario, driven mainly by an increase in the temperature of air masses, with little change in the strength of the systems. The regional climate model projected increase in frequency is smaller than for the global climate models used as input, with a large model range and natural variability. We also demonstrate how a blocking Foehn effect and topographic channelling contributed to the extreme conditions during an extreme fire weather day in Tasmania in January 2013. Effects such as these are likely to contribute to high fire danger throughout the century. Regional climate models are useful tools that enable various meteorological drivers of fire danger to be considered in projections of future fire danger.     

Trends in Abundance Indices of Fishes in Maryland’s Coastal Bays During 1972–2009

Maryland’s coastal bays provide habitat for juveniles of many commercially and recreationally important species of shellfish and finfish. Since 1972, the Maryland Department of Natural Resources has conducted the Maryland Coastal Bays Trawl and Seine Survey to monitor the populations of key species. The survey has undergone substantial spatial and methodological changes affecting the interpretation of simple indices of abundance. We developed generalized linear models to standardize the indices of abundance of five commonly caught fish species: Atlantic menhaden Brevoortia tyrannus, weakfish Cynoscion regalis, spot Leiostomus xanthurus, bay anchovy Anchoa mitchilli, and summer flounder Paralichthys dentatus. Density declined significantly since 1972 for menhaden, bay anchovy, and spot in at least one region within the coastal bays. The northern bays had significantly higher densities than the southern bays for all species. Changes in abundance indices of the five species examined were not related to sea grass coverage, temperature, salinity, nitrogen-to-phosphorus ratios, and other habitat variables but were likely a result of stock-wide recruitment processes.