The effect of water on accessory phase solubility in subaluminous and peralkaline granitic melts

The solubilities of columbite, tantalite, wolframite, rutile, zircon and hafnon were determined as a function of the water contents in peralkaline and subaluminous granite melts. All experiments were conducted at 1035 °C and 2 kbar and the water contents of the melts ranged from nominally dry to approximately 6 wt.% H₂O. Accessory phase solubilities are not affected by the water content of the peralkaline melt. By contrast, solubilities are affected by the water content of the subaluminous melt, where the solubilities of all the accessory phases examined increase with the water content of the melt, up to 2 wt.% H₂O. At higher water contents, solubilities are nearly constant. It can be concluded that water is not an important control of accessory phase solubility, although the water content will affect diffusivities of components in the melt, thus whether or not accessory phases will be present as restite material. The solubility behaviour in the subaluminous and peralkaline melts supports previous spectroscopic studies, which have observed differences in the coordination of high field strength elements in dry vs. wet subaluminous granitic glasses, but not for peralkaline granitic glasses. Lastly, the fact that wolframite solubility increases with increasing water content in the subaluminous melt suggests that tungsten dissolved as a hexavalent species.     

Tidally driven stress accumulation and shear failure of Enceladus's tiger stripes

Straddling the south polar region of Saturn's moon Enceladus, the four principal “tiger stripe” fractures are a likely source of tectonic activity and plume generation. Here we investigate tidally driven stress conditions at the tiger stripe fractures through a combined analysis of shear and normal diurnal tidal stresses and accounting for additional stress at depth due to the overburden pressure. We compute Coulomb failure conditions to assess failure location, timing, and direction (right- vs left-lateral slip) throughout the Enceladus orbital cycle and explore a suite of model parameters that inhibit or promote shear failure at the tiger stripes. We find that low coefficients of friction (μf=0.1-0.2) and shallow overburden depths (z=2-4 km) permit shear failure along the tiger stripe faults, and that right- and/or left-lateral slip responses are possible. We integrate these conditions into a 3D time-dependent fault dislocation model to evaluate tectonic displacements and stress variations at depth during a tiger stripe orbital cycle. Depending on the sequence of stress accumulation and subsequent fault slip, which varies as a function of fault location and orientation, frictional coefficient, and fault depth, we estimate resolved shear stress accumulation of ∼70 kPa prior to fault failure, which produces modeled strike-slip displacements on the order of ∼0.5 m in the horizontal direction and ∼5 mm in the vertical direction per slip event. Our models also indicate that net displacements on the order of 0.1 m per orbital cycle, in both right- and left-lateral directions, are possible for particular fault geometries and frictional parameters. Tectonic activity inferred from these analyses correlates with observed plume activity and temperature anomalies at Enceladus's south polar region. Moreover, these analyses provide important details of stress accumulation and the faulting cycle for icy satellites subjected to diurnal tidal stress.     

High spatial resolution hyperspectral mapping of in-stream habitats, depths, and woody debris in mountain streams

This article evaluates the potential of 1-m resolution, 128-band hyperspectral imagery for mapping in-stream habitats, depths, and woody debris in third- to fifth-order streams in the northern Yellowstone region. Maximum likelihood supervised classification using principal component images provided overall classification accuracies for in-stream habitats (glides, riffles, pools, and eddy drop zones) ranging from 69% for third-order streams to 86% for fifth-order streams. This scale dependency of classification accuracy was probably driven by the greater proportion of transitional boundary areas in the smaller streams. Multiple regressions of measured depths (y) versus principal component scores (x₁, x₂,…, x_n) generated R² values ranging from 67% for high-gradient riffles to 99% for glides in a fifth-order reach. R² values were lower in third-order reaches, ranging from 28% for runs and glides to 94% for pools. The less accurate depth estimates obtained for smaller streams probably resulted from the relative increase in the number of mixed pixels, where a wide range of depths and surface turbulence occurred within a single pixel. Matched filter (MF) mapping of woody debris generated overall accuracies of 83% in the fifth-order Lamar River. Accuracy figures for the in-stream habitat and wood mapping may have been misleadingly low because the fine-resolution imagery captured fine-scale variations not mapped by field teams, which in turn generated false “misclassifications” when the image and field maps were compared.The use of high spatial resolution hyperspectral (HSRH) imagery for stream mapping is limited by the need for clear water to measure depth, by any tree cover obscuring the stream, and by the limited availability of airborne hyperspectral sensors. Nonetheless, the high accuracies achieved in northern Yellowstone streams indicate that HSRH imagery can be a powerful tool for watershed-wide mapping, monitoring, and modeling of streams.     

Seasonal connections between meteoric water and streamflow generation along a mountain headwater stream

In snowmelt-driven mountain watersheds, the hydrologic connectivity between meteoric waters and stream flow generation varies strongly with the season, reflecting variable connection to soil and groundwater storage within the watershed. This variable connectivity regulates how streamflow generation mechanisms transform the seasonal and elevational variation in oxygen and hydrogen isotopic composition (δ¹⁸O and δD) of meteoric precipitation. Thus, water isotopes in stream flow can signal immediate connectivity or more prolonged mixing, especially in high-relief mountainous catchments. We characterized δ¹⁸O and δD values in stream water along an elevational gradient in a mountain headwater catchment in southwestern Montana. Stream water isotopic compositions related most strongly to elevation between February and March, exhibiting higher δ¹⁸O and δD values with decreasing elevation. These elevational isotopic lapse rates likely reflect increased connection between stream flow and proximal snow-derived water sources heavily subject to elevational isotopic effects. These patterns disappeared during summer sampling, when consistently lower δ¹⁸O and δD values of stream water reflected contributions from snowmelt or colder rainfall, despite much higher δ¹⁸O and δD values expected in warmer seasonal rainfall. The consistently low isotopic values and absence of a trend with elevation during summer suggest lower connectivity between summer precipitation and stream flow generation as a consequence of drier soils and greater transpiration. As further evidence of intermittent seasonal connectivity between the stream and adjacent groundwaters, we observed a late-winter flush of nitrate into the stream at higher elevations, consistent with increased connection to accumulating mineralized nitrogen in riparian wetlands. This pattern was distinct from mid-summer patterns of nitrate loading at lower elevations that suggested heightened human recreational activity along the stream corridor. These observations provide insights linking stream flow generation and seasonal water storage in high elevation mountainous watersheds. Greater understanding of the connections between surface water, soil water and groundwater in these environments will help predict how the quality and quantity of mountain runoff will respond to changing climate and allow better informed water management decisions.     

Magnetostratigraphic evidence from the Cold Creek bar for onset of ice-age cataclysmic floods in eastern Washington during the Early Pleistocene

This study provides a detailed magnetostratigraphy of sediments composing the Cold Creek cataclysmic flood bar in the Pasco Basin, Washington. Our interpretation suggests onset of Missoula floods or similar events prior to 1.1 myr, later than previously suggested by Bjornstad et al. [Bjornstad, B.N., Fecht, K.R., Pluhar, C.J., 2001. Long history of pre-Wisconsin, Ice Age cataclysmic floods: evidence from southeastern Washington State. Journal of Geology 109 (6), 695-713]. Nonetheless these data suggest that Channeled Scabland features formed over a much longer timespan than commonly cited, that continental ice sheets of the early Pleistocene reached as far south as those of the late Pleistocene, and that similar physiography existed in eastern Washington and perhaps Montana to both generate and route Missoula-flood-like events. This study adds paleomagnetic polarity results from 213 new samples of silts and sands derived from nine new drill cores penetrating the Cold Creek cataclysmic flood bar to our previous database of 53 samples from four boreholes, resulting in a much more robust and detailed magnetostratigraphy. Rock magnetic studies on these sediments show pure magnetite to be the predominant remanence-carrying magnetic mineral, ruling out widespread remagnetization by secondary mineralization. The magnetostratigraphy at eastern Cold Creek bar is characterized by a normal polarity interval bracketed by reversed polarities. Equating the normal zone with the Jaramillo subchron (0.99-1.07 myr) affords the simplest correlation to the magnetic polarity timescale. Western Cold Creek bar was likely deposited during the Brunhes chron (0-0.78 myr) since it exhibits mainly normal polarities with only two thin reversed-polarity horizons that we interpret as magnetic excursions during the Brunhes.     

Reaction mechanism for the replacement of calcite by dolomite and siderite: implications for geochemistry, microstructure and porosity evolution during hydrothermal mineralisation

Carbonate reactions are common in mineral deposits due to CO₂-rich mineralising fluids. This study presents the first in-depth, integrated analysis of microstructure and microchemistry of fluid-mediated carbonate reaction textures at hydrothermal conditions. In doing so, we describe the mechanisms by which carbonate phases replace one another, and the implications for the evolution of geochemistry, rock microstructures and porosity. The sample from the 1.95 Moz Junction gold deposit, Western Australia, contains calcite derived from carbonation of a metamorphic amphibole—plagioclase assemblage that has further altered to siderite and dolomite. The calcite is porous and contains iron-rich calcite blebs interpreted to have resulted from fluid-mediated replacement of compositionally heterogeneous amphiboles. The siderite is polycrystalline but nucleates topotactically on the calcite. As a result, the boundaries between adjacent grains are low-angle boundaries (<10°), which are geometrically similar to those formed by crystal–plastic deformation and recovery. Growth zoning within individual siderite grains shows that the low-angle boundaries are growth features and not due to deformation. Low-angle boundaries develop due to the propagation of defects at grain faces and zone boundaries and by impingement of grains that nucleated with small misorientations relative to each other during grain growth. The cores of siderite grains are aligned with the twin planes in the parent calcite crystal showing that the reactant Fe entered the crystal along the twin boundaries. Dolomite grains, many of which appear to in-fill space generated by the siderite replacement, also show alignment of cores along the calcite twin planes, suggesting that they did not grow into space but replaced the calcite. Where dolomite is seen directly replacing calcite, it nucleates on the Fe-rich calcite due to the increased compatibility of the Fe-bearing calcite lattice relative to the pure calcite. Both reactions are interpreted as fluid-mediated replacement reactions which use the crystallography and elemental chemistry of the calcite. Experiments of fluid-mediated replacement reactions show that they proceed much faster than diffusion-based reactions. This is important when considering the rates of reactions relative to fluid flow in mineralising systems.     

Eclogite thermobarometry: The consistency between conventional thermobarometry and forward phase-equilibrium modelling

Eclogite thermobarometry is crucial for constraining the depths and temperatures to which oceanic and continental crust subduct. However, obtaining the pressure and temperature (P–T) conditions of eclogites is complex as they commonly display high-variance mineral assemblages, and the mineral compositions only vary slightly with P–T. In this contribution, we present a comparison between two independent and commonly used thermobarometric approaches for eclogites: conventional thermobarometry and forward phase-equilibrium modelling. We assess how consistent the thermobarometric calculations are using the garnet–clinopyroxene–phengite barometer and garnet–clinopyroxene thermometer with predictions from forward modelling (i.e. comparing the relative differences between approaches). Our results show that the overall mismatch in methods is typically ±0.2–0.3 GPa and ±29–42°C although differences as large as 80°C and 0.7 GPa are possible for a few narrow ranges of P–T conditions in the forward models. Such mismatch is interpreted as the relative differences among methods, and not as absolute uncertainties or accuracies for either method. For most of the investigated P–T conditions, the relatively minor differences between methods means that the choice in thermobarometric method itself is less important for geological interpretation than careful sample characterization and petrographic interpretation for deriving P–T from eclogites. Although thermobarometry is known to be sensitive to the assumed X_Fe³⁺ of a rock (or mineral), the relative differences between methods are not particularly sensitive to the choice of bulk-rock X_Fe³⁺, except at high temperatures (>650°C, amphibole absent) and for very large differences in assumed X_Fe³⁺ (0–0.5). We find that the most important difference between approaches is the activity–composition (a–x) relations, as opposed to the end-member thermodynamic data or other aspects of experimental calibration. When equivalent a–x relations are used in the conventional barometer, P calculations are nearly identical to phase-equilibrium models (ΔP < 0.1). To further assess the implications of these results for real rocks, we also evaluate common mathematical optimizations of reaction constants used for obtaining the maximum P–T with conventional thermobarometric approaches (e.g. using the highest aGrs² × aPrp in garnet and Si content in phengite, and the lowest aDi in clinopyroxene). These approaches should be used with caution, because they may not represent the compositions of equilibrium mineral assemblages at eclogite facies conditions and therefore systematically bias P–T calculations. Assuming method accuracy, geological meaningful P_max at a typical eclogite facies temperature of ~660°C will be obtained by using the greatest aDi, aCel, and aPrp and lowest aGrs and aMs; garnet and clinopyroxene with the lowest Fe²⁺/Mg ratios may yield geological meaningful T_max at a typical eclogite facies pressure of 2.5 GPa.     

Woodleigh impact structure,Australia: Shock petrography and geochemical studies

Abstract— The large, complex Woodleigh structure in the Carnarvon basin of Western Australia has recently been added to the terrestrial impact crater record. Many aspects of this structure are, however, still uncertain. This work provides a detailed petrographic assessment of a suite of representative drill core samples from the borehole Woodleigh 1 that penetrated uplifted basement rocks of the central part of this structure. Fundamental rock and mineral deformation data and high‐precision chemical data, including results of PGE and oxygen isotopic analysis, are presented. The sampled interval displays likely impact‐produced macrodeformation in the form of fracturing and breccia veining at the microscopic scale. Contrary to earlier reports that these breccias represent pseudotachylite (friction melt) or even shock/shear‐produced pseudotachylitic melt breccia cannot be confirmed due to pervasive post‐impact alteration. Abundant planar deformation features (PDFs) in quartz, in addition to diaplectic glass and partial isotropization, are the main shock deformation effects observed, confirming that Woodleigh is of impact origin. Over the investigated depth interval, the statistics of quartz grains with a variable number of sets of PDFs does not change significantly, and the patterns of crystallographic orientations of PDFs in randomly selected quartz grains does not indicate a change in absolute shock pressure with depth either. The value of oxygen isotopes for the recognition of meteoritic contamination, as proposed by earlier Woodleigh workers, is critically assessed. Neither INA nor PGE analyses of our samples support the presence of a meteoritic component within this basement section, as had been claimed in earlier work.     

Tracing the exhumation of the Eclogite Zone (Tauern Window,Eastern Alps) by <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar dating of white mica in eclogites

New radiometric ages from the Subpenninic nappes (Eclogite Zone and Rote Wand – Modereck Nappe, Tauern Window) show that phengites formed under eclogite-facies metamorphic conditions retain their initial isotopic signature, even when associated lithologies were overprinted by greenschist- to amphibolite-facies metamorphism. Different stages of the eclogite-facies evolution can be dated provided ⁴⁰Ar/³⁹Ar dating is combined with micro-structural analyses. An age of 39 Ma from the Rote Wand – Modereck Nappe is interpreted to be close to the burial age of this unit. Eclogite deformation within the Eclogite Zone started at the pressure peak along distinct shear zones, and prevailed along the exhumation path. An age of ca. 38 Ma is only observed for eclogites not affected by subsequent deformation and is interpreted as maximum age due to the possible influence of homogenously distributed excess argon. During exhumation deformation was localised along distinct mylonitic shear zones. This stage is mainly characterised by the formation of dynamically recrystallized omphacite2 and phengite. Deformation resulted in the resetting of the Ar isotopic system within the recrystallized white mica. Flat argon release spectra showing ages of 32 Ma within mylonites record the timing of cooling along the exhumation path, and the emplacement onto the Venediger Nappe. Ar-release patterns and ³⁶Ar/⁴⁰Ar vs.³⁹Ar/⁴⁰Ar isotope correlation analyses indicate no significant ⁴⁰Ar-loss after initial closure, and only a negligible incorporation of excess argon. From the pressure peak onwards, eclogitic conditions prevailed for almost 8–10 Ma.     

Luminescence chronology of river adjustment and incision of Quaternary sediments in the alluvial plain of the Sabarmati River, north Gujarat, India

River adjustment and incision in the Sabarmati basin, Gujarat, India have been examined at a site near Mahudi. Towards this, the morphostratigraphy and depositional chronometry of the middle alluvial plains were investigated. The upper fluvial sequence, along with the overlying aeolian sand and riverbed scroll plains, provide clues to the evolution of the present Sabarmati River. Sedimentological analyses of the upper fluvial sequence indicate its deposition by a meandering river system during what is believed to be a persistent wetter phase. Luminescence chronology bracketed this sequence to between 54 and 30 ka, which corresponds to Oxygen Isotope Stage-3, during which the SW monsoon was enhanced. The overlying aeolian sand has been dated to 12 ka, indicating that dune accretion occurred simultaneously with the strengthening of the SW monsoon during the Early Holocene. Adjustment of the Sabarmati along a N–S transect is placed around 12 ka and the incision is bracketed between 12 and 4.5 ka. River adjustment could have been tectonic; however, the incision was facilitated by the availability of continuous flow in the river caused by the SW monsoon. The basin experienced two tectonic events at about 3 and 0.3 ka, as demonstrated by the morphology of the scroll plains.