Environmental change and oyster colonization within the Hudson River estuary linked to Holocene climate

Geophysical mapping and sampling data provide a record of changing environmental and faunal conditions within the Hudson River estuary during the mid- to late Holocene. On the shallow, broad marginal flats of the mesohaline Hudson, fossil oyster beds (Crassostrea virginica) are found exposed on the river bottom and buried by sediment. The shallowest beds are well imaged in chirp sub-bottom and side-scan sonar data and form discrete flow-perpendicular bands, 0.6–1.0 km wide and up to 3 km long, which cover 30% of the river bottom. Radiocarbon-dated sediment cores indicate oysters thrived within two time periods from ~500–2,400 and ~5,600–6,100 cal. years b.p. Sediment and physical property data indicate a changing depositional regime consistent with the oyster chronology. Similar changes in oyster presence are found in local shell midden sites of the Lower Hudson Valley as well as elsewhere along the Atlantic coast, and may reflect climatic controls associated with warm–cool cycles during the Holocene. Oysters flourished during the mid-Holocene warm period, disappeared with the onset of cooler climate at 4,000–5,000 cal. years b.p., and returned during warmer conditions of the late Holocene. The most recent demise of oysters within the Hudson at 500–900 cal. years b.p. may have accompanied the Little Ice Age.     

Fluid inclusion studies of the Abawso gold prospect,near the Ashanti Belt,Ghana

The major Ghanaian lode gold deposits are preferentially aligned along the western and eastern contacts of the Kumasi Basin with the Ashanti and Sefwi Belts, respectively. The investigated area of the Abawso small-scale concession, covering the workings of the old Ettadom mine, is situated 3 km west of the lithological contact of the Birimian metavolcanic rocks of the Akropong Belt in the east with the Birimian metasedimentary rocks of the Kumasi Basin in the west. The rocks of the Abawso concession represent a steeply NW-dipping limb of a SE-verging anticline with an axis plunging to the SW. Quartz veining occurs predominantly in the form of en échelon dilatational veins along NNE–SSW-striking shear zones of a few metres width and shows evidence of brittle and ductile deformation. Also stockwork-style quartz veining occurs in the vicinity of the main shaft of the old Ettadom mine. Hydrothermal alteration includes sericitisation, sulphidation and locally carbonatisation. The auriferous quartz veins mainly follow the trend of brittle to ductile deformed quartz veins; however, some occur in stockwork. Fluid inclusion studies reveal a large number of H₂O inclusions along intragranular trails in auriferous quartz vein samples, as well as an overall dominance of H₂O and H₂O-CO₂ inclusions over CO₂ inclusions. Textural observations and physico-chemical fluid inclusion properties indicate post-entrapment modifications for all quartz vein samples due to grain boundary migration recrystallisation. This process is interpreted to be responsible for the generation of the CO₂ inclusions from a H₂O-CO₂ parent fluid. In comparison with mineralisation at the Ashanti and Prestea deposits, which are characterised by CO₂±N₂ inclusions, the observed inclusion assemblage may be due to a shallower crustal level of mineralisation, or different degrees and styles of recrystallisation, or a less pronounced development of laminated quartz veins due to comparably restricted pressure fluctuations. Furthermore, the microthermometric observations allow the reconstruction of a possible retrograde P-T path, depicting near-isothermal decompression in the P-T range of the brittle/ductile transition.Editorial handling: E. Frimmel     

A modeling of the structure and favorable H-docking sites and defects for the high-pressure silica polymorph stishovite

Employing first-principles methods, the docking sites for H were determined and H, Al, and vacancy defects were modeled with an infinite periodic array of super unit cells each consisting of 27 contiguous symmetry nonequivalent unit cells of the crystal structure of stishovite. A geometry optimization of the super-cell structure reproduces the observed bulk structure within the experimental error when P1 translational symmetry was assumed and an array of infinite extent was generated. A mapping of the valence electrons for the structure displays mushroom-shaped isosurfaces on the O atom, one on each side of the plane of the OSi₃ triangle in the nonbonded region. An H atom, placed in a cell near the center of the super cell, was found to dock upon geometry optimization at a distance of 1.69 Å from the O atom with the OH vector oriented nearly perpendicular to the plane of the triangle such that the OH vector makes a angle of 91° with respect to [001]. However, an optimization of a super cell with an Al atom replacing Si and an H atom placed nearby in a centrally located cell resulted in an OH distance of 1.02 Å with the OH vector oriented perpendicular to [001] as observed in infrared studies. The geometry-optimized position of the H atom was found to be in close agreement with that (0.44, 0.12, 0.0) determined in an earlier study of the theoretical electron density distribution. The docking of the H atom at this site was found to be 330 kJ mol^–1 more stable than a docking of the atom just off the shared OO edge of the octahedra as determined for rutile. A geometry optimization of a super cell with a missing Si generated a vacant octahedra that is 20% larger than that of the SiO₆ octahedra. The valence electron density distribution displayed by the two-coordinate O atoms that coordinate the vacant octahedral site is very similar to those displayed by the bent SiOSi angles in coesite. The internal distortions induced by the defect were found to diminish rather rapidly with distance, with the structure annealing to that observed in the bulk crystal to within about three coordination spheres.     

A reappraisal of episodic burial metamorphism in the Andes of central Chile

The Jurassic to Miocene sequences of the central Andes, east of Santiago, reputedly show repeated cycles of episodic sub-greenschist facies, burial metamorphism that are identified by sharp breaks in metamorphic grade at major stratigraphic boundaries. This paper presents the first detailed petrochemical analysis of these low-grade metamorphic sequences by examining the progressive development of secondary minerals, reaction progress in mafic phyllosilicates, and topological variations in the low-grade assemblages as a means of testing this model. The results indicate a progressive increase from zeolite facies through to close to the onset of greenschist facies from Miocene to Jurassic rocks. Combined analysis of reaction progress in mafic phyllosilicates and petrochemical relationships of chlorite–pumpellyite–actinolite in metabasites provides no evidence for sharp metamorphic breaks at major stratigraphic boundaries. Integrating the results presented here with the most recent models of stratigraphic/tectonic development of the central Andes shows that the metamorphism took place in two episodes, and was not episodic on a 40-million-year cycle. An absence of sharp breaks in metamorphic grade in any part of the succession, as demonstrated here, shows that the original petrographic establishment of low-grade facies provided insufficient resolution of changes in metamorphic conditions to establish definitive evidence of such breaks. Accordingly, this study suggests that re-assessment of metamorphic breaks reputedly developed in other areas of the Andean Cordillera is imperative in order to resolve the questions raised here about the origin of the low-grade metamorphism.Editorial responsibility: B. Collins     

Dauphiné twinning and texture memory in polycrystalline quartz. Part 3: texture memory during phase transformation

Samples of quartz-bearing rocks were heated above the α (trigonal)–β (hexagonal) phase transformation of quartz (625–950°C) to explore changes in preferred orientation patterns. Textures were measured both in situ and ex situ with neutron, synchrotron X-ray and electron backscatter diffraction. The trigonal–hexagonal phase transformation does not change the orientation of c- and a-axes, but positive and negative rhombs become equal in the hexagonal β-phase. In naturally deformed quartzites measured by neutron diffraction a perfect texture memory was observed, i.e. crystals returned to the same trigonal orientation they started from, with no evidence of twin boundaries. Samples measured by electron back-scattered diffraction on surfaces show considerable twinning and memory loss after the phase transformation. In experimentally deformed quartz rocks, where twinning was induced mechanically before heating, the orientation memory is lost. A mechanical model can explain the memory loss but so far it does not account for the persistence of the memory in quartzites. Stresses imposed by neighboring grains remain a likely cause of texture memory in this mineral with a very high elastic anisotropy. If stresses are imposed experimentally the internal stresses are released during the phase transformation and the material returns to its original state prior to deformation. Similarly, on surfaces there are no tractions and thus texture memory is partially lost.     

Evaluation of a Photosynthesis-Based Canopy Resistance Formulation in the Noah Land-Surface Model

Accurately representing complex land-surface processes balancing complexity and realism remains one challenge that the weather modelling community is facing nowadays. In this study, a photosynthesis-based Gas-exchange Evapotranspiration Model (GEM) is integrated into the Noah land-surface model replacing the traditional Jarvis scheme for estimating the canopy resistance and transpiration. Using 18-month simulations from the High Resolution Land Data Assimilation System (HRLDAS), the impact of the photosynthesis-based approach on the simulated canopy resistance, surface heat fluxes, soil moisture, and soil temperature over different vegetation types is evaluated using data from the Atmospheric Radiation Measurement (ARM) site, Oklahoma Mesonet, 2002 International H₂O Project (IHOP_2002), and three Ameriflux sites. Incorporation of GEM into Noah improves the surface energy fluxes as well as the associated diurnal cycle of soil moisture and soil temperature during both wet and dry periods. An analysis of midday, average canopy resistance shows similar day-to-day trends in the model fields as seen in observed patterns. Bias and standard deviation analyses for soil temperature and surface fluxes show that GEM responds somewhat better than the Jarvis scheme, mainly because the Jarvis approach relies on a parametrised minimum canopy resistance and meteorological variables such as air temperature and incident radiation. The analyses suggest that adding a photosynthesis-based transpiration scheme such as GEM improves the ability of the land-data assimilation system to simulate evaporation and transpiration under a range of soil and vegetation conditions.     

Climate Suitability for Stable Malaria Transmission in Zimbabwe Under Different Climate Change Scenarios

Climate is one factor that determines the potential range of malaria. As such, climate change may work with or against efforts to bring malaria under control. We developed a model of future climate suitability for stable Plasmodium falciparum malaria transmission in Zimbabwe. Current climate suitability for stable malaria transmission was based on the MARA/ARMA model of climatic constraints on the survival and development of the Anopheles vector and the Plasmodium falciparum malaria parasite. We explored potential future geographic distributions of malaria using 16 projections of climate in 2100. The results suggest that, assuming no future human-imposed constraints on malaria transmission, changes in temperature and precipitation could alter the geographic distribution of malaria in Zimbabwe, with previously unsuitable areas of dense human population becoming suitable for transmission. Among all scenarios, the highlands become more suitable for transmission, while the lowveld and areas with low precipitation show varying degrees of change, depending on climate sensitivity and greenhouse gas emission stabilization scenarios, and depending on the general circulation model used. The methods employed can be used within or across other African countries.     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

Numerical evaluation of mean-field homogenisation methods for predicting shale elastic response

Homogenisation techniques have been successfully used to estimate the mechanical response of synthetic composite materials, due to their ability to relate the macroscopic mechanical response to the material microstructure. The adoption of these mean-field techniques in geo-composites such as shales is attractive, partly because of the practical difficulties associated with the experimental characterisation of these highly heterogeneous materials. In this paper, numerical modelling has been undertaken to investigate the applicability of homogenisation methods in predicting the macroscopic, elastic response of clayey rocks. The rocks are considered as two-level composites consisting of a porous clay matrix at the first level and a matrix-inclusion morphology at the second level. The simulated microstructures ranged from a simple system of one inclusion/void embedded in a matrix to complex, random microstructures. The effectiveness and limitations of the different homogenisation schemes were demonstrated through a comparative evaluation of the macroscopic elastic response, illustrating the appropriate schemes for upscaling the microstructure of shales. Based on the numerical simulations and existing experimental observations, a randomly distributed pore system for the micro-structure of porous clay matrix has been proposed which can be used for the subsequent development and validation of shale constitutive models. Finally, the homogenisation techniques were used to predict the experimental measurements of elastic response of shale core samples. The developed methodology is proved to be a valuable tool for verifying the accuracy and performance of the homogenisation techniques.