Monte Carlo and molecular dynamics simulations of methane in potassium montmorillonite clay hydrates at elevated pressures and temperatures

The structure and dynamics of methane in hydrated potassium montmorillonite clay have been studied under conditions encountered in sedimentary basin and compared to those of hydrated sodium montmorillonite clay using computer simulation techniques. The simulated systems contain two molecular layers of water and followed gradients of 150 barkm-1 and 30 Kkm-1 up to a maximum burial depth of 6 km. Methane particle is coordinated to about 19 oxygen atoms, with 6 of these coming from the clay surface oxygen. Potassium ions tend to move away from the center towards the clay surface, in contrast to the behavior observed with the hydrated sodium form. The clay surface affinity for methane was found to be higher in the hydrated K-form. Methane diffusion in the two-layer hydrated K-montmorillonite increases from 0.39×10-9 m2s-1 at 280 K to 3.27×10-9 m2s-1 at 460 K compared to 0.36×10-9 m2s-1 at 280 K to 4.26×10-9 m2s-1 at 460 K in Na-montmorillonite hydrate. The distributions of the potassium ions were found to vary in the hydrates when compared to those of sodium form. Water molecules were also found to be very mobile in the potassium clay hydrates compared to sodium clay hydrates. © 2004 Elsevier Inc. All All rights reserved.

Molecular dynamics simulation of methane in sodium montmorillonite clay hydrates at elevated pressures and temperatures

Computer simulation has been used to study the structure and dynamics of methane in hydrated sodium montmorillonite clays under conditions encountered in sedimentary basins. Systems containing approximately one, two, three and four molecular layers of water have followed gradients of 150 bar km-1 and 30Kkm-1, to a maximum burial depth of 6 km (900 bar and 460 K). Methane is coordinated to approximately 19 oxygen atoms, of which typically 6 are provided by the clay surface. Only in the three- and four-layer hydrates is methane able to leave the clay surface. Diffusion depends strongly on the porosity (water content) and burial depth: self-diffusion coefficients are in the range 0.12 × 10-9m2s-1 for water and 0.04 × 10−9m2s−1 < D < 8.64 × 10−9m2s−1 for methane. Bearing in mind that porosity decreases with burial depth, it is estimated that maximum diffusion occurs at around 3 km. This is in good agreement with the known location of methane reservoirs in sedimentary basins.

The study of silicates and refractory materials at high pressures and temperatures

To better understand high pressure behavior of solids, both silicates and oxides have been investigated to clarify the high pressure melting, phase transformations and thermal parameters as well as their size dependences, both theoretically and experimentally. ^ To judge the precision of data determined experimentally, the reliabilities of different high pressure techniques have been discussed. A thermodynamic model has been developed and demonstrated to be able to closely reproduce the melting of solids by comparison between results calculated and data obtained experimentally, including metals (Al, Ni and Pt), Silicates (Mg3Al 2Si3O12 and CaMgSi2O6), Halides (NaCl, CsCl and LiF) and Oxides (MgO, FeO and Al2O3). The melting data obtained have been discussed to address the dynamics of the Earth's interior. ^ Results obtained with Raman spectroscopy and x-ray diffraction show that solids including silicates (andradite and pyrope) and oxides (CeO2 and TiO2) undergo a series of pressure-induced phase transformations. The effects of particle size under high pressures have been investigated. The results obtained indicate that the reduction of particle size leads to the enhancement of the bulk modulus and a significant decrease of transition pressure in TiO2 (rutile) and CeO2. The pressure-induced amorphization in anatase also results from the size effects. ^ Combining the data obtained with global seismic tomography, the physics and chemistry of the Earth's mantle and the dynamics of the core-mantle interaction have been discussed. The high pressure phases of Al3+- and Fe3+-bearing minerals play important roles in the dynamics of the lower mantle. ^

Exploring thermal and mechanical properties of selected transition elements under extreme conditions: Experiments at high pressures and high temperatures

Transition metals (Ti, Zr, Hf, Mo, W, V, Nb, Ta, Pd, Pt, Cu, Ag, and Au) are essential building units of many materials and have important industrial applications. Therefore, it is important to understand their thermal and physical behavior when they are subjected to extreme conditions of pressure and temperature. This dissertation presents: • An improved experimental technique to use lasers for the measurement of thermal conductivity of materials under conditions of very high pressure (P, up to 50 GPa) and temperature (T up to 2500 K). • An experimental study of the phase relationship and physical properties of selected transition metals, which revealed new and unexpected physical effects of thermal conductivity in Zr, and Hf under high P-T. • New phase diagrams created for Hf, Ti and Zr from experimental data. • P-T dependence of the lattice parameters in α-hafnium. Contrary to prior reports, the α-ω phase transition in hafnium has a negative dT/dP slope. • New data on thermodynamic and physical properties of several transition metals and their respective high P-T phase diagrams. • First complete thermodynamic database for solid phases of 13 common transition metals was created. This database has: All the thermochemical data on these elements in their standard state (mostly available and compiled); All the equations of state (EoS) formulated from pressure-volume-temperature data (measured as a part of this study and from literature); Complete thermodynamic data for selected elements from standard to extreme conditions. The thermodynamic database provided by this study can be used with available thermodynamic software to calculate all thermophysical properties and phase diagrams at high P-T conditions. For readers who do not have access to this software, tabulated values of all thermodynamic and volume data for the 13 metals at high P-T are included in the APPENDIX. In the APPENDIX, a description of several other high-pressure studies of selected oxide systems is also included. Thermophysical properties (Cp, H, S, G) of the high P-T ω-phase of Ti, Zr and Hf were determined during the optimization of the EoS parameters and are presented in this study for the first time. These results should have important implications in understanding hexagonal-close-packed to simple-hexagonal phase transitions in transition metals and other materials.

Impact of high pCO2 and warmer temperatures on the process of silica biomineralization in the sponge Mycale grandis

Siliceous sponges have survived pre-historical mass extinction events caused by ocean acidification and recent studies suggest that siliceous sponges will continue to resist predicted increases in ocean acidity. In this study, we monitored silica biomineralization in the Hawaiian sponge Mycale grandis under predicted pCO2 and sea surface temperature scenarios for 2100. Our goal was to determine if spicule biomineralization was enhanced or repressed by ocean acidification and thermal stress by monitoring silica uptake rates during short-term (48 h) experiments and comparing biomineralized tissue ratios before and after a long-term (26 d) experiment. In the short-term experiment, we found that silica uptake rates were not impacted by high pCO2 (1050 µatm), warmer temperatures (27°C), or combined high pCO2 with warmer temperature (1119 µatm; 27°C) treatments. The long-term exposure experiments revealed no effect on survival or growth rates of M. grandis to high pCO2 (1198 µatm), warmer temperatures (25.6°C), or combined high pCO2 with warmer temperature (1225 µatm, 25.7°C) treatments, indicating that M. grandis will continue to prosper under predicted increases in pCO2 and sea surface temperature. However, ash-free dry weight to dry weight ratios, subtylostyle lengths, and silicified weight to dry weight ratios decreased under conditions of high pCO2 and combined pCO2 warmer temperature treatments. Our results show that rising ocean acidity and temperature have marginal negative effects on spicule biomineralization and will not affect sponge survival rates of M. grandis.

Fertilisation and larval development in an Antarctic bivalve, Laternula elliptica, under reduced pH and elevated temperatures

Elevated temperatures associated with ocean warming and acidification can influence development and, ultimately, success of larval molluscs. The effect of projected oceanic changes on fertilisation and larval development in an Antarctic bivalve, Laternula elliptica, was investigated through successive larval stages at ambient temperature and pH conditions (-1.6°C and pH 7.98) and conditions representative of projections through to 2100 (-0.5°C to +0.4°C and pH 7.80 to pH 7.65). Where significant effects were detected, increased temperature had a consistently positive influence on larval development, regardless of pH level, while effects of reduced pH varied with larval stage and incubation temperature. Fertilisation was high and largely independent of stressors, with no loss of gamete viability. Mortality was unaffected at all development stages under experimental conditions. Elevated temperatures reduced occurrences of abnormalities in D-larvae and accelerated larval development through late veliger and D-larval stages, with D-larvae occurring 5 d sooner at 0.4°C compared to ambient temperature. Reduced pH did not affect occurrences of abnormalities in larvae, but it slowed the development of calcifying stages. More work is required to investigate the effects of developmental delays of the magnitude seen here in order to better determine the ecological relevance of these changes on longer term larval and juvenile success.

Estimation of sea surface temperatures and salinity of ODP Site 177-1089 (Appendix A)

ODP Site 1089 is optimally located in order to monitor the occurrence of maxima in Agulhas heat and salt spillage from the Indian to the Atlantic Ocean. Radiolarian-based paleotemperature transfer functions allowed to reconstruct the climatic history for the last 450 kyr at this location. A warm sea surface temperature anomaly during Marine Isotope Stage (MIS) 10 was recognized and traced to other oceanic records along the surface branch of the global thermohaline (THC) circulation system, and is particularly marked at locations where a strong interaction between oceanic and atmospheric overturning cells and fronts occurs. This anomaly is absent in the Vostok ice core deuterium, and in oceanic records from the Antarctic Zone. However, it is present in the deuterium excess record from the Vostok ice core, interpreted as reflecting the temperature at the moisture source site for the snow precipitated at Vostok Station. As atmospheric models predict a subtropical Indian source for such moisture, this provides the necessary teleconnection between East Antarctica and ODP Site 1089, as the subtropical Indian is also the source area of the Agulhas Current, the main climate agent at our study location. The presence of the MIS 10 anomaly in the delta13C foraminiferal records from the same core supports its connection to oceanic mechanisms, linking stronger Agulhas spillover intensity to increased productivity in the study area. We suggest, in analogy to modern oceanographic observations, this to be a consequence of a shallow nutricline, induced by eddy mixing and baroclinic tide generation, which are in turn connected to the flow geometry, and intensity, of the Agulhas Current as it flows past the Agulhas Bank. We interpret the intensified inflow of Agulhas Current to the South Atlantic as responding to the switch between lower and higher amplitude in the insolation forcing in the Agulhas Current source area. This would result in higher SSTs in the Cape Basin during the glacial MIS 10, due to the release into the South Atlantic of the heat previously accumulating in the subtropical and equatorial Indian and Pacific Ocean. If our explanation for the MIS 10 anomaly in terms of an insolation variability switch is correct, we might expect that a future Agulhas SSST anomaly event will further delay the onset of next glacial age. In fact, the insolation forcing conditions for the Holocene (the current interglacial) are very similar to those present during MIS 11 (the interglacial preceding MIS 10), as both periods are characterized by a low insolation variability for the Agulhas Current source area. Natural climatic variability will force the Earth system in the same direction as the anthropogenic global warming trend, and will thus lead to even warmer than expected global temperatures in the near future.

Planktonic foraminiferal d18O and d13C values, derived sea-surface temperatures from d18O and selected planktonic foraminiferal species for core P69 of Hawke Bay off eastern North Island

Recent evidence suggests that the Subtropical Convergence (STC) zone east of New Zealand shifted little from its modern position along Chatham Rise during the last glaciation, and that offshore surface waters north of the STC zone cooled only slightly. However, at nearshore core site P69 (2195 m depth), 115 km off the east coast of North Island and ca 300 km north of the modern STC zone, planktonic foraminiferal species, transfer function data and stable oxygen and carbon isotope records suggest that surface waters were colder by up to 6°C during the late last glacial period compared to the Holocene, and included a strong upwelling signature. Presently site P69 is bathed by south-flowing subtropical waters in the East Cape Current. The nearshore western end of Chatham Rise supports a major bathymetric depression, the Mernoo Saddle, through which some exchange between northern subtropical and southern subantarctic water presently occurs. It is proposed that as a result of much intensified current flows south of the Rise during the last glaciation, a consequence of more compressed subantarctic water masses, lowered sea level, and an expanded and stronger Westerly Wind system, there was accelerated leakage northwards of both Australasian Subantarctic Water and upwelled Antarctic Intermediate Water over Mernoo Saddle in a modified and intensified Southland Current. The expanded cold water masses displaced the south-flowing warm East Cape Current off southeastern North Island, and offshore divergence was accompanied by wind-assisted upwelling of nutrient-rich waters in the vicinity of P69. A comparable kind of inshore cold water jetting possibly characterised most glacial periods since the latest Miocene, and may account for the occasional occurrence of subantarctic marine fossils in onland late Cenozoic deposits north of the STC zone, rather than invoking wholesale major oscillations of the oceanic STC itself.

Estimates of annual mean surface temperatures in the late Quaternary of the tropical Atlantic

At least two modes of glacial-interglacial climate change have existed within the tropical Atlantic Ocean during the last 20,000 years. The first mode (defined by cold glacial and warm interglacial conditions) occurred symmetrically north and south of the equator and dominated the eastern boundary currents and tropical upwelling areas. This pattern suggests that mode 1 is driven by a glacial modification of surface winds in both hemispheres. The second mode of oceanic climate change, defined by temperature extremes centered on the deglaciation, was hemispherically asymmetrical, with the northern tropical Atlantic relatively cold and the southern tropical Atlantic relatively warm during deglaciation. A likely cause for this pattern of variation is a reduction of the presently northward cross-equatorial heat flux during deglaciation. No single mechanism accounts for all the data. Potential contributors to oceanic climate changes are linkage to high-latitude climates, modification of monsoonal winds by ice sheet and/or insolation changes, atmospheric CO2 and greenhouse effects, indirect effects of glacial meltwater, and variations in thermohaline overturn of the oceans.