Testing the fidelity of thermometers at ultrahigh temperatures

A highly residual granulite facies rock (sample RG07‐21) from Lunnyj Island in the Rauer Group, East Antarctica, presents an opportunity to compare different approaches to constraining peak temperature in high‐grade metamorphic rocks. Sample RG07‐21 is a coarse‐grained pelitic migmatite composed of abundant garnet and orthopyroxene along with quartz, biotite, cordierite, and plagioclase with accessory rutile, ilmenite, zircon, and monazite. The inferred sequence of mineral growth is consistent with a clockwise pressure–temperature (P–T) evolution when compared with a forward model (P–T pseudosection) for the whole‐rock chemical composition. Peak metamorphic conditions are estimated at 9 ± 0.5 kbar and 910 ± 50°C based on conventional Al‐in‐orthopyroxene thermobarometry, Zr‐in‐rutile thermometry, and calculated compositional isopleths. U–Pb ages from zircon rims and neocrystallized monazite grains yield ages of c. 514 Ma, suggesting that crystallization of both minerals occurred towards the end of the youngest pervasive metamorphic episode in the region known as the Prydz Tectonic Event. The rare earth element compositions of zircon and garnet are consistent with equilibrium growth of these minerals in the presence of melt. When comparing the thermometry methods used in this study, it is apparent that the Al‐in‐orthopyroxene thermobarometer provides the most reliable estimate of peak conditions. There is a strong textural correlation between the temperatures obtained using the Zr‐in‐rutile thermometer––maximum temperatures are recorded by a single rutile grain included within orthopyroxene, whereas other grains included in garnet, orthopyroxene, quartz, and biotite yield a range of temperatures down to 820°C. Ti‐in‐zircon thermometry returns significantly lower temperature estimates of 678–841°C. Estimates at the upper end of this range are consistent with growth of zircon from crystallizing melt at temperatures close to the elevated (H₂O undersaturated) solidus. Those estimates, significantly lower than the calculated temperature of this residual solidus, may reflect isolation of rutile from the effective equilibration volume leading to an activity of TiO₂ that is lower than the assumed value of unity.     

Titanite dates crystallization: Slow Pb diffusion during super‐solidus re‐equilibration

Titanite can be found in rocks of wide compositional range, is reactive, growing or regrowing during metamorphic and hydrothermal events, and is generally amenable to U–Pb geochronology. Experimental evidence suggest that titanite has a closure temperature for Pb ranging from 550 to 650°C, and thus titanite dates are commonly interpreted as cooling ages. However, this view has been challenged in recent years by evidence from natural titanite which suggests the closure temperature may be significantly higher (up to 800°C). Here, we investigate titanite in an enclave of migmatitic gneiss included within a granite intrusion. The titanite crystals exhibit textural features characteristic of fluid‐mediated mass transfer processes on length scales of <100 µm. These textural features are associated with variation in both Pb concentrations and distinct U–Pb isotopic compositions. Zr‐in‐titanite thermometry indicates that modification of the titanite occurred at temperatures in excess of 840°C, in the presence of a high‐T silicate melt. The Pb concentration gradients preserved in these titanite crystals are used to determine the diffusivity of Pb in titanite under high‐T conditions. We estimate diffusivities ranging from 2 × 10^?22 to 5 × 10^?25 m²/s. These results are significantly lower than experimental data predict yet are consistent with other empirical data on natural titanites, suggesting that Pb diffusivity is similar to that of Sr. Thus our data challenge the wide‐held assumption that U–Pb titanite dates only reflect cooling ages.     

Petrochronology of oxidized granulites from southern Peru

Sapphirine–quartz granulites from the Cocachacra region of the Arequipa Massif in southern Peru record early Neoproterozoic ultrahigh‐temperature metamorphism. Phase equilibrium modelling and zircon petrochronology are used to quantify timing and pressure–temperature (P–T) conditions of metamorphism. Modelling of three magnetite‐bearing sapphirine–quartz samples indicates peak temperatures of >950°C at ~0.7 GPa and a clockwise P–T evolution. Elevated concentrations of Al in orthopyroxene are also consistent with ultrahigh‐temperature conditions. Neoblastic zircon records ages of c. 1.0–0.9 Ga that are interpreted to record protracted ultrahigh‐temperature metamorphism. Th/U ratios of zircon of up to 100 reflect U‐depleted whole‐rock compositions. Concentrations of heavy rare earth elements in zircon do not show systematic trends with U–Pb age but do correlate with variable whole‐rock compositions. Very large positive Ce anomalies in zircon from two samples probably relate to strongly oxidizing conditions during neoblastic zircon crystallization. Low concentrations of Ti‐in‐zircon (<10 ppm) are interpreted to result from reduced titania activities due to the strongly oxidized nature of the granulites and the sequestration of titanium‐rich minerals away from the reaction volume. Whole‐rock compositions and oxidation state have a strong influence on the trace element composition of metamorphic zircon, which has implications for interpreting the geological significance of ages retrieved from zircon in oxidized metamorphic rocks.     

Late-glacial and Holocene vegetation and climatic history of the Cass Basin, central South Island, New Zealand

Lithology, pollen, macrofossils, and stable carbon isotopes from an intermontane basin bog site in southern New Zealand provide a detailed late-glacial and early Holocene vegetation and climate record. Glacial retreat occurred before 17,000 cal yr B.P., and tundra-like grassland–shrubland occupied the basin shortly after. Between 16,500 and 14,600 cal yr B.P., a minor regional expansion of forest patches occurred in response to warming, but the basin remained in shrubland. Forest retreated between 14,600 and 13,600 cal yr B.P., at about the time of the Antarctic Cold Reversal. At 13,600 cal yr B.P., a steady progression from shrubland to tall podocarp forest began as the climate ameliorated. Tall, temperate podocarp trees replaced stress-tolerant shrubs and trees between 12,800 and 11,300 cal yr B.P., indicating sustained warming during the Younger Dryas Chronozone (YDC). Stable isotopes suggest increasing atmospheric humidity from 11,800 to 9300 cal yr B.P. Mild (annual temperatures at least 1°C higher than present), and moist conditions prevailed from 11,000 to 10,350 cal yr B.P. Cooler, more variable conditions followed, and podocarp forest was completely replaced by montane Nothofagus forest at around 7500 cal yr B.P. with the onset of the modern climate regime. The Cass Basin late-glacial climate record closely matches the Antarctic ice core records and is in approximate antiphase with the North Atlantic.     

Improved decentralized multi-sensor navigation system for airborne applications

The integration of Global Navigation Satellite System (GNSS) and Inertial Navigation System (INS) technologies is a very useful navigation option for high-accuracy positioning in many applications. However, its performance is still limited by GNSS satellite availability and satellite geometry. To address such limitations, a non-GNSS-based positioning technology known as “Locata” is used to augment a standard GNSS/INS system. The conventional methods for multi-sensor integration can be classified as being either in the form of centralized Kalman filtering (CKF), or decentralized Kalman filtering. However, these two filtering architectures are not always ideal for real-world applications. To satisfy both accuracy and reliability requirements, these three integration algorithms—CKF, federated Kalman filtering (FKF) and an improved decentralized filtering, known as global optimal filtering (GOF)—are investigated. In principle, the GOF is derived from more information resources than the CKF and FKF algorithms. These three algorithms are implemented in a GPS/Locata/INS integrated navigation system and evaluated using data obtained from a flight test. The experimental results show that the position, velocity and attitude solution derived from the GOF-based system indicate improvements of 30, 18.4 and 20.8% over the CKF- and FKF-based systems, respectively.     

Modelling tropical cyclone risks for present and future climate change scenarios using geospatial techniques

Tropical cyclones and their devastating impacts are of great concern to coastal communities globally. An appropriate approach integrating climate change scenarios at local scales is essential for producing detailed risk models to support cyclone mitigation measures. This study developed a simple cyclone risk-modelling approach under present and future climate change scenarios using geospatial techniques at local scales, and tested using a case study in Sarankhola Upazila from coastal Bangladesh. Linear storm-surge models were developed up to 100-year return periods. A local sea level rise scenario of 0.34?m for the year 2050 was integrated with surge models to assess the climate change impact. The resultant storm-surge models were used in the risk-modelling procedures. The developed risk models successfully identified the spatial extent and levels of risk that match with actual extent and levels within an acceptable limit of deviation. The result showed that cyclone risk areas increased with the increase of return period. The study also revealed that climate change scenario intensified the cyclone risk area by 5–10% in every return period. The findings indicate this approach has the potential to model cyclone risk in other similar coastal environments for developing mitigation plans and strategies.     

Modelling ecohydrological feedbacks in forest and grassland plots under a prolonged drought anomaly in Central Europe 2018–2020

Recent studies have highlighted the importance of understanding ecohydrological drought feedbacks to secure water resources under a changing climate and increasing anthropogenic impacts. In this study, we monitored and modelled feedbacks in the soil–plant-atmosphere continuum to the European drought summer 2018 and the following 2 years. The physically based, isotope-aided model EcH₂O-iso was applied to generic vegetation plots (forest and grassland) in the lowland, groundwater-dominated research catchment Demnitzer Millcreek (NE Germany; 66 km²). We included, inter alia, soil water isotope data in the model calibration and quantified changing “blue” (groundwater recharge) and “green” (evapotranspiration) water fluxes and ages under each land use as the drought progressed. Novel plant xylem isotope data were excluded from calibration but were compared with simulated root uptake signatures in model validation. Results indicated inter-site differences in the dynamics of soil water storage and fluxes with contrasting water age both during the drought and the subsequent 2 years. Forest vegetation consistently showed a greater moisture stress, more rapid recovery and higher variability in root water uptake depths from a generally younger soil water storage. In contrast, the grassland site, which had more water-retentive soils, showed higher and older soil water storage and groundwater recharge fluxes. The damped storage and flux dynamics under grassland led to a slower return to younger water ages at depth. Such evidence-based and quantitative differences in ecohydrological feedbacks to drought stress in contrasting soil-vegetation units provide important insights into Critical Zone water cycling. This can help inform future progress in the monitoring, modelling and development of climate mitigation strategies in drought-sensitive lowlands.     

How well do integrated assessment models simulate climate change?

Integrated assessment models (IAMs) are regularly used to evaluate different policies of future emissions reductions. Since the global costs associated with these policies are immense, it is vital that the uncertainties in IAMs are quantified and understood. We first demonstrate the significant spread in the climate system and carbon cycle components of several contemporary IAMs. We then examine these components in more detail to understand the causes of differences, comparing the results with more complex climate models and earth system models (ESMs), where available. Our results show that in most cases the outcomes of IAMs are within the range of the outcomes of complex models, but differences are large enough to matter for policy advice. There are areas where IAMs would benefit from improvements (e.g. climate sensitivity, inertia in climate response, carbon cycle feedbacks). In some cases, additional climate model experiments are needed to be able to tune some of these improvements. This will require better communication between the IAM and ESM development communities.