Direct visualization of solute locations in laboratory ice samples, links to videos

Many important chemical reactions occur in polar snow, where solutes may be present in several reservoirs, including at the air-ice interface and in liquid-like regions within the ice matrix. Some recent laboratory studies suggest chemical reaction rates may differ in these two reservoirs. While investigations have examined where solutes are found in natural snow and ice, similar research has not identified solute locations in laboratory samples, nor the possible factors controlling solute segregation. To address this, we examined solute locations in ice samples prepared from either aqueous cesium chloride (CsCl) or Rose Bengal solutions that were frozen using several different methods. Samples frozen in a laboratory freezer had the largest liquid-like inclusions and air bubbles, while samples frozen in a custom freeze chamber had somewhat smaller air bubbles and inclusions; in contrast, samples frozen in liquid nitrogen showed much smaller concentrated inclusions and air bubbles, only slightly larger than the resolution limit of our images (~2 µm). Freezing solutions in plastic versus glass vials had significant impacts on the sample structure, perhaps because the poor heat conductivity of plastic vials changes how heat is removed from the sample as it cools. Similarly, the choice of solute had a significant impact on sample structure, with Rose Bengal solutions yielding smaller inclusions and air bubbles compared to CsCl solutions frozen using the same method. Additional experiments using higher-resolution imaging of an ice sample show that CsCl moves in a thermal gradient, supporting the idea that the solutes in ice are present in liquid-like regions. Our work shows that the structure of laboratory ice samples, including the location of solutes, is sensitive to freezing method, sample container, and solute characteristics, requiring careful experimental design and interpretation of results.

Shear experiments from IODP Expedition 322 and 333 samples

The mechanical behavior of the plate boundary fault zone is of paramount importance in subduction zones, because it controls megathrust earthquake nucleation and propagation as well as the structural style of the forearc. In the Nankai area along the NanTroSEIZE (Kumano) drilling transect offshore SW Japan, a heterogeneous sedimentary sequence overlying the oceanic crust enters the subduction zone. In order to predict how variations in lithology, and thus mechanical properties, affect the formation and evolution of the plate boundary fault, we conducted laboratory tests measuring the shear strengths of sediments approaching the trench covering each major lithological sedimentary unit. We observe that shear strength increases nonlinearly with depth, such that the (apparent) coefficient of friction decreases. In combination with a critical taper analysis, the results imply that the plate boundary position is located on the main frontal thrust. Further landward, the plate boundary is expected to step down into progressively lower stratigraphic units, assisted by moderately elevated pore pressures. As seismogenic depths are approached, the décollement may further step down to lower volcaniclastic or pelagic strata but this requires specific overpressure conditions. High-taper angle and elevated strengths in the toe region may be local features restricted to the Kumano transect.

Table 1. Pressure, temperature and isotope ratios between H2O and CO2 for all experiments

Traditionally, the application of stable isotopes in Carbon Capture and Storage (CCS) projects has focused on d13C values of CO2 to trace the migration of injected CO2 in the subsurface. More recently the use of d18O values of both CO2 and reservoir fluids has been proposed as a method for quantifying in situ CO2 reservoir saturations due to O isotope exchange between CO2 and H2O and subsequent changes in d18OH2O values in the presence of high concentrations of CO2. To verify that O isotope exchange between CO2 and H2O reaches equilibrium within days, and that d18OH2O values indeed change predictably due to the presence of CO2, a laboratory study was conducted during which the isotope composition of H2O, CO2, and dissolved inorganic C (DIC) was determined at representative reservoir conditions (50°C and up to 19 MPa) and varying CO2 pressures. Conditions typical for the Pembina Cardium CO2 Monitoring Pilot in Alberta (Canada) were chosen for the experiments. Results obtained showed that d18O values of CO2 were on average 36.4±2.2 per mil (1 sigma, n=15) higher than those of water at all pressures up to and including reservoir pressure (19 MPa), in excellent agreement with the theoretically predicted isotope enrichment factor of 35.5 per mil for the experimental temperatures of 50°C. By using 18O enriched water for the experiments it was demonstrated that changes in the d18O values of water were predictably related to the fraction of O in the system sourced from CO2 in excellent agreement with theoretical predictions. Since the fraction of O sourced from CO2 is related to the total volumetric saturation of CO2 and water as a fraction of the total volume of the system, it is concluded that changes in d18O values of reservoir fluids can be used to calculate reservoir saturations of CO2 in CCS settings given that the d18O values of CO2 and water are sufficiently distinct.