Basal slip along glaciers and ice streams can be significantly modified by external time-dependent forcing, although it is not clear why some systems are more sensitive to tidal stresses. We have conducted a series of laboratory experiments to explore the effect of time varying load point velocity on ice-on-rock friction. Varying the load point velocity induces shear stress forcing, making this an analogous simulation of aspects of ice stream tidal modulation. Ambient pressure, double-direct shear experiments were conducted in a cryogenic servo-controlled biaxial deformation apparatus at temperatures between −2°C and −16°C. In addition to a background, median velocity (1 and 10 μm/s), a sinusoidal velocity was applied to the central sliding sample over a range of periods and amplitudes. Normal stress was held constant over each run (0.1, 0.5 or 1 MPa) and the shear stress was measured. Over the range of parameters studied, the full spectrum of slip behavior from creeping to slow-slip to stick-slip was observed, similar to the diversity of sliding styles observed in Antarctic and Greenland ice streams. Under conditions in which the amplitude of oscillation is equal to the median velocity, significant healing occurs as velocity approaches zero, causing a high-amplitude change in friction. The amplitude of the event increases with increasing period (i.e. hold time). At high normal stress, velocity oscillations force an otherwise stable system to behave unstably, with consistently-timed events during every cycle. Rate-state friction parameters determined from velocity steps show that the ice-rock interface is velocity strengthening. A companion paper describes a method of analyzing the oscillatory data directly. Forward modeling of a sinusoidally-driven slider block, using rate-and-state dependent friction formulation and experimentally derived parameters, successfully predicts the experimental output in all but a few cases.
more »
« less
The Bridgmanite–Akimotoite–Majorite Triple Point Determined in Large Volume Press and Laser-Heated Diamond Anvil Cell
The bridgmanite–akimotoite–majorite (Bm–Ak–Mj or BAM) triple point in MgSiO3 has been measured in large-volume press (LVP; COMPRES 8/3 assembly) and laser-heated diamond anvil cell (LHDAC). For the LVP data, we calculated pressures from the calibration provided for the assembly. For the LHDAC data, we conducted in situ determination of pressure at high temperature using the Pt scale at synchrotron. The measured temperatures of the triple point are in good agreement between LVP and LHDAC at 1990–2000 K. However, the pressure for the triple point determined from the LVP is 3.9 ± 0.6 GPa lower than that from the LHDAC dataset. The BAM triple point determined through these experiments will provide an important reference point in the pressure–temperature space for future high-pressure experiments and will allow mineral physicists to compare the pressure–temperature conditions measured in these two different experimental methods.
more »
« less
- Award ID(s):
- 1725094
- NSF-PAR ID:
- 10201045
- Date Published:
- Journal Name:
- Minerals
- Volume:
- 10
- Issue:
- 1
- ISSN:
- 2075-163X
- Page Range / eLocation ID:
- 67
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Experimental determination of quartz solubility in H2O-CaCl2 solutions at 600–900 °C and 0.6–1.4 GPamore » « less
Abstract Fluid-mediated calcium metasomatism is often associated with strong silica mobility and the presence of chlorides in solution. To help quantify mass transfer at lower crustal and upper mantle conditions, we measured quartz solubility in H2O-CaCl2 solutions at 0.6–1.4 GPa, 600–900 °C, and salt concentrations to 50 mol%. Solubility was determined by weight loss of single-crystals using hydrothermal piston-cylinder methods. All experiments were conducted at salinity lower than salt saturation. Quartz solubility declines exponentially with added CaCl2 at all conditions investigated, with no evidence for complexing between silica and Ca. The decline in solubility is similar to that in H2O-CO2 but substantially greater than that in H2O-NaCl at the same pressure and temperature. At each temperature, quartz solubility at low salinity (XCaCl2 < 0.1) depends strongly on pressure, whereas at higher XCaCl2 it is nearly pressure independent. This behavior is consistent with a transition from an aqueous solvent to a molten salt near XCaCl2 ~0.1. The solubility data were used to develop a thermodynamic model of H2O-CaCl2 fluids. Assuming ideal molten-salt behavior and utilizing previous models for polymerization of hydrous silica, we derived values for the activity of H2O (aH2O), and for the CaCl2 dissociation factor (α), which may vary from 0 (fully associated) to 2 (fully dissociated). The model accurately reproduces our data along with those of previous work and implies that, at conditions of this study, CaCl2 is largely associated (<0.2) at H2O density <0.85 g/cm3. Dissociation rises isothermally with increasing density, reaching ~1.4 at 600 °C, 1.4 GPa. The variation in silica molality with aH2O in H2O-CaCl2 is nearly identical to that in H2O-CO2 solutions at 800 °C and 1.0 GPa, consistent with the absence of Ca-silicate complexing. The results suggest that the ionization state of the salt solution is an important determinant of aH2O, and that H2O-CaCl2 fluids exhibit nearly ideal molecular mixing over a wider range of conditions than implied by previous modeling. The new data help interpret natural examples of large-scale Ca-metasomatism in a wide range of lower crustal and upper mantle settings.
-
Titanite has the ability to incorporate significant amounts of common Pb, which leads to uncertainty when applying the U-Pb decay series for geochronology. The isobaric interference of 204Hg on 204Pb poses an additional complexity in applying common Pb corrections. Here we investigate the removal of 204Hg interferences during titanite U-Pb dating using reaction cell gas chemistry via triple quadrupole mass spectrometry. U-Pb dates were determined for the natural titanite reference materials MKED-1 and BLR1 using an ESI NWR193UC excimer laser coupled to an Agilent 8900 ‘triple quad’ mass spectrometer. The 8900 is equipped with an octopole collision/reaction cell, which enables online interference removal. Two experiments were run, one in which we collected data in NoGas mode, and one in which NH3 was used as a reaction cell gas in MS/MS mode, in order to assess the feasibility of determining U/Pb ratios with mass shifted isotopes. In all experiments, a signal smoothing device was placed inline just before the ICP-MS interface, downstream from the addition of the Ar nebulizer gas to the He carrier gas stream. For the NoGas experiment, titanite was ablated using a 25 µm spot, with a beam energy density of 3 J/cm2, and a pulse rate of 4 Hz. In NoGas mode, signal intensities for the isotopes 201Hg, 202Hg, 204Pb, 206Pb, 207Pb, 232Th, 235U, and 238U were counted. In MS/MS mode, titanite was ablated using a 40 µm spot, with a beam energy density of 5 J/cm2, and a pulse rate of 4 Hz. A larger spot size in this experiment was used to counteract the decrease in signal intensity due to use of the reaction cell. In MS/MS mode, NH3 was flowed through the reaction cell in order to enable a charge transfer reaction between NH3 and Hg+, effectively neutralizing Hg. The isotopes 201Hg, 202Hg, 204Pb, 206Pb, and 207Pb were measured on-mass, as the isotopes of Pb are not affected by the NH3 gas. Uranium and Th both exhibit partial reaction with NH3 gas; therefore, the isotopes 232Th, 235U, and 238U were measured mass-shifted up 15 mass units, at masses 247, 250, and 253 respectively. Ratios of 207Pb/235U, 206Pb/238U, and 207Pb/206Pb were determined using the UPbGeochron4 DRS in Iolite (v.3.71) with MKED-1 as the primary reference material. Dates were calculated using IsoplotR by applying the Stacey-Kramers correction for common Pb. All isotopes of Hg were effectively neutralized by the NH3 charge transfer reaction in MS/MS mode; zero counts were detected for Hg isotopes. Dates for the BLR-1 titanite were 1050.55 ± 2.72 (2σ, n=12) Ma in NoGas mode, and 1048 ± 1.88 (2σ, n=15) Ma in MS/MS mode. These dates are in excellent agreement with the TIMS 206Pb/238U date for the BLR-1 titanite of 1047.1 ± 0.4 Ma. This method has the potential to enable measurement of 204Pb without needing to correct for Hg interferences.more » « less
-
more » « less
Abstract Recent advances in remote sensing and the upcoming launch of the joint NASA/CNES/CSA/UKSA Surface Water and Ocean Topography (SWOT) satellite point toward improved river discharge estimates in ungauged basins. Existing discharge methods rely on “prior river knowledge” to infer parameters not directly measured from space. Here, we show that discharge estimation is improved by classifying and parameterizing rivers based on their unique geomorphology and hydraulics. Using over 370,000 in situ hydraulic observations as training data, we test unsupervised learning and an “expert” method to assign these hydraulics and geomorphology to rivers via remote sensing. This intervention, along with updates to model physics, constitutes a new method we term “geoBAM,” an update of the Bayesian At‐many‐stations hydraulic geometry‐Manning's (BAM) algorithm. We tested geoBAM on Landsat imagery over more than 7,500 rivers (108 are gauged) in Canada's Mackenzie River basin and on simulated hydraulic data for 19 rivers that mimic SWOT observations without measurement error. geoBAM yielded considerable improvement over BAM, improving the median Nash‐Sutcliffe efficiency (NSE) for the Mackenzie River from −0.05 to 0.26 and from 0.16 to 0.46 for the SWOT rivers. Further, NSE improved by at least 0.10 in 78/108 gauged Mackenzie rivers and 8/19 SWOT rivers. We attribute geoBAM improvement to parameterizing rivers by type rather than globally, but prediction accuracy worsens if parameters are misassigned. This method is easily mapped to rivers at the global scale and paves the way for improving future discharge estimates, especially when coupled with hydrologic models.
-
Abstract We perform ab initio simulations of beryllium (Be) and magnesium oxide (MgO) at megabar pressures and compare their structural and thermodynamic properties. We make a detailed comparison of our two recently derived phase diagrams of Be (Wu et al 2021 Phys. Rev. B 104 014103) and MgO (Soubiran and Militzer 2020 Phys. Rev. Lett. 125 175701) using the thermodynamic integration technique, as they exhibit striking similarities regarding their shape. We explore whether the Lindemann criterion can explain the melting temperatures of these materials through the calculation of the Debye temperature at high pressure. From our free energy calculations, we find that the melting line of both materials is well represented by the Simon–Glazel fit T m ( P ) = T 0 (1 + P / a ) 1/ c , where T 0 = 1564 K, a = 15.8037 GPa and c = 2.4154 for Be, while T 0 = 3010 K, a = 10.5797 GPa and c = 2.8683 for the MgO in the B1. For the B2 phase, we use the values a = 26.1163 GPa and c = 2.2426. Both materials exhibit negative Clapeyron slopes on the boundaries between the two solid phases that are strongly affected by anharmonic effects, which also influence the location of the solid–solid–liquid triple point. We find that the quasi-harmonic approximation underestimates the stability range of the low-pressure phases, namely hcp for Be and B1 for MgO. We also compute the phonon dispersion relations at low and high pressure for each of the phases of these materials, and also explore how the phonon density of states is modified by temperature. Finally, we derive secondary shock Hugoniot curves in addition to the principal Hugoniot curve for both materials, and study their offsets in pressure between solid and liquid branches.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/min10010067
