Heat-producing crust regulation of subsurface temperatures : a stochastic model re-evaluation of the geothermal potential in western Queensland, Australia

A large subsurface, elevated temperature anomaly is well documented in Central Australia. High Heat Producing Granites (HHPGs) intersected by drilling at Innamincka are often assumed to be the dominant cause of the elevated subsurface temperatures, although their presence in other parts of the temperature anomaly has not been confirmed. Geological controls on the temperature anomaly remain poorly understood. Additionally, methods previously used to predict temperature at 5 km depth in this area are simplistic and possibly do not give an accurate representation of the true distribution and magnitude of the temperature anomaly. Here we re-evaluate the geological controls on geothermal potential in the Queensland part of the temperature anomaly using a stochastic thermal model. The results illustrate that the temperature distribution is most sensitive to the thermal conductivity structure of the top 5 km. Furthermore, the results indicate the presence of silicic crust enriched in heat producing elements between and 40 km.

Cyclic Salinity Model for deposition of large-scale coal deposits

Existing field data for Rangal coals (Late Permian) of the Bowen Basin, Queensland, Australia, are inconsistent with the depositional model generally accepted in the current geological literature to explain coal deposition. Given the apparent unsuitability of the current depositional model to the Bowen Basin coal data, a new depositional model, here named the Cyclic Salinity Model, is proposed and tested in this study.

A general geomorphological recession flow model for river basins

Recession flows in a basin are controlled by the temporal evolution of its active drainage network (ADN). The geomorphological recession flow model (GRFM) assumes that both the rate of flow generation per unit ADN length (q) and the speed at which ADN heads move downstream (c) remain constant during a recession event. Thereby, it connects the power law exponent of -dQ/dt versus Q (discharge at the outlet at time t) curve, , with the structure of the drainage network, a fixed entity. In this study, we first reformulate the GRFM for Horton-Strahler networks and show that the geomorphic ((g)) is equal to D/(D-1), where D is the fractal dimension of the drainage network. We then propose a more general recession flow model by expressing both q and c as functions of Horton-Strahler stream order. We show that it is possible to have = (g) for a recession event even when q and c do not remain constant. The modified GRFM suggests that is controlled by the spatial distribution of subsurface storage within the basin. By analyzing streamflow data from 39 U.S. Geological Survey basins, we show that is having a power law relationship with recession curve peak, which indicates that the spatial distribution of subsurface storage varies across recession events. Key Points The GRFM is reformulated for Horton-Strahler networks. The GRFM is modified by allowing its parameters to vary along streams. Sub-surface storage distribution controls recession flow characteristics.

Source Model for Blasting Vibration

By analyzing and comparing the experimental data, the point source moment theory and the cavity theory, it is concluded that the vibrating signals away from the blasting explosive come mainly from the natural vibrations of the geological structures near the broken blasting area. The source impulses are not spread mainly by the inelastic properties (such as through media damping, as believed to be the case by many researchers) of the medium in the propagation pass, but by this structure. Then an equivalent source model for the blasting vibrations of a fragmenting blasting is proposed, which shows the important role of the impulse of the source's time function under certain conditions. For the purpose of numerical simulation, the model is realized in FEM, The finite element results are in good agreement with the experimental data.

New insights into the formation and modification of carbonate-bearing minerals and methane gas in geological systems using multiply substituted isotopologues

This thesis describes the use of multiply-substituted stable isotopologues of carbonate minerals and methane gas to better understand how these environmentally significant minerals and gases form and are modified throughout their geological histories. Stable isotopes have a long tradition in earth science as a tool for providing quantitative constraints on how molecules, in or on the earth, formed in both the present and past. Nearly all studies, until recently, have only measured the bulk concentrations of stable isotopes in a phase or species. However, the abundance of various isotopologues within a phase, for example the concentration of isotopologues with multiple rare isotopes (multiply substituted or 'clumped' isotopologues) also carries potentially useful information. Specifically, the abundances of clumped isotopologues in an equilibrated system are a function of temperature and thus knowledge of their abundances can be used to calculate a sample’s formation temperature. In this thesis, measurements of clumped isotopologues are made on both carbonate-bearing minerals and methane gas in order to better constrain the environmental and geological histories of various samples.

Clumped-isotope-based measurements of ancient carbonate-bearing minerals, including apatites, have opened up paleotemperature reconstructions to a variety of systems and time periods. However, a critical issue when using clumped-isotope based measurements to reconstruct ancient mineral formation temperatures is whether the samples being measured have faithfully recorded their original internal isotopic distributions. These original distributions can be altered, for example, by diffusion of atoms in the mineral lattice or through diagenetic reactions. Understanding these processes quantitatively is critical for the use of clumped isotopes to reconstruct past temperatures, quantify diagenesis, and calculate time-temperature burial histories of carbonate minerals. In order to help orient this part of the thesis, Chapter 2 provides a broad overview and history of clumped-isotope based measurements in carbonate minerals.

In Chapter 3, the effects of elevated temperatures on a sample’s clumped-isotope composition are probed in both natural and experimental apatites (which contain structural carbonate groups) and calcites. A quantitative model is created that is calibrated by the experiments and consistent with the natural samples. The model allows for calculations of the change in a sample’s clumped isotope abundances as a function of any time-temperature history.

In Chapter 4, the effects of diagenesis on the stable isotopic compositions of apatites are explored on samples from a variety of sedimentary phosphorite deposits. Clumped isotope temperatures and bulk isotopic measurements from carbonate and phosphate groups are compared for all samples. These results demonstrate that samples have experienced isotopic exchange of oxygen atoms in both the carbonate and phosphate groups. A kinetic model is developed that allows for the calculation of the amount of diagenesis each sample has experienced and yields insight into the physical and chemical processes of diagenesis.

The thesis then switches gear and turns its attention to clumped isotope measurements of methane. Methane is critical greenhouse gas, energy resource, and microbial metabolic product and substrate. Despite its importance both environmentally and economically, much about methane’s formational mechanisms and the relative sources of methane to various environments remains poorly constrained. In order to add new constraints to our understanding of the formation of methane in nature, I describe the development and application of methane clumped isotope measurements to environmental deposits of methane. To help orient the reader, a brief overview of the formation of methane in both high and low temperature settings is given in Chapter 5.

In Chapter 6, a method for the measurement of methane clumped isotopologues via mass spectrometry is described. This chapter demonstrates that the measurement is precise and accurate. Additionally, the measurement is calibrated experimentally such that measurements of methane clumped isotope abundances can be converted into equivalent formational temperatures. This study represents the first time that methane clumped isotope abundances have been measured at useful precisions.

In Chapter 7, the methane clumped isotope method is applied to natural samples from a variety of settings. These settings include thermogenic gases formed and reservoired in shales, migrated thermogenic gases, biogenic gases, mixed biogenic and thermogenic gas deposits, and experimentally generated gases. In all cases, calculated clumped isotope temperatures make geological sense as formation temperatures or mixtures of high and low temperature gases. Based on these observations, we propose that the clumped isotope temperature of an unmixed gas represents its formation temperature — this was neither an obvious nor expected result and has important implications for how methane forms in nature. Additionally, these results demonstrate that methane-clumped isotope compositions provided valuable additional constraints to studying natural methane deposits.

Recycle of transuranium elements in light water reactors for reduction of geological storage requirements

This paper presents results of a feasibility study aimed at developing a zero-transuranic-discharge fuel cycle based on the U-Th-TRU ternary cycle. The design objective is to find a fuel composition (mixture of thorium, enriched uranium, and recycled transuranic components) and fuel management strategy resulting in an equilibrium charge-discharge mass flow. In such a fuel cycle scheme, the quantity and isotopic vector of the transuranium (TRU) component is identical at the charge and discharge time points, thus allowing the whole amount of the TRU at the end of the fuel irradiation period to be separated and reloaded into the following cycle. The TRU reprocessing activity losses are the only waste stream that will require permanent geological storage, virtually eliminating the long-term radiological waste of the commercial nuclear fuel cycle. A detailed three-dimensional full pressurized water reactor (PWR) core model was used to analyze the proposed fuel composition and management strategy. The results demonstrate the neutronic feasibility of the fuel cycle with zero-TRU discharge. The amount of TRU and enriched uranium loaded reach equilibrium after about four TRU recycles. The reactivity coefficients were found to be within a range typical for a reference PWR core. The soluble boron worth is reduced by a factor of ∼2 from a typical PWR value. Nevertheless, the results indicate the feasibility of an 18-month fuel cycle design with an acceptable beginning-of-cycle soluble boron concentration even without application of burnable poisons.

Stratigraphic Model Predictions of Geoacoustic Properties

Geoacoustic properties of the seabed have a controlling role in the propagation and reverberation of sound in shallow-water environments. Several techniques are available to quantify the important properties but are usually unable to adequately sample the region of interest. In this paper, we explore the potential for obtaining geotechnical properties from a process-based stratigraphic model. Grain-size predictions from the stratigraphic model are combined with two acoustic models to estimate sound speed with distance across the New Jersey continental shelf and with depth below the seabed. Model predictions are compared to two independent sets of data: 1) Surficial sound speeds obtained through direct measurement using in situ compressional wave probes, and 2) sound speed as a function of depth obtained through inversion of seabed reflection measurements. In water depths less than 100 m, the model predictions produce a trend of decreasing grain-size and sound speed with increasing water depth as similarly observed in the measured surficial data. In water depths between 100 and 130 m, the model predictions exhibit an increase in sound speed that was not observed in the measured surficial data. A closer comparison indicates that the grain-sizes predicted for the surficial sediments are generally too small producing sound speeds that are too slow. The predicted sound speeds also tend to be too slow for sediments 0.5-20 m below the seabed in water depths greater than 100 m. However, in water depths less than 100 m, the sound speeds between 0.5-20-m subbottom depth are generally too fast. There are several reasons for the discrepancies including the stratigraphic model was limited to two dimensions, the model was unable to simulate biologic processes responsible for the high sound-speed shell material common in the model area, and incomplete geological records necessary to accurately predict grain-size

Genetic model of the hydrate system in the fine grain sediments in the northern continental slope of South China Sea

Gas hydrate samples were obtained firstly in China by drilling on the northern margin of South China Sea (SCS). To understand the formation mechanism of this unique accumulation system, this paper discusses the factors controlling the formation of the system by accurate geophysical interpretation and geological analysis, based on the high precision 2-D and 3-D multichannel seismic data in the drilling area. There are three key factors controlling the accumulation of the gas hydrate system in fine grain sediment: (1) large volume of fluid bearing methane gas Joins the formation of gas hydrate. Active fluid flow in the northern South China Sea makes both thermal gas and/or biogenic gas migrate into shallow strata and form hydrate in the gas hydrate stability zone (GHSZ). The fluid flow includes mud diapir and gas chimney structure. They are commonly characterized by positive topographic relief, acoustic turbidity and push-down, and low reflection intensity on seismic profiles. The gas chimneys can reach to GHSZ, which favors the development of BSRs. It means that the active fluid flow has a close relationship with the formation and accumulation of gas hydrate. (2) The episodic process of fracture plays an important role in the generation of gas hydrate. It may provide the passage along which thermogenic or biogenic gas migrated into gas hydrate stability zone (GHSZ) upward. And it increases the pore space for the growth of hydrate crystal. (3) Submarine landslide induced the anomalous overpressure activity and development of fracture in the GHSZ. The formation model of high concentration gas hydrate in the drilling sea area was proposed on the basis of above analysis.

Simultaneous determination of four phases of copper in geological samples

According to the specific property of respective phases in copper ore, optimal extracting condition for, tour phases of copper with five infusion was examined. After a series of experiments, for example, different weight Of sample, artificial mixed single phases sample and repeated experiment, the extracting ratio constants foil all phases in different infusions were obtained. The mathematical model which is on the basis of the stable extracting constant is established. Based on the purpose-control genetic algorithms, we only need to determine the total copper in the sample in different infusions under the selected condition, then the value of respective phase could be calculated. The-relative standard deviation are (%) free oxidized copper phase:1.5, conjunction oxidized copper phase: 17.6,, secondary copper sulfide phase: 1.9, primary copper sulfide phase : 2.7,total copper: 0.8.

Physical and economic potential of geological CO2 storage in saline aquifers.

Carbon sequestration in sandstone saline reservoirs holds great potential for mitigating climate change, but its storage potential and cost per ton of avoided CO2 emissions are uncertain. We develop a general model to determine the maximum theoretical constraints on both storage potential and injection rate and use it to characterize the economic viability of geosequestration in sandstone saline aquifers. When applied to a representative set of aquifer characteristics, the model yields results that compare favorably with pilot projects currently underway. Over a range of reservoir properties, maximum effective storage peaks at an optimal depth of 1600 m, at which point 0.18-0.31 metric tons can be stored per cubic meter of bulk volume of reservoir. Maximum modeled injection rates predict minima for storage costs in a typical basin in the range of $2-7/ ton CO2 (2005 U.S.$) depending on depth and basin characteristics in our base-case scenario. Because the properties of natural reservoirs in the United States vary substantially, storage costs could in some cases be lower or higher by orders of magnitude. We conclude that available geosequestration capacity exhibits a wide range of technological and economic attractiveness. Like traditional projects in the extractive industries, geosequestration capacity should be exploited starting with the low-cost storage options first then moving gradually up the supply curve.

Assessing the environmental consequences of CO2 leakage from geological CCS: Generating evidence to support environmental risk assessment

At the start of the industrial revolution (circa 1750) the atmospheric concentration of carbon dioxide (CO2) was around 280 ppm. Since that time the burning of fossil fuel, together with other industrial processes such as cement manufacture and changing land use, has increased this value to 400 ppm, for the first time in over 3 million years. With CO2 being a potent greenhouse gas, the consequence of this rise for global temperatures has been dramatic, and not only for air temperatures. Global Sea Surface Temperature (SST) has warmed by 0.4–0.8 °C during the last century, although regional differences are evident (IPCC, 2007). This rise in atmospheric CO2 levels and the resulting global warming to some extent has been ameliorated by the oceanic uptake of around one quarter of the anthropogenic CO2 emissions (Sabine et al., 2004). Initially this was thought to be having little or no impact on ocean chemistry due to the capacity of the ocean’s carbonate buffering system to neutralise the acidity caused when CO2 dissolves in seawater. However, this assumption was challenged by Caldeira and Wickett (2005) who used model predictions to show that the rate at which carbonate buffering can act was far too slow to moderate significant changes to oceanic chemistry over the next few centuries. Their model predicted that since pre-industrial times, ocean surface water pH had fallen by 0.1 pH unit, indicating a 30% increase in the concentration of H+ ions. Their model also showed that the pH of surface waters could fall by up to 0.4 units before 2100, driven by continued and unabated utilisation of fossil fuels. Alongside increasing levels of dissolved CO2 and H+ (reduced pH) an increase in bicarbonate ions together with a decrease in carbonate ions occurs. These chemical changes are now collectively recognised as “ocean acidification”. Concern now stems from the knowledge that concentrations of H+, CO2, bicarbonate and carbonate ions impact upon many important physiological processes vital to maintaining health and function in marine organisms. Additionally, species have evolved under conditions where the carbonate system has remained relatively stable for millions of years, rendering them with potentially reduced capacity to adapt to this rapid change. Evidence suggests that, whilst the impact of ocean acidification is complex, when considered alongside ocean warming the net effect on the health and productivity of the oceans will be detrimental.