New mineral activity–composition relations for thermodynamic calculations in metapelitic systems

New activity–composition (a–x) relations for minerals commonly occurring in metapelites are presented for use with the internally consistent thermodynamic dataset of Holland & Powell ( 2011 , Journal of Metamorphic Geology, 29 , 333–383). The a–x relations include a broader consideration of Fe₂O₃ in minerals, changes to the formalism of several phases and order–disorder in all ferromagnesian minerals where Fe–Mg mixing occurs on multiple sites. The a–x relations for chlorite, biotite, garnet, chloritoid, staurolite, cordierite, orthopyroxene, muscovite, paragonite and margarite have been substantially reparameterized using the approach outlined in the companion paper in this issue. For the first time, the entire set of a–x relations for the common ferromagnesian minerals in metapelitic rocks is parameterized simultaneously, with attention paid to ensuring that they can be used together to calculate phase diagrams of geologically appropriate topology. The a–x relations developed are for use in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCKFMASHTO) system for both subsolidus and suprasolidus conditions. Petrogenetic grids in KFMASH and KFMASHTO are similar in topology to those produced with earlier end‐member datasets and a–x relations, but with some notable differences. In particular, in subsolidus equilibria, the FeO/(FeO + MgO) of garnet is now greater than in coexisting staurolite, bringing a number of key staurolite‐bearing equilibria into better agreement with inferences from field and petrographic observations. Furthermore, the addition of Fe³⁺ and Ti to a number of silicate phases allows more plausible equilibria to be calculated in relevant systems. Pseudosections calculated with the new a–x relations are also topologically similar to equivalent diagrams using earlier a–x relations, although with many low variance fields shifting in P–T space to somewhat lower pressure conditions.     

Constraints on the timing and conditions of high‐grade metamorphism,charnockite formation and fluid–rock interaction in the Trivandrum Block,southern India

Incipient charnockites have been widely used as evidence for the infiltration of CO₂‐rich fluids driving dehydration of the lower crust. Rocks exposed at Kakkod quarry in the Trivandrum Block of southern India allow for a thorough investigation of the metamorphic evolution by preserving not only orthopyroxene‐bearing charnockite patches in a host garnet–biotite felsic gneiss, but also layers of garnet–sillimanite metapelite gneiss. Thermodynamic phase equilibria modelling of all three bulk compositions indicates consistent peak‐metamorphic conditions of 830–925 °C and 6–9 kbar with retrograde evolution involving suprasolidus decompression at high temperature. These models suggest that orthopyroxene was most likely stabilized close to the metamorphic peak as a result of small compositional heterogeneities in the host garnet–biotite gneiss. There is insufficient evidence to determine whether the heterogeneities were inherited from the protolith or introduced during syn‐metamorphic fluid flow. U–Pb geochronology of monazite and zircon from all three rock types constrains the peak of metamorphism and orthopyroxene growth to have occurred between the onset of high‐grade metamorphism at c. 590 Ma and the onset of melt crystallization at c. 540 Ma. The majority of metamorphic zircon growth occurred during protracted melt crystallization between c. 540 and 510 Ma. Melt crystallization was followed by the influx of aqueous, alkali‐rich fluids likely derived from melts crystallizing at depth. This late fluid flow led to retrogression of orthopyroxene, the observed outcrop pattern and to the textural and isotopic modification of monazite grains at c. 525–490 Ma.     

The Rock of Gibraltar re‐visited

Gibraltar, a 6km² peninsula jutting south from Spain at the western entrance to the Mediterranean Sea, is dominated by its 424 m‐high Rock: famous as a landmark to seafarers since ancient times. Twenty‐five years ago, an article in Geology Today (1991, v.7, pp.95–101) interpreted the Rock as a partly overturned mass of Early Jurassic dolomitic limestone, thrust into position during continent–continent collision about 15–20 million years before present and shaped finally by shoreline processes active during Quaternary tectonic uplift. Later articles featured the complex of tunnels and chambers excavated between 1782 and 1968 by British military engineers within its 2.6 km‐long main ridge (1992, v.8, pp.92–98), and the Neanderthals known to have used some of its caves (1997, v.13, pp.179–184). Significant aspects of the Rock's geological history have clarified since then.     

Recent progress on combining geomorphological and geochronological data with ice sheet modelling,demonstrated using the last British–Irish Ice Sheet

Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.     

A middle–late Ediacaran volcano‐sedimentary record from the eastern Arabian‐Nubian shield

The Ediacaran Jibalah Group comprises volcano‐sedimentary successions that filled small fault‐bound basins along the NW–SE‐trending Najd fault system in the eastern Arabian‐Nubian Shield. Like several other Jibalah basins, the Antaq basin contains exquisitely preserved sedimentary structures and felsic tuffs, and hence is an excellent candidate for calibrating late Ediacaran Earth history. Shallow‐marine strata from the upper Jibalah Group (Muraykhah Formation) contain a diversity of load structures and intimately related textured organic (microbial) surfaces, along with a fragment of a structure closely resembling an Ediacaran frond fossil and a possible specimen of Aspidella. Interspersed carbonate beds through the Muraykhah Formation record a positive δ¹³C shift from ?6 to 0‰. U‐Pb zircon geochronology indicates a maximum depositional age of ~570 Ma for the upper Jibalah Group, consistent with previous age estimates. Although this age overlaps with that of the upper Huqf Supergroup in nearby Oman, these sequences were deposited in contrasting tectonic settings on opposite sides of the final suture of the East African Orogen.     

Inferring air–water temperature relationships from river and catchment properties

Water temperature (Tw) is a key determinant of freshwater ecosystem status and cause for concern under a changing climate. Hence, there is growing interest in the feasibility of moderating rising Tw through management of riparian shade. The Loughborough University Temperature Network (LUTEN) is an array of 36 water and air temperature (Ta) monitoring sites in the English Peak District set‐up to explore the predictability of local Tw, given Ta, river reach, and catchment properties. Year 1 of monitoring shows that 84%–94% of variance in daily Tw is explained by Ta. However, site‐specific logistic regression parameters exhibit marked variation and dependency on upstream riparian shade. Perennial spring flows in the lower River Dove also affect regression model parameters and strongly buffer daily and seasonal mean Tw. The asymptote of the models (i.e. maximum expected Tw) is particularly sensitive to groundwater inputs. We conclude that reaches with spring flows potentially offer important thermal refuges for aquatic organisms against expected long‐term warming of rivers and should be afforded special protection. Copyright © 2013 John Wiley & Sons, Ltd.     

Engagement in the Third U.S. National Climate Assessment: commitment,capacity, and communication for impact

This paper examines the dynamics of energy investments and clean energy Research and Development (R&D) using a scenario-based modeling approach. Starting from the global scenarios proposed in the RoSE model ensemble experiment, we analyze the dynamics of investments under different assumptions regarding economic and population growth as well as availability of fossil fuel resources, in the absence of a climate policy. Our analysis indicates that economic growth and the speed of income convergence across countries matters for improvements in energy efficiency, both via dedicated R&D investments but mostly through capital-energy substitution. In contrast, fossil fuel prices, by changing the relative competitiveness of energy sources, create an economic opportunity for radical innovation in the energy sector. Indeed, our results suggest that fossil fuel availability is the key driver of investments in low carbon energy innovation. However, this innovation, by itself, is not sufficient to induce emission reductions compatible with climate stabilization objectives.     

Development of large‐scale seismites in Upper Cretaceous fluvial sandstones in a fault‐proximal setting

Large‐scale soft‐sediment deformation structures occur within fluvial sandstone bodies of the Upper Cretaceous Wahweap Formation in the Kaiparowits basin, southern Utah, USA. These structures represent an exceptional example of metre‐scale fault‐proximal, seismogenic load structures in nearly homogenous sandstones. The load structures consist of two types: large‐scale load casts and wedge‐shaped load structures. Large‐scale load casts penetrate up to 4·5 m into the underlying sandstone bed. Wedge‐shaped load structures include metre‐scale, parallel, sub‐vertical features and decimetre‐scale features along the periphery of the large‐scale load casts or other wedge‐shaped load structures. Wedge‐shaped load structures contain well‐developed, medial cataclastic shear deformation bands. All load structures contain pervasive well‐defined millimetre‐thick to centimetre‐thick internal laminae, oriented parallel to the outside form of the load structures and asymptotic to deformation bands. Both types of load structures formed because of an inverted density profile, earthquake‐triggered liquefaction and growth of irregularities (a Rayleigh–Taylor instability) on the sandstone–sandstone erosional contact. The internal laminae and deformation bands formed during deformation and clearly demonstrate polyphase deformation, recording a transition from liquefied to hydroplastic to brittle modes of deformation. Decimetre‐scale wedge‐shaped load structures on the edge of the large‐scale load casts probably formed towards the end of a seismic event after the sediment dewatered and increased the frictional contact of grains enough to impart strength to the sands. Metre‐scale wedge‐shaped load structures were created as the tips of downward foundering sediments were driven into fractures, which widened incrementally with seismic pulsation. With each widening of the fracture, gravity and a suction effect would draw additional sediment into the fracture. Superimposed laminae indicate a secondary syndeformational origin for internal laminae, probably by flow‐generated shearing and vibrofluidization mechanisms. Large‐scale and wedge‐shaped load structures, polyphase deformation and secondary laminae may characterize soft‐sediment deformation in certain fault‐proximal settings.     

Searching for Added Value in Simulating Climate Extremes with a High-Resolution Regional Climate Model over Western Canada

We evaluate the capacity of a regional climate model to represent observed extreme temperature and precipitation events and also examine the impact of increased resolution, in an effort to identify added value in this respect. Two climate simulations of western Canada (WCan) were conducted with the Canadian Regional Climate Model (version 4) at 15 (CRCM15) and 45?km (CRCM45) horizontal resolution driven at the lateral boundaries by data from the European Centre for Medium-range Weather Forecasts (ECMWF) 40-year Reanalysis (ERA-40) for the period 1973–1995. The simulations were evaluated using the spline-interpolated dataset ANUSPLIN, a daily observational gridded surface temperature and precipitation product with a nominal resolution of approximately 10?km. We examine a range of climate extremes, comprising the 10th and 90th percentiles of daily maximum (TX) and minimum (TN) temperatures, the 90th percentile of daily precipitation (PR90), and the 27 core Climate Daily Extremes (CLIMDEX) indices.

Both simulations exhibit cold biases compared with observations over WCan, with the bias exacerbated at higher resolution, suggesting little added value for temperature overall. There are instances, however, of regional improvement in the spatial pattern of temperature extremes at the higher resolution of CRCM15 (e.g., the CLIMDEX index for the annual number of days when TX?>?25°C). The high-resolution simulations also reveal similarly localized features in precipitation (e.g., rain shadows) that are not resolved at the 45?km resolution. With regard to precipitation extremes, although both simulations generally display wet biases, CRCM15 features a reduced bias in PR90 in all seasons except winter. This improvement occurs despite the fact that spatial and interannual variability of PR90 in CRCM15 is significantly overestimated relative to both CRCM45 and ANUSPLIN. We posit that these characteristics are the result of demonstrable differences between corresponding topographical datasets used in the gridded observations and CRCM, the resulting errors propagated to physical variables tied to elevation and the beneficial effect of subsequent spatial averaging. Because topographical input is often discordant between simulations and gridded observations, it is argued that a limited form of spatial averaging may contribute added value beyond that which has already been noted in previous studies with respect to small-scale climate variability.