The effect of Ca-Tschermaks component on trace element partitioning between clinopyroxene and silicate melt

We have studied the influence of Ca-Tschermaks (Calcium Tschermaks or CaTs) content of clinopyroxene on the partitioning of trace elements between this phase and silicate melt at fixed temperature and pressure. Ion probe analyses of experiments carried out in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂, at 0.1 MPa and 1218°C, produced crystal-melt partition coefficients (D) of 36 trace elements (Li, Cl, Sc, Ti, V, Cr, Fe, Co, Ge, Sr, Y, Zr, Nb, Mo, Ru, Rh, In, Sn, Sb, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta and W), for clinopyroxene compositions between 10 and 32 mol% CaTs. Partition coefficients for 2+ to 5+ cations show, for each charge, a near parabolic dependence of log D on ionic radius of the substituting cation, for partitioning into both the M1 and M2 sites of clinopyroxene. Fitting the results to the elastic strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal-melt partition coefficients from elastic moduli. Nature 372, 452–454] we obtain results for the strain-free partition coefficients of theoretical cations (D₀), with site radius r₀, and for the site's Young's Modulus (E).

In agreement with earlier data our results show that increasing ^ivAl concentration in cpx is matched by increasing D, E_M1, E_M2 and D₀ for tri-, tetra- and pentavalent cations. The degree of fractionation between chemically similar elements (i.e. Ta/Nb, Zr/Hf) also increases. In contrast, D values for mono-, di- and hexavalent cations decrease with increasing ^ivAl in the cpx. The large suite of trace elements used has allowed us to study the effects of cation charge on D₀, r₀ and E. We have found that D₀ and r₀ decrease with increasing cation charge, e.g. r₀=0.66 Å for 4+ cations and 0.59 Å for 5+ cations substituting into M1. Values of E_M1 and E_M2 increase with cation charge as well as with increasing ^ivAl content. The increase in E_M2 is linear and close to the trend set by Hazen and Finger [Hazen, R.M., Finger, L.W., 1979. Bulk modulus-volume relationship for cation–anion polyhedra. J. Geophys. Res. 84 (10) 6723–6728] for oxides. E_M1 values are much higher and do not fit the trend predicted by the Hazen and Finger relationship.     

The role of closed ecological systems in carbon cycle modelling

Acquiring a mechanistic understanding of the role of biotic feedbacks for the links between atmospheric CO₂ concentrations and temperature is essential for trustworthy climate predictions. Currently, computer-based simulations are the only available tool to estimate the global impact of biotic feedbacks on future atmospheric CO₂ and temperatures. Here we propose an alternative and complementary approach by using materially closed, energetically open analogue/physical models of the carbon cycle. We argue that there is unexplored potential in using a materially closed approach to improve our understanding of the magnitude and direction of many biotic carbon feedbacks and that recent technological advances make this feasible. We also suggest how such systems could be designed and discuss the advantages and limitations of establishing physical models of the global carbon cycle.     

Tourmaline mineralization in the Barberton greenstone belt,South Africa: early Archean metasomatism by evaporite-derived boron

Tourmaline-rich rocks are common in the lowgrade, interior portions of the Barberton greenstone belt of South Africa, where shallow-marine sediments and underlying altered basaltic and komatiitic lavas contain up to 50% tourmaline. The presence of tourmaline-bearing rip-up clasts, intraformational tourmalinite pebbles, and tourmaline-coated grains indicates that boron mineralization was a low-temperature, surficial process. The association of these lithologies with stromatolites, evaporites, and shallow-water sedimentary structures and the virtual absence of tourmaline in correlative deep-water facies rocks in the greenstone bels strengthens this model.Five tourmaline-bearing lithologic groups (basalts, komatiites, evaporite-bearing sediments, stromatolitic sediments, and quartz veins) are distinguished based on field, petrographic, and geochemical criteria. Individual tourmaline crystals within these lithologies show internal chemical and textural variations that reflect continued growth through intervals of change in bulk-rock and fluid composition accompanying one or more metasomatic events. Large single-crystal variations exist in Fe/Mg, Al/Fe, and alkali-site vacancies. A wide range in tourmaline composition exists in rocks altered from similar protoliths, but tourmalines in sediments and lavas have similar compositional variations. Boron-isotope analysis of the tourmalines suggest that the boron enrichment in these rocks has a major marine evaporitic component. Sediments with gypsum pseudomorphs and lavas altered at low temperatures by shallow-level brines have the highest ¹¹B values (+2.2 to-1.9); lower ¹¹B values of late quartz veins (-3.7 to-5.7) reflect intermediate temperature, hydrothermal remobilization of evaporitic boron. The ¹¹B values of tourmaline-rich stromatolitic sediments (-9.8 and-10.5) are consistent with two-stage boron enrichment, in which earlier marine evaporitic boron was hydrothermally remobilized and vented in shallow-marine or subaerial sites, mineralizing algal stromatolites. The stromatolite-forming algae preferentially may have lived near the sites of hydrothermal discharge in Archean times.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

α-Emitting mineral inclusions in apatite, their effect on (U-Th)/He ages, and how to reduce it

U-Th rich mineral inclusions in apatite are often held responsible for erroneously old (U-Th)/He ages, because they produce “parentless” He. Three aspects associated with this problem are discussed here. First, simple dimensional considerations indicate that for small mineral inclusions, the parentless helium problem might not be as serious as generally thought. For example, a mineral inclusion that is 10% the length, width and height of its host apatite needs to be a thousand times more concentrated in U and Th to produce an equal amount of He. Therefore, single isolated inclusions smaller than a few μm are unlikely to contribute significant helium. For larger or more abundant inclusions, the parentless helium problem can be solved by dissolution of the apatite and its inclusions in hot HF. Second, besides creating parentless helium, inclusions also complicate α-ejection corrections. Mathematical exploration of this latter problem for spherical geometries reveals that for randomly distributed inclusions, the probability distribution of single-grain ages is predicted to have a sharp mode at the mean age, with tails towards younger and older ages. Multiple-grain measurements will yield accurate and precise age estimates if 10 or more randomly distributed α-emitting mineral inclusions are present in a sample. Third, thermal modeling indicates that mineral inclusions have a non-trivial but minor (<5 °C) effect on the closure temperature. These predictions were tested on apatites from rapidly cooled migmatites of Naxos (Greece) which contain abundant U-rich zircon inclusions. Thirty-seven samples were subjected to two kinds of treatment. The “pooled” age (i.e., the synthetic multi-grain age computed from a number of single-grain analyses) of four inclusion-free samples (13 apatites), prepared in HNO₃ is 10.9 Ma, close to apatite and zircon fission-track ages from the same rock. (U-Th)/He ages of 14 inclusion-bearing samples dissolved in HNO₃ range between 9 and 45 Ma, with a pooled age of 22.6 Ma. The ages of 19 HF-treated samples range between 5 and 16 Ma, with 10 of 14 single-grain samples between 9 and 13 Ma and a pooled age of 10.9 Ma. These observations agree with the theoretical predictions and support the addition of HF-treated apatite (U-Th)/He dating to the thermochronological toolbox.     

The 2020 glacial lake outburst flood process chain at Lake Salkantaycocha (Cordillera Vilcabamba,Peru)

Glacial lakes represent a threat for the populations of the Andes and numerous disastrous glacial lake outburst floods (GLOFs) occurred as a result of sudden dam failures or dam overtoppings triggered by landslides such as rock/ice avalanches into the lake. This paper investigates a landslide-triggered GLOF process chain that occurred on February 23, 2020, in the Cordillera Vilcabamba in the Peruvian Andes. An initial slide at the SW slope of Nevado Salkantay evolved into a rock/ice avalanche. The frontal part of this avalanche impacted the moraine-dammed Lake Salkantaycocha, triggering a displacement wave which overtopped and surficially eroded the dam. Dam overtopping resulted in a far-reaching GLOF causing fatalities and people missing in the valley downstream. We analyze the situations before and after the event as well as the dynamics of the upper portion of the GLOF process chain, based on field investigations, remotely sensed data, meteorological data and a computer simulation with a two-phase flow model. Comparison of pre- and post-event field photographs helped us to estimate the initial landslide volume of 1–2 million m³. Meteorological data suggest rainfall and/or melting/thawing processes as possible causes of the landslide. The simulation reveals that the landslide into the lake created a displacement wave of 27 m height. The GLOF peak discharge at the dam reached almost 10,000 m³/s. However, due to the high freeboard, less than 10% of the lake volume drained, and the lake level increased by 10–15 m, since the volume of landslide material deposited in the lake (roughly 1.3 million m³) was much larger than the volume of released water (57,000 m³, according to the simulation). The model results show a good fit with the observations, including the travel time to the uppermost village. The findings of this study serve as a contribution to the understanding of landslide-triggered GLOFs in changing high-mountain regions.

     

The identification and interaction of network,intranetwork, and ephemeral-region magnetic fields

Network magnetic fields, ephemeral active regions, and intranetwork magnetic fields are illustrated and discussed in several contexts. First, they are presented in relation to the appearance and disappearance of magnetic flux. Second, their properties in common with all solar magnetic features are discussed. Third, their distinguishing characteristics are emphasized. Lastly, their interactions are illustrated.Network magnetic fields are no longer considered to be just the aged remnants of active regions. The network is the dynamic product of the merging and cancelling of intranetwork fields, ephemeral regions, and the remnants of active regions. Intranetwork fields are magnetic fields of mixed polarity that appear to originate continuously from localized source sites in between the network. The intranetwork magnetic fields are characterized by flow of successive fragments in approximately radial patterns away from their apparent source sites and by the relative weakness of their magnetic fields. Ephemeral active regions are small, new bipoles that grow as a unit or a succession of bipolar units and whose poles move in opposite directions from their apparent site of origin. Large ephemeral regions are not distinguishable from small active regions.Solar Cycle Workshop Paper.     

Particle composition and size distributions in and around a deep-pit swine operation,Ames, IA

The contribution of emissions from agricultural facilities is rapidly becoming a major concern for local and regional air quality. Characterization of particle properties such as physical size distribution and chemical composition can be valuable in understanding the processes contributing to emissions and ultimate fate of particulate matter from agricultural facilities. A measurement campaign was conducted at an Iowa, deep-pit, three-barn swine finishing facility to characterize near-source ambient particulate matter. Size-specific mass concentrations were determined using minivol samplers, with additional size distribution information obtain using optical particle counters. Particulate composition was determined via ion chromatographic analysis of the collected filters. A thermal-CO₂ elemental/organic carbon analyzer measured particulate carbon. The chemical composition and size distribution of sub-micron particles were determined via real-time aerosol mass spectrometry. Primary particulate was not found to be a major emission from the examined facility, with filter-based impactor samples showing average near-source increases (~15–50 m) in ambient PM₁₀ of 5.8 ± 2.9 μg m⁻³ above background levels. PM_2.5 also showed contribution attributable to the facility (1.7 ± 1.1 μg m⁻³). Optical particle counter analysis of the numerical size distributions showed bimodal distributions for both the upwind and downwind conditions, with maximums around 2.5 μm and below the minimum quantified diameter of 0.3 μm. The distributions showed increased numbers of coarse particles (PM₁₀) during periods when wind transport came from the barns, but the differences were not statistically significant at the 95% confidence level. The PM₁₀ aerosols showed statistically increased concentrations of sulfate, nitrate, ammonium, calcium, organic carbon, and elemental carbon when the samplers were downwind from the pig barns. Organic carbon was the major constituent of the barn-impacted particulate matter in both sub-micron (54%) and coarse size (20%) ranges. The AMS PM₁ chemical speciation showed similar species increases, with the exception of and Ca⁺², the latter not quantified by the AMS.