Overstepping the garnet isograd: a comparison of QuiG barometry and thermodynamic modeling

The consequences of overstepping the garnet isograd reaction have been investigated by comparing the composition of garnet formed at overstepped P–T conditions (the overstep or “OS” model) with the P–T conditions that would be inferred by assuming garnet nucleated in equilibrium with the matrix assemblage at the isograd (the equilibrium or “EQ” model). The garnet nucleus composition formed at overstepped conditions is calculated as the composition that produces the maximum decrease in Gibbs free energy from the equilibrated, garnet-absent, matrix assemblage for the bulk composition under study. Isopleths were then calculated for this garnet nucleus composition assuming equilibrium with the matrix assemblage (the EQ model). Comparison of the actual P–T conditions of nucleation (the OS model) with those inferred from the EQ model reveals considerable discrepancy between the two. In general, the inferred garnet nucleation P–T conditions (the EQ model) are at a lower temperature and higher or lower pressure (depending on the coexisting calcic phase(s)) than the actual (OS model) nucleation conditions. Moreover, the degree of discrepancy increases with the degree of overstepping. Independent estimates of the pressure of nucleation of garnet were made using the Raman shift of quartz inclusions in garnet (quartz-in-garnet or QuiG barometry). To test the validity of this method, an experimental synthesis of garnet containing quartz inclusions was made at 800 °C, 20 kbar, and the measured Raman shift reproduced the synthesis conditions to within 120 bars. Raman band shifts from three natural samples were then used to calculate an isochore along which garnet was presumed to have nucleated. Model calculations were made at several temperatures along this isochore (the OS model), and these P–T conditions were compared to those computed assuming equilibrium nucleation (the EQ model) to estimate the degree of overstepping displayed by these samples. A sample from the garnet isograd in eastern Vermont is consistent with overstepping of around 10 degrees and 0.6 kbar (affinities of around 2 kJ/mole garnet). A sample from the staurolite–kyanite zone in the same terrane requires overstepping of around 50 °C and 2–5 kbar (affinities of around 10–18 kJ/mole garnet). A similar amount of overstepping was inferred for a blueschist sample from Sifnos, Greece. These results indicate that overstepping of garnet nucleation reactions may be common and pronounced in regionally metamorphosed terranes, and that the P–T conditions and paths inferred from garnet zoning studies may be egregiously in error.     

Thermal and rheological controls on the formation of mafic enclaves or banded pumice

Magma mixing can occur in a fluid manner to produce banded pumice or in a brittle manner to form enclaves. We propose that the critical control on mixing style is a competition between developing networks of crystals in the intruding magma that impart a strength to the magma and melting and disrupting those networks in the host. X-ray computed tomography analysis demonstrates that banded pumice from the 1915 Mt. Lassen eruption lacks crystal networks. In contrast, rhyodacite hosts with mafic enclaves from Chaos Crags contain well-developed networks of large crystals. We present a one-dimensional conductive cooling model that predicts mixing style, either ductile or brittle, as a function of magma compositions, temperatures, and the size of the intruding dike. Our model relies on three assumptions: (1) Mixing is initiated by the injection of a hot dike into a cooler magma body with a yield strength; (2) when magma crystallinity exceeds a critical value, 13 vol% plagioclase, the magma develops a yield strength; and (3) when total crystallinity exceeds 40 vol%, the magma has a penetrative crystal network and is effectively solid. Importantly, because the two magmas are of different compositions, their crystallinities and viscosities do not have the same variations with temperature. As the intruding magma cools, it crystallizes from the outside in, while simultaneously, host magma temperature near the intruder rises. Mixing of the two magmas begins when the host magma is heated sufficiently to (1) disrupt the crystal network and (2) initiate convection. If the shear stress exerted by the convecting host magma on the dike is greater than the yield strength of the dike margin (and dike crystallinity does not exceed 40 %), then fluid mixing occurs, otherwise enclaves form by brittle deformation of the dike. Application of the model to magma compositions representative of Lassen and Chaos Crags shows that emplacement of dikes <1 m thick should produce enclaves, whereas thicker dikes should generate fluid mixing and form banded pumice within days to weeks of emplacement. Similar relationships apply to other modeled magmatic systems, including Pinatubo, Unzen, and Ksudach/Shtuybel’ volcanoes. For all studied systems, the absolute size of the intruding dike, not just its proportion relative to the host, influences mixing style.     

Orbit Determination of Binary Asteroids

In addition to the detection of an asteroid moon or a binary asteroid, the knowledge of the satellite’s true orbit is of high importance to derive fundamental physical parameters of the binary system such as its mass and to shed light on its possible formation history and dynamical evolution (prograde/retrograde orbit, large/small eccentricity or inclination, etc.). A new methodology for preliminary orbit determination of binary asteroids – and visual binaries in general – is proposed. It is based on Thiele–Innes method combined with a ‘trial and error’ Monte-Carlo technique. This method provides the full set of solutions (bundle of orbits, with the 7 orbital elements) even for a reduced number of observations. The mass is a direct by-product of this orbit determination, from which one can next infer the bulk-density and porosity. In addition to the bundle of orbits, the method provides the marginal probability densities of the foreseen parameters. Such error analysis – since it avoids linear approximation – can be of importance for the prediction of the satellite’s position in the plane-of-sky during future stellar occultations or subsequent observations, but also for the analysis of the orbit’s secular evolution. After briefly describing the method, we present the algorithm and its application to some practical cases, with particular emphasis on asteroids binaries and applications on orbital evolution.     

Activity-composition relations in the system CaCO3-MgCO3 predicted from static structure energy calculations and Monte Carlo simulations

Thermodynamic mixing properties and subsolidus phase relations of the rhombohedral carbonate system, (1 − x) · CaCO₃ − x · MgCO₃, were modelled in the temperature range of 623-2023 K with static structure energy calculations based on well-parameterised empirical interatomic potentials. Relaxed static structure energies of a large set of randomly varied structures in a 4 × 4 × 1 supercell of calcite (a = 19.952 Å, c = 17.061 Å) were calculated with the General Utility Lattice Program (GULP). These energies were cluster expanded in a basis set of 12 pair-wise effective interactions. Temperature-dependent enthalpies of mixing were calculated by the Monte Carlo method. Free energies of mixing were obtained by thermodynamic integration of the Monte Carlo results. The calculated phase diagram is in good agreement with experimental phase boundaries.     

Simulated effects of land immersion on regional arid climate: a case study of the pre-Saharan playa of Chott el-Jerid (south of Tunisia)

The potential effects on the regional climate induced by partially immersing the arid pre-Saharan playa basin of Chott el-Jerid (south of Tunisia) are investigated by comparing two multi-year (1991–2011) sets of numerical simulations each consisting of ten-member ensemble and performed using the WRF regional climate model. The first WRF ensemble is performed under current land use and land cover, while the second is carried out after introducing a virtual large and shallow surface water reservoir (a lake) in Chott el-Jerid. The most pronounced effects generated by the artificial lake are circumscribed over its surface and slightly spread downwind to the other parts of the Chott. The lake has a clear moderating effect on near-surface air temperatures by increasing (decreasing) the wintertime (summertime) air temperatures. Sensible heat fluxes are remarkably increased in winter and decreased in summer over the lake following the temperature gradient between the lake surface and the overlying atmosphere. Latent heat fluxes, moisture convergence, and water vapor mixing ratio are increased over the lake throughout the year, especially in winter. The lake also induces domain-wide decreased (increased) surface pressures and land (lake) breeze circulation in winter (summer). Simulated rainfall amounts are most increased over the lake in winter likely because of an enhanced atmospheric instability, while they slightly decrease in summer.     

Correction to: Simulated effects of land immersion on regional arid climate: a case study of the pre-Saharan playa of Chott el-Jerid (South of Tunisia)

Theoretical and Applied Climatology - In the original article, the last panel of Fig. 10 should be labeled “h”. Also some data and Table 2 caption were unfortunately incorrect.     

Hiatus of wetland methane emissions associated with recent La Niña episodes in the Asian monsoon region

The variability of methane emissions from wetlands in the tropics and northern temperate regions can explain more than 70% of the interannual variation in global wetland methane emissions, which are largely driven by climate variability. We use climate reanalysis, remote sensing wetland area dataset and simulations from 11 land models contributing to Global Methane Budget to investigate the interannual variation and anomalies of wetland methane emissions in the Asian Monsoon region. Methane emissions in this region steadily increased over 2000–2012. However, abnormally low methane emissions were found in equatorial fully humid (Af), warm temperate winter dry (Cw), and warm temperate fully humid (Cf) Asian Monsoon climate sub-regions in 2008, 2009 and 2011, respectively. These spatially-shifting low emissions occurred simultaneously with observed wetland area shrinkage due to abnormally low precipitation. Interannual variability of wetland methane emissions in Asian Monsoon region are primarily driven by South Asian Monsoon system. However, the abnormally low emissions are related to strong La Niña events, and its accompanying effect of weakened East Asian Monsoon system and eastward Western Pacific subtropical high, which drives the shifting pattern of rainfall, and thus the spatial pattern of methane emission anomalies.     

Heat waves in the United States: definitions, patterns and trends

High temperatures and heat waves are related but not synonymous concepts. Heat waves, generally understood to be acute periods of extreme warmth, are relevant to a wide range of stakeholders because of the impacts that these events have on human health and activities and on natural environments. Perhaps because of the diversity of communities engaged in heat wave monitoring and research, there is no single, standard definition of a heat wave. Experts differ in which threshold values (absolute versus relative), duration and ancillary variables to incorporate into heat wave definitions. While there is value in this diversity of perspectives, the lack of a unified index can cause confusion when discussing patterns, trends, and impacts. Here, we use data from the North American Land Data Assimilation System to examine patterns and trends in 15 previously published heat wave indices for the period 1979–2011 across the Continental United States. Over this period the Southeast region saw the highest number of heat wave days for the majority of indices considered. Positive trends (increases in number of heat wave days per year) were greatest in the Southeast and Great Plains regions, where more than 12 % of the land area experienced significant increases in the number of heat wave days per year for the majority of heat wave indices. Significant negative trends were relatively rare, but were found in portions of the Southwest, Northwest, and Great Plains.