SPH approach for simulating hydro-mechanical processes with large deformations and variable permeabilities

A simulation framework based on Smoothed Particle Hydrodynamics (SPH) is introduced to model problems involving the interaction between flowing water and soil deformation. Changes in soil porosity and associated permeability are automatically adjusted within this framework. The framework’s capabilities are presented and discussed for three geotechnical problems caused by flowing water. The comparison between simulation results and experiments shows that SPH with the proposed concept is capable of quantitatively simulating the hydro-mechanical processes beyond limit state with satisfactory agreement. To improve the computational stability, a correction procedure and a new algorithm for the selection of the optimal time step are introduced.     

Deformation history of the Naraku Batholith,Mt Isa Inlier,Australia: Implications for pluton ages and geometries from structural study of the Dipvale Granodiorite and Levian Granite

Plutons of the Naraku Batholith were emplaced into Proterozoic metasediments of the northern portion of the Eastern Fold Belt of the Mt Isa Inlier during two intrusive episodes approximately 200 million years apart. Structural relationships and geochronological data suggest that the older plutons (ca 1750 Ma) are contemporaneous with granites of the Wonga Batholith to the west. The Dipvale Granodiorite and the Levian Granite represent these older intrusive phases of the Naraku Batholith, and both contain an intense tectonic foliation, S₁, which is interpreted to have formed during the north‐south shortening associated with D₁ of the Isan Orogeny. The geometry of S₁ form surfaces at the southern end of the Dipvale Granodiorite, and of the previously unrecognised sheeted contact, defines a macroscopic, steeply south‐southwest‐plunging antiform, which was produced by the regional D₂ of the Isan Orogeny. S₁ form surfaces in the Levian Granite define open F₂ folds with wavelengths of several hundred metres. The structural age of emplacement of the Dipvale Granodiorite and the Levian Granite is interpreted to be pre‐ or syn‐ the regional D₁. An intense foliation present in some of the younger (ca 1505 Ma) granites that comprise the bulk of the Naraku Batholith is interpreted to represent S₃ of the Isan Orogeny. Foliations commonly have similar styles and orientations in both the pre‐D₁ and younger plutons. This emphasises the simplicity with which regional fabrics can be, and probably have been, miscorrelated in the Eastern Fold Belt, and that the classification of granites in general on the basis of structural and geometric criteria alone is fraught with danger.     

Pn arrivals and lateral variations of Moho geometry beneath the Kaapvaal craton

Pn arrivals from mining-induced earthquakes on the edge of the Witwatersrand basin show that the P wavespeeds in the uppermost mantle are almost constant throughout most of the Kaapvaal craton. The presence of only small wavespeed variations allows the use of a simple method of estimating crustal thicknesses below the stations of the Kaapvaal broad-band network using Pn times that has been compared with results from receiver functions. One thousand three hundred thirty-seven Pn arrivals were used to derive crustal thicknesses at 46 stations on the Kaapvaal craton. The average crustal thicknesses for 19 centrally located stations on each of the northern and southern regions of the craton that yielded well-constrained thicknesses were 50.52±0.88 km and 38.07±0.85 km, respectively. In contrast, the corresponding average thicknesses determined from receiver functions were 43.58±0.57 km and 37.58±0.70 km, respectively. The systematically lower values for receiver functions in the northern part of the Kaapvaal craton that was affected by the Bushveld magmatism at 2.05 Ga, suggest that the receiver functions do not enable the petrological crust mantle boundary to be reliably resolved due to variations in composition and metamorphic grade in a mafic lower crust. The Pn times also suggest pervasive azimuthal anisotropy with maximum wavespeeds of about 8.40 km/s at azimuths of about 15° and 217° in the northern and southern regions of the craton, respectively, and minimum wavespeeds of about 8.25 km/s.     

Gold partitioning in melt-vapor-brine systems

We used laser-ablation inductively coupled plasma mass spectrometry to measure the solubility of gold in synthetic sulfur-free vapor and brine fluid inclusions in a vapor + brine + haplogranite + magnetite + gold metal assemblage. Experiments were conducted at 800°C, oxygen fugacity buffered at Ni-NiO (NNO), and pressures ranging from 110 to 145 MPa. The wt% NaCl eq. of vapor increases from 2.3 to 19 and that of brine decreases from 57 to 35 with increasing pressure. The composition of the vapors and brines are dominated by NaCl + KCl + FeCl₂ + H₂O. Gold concentrations in vapor and brine decrease from 36 to 5 and 50 to 28 μg/g, respectively, and the calculated vapor:brine partition coefficients for gold decrease from 0.72 to 0.17 as pressure decreases from 145 to 110 MPa. These data are consistent with the thermodynamic boundary condition that the concentration of gold in the vapor and brine must approach a common value as the critical pressure is approached along the 800°C isotherm in the NaCl-KCl-FeCl₂-HCl-H₂O system.We use the equilibrium constant for gold dissolution as AuOH⁰, extrapolated from lower temperature and overlapping pressure range, to calculate expected concentrations of AuOH⁰ in our experimental vapors. These calculations suggest that a significant quantity of gold in our experimental vapors is present as a non-hydroxide species. Possible chloridogold(I) species are hypothesized based on the positively correlated gold and chloride concentrations in our experimental vapors. The absolute concentration of gold in our synthetic vapor, brine, and melt and calculated mass partition coefficients for gold between these physicochemically distinct magmatic phases suggests that gold solubility in aqueous fluids is a function of aqueous phase salinity, specifically total chloride concentration, at magmatic conditions. However, though we highlight here the effect of salinity, the combination of our data with data sets from lower temperatures evinces a significant decrease in gold solubility as temperature drops from 800°C to 600°C. This decrease in solubility has implications for gold deposition from ascending magmatic fluids.     

The geochemistry of primary and weathered oil shale and coquina across the Julia Creek vanadium deposit (Queensland, Australia)

A significant resource of vanadium and molybdenum exists near Julia Creek, Australia, where the middle Cretaceous organic-rich Toolebuc Formation lies between 0 and 25 m of the surface. We present and discuss a comprehensive geochemical study of the Toolebuc Formation and its enclosing stratigraphy near Julia Creek to understand this ore deposit. V and Mo contents in fresh facies are strongly associated with total organic carbon (TOC) contents, but not with Al or CaCO₃; this suggests that V and Mo were originally concentrated in the organic fraction. However, chemical extractions using H₂O₂ indicate that Mo was originally concentrated in pyrite. The data also suggest that V was mobilised from organic matter during early diagenesis and became associated with clays as little V was extracted by H₂O₂ in the fresh samples. TOC contents in the Toolebuc Formation were removed during weathering, residually enriching trace metals including V and Mo, and as a result, the TOC relationship with V and Mo disintegrates. With weathering, both V and Mo predominantly became associated with iron oxide/hydroxide phases (and possibly other unidentified phases) as these elements in the weathered facies were highly soluble in the sodium citrate–sodium dithionite digestion. Large shale-hosted V and Mo deposits such as Julia Creek offer a potentially viable alternative to the currently mined magnetite-hosted deposits. A thorough understanding of the formation and host mineral phases for V and Mo of these shale deposits, however, is critical to ensure that these valuable metals can be feasibly extracted.     

Formation of refractory inclusions by evaporation of condensate precursors

Berman’s (1983) activity-composition model for CaO-MgO-Al₂O₃-SiO₂ liquids is used to calculate the change in bulk chemical and isotopic composition during simultaneous cooling, evaporation, and crystallization of droplets having the compositions of reasonable condensate precursors of Types A and B refractory inclusions in CV3 chondrites. The degree of evaporation of MgO and SiO₂, calculated to be faithfully recorded in chemical and isotopic zoning of individual melilite crystals, is directly proportional to evaporation rate, which is a sensitive function of P_H2, and inversely proportional to the droplet radius and cooling rate. When the precursors are partially melted in pure hydrogen at peak temperatures in the vicinity of the initial crystallization temperature of melilite, their bulk chemical compositions evolve into the composition fields of refractory inclusions, mass-fractionated isotopic compositions of Mg, Si, and O are produced that are in the range of the isotopic compositions of natural inclusions, and melilite zoning profiles result that are similar to those observed in real inclusions. For droplets of radius 0.25 cm evaporating at P_H2 = 10⁻⁶ bar, precursors containing 8 to 13 wt.% MgO and 20 to 23% SiO₂ evolve into objects similar to compact Type A inclusions at cooling rates of 2 to 12 K/h, depending on the precise starting composition. Precursors containing 13 to 14 wt.% MgO and 23 to 26% SiO₂ evolve into objects with the characteristics of Type B1 inclusions at cooling rates of 1.5 to 3 K/h. The relatively SiO₂-poor members of the Type B2 group can be produced from precursors containing 14 to 16 wt.% MgO and 27 to 33% SiO₂ at cooling rates of <1 K/h. Type B2’s containing 27 to 35 wt.% SiO₂ and <12% MgO require precursors with higher SiO₂/MgO ratios at MgO > 15% than are found on any condensation curve. The characteristics of fluffy Type A inclusions, including their reversely zoned melilite, can only be understood in the context of this model if they contain relict melilite.     

Geochemical studies in the drainage basin of the Rio Vouga,Portugal. IV. Impact of land use on the hydrogeochemistry of natural waters in the Vouzela region

A hydrogeochemical study of the impact of land use on the composition of natural waters in the Vouzela region in northern Portugal was carried out during the summer of 1983. Water samples were collected from spring and major streams in the area and analyzed for major constituents and some trace elements. Analysis of variance and subsequent pairwise contrast tests demonstrated that waters from agricultural areas are significantly enriched in constiluents such as Na, K, Ca, Mg, Sr, Cl, SO₄, and NO₃. High concentrations of these ions in agricultural areas are the result of both the application of fertilizers and enhanced evaporation by the intensive irrigation of agricultural lands. Constitutents such as H₄SiO₄, HCO₃, F, and probably Li, which are typically related to mineral weathering, were not affected by land use. In waters from forested areas the concentration of most constituents was about 20% higher than in waters from uncultivated areas with mainly a grass cover. This reflects the difference in the evapotranspiration of these two vegetation types. The chemistry of the Rio Zela clearly reflects differences in land use in the Rio Zela valley.     

Diversity of recent planktonic foraminifera in the southern Indian Ocean and Late Pleistocene paleotemperatures

Planktonic foraminiferal assemblages have been examined in 25 trigger core top samples and 51 piston core top samples collected between latitudes 28° S and 55° S and longitudes 79° E and 120° E from the southern Indian Ocean during cruises of the U.S.N.S. Eltanin. Samples taken from water depths exceeding 4000 m and/or showing evidence of calcium carbonate dissolution were eliminated from further analysis. The final piston core data set consists of 34 samples; the trigger core data set containing 21 samples. A close relationship exists between changes in the planktonic foraminiferal assemblages in the surface sediments and surface water temperatures. Species diversity values were computed for each of the core top assemblages using the Shannon-Wiener Index and the Brillouin Index, each of which takes into consideration the number of species and the proportionment of individuals among the species. The Shannon and Brillouin diversity values for all samples are positively correlated (correlation coefficient (r) = +.999). Regression analysis of latitude versus Shannon diversity values in the trigger core samples clearly shows a decrease in diversity with increasing latitude (r = ?.979). Furthermore, a strong correlation (r = +.977) exists between decreasing species diversity (Shannon) and decreasing average summer-winter temperature of the overlying surface waters. A paleotemperature equation derived from the relationship of diversity in trigger core samples and surface water temperature was used to generate paleotemperature curves for five trigger cores and a 6 m piston core of Late Pleistocene age, located beneath the present position of the Subtropical Convergence. A 7–8° C temperature range is suggested between the interglacial and glacial episodes in this Late Pleistocene sequence, and probably reflects latitudinal shifts of the Subtropical Convergence and Australasian Front during the Late Pleistocene.     

Cation disorder in garnets along the Mg3Al2Si3O12-Mg4Si4O12 join: an infrared,Raman and NMR study

We have obtained infrared and Raman spectra for garnets synthesized at high (static) pressures and temperatures along the join Mg₃Al₂Si₃O₁₂ (pyrope) — Mg₄Si₄O₁₂ (magnesium majorite). The vibrational spectra of Mg-majorite show a large number of additional weak peaks compared with the spectra of cubic pyrope garnet, consistent with tetragonal symmetry for the MgSiO₃ garnet phase. The Raman bands for this phase show no evidence for line broadening, suggesting that Mg and Si are ordered on octahedral sites in the garnet. The bands for the intermediate garnet compositions are significantly broadened compared with the end-members pyrope and Mg-majorite, indicating cation disorder in the intermediate phases. Solid state ²⁷Al NMR spectroscopy for pyrope and two intermediate compositions show that Al is present only on octahedral sites, so the cation disorder is most likely confined to Mg-Al-Si mixing on the octahedral sites. We have also obtained a Raman spectrum for a natural, shock-produced (Fe,Mg) majorite garnet. The sharp Raman peaks suggest little or no cation disorder in this sample.