Matrix Diffusion Effects on Contaminant Migration from an Injection Well in Fractured Sandstone

Deep-well injection into fractured sandstone is an option for the disposal of contaminated mine dewatering discharge from an open pit uranium mine. As part of the assessment of potential contaminant migration from deep-well injection, the effect of matrix diffusion was evaluated. An analytical mathematical model was developed for the simulation of the radial movement of a contaminant front away from an injection point under steady flow conditions in a planar fracture with uniform properties. The model includes the effects of advection in the fracture, diffusion of contaminants from the fracture into the rock matrix, and equilibrium adsorption on the fracture surface as well as in the rock matrix. Effective diffusion coefficients obtained from laboratory experiments on 11 intact core samples varied from 3.4 × 10⁻⁸ to 3.2 × 10⁻⁷ cm²/s. Model simulations were made with diffusion coefficient values in this range and with single-fracture injection rates estimated from fracture frequencies in boreholes, and from bulk hydraulic conductivity values obtained from field tests. Because of matrix diffusion, the rate of outward movement of the front of the nonreactive contaminants from the injection well is much slower than the rate of water flow in the fractures. Simulations of the movement of contaminants that undergo adsorption indicate that even a small distribution coefficient for the rock matrix causes the contaminants to remain very close to the injection well during the one-year period. The results of the simplified model demonstrate that matrix diffusion is an important process that cannot be neglected in the assessment of a waste disposal scheme located in fractured porous rock. However, in order to make a definitive assessment of the capability of matrix diffusion and associated matrix adsorption to significantly limit the extent of contaminant migration around injection wells, it would be necessary to conduct field tests such as a preliminary or experimental injection.     

The Ricker Hills Tillite provides evidence of Oligocene warm-based glaciation in Victoria Land, Antarctica

The relationship between the Ricker Hills Tillite (RHT), which represents the northernmost outcrop of lithified continental glacial deposits in Victoria Land, is discussed with respect to the glacial landscape assemblage of the Ricker Hills, a nunatak at the internal border of the Transantarctic Mountains. A warm-based ice sheet deposited the tillite and induced syn- to post-depositional glacial deformation under wet conditions both of the tillite and of the bedrock. The thickness of the ice sheet on the nunatak is estimated to have been 600 m, at most. The area had been deeply eroded before deposition of the RHT as documented by the low elevation of tillite outcrops located in overdeepened depressions of the nunatak. Micropaleontological analysis evidences only the presence of Permian to Jurassic palynomorphs. X-ray diffraction and SEM–EDS analyses of clay minerals in the RHT indicate continental chemical weathering under wet conditions after the RHT deposition. As documented by clay mineral assemblage variation in CRP drillholes, the progressive cooling of the Antarctic continent allowed chemical weathering in “warm” conditions until the Late Oligocene period in southern Victoria Land, leading to a chronological constrain for RHT deposition. Conservatively estimating the sea level to have been between the tillite outcrops and the erosional trimline limiting horns in the Ricker Hills, at the time of RHT deposition, we suggest that the maximum uplift of the area would not have exceeded 900–1500 m since at least Late Oligocene.     

Determination of Rare Earth Elements and Yttrium in Nine Geochemical Reference Samples Using a Novel Group Separation Procedure Involving Mixed-Acid Elution Ion-Exchange Chromatography

Rare earth data for nine well-known geochemical reference samples were obtained using a rapid and efficient cation exchange chromatographic separation procedure. It involves the use of a mixed-acid eluent (HCI and HNO3) on a 15 x 1 cm column of Bio-Rad AG50–X8 (200–400 mesh) to remove matrix and many trace elements from rock solutions. The rare earth elements are removed from the column using 60 ml 7M HNO3. The pre-concentration procedure described is routinely used in this laboratory for the determination of the rare earth elements by mass spectrometry-isotope dilution, and has been satisfactorily used for their determination using inductively-coupled plasma-source atomic emission spectrometry.     

Numerical simulation of mud-flows impacting structures

The study of the interaction of mud-flows with obstacles is important to define inundation zones in urban areas and to design the possible structural countermeasures. The paper numerically investigates the impact of a mud-flow on rigid obstacles to evaluate the force acting on them using two different depth-integrated theoretical models, Single-Phase Model(SPM) and Two-Phase Model(TPM), to compare their performance and limits. In the first one the water-sediment mixture is represented as a homogeneous continuum described by a shearthinning power-law rheology. Alternatively, the twophase model proposed by Di Cristo et al in 2016 is used, which separately accounts for the liquid and solid phases. The considered test cases are represented by a 1D landslide flowing on a steep slope impacting on a rigid wall and a 2D mud dam-break flowing on a horizontal bottom in presence of single and multiple rigid obstacles. In the 1D test case, characterized by a very steep slope, the Two-Phase Model predicts the separation between the two phases with a significant longitudinal variation of the solid concentration. In this case the results indicate appreciable differences between the two models in the estimation of both the wave celerity and the magnitude of the impact, with an overestimation of the peak force when using the Single-Phase Model. In the 2D test-cases, where the liquid and solid phases remain mixed, even if the flow fields predicted by the two models present some differences, the essential features of the interaction with the obstacles, along with the maximum impact force, are comparable.     

Trace and minor elements in sphalerite: A LA-ICPMS study

Sphalerite is an important host mineral for a wide range of minor and trace elements. We have used laser-ablation inductively coupled mass spectroscopy (LA-ICPMS) techniques to investigate the distribution of Ag, As, Bi, Cd, Co, Cu, Fe, Ga, Ge, In, Mn, Mo, Ni, Pb, Sb, Se, Sn and Tl in samples from 26 ore deposits, including specimens with wt.% levels of Mn, Cd, In, Sn and Hg. This technique provides accurate trace element data, confirming that Cd, Co, Ga, Ge, In, Mn, Sn, As and Tl are present in solid solution. The concentrations of most elements vary over several orders of magnitude between deposits and in some cases between single samples from a given deposit. Sphalerite is characterized by a specific range of Cd (typically 0.2-1.0 wt.%) in each deposit. Higher Cd concentrations are rare; spot analyses on samples from skarn at Baisoara (Romania) show up to 13.2 wt.% (Cd²⁺ ↔ Zn²⁺ substitution). The LA-ICPMS technique also allows for identification of other elements, notably Pb, Sb and Bi, mostly as micro-inclusions of minerals carrying those elements, and not as solid solution. Silver may occur both as solid solution and as micro-inclusions. Sphalerite can also incorporate minor amounts of As and Se, and possibly Au (e.g., Magura epithermal Au, Romania). Manganese enrichment (up to ∼4 wt.%) does not appear to enhance incorporation of other elements. Sphalerite from Toyoha (Japan) features superimposed zoning. Indium-sphalerite (up to 6.7 wt.% In) coexists with Sn-sphalerite (up to 2.3 wt.%). Indium concentration correlates with Cu, corroborating coupled (Cu⁺In³⁺) ↔ 2Zn²⁺ substitution. Tin, however, correlates with Ag, suggesting (2Ag⁺Sn⁴⁺) ↔ 3Zn²⁺ coupled substitution. Germanium-bearing sphalerite from Tres Marias (Mexico) contains several hundred ppm Ge, correlating with Fe. We see no evidence of coupled substitution for incorporation of Ge. Accordingly, we postulate that Ge may be present as Ge²⁺ rather than Ge⁴⁺. Trace element concentrations in different deposit types vary because fractionation of a given element into sphalerite is influenced by crystallization temperature, metal source and the amount of sphalerite in the ore. Epithermal and some skarn deposits have higher concentrations of most elements in solid solution. The presence of discrete minerals containing In, Ga, Ge, etc. also contribute to the observed variance in measured concentrations within sphalerite.     

Site classification using equivalent soil profiles for building-liquefaction interaction

The seismic behaviour of a building on a liquefiable deposit is a complex interaction which involves quantifying both shaking induced damage and permanent ground deformation-related damage. In this paper the key parameters that influence both surface shaking and foundation settlements have been identified as the depth, thickness and liquefaction resistance of an equivalent liquefiable layer. These parameters can be used to develop an ‘equivalent soil profile’ that is analogous to the equivalent single degree-of-freedom that reduces the complexity of the dynamic response of a building into comparable and easily understood quantities. The equivalent soil profile is quantified independent of the seismic hazard, making it compatible with performance based design and assessment frameworks such that the building and soil profile can be directly assessed at different levels of seismic hazard. Several numerical studies are presented that demonstrate the influence of these key parameters on the ground surface shaking and foundation settlement. A set of criteria are proposed for classifying soil profiles into 22 different soil classes for regional loss assessment. An algorithm was developed for automatically fitting the equivalent soil profile to a cone penetration test trace and issues with the fitting are discussed. Field reconnaissance was undertaken to collect additional data to support existing datasets on the performance of buildings in Adapazari, during the 1999 Kocaeli, Turkey, earthquake (Mw = 7.4). The field case history data was used to investigate the correlation between the depth, thickness and liquefaction resistance of an equivalent liquefiable layer, on the extent of foundation permanent deformation. The case history data showed that in general a shallow, thick and weak liquefiable layer near the surface results in significant settlement but a lack of data for buildings on non-liquefiable deposits and the additional complexities involved with real buildings and soil deposits, meant that the trends observed in the idealised numerical models could not identified in the field case history data set.

     

Late Holocene evolution of Tiber river delta and geoarchaeology of Claudius and Trajan Harbor,Rome

This study provides new information on the evolution of the Tiber delta during the late Holocene, and describes the record of sedimentation in the Port of Claudius and Trajan. The Tiber flowed into the Tyrrhenian Sea in the area of the Port of Claudius and Trajan between the 8th and 5th centuries B.C. The harbor basins excavated by the Romans during the 1st and 2nd centuries A.D. were filled with sediments containing saltwater fauna, then brackish‐water fauna, and finally by alluvial and marsh deposits. The Port of Claudius and Trajan was silted in as a consequence of episodes of alluviation that correlate with environmental events, possibly triggered by global climate changes. © 2009 Wiley Periodicals, Inc.