Diffusion of Ar in muscovite

Hydrothermal treatment of closely sized muscovite aggregates in a piston-cylinder apparatus induced ⁴⁰Ar^∗ loss that is revealed in ⁴⁰Ar/³⁹Ar step heating spectra. Age spectra and Arrhenius data, however, differ from that expected from a single diffusion length scale. A numerical model of episodic loss assuming the presence of multiple diffusion domains yields excellent fits between synthetic and actual degassing spectra. We used this model to isolate ⁴⁰Ar^∗ loss from the grains that remained intact during hydrothermal treatment at 10 kbar permitting calculation of diffusion coefficients in the temperature range 730-600 °C. Diffusion data generated in this manner yield an activation energy (E) of 63 ± 7 kcal/mol and frequency factor (D_o) of 2.3  cm²/s. Experiments at 20 kbar yield diffusivities lower by about an order of magnitude and correspond to an activation volume of ∼14 cm³/mol. Together, these parameters predict substantially greater retentivity of Ar in muscovite than previously assumed and correspond to a closure temperature (T_c) of 425 °C for a 100 μm radius grain cooling at 10 °C/Ma at 10 kbar (T_c = 405 °C at 5 kbar. Age and log (r/r_o) spectra for the run products show strong correlations indicating that muscovites can retain Ar diffusion boundaries and mechanisms that define their natural retentivity during vacuum step heating. This may permit the application of high resolution, continuous ⁴⁰Ar/³⁹Ar thermochronology to low grade, regionally metamorphosed terranes.     

An overview of Task 6 of the Äspö Task Force: modelling groundwater and solute transport: improved understanding of radionuclide transport in fractured rock

An overview is presented of a 4-year study by the Äspö Task Force on Modelling of Groundwater Flow and Transport of Solutes, whose primary aim was to build a bridge between the approaches used for site characterisation (SC) and performance assessment (PA) associated with nuclear waste repositories. Eleven modelling teams representing six national radioactive waste organisations participated in eight modelling exercises whose objectives were: to assess simplifications used in PA models; to determine how, and to what extent, experimental tracer and flow experiments can constrain the range of parameters used in PA models; to support the design of SC programmes to assure that the results have optimal value for PA calculations; and to improve the understanding of site-specific flow and transport behaviour at different scales using SC models. The modelling tasks were concerned with flow and transport through single and multiple near-planar features on SC and PA timescales, including the diffusion of solutes into multiple immobile zones adjacent to fracture surfaces. In general, tracer tests provide only limited quantitative constraints on retention parameter values relevant to PA but nevertheless provide insight about the flow and transport processes, which is a key element of the bridge between SC and PA.     

A 2000-yr reconstruction of air temperature in the Great Basin of the United States with specific reference to the Medieval Climatic Anomaly

A sediment core representing the past two millennia was recovered from Stella Lake in the Snake Range of the central Great Basin in Nevada. The core was analyzed for sub-fossil chironomids and sediment organic content. A quantitative reconstruction of mean July air temperature (MJAT) was developed using a regional training set and a chironomid-based WA-PLS inference model (r²_jack = 0.55, RMSEP = 0.9°C). The chironomid-based MJAT reconstruction suggests that the interval between AD 900 and AD 1300, corresponding to the Medieval Climate Anomaly (MCA), was characterized by MJAT elevated 1.0°C above the subsequent Little Ice Age (LIA), but likely not as warm as recent conditions. Comparison of the Stella Lake temperature reconstruction to previously published paleoclimate records from this region indicates that the temperature fluctuations inferred to have occurred at Stella Lake between AD 900 and AD 1300 correspond to regional records documenting hydroclimate variability during the MCA interval. The Stella Lake record provides evidence that elevated summer temperature contributed to the increased aridity that characterized the western United States during the MCA.     

Anisotropic dual-continuum representations for multiscale poroelastic materials: Development and numerical modelling

Dual-continuum (DC) models can be tractable alternatives to explicit approaches for the numerical modelling of multiscale materials with multiphysics behaviours. This work concerns the conceptual and numerical modelling of poroelastically coupled dual-scale materials such as naturally fractured rock. Apart from a few exceptions, previous poroelastic DC models have assumed isotropy of the constituents and the dual-material. Additionally, it is common to assume that only one continuum has intrinsic stiffness properties. Finally, little has been done into validating whether the DC paradigm can capture the global poroelastic behaviours of explicit numerical representations at the DC modelling scale. We address the aforementioned knowledge gaps in two steps. First, we utilise a homogenisation approach based on Levin's theorem to develop a previously derived anisotropic poroelastic constitutive model. Our development incorporates anisotropic intrinsic stiffness properties of both continua. This addition is in analogy to anisotropic fractured rock masses with stiff fractures. Second, we perform numerical modelling to test the DC model against fine-scale explicit equivalents. In doing, we present our hybrid numerical framework, as well as the conditions required for interpretation of the numerical results. The tests themselves progress from materials with isotropic to anisotropic mechanical and flow properties. The fine-scale simulations show that anisotropy can have noticeable effects on deformation and flow behaviour. However, our numerical experiments show that the DC approach can capture the global poroelastic behaviours of both isotropic and anisotropic fine-scale representations.     

Evaluating the impact of caprock and reservoir properties on potential risk of CO2 leakage after injection

Numerical models are essential tools in fully understanding the fate of injected CO₂ for commercial-scale sequestration projects and should be included in the life cycle of a project. Common practice involves modeling the behavior of CO₂ during and after injection using site-specific reservoir and caprock properties. Little has been done to systematically evaluate and compare the effects of a broad but realistic range of reservoir and caprock properties on potential CO₂ leakage through caprocks. This effort requires sampling the physically measurable range of caprock and reservoir properties, and performing numerical simulations of CO₂ migration and leakage. In this study, factors affecting CO₂ leakage through intact caprocks are identified. Their physical ranges are determined from the literature from various field sites. A quasi-Monte Carlo sampling approach is used such that the full range of caprock and reservoir properties can be evaluated without bias and redundant simulations. For each set of sampled properties, the migration of injected CO₂ is simulated for up to 200 years using the water–salt–CO₂ operational mode of the STOMP simulator. Preliminary results show that critical factors determining CO₂ leakage rate through caprocks are, in decreasing order of significance, the caprock thickness, caprock permeability, reservoir permeability, caprock porosity, and reservoir porosity. This study provides a function for prediction of potential CO₂ leakage risk due to permeation of intact caprock and identifies a range of acceptable seal thicknesses and permeability for sequestration projects. The study includes an evaluation of the dependence of CO₂ injectivity on reservoir properties.     

A Cambrian mixed carbonate–siliciclastic platform in SW Gondwana: evidence from the Western Sierras Pampeanas (Argentina) and implications for the early Paleozoic paleogeography of the proto-Andean margin

International Journal of Earth Sciences - The Western Sierras Pampeanas (WSP) of Argentina record a protracted geological history from the Mesoproterozoic assembly of the Rodinia supercontinent to...     

Occurrence of juvenile Atlantic sturgeon,Acipenser oxyrhynchus,in the upper tidal Delaware River

Eighty-nine juvenile Atlantic sturgeon,Acipenser oxyrhynchus oxyrhynchus, including 10 recaptures, were collected, tagged and released in the upper tidal Delaware River between July 1981 and December 1984. All were captured in the river channel between Roebling and Trenton, New Jersey, using bottom-set experimental gill nets during daylight hours. The species was present from July through December and collected in increased numbers in the past two years. The fork lengths ranged from 284 to 862 mm (mean±SD = 516 mm±106 mm) and the weights from 140 to 4,250 g (mean±SD = 1,369 g±815 g). Ten juvenile Atlantic sturgeon were recaptured between September 1983 and December 1984 and were at large from one to 418 d. Recapture data suggest that these sturgeon utilize this area annually from July through December, possibly as a nursery. Juvenile Atlantic sturgeon appear to utilize the upper tidal portion of the Delaware River for a much longer period of time and at lower temperatures than in other river systems.     

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.