Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.     

Lithium isotopes in the system Qz-Ms-fluid:An experimental study

The fractionation of lithium isotopes among quartz, muscovite, and a chloride-bearing aqueous fluid has been investigated experimentally at 400°-500°C and 50-100 MPa. Experiments were performed for 15-60 days in cold seal vessels with natural mineral specimens. Lithium was introduced primarily through the fluid, which also contained KCl and HCl. In most runs, the fluid was prepared with the L-SVEC standard (δ⁷Li = 0) and was 1 M in total chloride with K/Li/H = 100/10/1. In two experiments, a ⁶Li spike was employed. The experiments demonstrate that quartz and muscovite are susceptible to pronounced, rapid shifts in Li isotopic composition by diffusion through interaction with a Li-bearing fluid, particularly at 500°C. At 500°C, fractionation factors were determined to be Δ_{quartz-aqueous} ≅ +8 to + 12‰ and Δ_{muscovite-aqueous} ≅ +18 to + 20‰. An intermineral fractionation factor is given by Δ_{muscovite-quartz} ≈ +9‰. At 400°C, the results suggest Δ_{quartz-aqueous} ≈ +4 to + 6‰. The study provides evidence of systematic fractionation in lithium isotopes at the temperatures of some magmatic processes, such as those associated with porphyry-type ore systems and pegmatites.     

Importance of climate,forest fires and human population size in the Holocene boreal forest composition change in northern Europe

The relative importance of climate, forest fires and human population size on long‐term boreal forest composition were statistically investigated at regional and local scales in Fennoscandia. We employ pollen data from lakes, reflecting regional vegetation, and small forest hollows, reflecting local vegetation, from Russia, Finland and Sweden to reconstruct the long‐term forest composition. As potential drivers of the Holocene forest dynamics we consider climate, generated from a climate model and oxygen isotope data, past forest fires generated from sedimentary charcoal data and human population size derived from radiocarbon dated archaeological findings. We apply the statistical method of variation partitioning to assess the relative importance of these environmental variables on long‐term boreal forest composition. The results show that climate is the main driver of the changes in Holocene boreal forest composition at the regional scale. However, at the local scale the role of climate is relatively small. In general, the importance of forest fires is low both at regional and local scales. The fact that both climate and forest fires explain relatively small proportions of variation in long‐term boreal vegetation in small forest hollow records demonstrates the complexity of factors affecting stand‐scale forest dynamics. The relative importance of human population size was low in both the prehistorical and the historical time periods. However, this is the first time that this type of data has been used to statistically assess the importance of human population size on boreal vegetation and the spatial representativeness of the data may cause bias to the analysis.     

Landscape controls on long‐term runoff in subhumid heterogeneous Boreal Plains catchments

We compared median runoff (R) and precipitation (P) relationships over 25 years from 20 mesoscale (50 to 5,000 km²) catchments on the Boreal Plains, Alberta, Canada, to understand controls on water sink and source dynamics in water‐limited, low‐relief northern environments. Long‐term catchment R and runoff efficiency (RP^?1) were low and varied spatially by over an order of magnitude (3 to 119 mm/year, 1 to 27%). Intercatchment differences were not associated with small variations in climate. The partitioning of P into evapotranspiration (ET) and R instead reflected the interplay between underlying glacial deposit texture, overlying soil‐vegetation land cover, and regional slope. Correlation and principal component analyses results show that peatland‐swamp wetlands were the major source areas of water. The lowest estimates of median annual catchment ET (321 to 395 mm) and greatest R (60 to 119 mm, 13 to 27% of P) were observed in low‐relief, peatland‐swamp dominated catchments, within both fine‐textured clay‐plain and coarse‐textured glacial deposits. In contrast, open‐water wetlands and deciduous‐mixedwood forest land covers acted as water sinks, and less catchment R was observed with increases in proportional coverage of these land covers. In catchments dominated by hummocky moraines, long‐term runoff was restricted to 10 mm/year, or 2% of P. This reflects the poor surface‐drainage networks and slightly greater regional slope of the fine‐textured glacial deposit, coupled with the large soil‐water and depression storage and higher actual ET of associated shallow open‐water marsh wetland and deciduous‐forest land covers. This intercatchment study enhances current conceptual frameworks for predicting water yield in the Boreal Plains based on the sink and source functions of glacial landforms and soil‐vegetation land covers. It offers the capability within this hydro‐geoclimatic region to design reclaimed catchments with desired hydrological functionality and associated tolerances to climate or land‐use changes and inform land management decisions based on effective catchment‐scale conceptual understanding.     

Middle Paleolithic occupation on a Marine Isotope Stage 5 lakeshore in the Nefud Desert,Saudi Arabia

Major hydrological variations associated with glacial and interglacial climates in North Africa and the Levant have been related to Middle Paleolithic occupations and dispersals, but suitable archaeological sites to explore such relationships are rare on the Arabian Peninsula. Here we report the discovery of Middle Paleolithic assemblages in the Nefud Desert of northern Arabia associated with stratified deposits dated to 75,000 years ago. The site is located in close proximity to a substantial relict lake and indicates that Middle Paleolithic hominins penetrated deeply into the Arabian Peninsula to inhabit landscapes vegetated by grasses and some trees. Our discovery supports the hypothesis of range expansion by Middle Paleolithic populations into Arabia during the final humid phase of Marine Isotope Stage 5, when environmental conditions were still favorable.     

A framework of Holocene and Late Pleistocene environmental change in eastern Iran inferred from the dating of periods of alluvial fan abandonment,river terracing,and lake deposition

We review studies of the Holocene and Late Pleistocene stratigraphy of eastern Iran to infer past changes in the environment within this presently arid region. We build a scenario of widespread, and presumably climatically driven, evolution of the landscape through the Holocene. Six sites, covering a 10° range in latitude, indicate a regional abandonment of alluvial fan surfaces at ～10 ± 3 ka, with the younger (～9 ka) end of this age range supported by several of the best-constrained studies. Incision of rivers into the fan surfaces has occurred in discrete stages in the early to mid-Holocene (～9–7 ka) leading to the formation of flights of river terraces. Detailed records of lakebed deposition in the presently arid interior of Iran are rare, though the available data indicate lake highstand conditions at <7.8 ka at South Golbaf in SE Iran and at < 8.7 ± 1.1 ka at the Nimbluk plain in NE Iran. The major periods of Holocene landscape development hence correlate with a period of time where water was more abundant than at present, with incision of rivers into thick alluvial deposits possibly occurring due to a combination of decreased sediment supply and high levels of precipitation, and with the formation of inset river terraces possibly responding to century-scale fluctuations in precipitation. No major geomorphic changes are identified within the later part of the Holocene, from which we infer that increased aridity has slowed evolution of the landscape. A decrease in precipitation in the mid-Holocene may have had a detrimental effect on bronze age societies in eastern Iran as has been inferred elsewhere in the eastern Mediterranean region. The pre-Holocene environmental changes in eastern Iran are less well constrained, though there are suggestions of alluvial fan abandonment at 40–60 ka, at ～80 ka, and at ～120 ka.     

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 probabilistic approach for computing water retention of particulate systems from statistics of grain size and tessellated pore network

The paper offers an analytical determination of the hydraulic properties of an unsaturated soil with reference to its retention curve, which describes the relationship between the volumetric water content and capillarity through matric suction. The analysis combines a particulate approach focused on the physics at the pore scale, including microstructural aspects, with a probabilistic approach where the void space and grain size are considered as random variables. In the end, the soil water characteristic curve of an unsaturated granular medium along a drying path can be derived analytically based on the sole information of particle size distribution. The analysis hinges on the tessellation of a wet granular system into an assemblage of tetrahedral unit cells revealing a pore network upon which capillary physics are computed with respect to pore throat invasion by a non-wetting fluid with evolving pendular capillary bridges. The crux of the paper is to pass from particle size probability distribution to a matching void space distribution to eventually reveal key information such as void cell and solid volume statistics. Making reasonable statistically based assumptions to render calculations tractable, the water retention curve can be readily constructed. Model predictions compare quite favourably with experimental data available for actual soils, especially in the high saturation range. Having a sound scientific basis, the model can be made amenable to address a variety of soils with a wider range of particle sizes.