Be-daughter minerals in fluid and melt inclusions: implications for the enrichment of Be in granite–pegmatite systems

The study of re-homogenized melt inclusions in the same growth planes of quartz of pegmatites genetically linked to the Variscan granite of the Ehrenfriedersdorf complex, Erzgebirge, Germany, by ion microprobe analyses has determined high concentrations of Be, up to 10,000 ppm, in one type of melt inclusion, as well as moderate concentrations in the 100 ppm range in a second type of melt inclusion. Generally, the high Be concentrations are associated with the H₂O- and other volatile-rich type-B melt inclusions, and the lower Be concentration levels are connected to H₂O-poor type-A melt inclusions. Both inclusion types, representing conjugate melt pairs, are formed by a liquid–liquid immiscibility separation process. This extremely strong and very systematic scattering in Be provides insights into the origin of Be concentration and transport mechanisms in pegmatite-forming melts. In this contribution, we present more than 250 new analytical data and show with ion microprobe and fs-LA-ICPMS studies on quenched glasses, as well as with confocal Raman spectroscopy of daughter minerals in unheated melt inclusions, that the concentrations of Be may achieve such extreme levels during melt–melt immiscibility of H₂O-, B-, F-, P-, ± Li-enriched pegmatite-forming magmas. Starting from host granite with about 10 ppm Be, melt inclusions with 10,000 ppm Be correspond to enrichment by a factor of over 1,000. This strong enrichment of Be is the result of processes of fractional crystallization and further enrichment in melt patches of pegmatite bodies due to melt–melt immiscibility at fluid saturation. We also draw additional conclusions regarding the speciation of Be in pegmatite-forming melt systems from investigation of the Be-bearing daughter mineral phases in the most H₂O-rich melt inclusions. In the case of evolved volatile and H₂O-rich pegmatite systems, B, P, and carbonates are important for the enrichment and formation of stable Be complexes.     

Automated AF-demagnetization on the 2G-Enterprises through-bore,cryogenic magnetometer

Since the early sixties, alternating field demagnetization (AFD) has been a standard laboratory technique for demagnetizing rocks to expose the multicomponent structure of their natural remanent magnetization (NRM). In the majority of AFD implementations, however, the procedure remains as labour-intensive as ever. The implementation that we have developed at the Australian Geological Survey Organisation, automates the procedure for AFD based on the static method, and results in significant productivity and efficiency gains without compromising data quality. A properly formulated procedure for static AFD may be the only method of retrieving higher-coercivity components of natural remanence in samples prone to developing gyroremanence at higher alternating fields (AFs). Our AFD environment comprises: a 2G-Enterprises through-bore, cryogenic magnetometer; 2G AF-coils and control equipment; and personal computer software, developed by us, to control all procedural aspects for a complete AFD of a sample including, importantly, a counteracting procedure to neutralize the effects of gyroremanence build-up at higher AFs. With our system, AFD of 8 samples/day, each of 20+ steps, requires only 20 min of user attention compared with a full day for conventional systems.     

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.     

High-temperature emplacement of the Cerro Galán and Toconquis Group ignimbrites (Puna plateau,NW Argentina) determined by TRM analyses

Estimates of pyroclastic flow emplacement temperatures in the Cerro Galán ignimbrite and Toconquis Group ignimbrites were determined using thermal remanent magnetization of lithic clasts embedded within the deposits. These ignimbrites belong to the Cerro Galán volcanic system, one of the largest calderas in the world, in the Puna plateau, NW Argentina. Temperature estimates for the 2.08-Ma Cerro Galán ignimbrite are retrieved from 40 sites in 14 localities (176 measured clasts), distributed at different distances from the caldera and different stratigraphic heights. Additionally, temperature estimates were obtained from 27 sample sites (125 measured clasts) from seven ignimbrite units forming the older Toconquis Group (5.60–4.51 Ma), mainly outcropping along a type section at Rio Las Pitas, Vega Real Grande. The paleomagnetic data obtained by progressive thermal demagnetization show that the clasts of the Cerro Galán ignimbrite have one single magnetic component, oriented close to the expected geomagnetic field at the time of emplacement. Results show therefore that most of the clasts acquired a new magnetization oriented parallel to the magnetic field at the moment of the ignimbrite deposition, suggesting that the clasts were heated up to or above the highest blocking temperature (T _b) of the magnetic minerals (T _b = 580°C for magnetite; T _b = 600–630°C for hematite). We obtained similar emplacement temperature estimations for six out of the seven volcanic units belonging to the Toconquis Group, with the exception of one unit (Lower Merihuaca), where we found two distinct magnetic components. The estimation of emplacement temperatures in this latter case is constrained at 580–610°C, which are lower than the other ignimbrites. These estimations are also in agreement with the lowest pre-eruptive magma temperatures calculated for the same unit (i.e., 790°C; hornblende–plagioclase thermometer; Folkes et al. 2011b). We conclude that the Cerro Galán ignimbrite and Toconquis Group ignimbrites were emplaced at temperatures equal to or higher than 620°C, except for Lower Merihuaca unit emplaced at lower temperatures. The homogeneity of high temperatures from proximal to distal facies in the Cerro Galán ignimbrite provides constraints for the emplacement model, marked by a relatively low eruption column, low levels of turbulence, air entrainment, surface–water interaction, and a high level of topographic confinement, all ensuring minimal heat loss.     

The Entrainment Rate for Buoyant Plumes in a Crossflow

We consider large-eddy simulations (LES) of buoyant plumes from a circular source with initial buoyancy flux F ₀ released into a stratified environment with constant buoyancy frequency N and a uniform crossflow with velocity U. We make a systematic comparison of the LES results with the mathematical theory of plumes in a crossflow. We pay particular attention to the limits [(U)\tilde] << 1{\tilde{U}\ll1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}, where [(U)\tilde]=U/(F₀ N)^1/4{\tilde{U}=U/(F_0 N)^{1/4}}, for which analytical results are possible. For [(U)\tilde] >> 1{\tilde{U}\gg 1}, the LES results show good agreement with the well-known two-thirds law for the rise in height of the plume. Sufficiently far above the source, the centreline vertical velocity of the LES plumes is consistent with the analytical z ^−1/3 and z ^−1/2 scalings for respectively [(U)\tilde] << 1{\tilde{U}\ll 1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}. In the general case, where the entrainment is assumed to be the sum of the contributions from the horizontal and vertical velocity components, we find that the discrepancy between the LES data and numerical solutions of the plume equations is largest for [(U)\tilde]=O(1){\tilde{U}=O(1)}. We propose a modified additive entrainment assumption in which the contributions from the horizontal and vertical velocity components are not equally weighted. We test this against observations of the plume generated by the Buncefield fire in the U.K. in December 2005 and find that the results compare favourably. We also show that the oscillations of the plume as it settles down to its final rise height may be attenuated by the radiation of gravity waves. For [(U)\tilde] << 1{\tilde{U}\ll 1} the oscillations decay rapidly due to the transport of energy away from the plume by gravity waves. For ${\tilde{U}>rsim 1}${\tilde{U}>rsim 1} the gravity waves travel in the same direction and at the same speed as the flow. In this case, the oscillations of the plume do not decay greatly by radiation of gravity waves.     

Non-uniform interhemispheric temperature trends over the past 550 years

The warming trend over the last century in the northern hemisphere (NH) was interrupted by cooling from ad 1940 to 1975, a period during which the southern hemisphere experienced pronounced warming. The cause of these departures from steady warming at multidecadal timescales are unclear; the prevailing explanation is that they are driven by non-uniformity in external forcings but recent models suggest internal climate drivers may play a key role. Paleoclimate datasets can help provide a long-term perspective. Here we use tree-rings to reconstruct New Zealand mean annual temperature over the last 550 years and demonstrate that this has frequently cycled out-of-phase with NH mean annual temperature at a periodicity of around 30–60 years. Hence, observed multidecadal fluctuations around the recent warming trend have precedents in the past, strongly implicating natural climate variation as their cause. We consider the implications of these changes in understanding and modelling future climate change.     

Evolution and modulation of tropical heating from the last glacial maximum through the twenty-first century

Twentieth century observations show that during the last 50?years the sea-surface temperature (SST) of the tropical oceans has increased by ~0.5°C and the area of SST >26.5 and 28°C (arbitrarily referred to as the oceanic warm pool: OWP) by 15 and 50% respectively in association with an increase in green house gas concentrations, with non-understood natural variability or a combination of both. Based on CMIP3 projections the OWP is projected to double during twenty-first century in a moderate CO₂ forcing scenario (IPCC A1B scenario). However, during the observational period the area of positive atmospheric heating (referred to as the dynamic warm pool, DWP), has remained constant. The threshold SST (T _H ), which demarks the region of net heating and cooling, has increased from 26.6°C in the 1950s to 27.1°C in the last decade and it is projected to increase to ~28.5°C by 2100. Based on climate model simulations, the area of the DWP is projected to remain constant during the twenty-first century. Analysis of the paleoclimate model intercomparison project (PMIP I and II) simulations for the Last Glacial maximum and the Mid-Holocene periods show a very similar behaviour, with a larger OWP in periods of elevated tropical SST, and an almost constant DWP associated with a varying T _H . The constancy of the DWP area, despite shifts in the background SST, is shown to be the result of a near exact matching between increases in the integrated convective heating within the DWP and the integrated radiative cooling outside the DWP as SST changes. Although the area of the DWP remains constant, the total tropical atmospheric heating is a strong function of the SST. For example the net heating has increased by about 10% from 1950 to 2000 and it is projected to increase by a further 20% by 2100. Such changes must be compensated by a more vigorous atmospheric circulation, with growth in convective heating within the warm pool, and an increase of subsiding air and stability outside the convective warm pool and an increase of vertical shear at the DWP boundaries. This finding is contrary to some conclusions from other studies but in accord with others. We discuss the similarities and differences at length.