Mississippian coral assemblages from the Khenifra area,Central Morocco: Biostratigraphy,biofacies, palaeoecology and palaeobiogeography

Analysis of Mississippian coral assemblages from the Khenifra region of Central Morocco together with data from foraminiferal/algal microfossils has established new age dating of 5 localities within the Azrou–Khenifra Basin: Souk El Had and Sidi Lamine, where corals occur mainly in biostromes protected by oolitic shoals, Tabainout, where corals have been recorded in different environments related to microbial mounds, Alhajra Almatkouba, where corals occur in biostromal reworked beds and Tiouinine, where corals occur in a well structured, fringing reef. This study demonstrates the presence of richer more diverse coral assemblages than previously recorded, in a variety of environmental settings. These coral assemblages strengthen correlations with the Adarouch area in the NE part of the Azrou–Khenifra Basin. It is emphasised that in the upper Viséan there are close similarities with rugose coral assemblages in other parts of the Western Palaeotethys including North Africa, SW Spain and NW Europe, and that all belong to the same biogeographic province.     

Mixing and differentiation in the Oruanui rhyolitic magma, Taupo, New Zealand: evidence from volatiles and trace elements in melt inclusions

Large pyroclastic rhyolites are snapshots of evolving magma bodies, and preserved in their eruptive pyroclasts is a record of evolution up to the time of eruption. Here we focus on the conditions and processes in the Oruanui magma that erupted at 26.5 ka from Taupo Volcano, New Zealand. The 530 km³ (void-free) of material erupted in the Oruanui event is comparable in size to the Bishop Tuff in California, but differs in that rhyolitic pumice and glass compositions, although variable, did not change systematically with eruption order. We measured the concentrations of H₂O, CO₂ and major and trace elements in zoned phenocrysts and melt inclusions from individual pumice clasts covering the range from early to late erupted units. We also used cathodoluminescence imaging to infer growth histories of quartz phenocrysts. For quartz-hosted inclusions, we studied both fully enclosed melt inclusions and reentrants (connecting to host melt through a small opening). The textures and compositions of inclusions and phenocrysts reflect complex pre-eruptive processes of incomplete assimilation/partial melting, crystallization differentiation, magma mixing and gas saturation. ‘Restitic’ quartz occurs in seven of eight pumice clasts studied. Variations in dissolved H₂O and CO₂ in quartz-hosted melt inclusions reflect gas saturation in the Oruanui magma and crystallization depths of ∼3.5–7 km. Based on variations of dissolved H₂O and CO₂ in reentrants, the amount of exsolved gas at the beginning of eruption increased with depth, corresponding to decreasing density with depth. Pre-eruptive mixing of magma with varying gas content implies variations in magma bulk density that would have driven convective mixing. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Shape and distribution analysis of Merensky Reef potholing, Northam Platinum Mine, western Bushveld Complex: implications for pothole formation and growth

Syn-magmatic removal of the cumulate pile during the formation of the Bushveld Complex resulted in “potholes”. Erosion progressed downward in the cumulate pile, resulting in a series of steep, transgressive contacts between locally conformable potholed reefs in the regional pothole sub-facies of the Swartklip Facies in the western limb of the Bushveld Complex. The deepest of these potholes, “third-order” or “FWP2” potholing, occurs where the base of the Merensky Cyclic Unit transgresses the Upper Pseudo-Reef Chromitite marker horizon. The base of a FWP2 pothole on Northam Platinum Mine consists of an unconformable stringer Merensky Chromitite overlain by a medium-grained, poikilitic orthopyroxenite and underlain by either a pegmatitic harzburgite or the medium-grained Lower Pseudo-Reef Anorthosite. Detailed shape and distribution analysis of FWP2 potholes reveals underlying patterns in their shape and distribution which, in turn, suggest a structural control. The ratio between pothole short vs long axes is 0.624 (N=1,385), although the ratio increases from 0.48 to 0.61 in the long axis range 10 to 60 m, then decreases from 0.61 to 0.57 from 61 to 100 m, increasing again from 0.57 to 0.61 from 101 to 400 m, suggesting that there is not a simple relationship between pothole shape and size. Shape (circularity, eccentricity, and dendricity) analysis of a subset of 638 potholes indicates that potholes with long axes <100 m have an elliptical, average normalized shape, elongate on a 120–150° orientation. Potholes with long axis lengths >100 m have an average normalized shape that is bilobate and elongate on a 120° orientation. The average aspect ratio (short axis length divided by long axis length) of potholes is highest for potholes with long axis lengths >100 m and lowest for potholes with long axis lengths between 35 and 60 m. The most common long axis orientation for potholes with long axis lengths <100 m is 150° but 120° for long axis lengths >100 m. Fractal analysis indicates that the distribution of pothole centers is controlled neither by a single nor several interacting fractal dimensions. Autocorrelation (Fry) analysis of the distribution of pothole centers shows recurring pothole distribution trends at 038, 070, and 110° for potholes over the full range of long axis lengths, while the trends of 008 and 152° occur in potholes with long axes lengths between 60 and 100 m. Chi-squared (X ²) analysis of the locations of pothole centers suggests that the distribution of small potholes is highly non-uniform but becomes exponentially more uniform with increasing pothole size. The model which best fits the observed shape and distribution analysis is a combination of protracted independent growth and “nearest neighbor” merging along specific orientations. For instance, the clustered distribution of original pothole centers resulted in merged potholes with long axes lengths of up to 60 m, exhibiting short vs long axes ratios of 0.61, preferred orientations of 150°, and alignment along 010 and 150° trends. Further independent growth allowed for merging of similar-sized (and smaller) neighboring potholes, generating potholes with long axes of up to 100 m in length, a preferred long axis orientation of 150°, and alignment along 010, 040, 075, and 150°. Subsequent preferential merging occurred along a 120° trend, thereby preserving a bilobate form. This implies that while pothole initiation and enlargement may be driven by a “top-down” (i.e., possibly thermomechanical) process, an underlying linear or structural catalyst/control is revealed in changes in pothole shape during enlargement and, furthermore, in the preferred trends along which potholes merged over a considerable period, possibly concomitant with adjustment of major structures in the footwall to the Bushveld Complex and pulses into the magma chamber.     

International Year of Astronomy 2009

IYA2009 will be a great opportunity for all of us working in astronomy (and space) to make an impact; certainly locally and nationally and, we hope, globally. So what is going on and how can you participate? Ian Robson has the answers.     

Large-scale impacts of bottom trawling on shelf primary productivity

Disturbance of the seabed resulting from bottom trawling affects ecosystem processes, such as the rate and magnitude of nutrient regeneration. The potential responses of the plankton community arising from such effects can be modelled, provided that reliable data on the effects on nutrient fluxes are available. In a north Cretan outer continental shelf and upper slope fishing ground (Heraklion Bay, Crete, Eastern Mediterranean) we applied a new field instrument which can simulate the passage of trawl groundropes across the sea floor and made direct seasonal measurements of the rate of dissolved and particulate nutrient releases resulting from seabed disturbance. These observational data were then integrated in a 3D ecosystem model. Results revealed that bottom trawling may trigger off considerable productivity pulses, in addition to pulses from the natural seasonal cycle.     

Deuterium tracing for the estimation of transpiration from trees Part 1. Field calibration

Flow rates calculated using the deuterium tracing method were compared with measurements of direct water uptake of excavated trees which had been inserted into water containers. Taking into account that the tracing method provides an estimate of the weighted mean flow rate over the period that the tracer pulse is present (the weighting factor is the relative concentration of the tracer as measured at the sampling point) agreement was obtained within the experimental error (about 10%) associated with the tracing method.     

Calcrete profile development in Quaternary alluvial sequences,southeast Spain: implications for using calcretes as a basis for landform chronologies

A detailed study of the morphology and micro‐morphology of Quaternary alluvial calcrete profiles from the Sorbas Basin shows that calcretes may be morphologically simple or complex. The ‘simple’ profiles reflect pedogenesis occurring after alluvial terrace formation and consist of a single pedogenic horizon near the land surface. The ‘complex’ profiles reflect the occurrence of multiple calcrete events during terrace sediment aggradation and further periods of pedogenesis after terrace formation. These ‘complex’ calcrete profiles are consequently described as composite profiles. The exact morphology of the composite profiles depends upon: (1) the number of calcrete‐forming events occurring during terrace sediment aggradation; (2) the amount of sediment accretion that occurs between each period of calcrete formation; and (3) the degree of pedogenesis after terrace formation. Simple calcrete profiles are most useful in establishing landform chronologies because they represent a single phase of pedogenesis after terrace formation. Composite profiles are more problematic. Pedogenic calcretes that form within them may inherit carbonate from calcrete horizons occurring lower down in the terrace sediments. In addition erosion may lead to the exhumation of older calcretes within the terrace sediment. Calcrete ‘inheritance’ may make pedogenic horizons appear more mature than they actually are and produce horizons containing carbonate embracing a range of ages. Calcrete exhumation exposes calcrete horizons whose morphology and radiometric ages are wholly unrelated to terrace surface age. Composite profiles are, therefore, only suitable for chronological studies if the pedogenic horizon capping the terrace sequence can be clearly distinguished from earlier calcrete‐forming events. Thus, a detailed morphological/micro‐morphological study is required before any chronological study is undertaken. This is the only way to establish whether particular calcrete profiles are suitable for dating purposes. Copyright © 2003 John Wiley & Sons, Ltd.     

Variations in the geomagnetic dipole moment during the Holocene and the past 50 kyr

All absolute paleointensity data published in peer-reviewed journals were recently compiled in the GEOMAGIA50 database. Based on the information in GEOMAGIA50, we reconstruct variations in the geomagnetic dipole moment over the past 50  kyr, with a focus on the Holocene period. A running-window approach is used to determine the axial dipole moment that provides the optimal least-squares fit to the paleointensity data, whereas associated error estimates are constrained using a bootstrap procedure. We subsequently compare the reconstruction from this study with previous reconstructions of the geomagnetic dipole moment, including those based on cosmogenic radionuclides (¹⁰Be and ¹⁴C). This comparison generally lends support to the axial dipole moments obtained in this study. Our reconstruction shows that the evolution of the dipole moment was highly dynamic, and the recently observed rates of change (5% per century) do not appear unique. We observe no apparent link between the occurrence of archeomagnetic jerks and changes in the geomagnetic dipole moment, suggesting that archeomagnetic jerks most likely represent drastic changes in the orientation of the geomagnetic dipole axis or periods characterized by large secular variation of the non-dipole field. This study also shows that the Holocene geomagnetic dipole moment was high compared to that of the preceding  40  kyr, and that  4 · 10²²  Am² appears to represent a critical threshold below which geomagnetic excursions and reversals occur.     

Quantifying volcanic ash fall hazard to electricity infrastructure

Quantifying the potential ash fall hazards from re-awakening volcanoes is a topic of great interest. While methods for calculating the probability of eruptions, and for numerical simulation of tephra dispersal and fallout exist, event records at most volcanoes that re-awaken sporadically on decadal to millennial cycles are inadequate to develop rigorous forecasts of occurrence, much less eruptive volume. Here we demonstrate a method by which eruption records from radiocarbon-dated sediment cores can be used to derive forecasting models for ash fall impacts on electrical infrastructure. Our method is illustrated by an example from the Taranaki region of New Zealand. Radiocarbon dates, expressed as years before present (B.P.), are used to define an age-depth model, classifying eruption ages (with associated errors) for a circa 1500–10 500 year B.P. record at Mt. Taranaki (New Zealand). In addition, data describing the youngest 1500 years of eruption activity is obtained from directly dated proximal deposits. Absence of trend and apparent independence in eruption intervals is consistent with a renewal model using a mix of Weibulls distributions, which was used to generate probabilistic forecasts of eruption recurrence. After establishing that interval length and tephra thickness were independent in the record, a thickness–volume relationship (from [Rhoades, D.A., Dowrick, D.J., Wilson, C.J.N., 2002. Volcanic hazard in New Zealand: Scaling and attenuation relations for tephra fall deposits from Taupo volcano. Nat. Hazards, 26:147–174]) was inverted to provide a frequency–volume relationship for eruptions. Monte Carlo simulation of the thickness–volume relationship was then used to produce probable ash fall thicknesses at any chosen site. Several critical electrical infrastructure sites in the Taranaki Region were analysed. This region, being the only gas and condensate-producing area in New Zealand, is of national economic importance, with activities in and around the area depending on uninterrupted power supplies. Forecasts of critical ash thicknesses (1 mm wet and 2 mm dry) that may cause short-circuiting, surges or power shutdowns in substations show that the annual probabilities of serious impact are between ~ 0.5% and 27% over a 50 year period. It was also found that while large eruptions with high ash plumes tend to affect “expected” areas in relation to prevailing winds, the direction impacts of small ash falls are far less predictable. In the Taranaki case study, areas out of normal downwind directions, but close to the volcano, have probabilities of impact for critical thicknesses of 1–2 mm of around half to 60% of those in downwind directions and therefore should not be overlooked in hazard analysis. Through this method we are able to definitively show that the potential ash fall hazard to electrical infrastructure in this area is low in comparison to other natural threats, and provide a quantitative measure for use in risk analysis and budget prioritisation for hazard mitigation measures.     

Passive L-band microwave soil moisture retrieval error arising from topography in otherwise uniform scenes

Soil moisture is an important component of the water cycle and will be measured for the first time on a global scale by a dedicated passive L-band microwave radiometer that is planned for launch in 2008. Here, the contribution of topography to the error budget is examined for a vegetated scene with uniform microwave emission. Dual-polarization brightness temperature curves were generated over a range of look angles for 1-D scenes with simple geometrical features, and the soil moisture was retrieved assuming a flat surface. The errors were small for the scenarios considered. Theoretical errors were tested for realistic topography with a DEM transect of a mountainous region, and were found to be comparable. Knowledge of the mean slope from high-resolution DEM data can be used to improve the accuracy of the retrieval.