Spatial Modeling of Geometallurgical Properties: Techniques and a Case Study

High-resolution spatial numerical models of metallurgical properties constrained by geological controls and more extensively by measured grade and geomechanical properties constitute an important part of geometallurgy. Geostatistical and other numerical techniques are adapted and developed to construct these high-resolution models accounting for all available data. Important issues that must be addressed include unequal sampling of the metallurgical properties versus grade assays, measurements at different scale, and complex nonlinear averaging of many metallurgical parameters. This paper establishes techniques to address each of these issues with the required implementation details and also demonstrates geometallurgical mineral deposit characterization for a copper–molybdenum deposit in South America. High-resolution models of grades and comminution indices are constructed, checked, and are rigorously validated. The workflow demonstrated in this case study is applicable to many other deposit types.     

Dissolved organic matter variations in coastal plain wetland watersheds: The integrated role of hydrological connectivity,land use,and seasonality

Dissolved organic matter (DOM) is integral to fluvial biogeochemical functions, and wetlands are broadly recognized as substantial sources of aromatic DOM to fluvial networks. Yet how land use change alters biogeochemical connectivity of upland wetlands to streams remains unclear. We studied depressional geographically isolated wetlands on the Delmarva Peninsula (USA) that are seasonally connected to downstream perennial waters via temporary channels. Composition and quantity of DOM from 4 forested, 4 agricultural, and 4 restored wetlands were assessed. Twenty perennial streams with watersheds containing wetlands were also sampled for DOM during times when surface connections were present versus absent. Perennial watersheds had varying amounts of forested wetland (0.4–82%) and agricultural (1–89%) cover. DOM was analysed with ultraviolet–visible spectroscopy, fluorescence spectroscopy, dissolved organic carbon (DOC) concentration, and bioassays. Forested wetlands exported more DOM that was more aromatic‐rich compared with agricultural and restored wetlands. DOM from the latter two could not be distinguished suggesting limited recovery of restored wetlands; DOM from both was more protein‐like than forested wetland DOM. Perennial streams with the highest wetland watershed cover had the highest DOC levels during all seasons; however, in fall and winter when temporary streams connect forested wetlands to perennial channels, perennial DOC concentrations peaked, and composition was linked to forested wetlands. In summer, when temporary stream connections were dry, perennial DOC concentrations were the lowest and protein‐like DOM levels the highest. Overall, DOC levels in perennial streams were linked to total wetland land cover, but the timing of peak fluxes of DOM was driven by wetland connectivity to perennial streams. Bioassays showed that DOM linked to wetlands was less available for microbial use than protein‐like DOM linked to agricultural land use. Together, this evidence indicates that geographically isolated wetlands have a significant impact on downstream water quality and ecosystem function mediated by temporary stream surface connections.     

Response of the Atlantic meridional overturning circulation to a reversal of greenhouse gas increases

The reversibility of the Atlantic meridional overturning circulation (AMOC) is investigated in multi-model experiments using global climate models (GCMs) where CO₂ concentrations are increased by 1 or 2 % per annum to 2× or 4× preindustrial conditions. After a period of stabilisation the CO₂ is decreased back to preindustrial conditions. In most experiments when the CO₂ decreases, the AMOC recovers before becoming anomalously strong. This "overshoot" is up to an extra 18.2Sv or 104 % of its preindustrial strength, and the period with an anomalously strong AMOC can last for several hundred years. The magnitude of this overshoot is shown to be related to the build up of salinity in the subtropical Atlantic during the previous period of high CO₂ levels. The magnitude of this build up is partly related to anthropogenic changes in the hydrological cycle. The mechanisms linking the subtropical salinity increase to the subsequent overshoot are analysed, supporting the relationship found. This understanding is used to explain differences seen in some models and scenarios. In one experiment there is no overshoot because there is little salinity build up, partly as a result of model differences in the hydrological cycle response to increased CO₂ levels and partly because of a less aggressive scenario. Another experiment has a delayed overshoot, possibly as a result of a very weak AMOC in that GCM when CO₂ is high. This study identifies aspects of overshoot behaviour that are robust across a multi-model and multi-scenario ensemble, and those that differ between experiments. These results could inform an assessment of the real-world AMOC response to decreasing CO₂.     

Short-term residence and movement patterns of the annular seabream Diplodus annularis in a temperate marine reserve

The short-term movements of a small temperate fish, the annular seabream Diplodus annularis (Linnaeus 1758), were examined using standard tag-recapture and passive acoustic telemetry in Palma Bay (NW Mediterranean), a marine protected area (MPA). The study aimed to provide valuable information for assessing the recreational fishery and its results suggest that MPAs can be used to protect the adult stock of D. annularis. All the fish tagged with standard tags were recaptured near the release locations, with a maximum distance of ∼300 m. The maximum time between release and recapture was 185 d. Two different arrays of acoustic receivers were deployed, one in 2008 and another in 2009, within the MPA. Twenty adults were surgically tagged with acoustic transmitters. Fish monitored in 2008 (n = 12) were translocated from the point of capture to analyse the movement behaviour after artificial displacement. Upon release at displaced locations, 67% of the fish moved towards the original capture location using a time of return that ranged from 0.75 to 15.25 h. Fish monitored in 2009 (n = 8) were released at the point of capture. They showed high site fidelity with a maximum period of 27 d between the first time and the last time they were detected.     

Flow directions in ash-flow tuffs: a comparison of geological and magnetic susceptibility measurements, Tshirege member (upper Bandelier Tuff), Valles caldera, New Mexico, USA

This study concludes that the elongation axis (K ₁) of the ellipsoid of anisotropic magnetic susceptibility (AMS) is a suitable proxy for flow axis in ashflow tuffs. 153 oriented samples (176 specimens) were studied from 18 sites in the 1.1 Ma Tshirege member of the Bandelier Tuff. These sites are distributed around the Valles caldera at distances of 5–25 km outside of the rim.K ₁ axes correlate well with postulated radial flow axes at 13 sites.K ₁ also agrees with measured geological flow indicators, mainly imbricated larger clasts, at 7 sites. At 2 of the 5 sites where significant disagreement is seen between theoretical radial flow directions and measuredK ₁ axes, theK ₁ axes correspond well with geological flow indicators, indicating that the divergence of flow from the predicted radial flow pattern is real. Two major topographic buttresses are suggested as the cause of flow divergence for the Tshirege ash flows: the San Pedro buttress northwest of the caldera, and the San Miguel buttress in the southeast. In situK ₁ axes plunge about 7° toward the source at two-thirds of the sites; therefore the plunge ofK ₁ is a plausible in situ indicator for thedirection of flow. Multiple flow zones in sections of several meters thickness indicate changes of flow direction that are both rapid and large during ash-flow emplacement. These observations raisre the question of how best to represent mean flow directions in ash-flow sheets: by eigenvector methods, by vector-sum methods, or by modes. A method for measuring imbrication of larger clasts using apparent dips in vertical joints is outlined. Imbrication, determined in this way at one-third of the sites, dips toward the source, i.e., up-flow. The minimum (K ₃) axis of the AMS ellipsoid correlates with the flow foliation rather than with the larger clast imbrication. The flow axes of ash flows correspond with theK ₁ axes, not with the declination ofK ₃ axes as suggested by some authors. Initial dip of the sampled ash flows is not large and does not affect the paleomagnetic remanence direction, which is reversed with a mean ofD=173.5°,I=-38.4°, ₉₅=3.4°N=18. This mean is not different at the 95% confidence level from that of earlier workers. The mean pole, at 098.0°E, 74.8°N,A ₉₅=3.3°,N=18, is about 15° far-sided relative to the expected time-averaged geomagnetic pole, suggesting a history of emplacement too short to adequately average secular variation.     

Testing solar forcing of pervasive Holocene climate cycles

The temporal and spatial extent of Holocene climate change is an area of considerable uncertainty, with solar forcing recently proposed to be the origin of cycles identified in the North Atlantic region. To address these issues we have developed an annually resolved record of changes in Irish bog tree populations over the last 7468 years which, together with radiocarbon‐dated bog and lake‐edge populations, extend the dataset back to ～9000 yr ago. The Irish trees underpin the internationally accepted radiocarbon calibration curve, used to derive a proxy of solar activity, and allow us to test solar forcing of Holocene climate change. Tree populations and age structures provide unambiguous evidence of major shifts in Holocene surface moisture, with a dominant cyclicity of 800 yr, similar to marine cycles in the North Atlantic, indicating significant changes in the latitude and intensity of zonal atmospheric circulation across the region. The cycles, however, are not coherent with changes in solar activity (both being on the same absolute timescale), indicating that Holocene North Atlantic climate variability at the millennial and centennial scale is not driven by a linear response to changes in solar activity. Copyright © 2005 John Wiley & Sons, Ltd.     

Variations in dissolved CO2 and CH4 in a first‐order stream and catchment: an investigation of soil–stream linkages

Spatial and seasonal variations in CO₂ and CH₄ concentrations in streamwater and adjacent soils were studied at three sites on Brocky Burn, a headwater stream draining a peatland catchment in upland Britain. Concentrations of both gases in the soil atmosphere were significantly higher in peat and riparian soils than in mineral soils. Peat and riparian soil CO₂ concentrations varied seasonally, showing a positive correlation with air and soil temperature. Streamwater CO₂ concentrations at the upper sampling site, which mostly drained deep peats, varied from 2·8 to 9·8 mg l^?1 (2·5 to 11·9 times atmospheric saturation) and decreased markedly downstream. Temperature‐related seasonal variations in peat and riparian soil CO₂ were reflected in the stream at the upper site, where 77% of biweekly variation was explained by an autoregressive model based on: (i) a negative log‐linear relationship with stream flow; (ii) a positive linear relationship with soil CO₂ concentrations in the shallow riparian wells; and (iii) a negative linear relationship with soil CO₂ concentrations in the shallow peat wells, with a significant 2‐week lag term. These relationships changed markedly downstream, with an apparent decrease in the soil–stream linkage and a switch to a positive relationship between stream flow and stream CO₂. Streamwater CH₄ concentrations also declined sharply downstream, but were much lower (<0·01 to 0·12 mg l^?1) than those of CO₂ and showed no seasonal variation, nor any relationship with soil atmospheric CH₄ concentrations. However, stream CH₄ was significantly correlated with stream flow at the upper site, which explained 57% of biweekly variations in dissolved concentrations. We conclude that stream CO₂ can be a useful integrative measure of whole catchment respiration, but only at sites where the soil–stream linkage is strong. Copyright © 2004 John Wiley & Sons, Ltd.     

Dissolution/precipitation kinetics of boehmite and gibbsite: Application of a pH-relaxation technique to study near-equilibrium rates

The dissolution and precipitation rates of boehmite, AlOOH, at 100.3 °C and limited precipitation kinetics of gibbsite, Al(OH)₃, at 50.0 °C were measured in neutral to basic solutions at 0.1 molal ionic strength (NaCl + NaOH + NaAl(OH)₄) near-equilibrium using a pH-jump technique with a hydrogen-electrode concentration cell. This approach allowed relatively rapid reactions to be studied from under- and over-saturation by continuous in situ pH monitoring after addition of basic or acidic titrant, respectively, to a pre-equilibrated, well-stirred suspension of the solid powder. The magnitude of each perturbation was kept small to maintain near-equilibrium conditions. For the case of boehmite, multiple pH-jumps at different starting pHs from over- and under-saturated solutions gave the same observed, first order rate constant consistent with the simple or elementary reaction: .

This relaxation technique allowed us to apply a steady-state approximation to the change in aluminum concentration within the overall principle of detailed balancing and gave a resulting mean rate constant, (2.2 ± 0.3) × 10⁻⁵ kg m⁻² s⁻¹, corresponding to a 1σ uncertainty of 15%, in good agreement with those obtained from the traditional approach of considering the rate of reaction as a function of saturation index. Using the more traditional treatment, all dissolution and precipitation data for boehmite at 100.3 °C were found to follow closely the simple rate expression:

R_net,boehmite=10^-5.485{m_OH^-}{1-exp(ΔG_r/RT)}, with R_net in units of mol m⁻² s⁻¹. This is consistent with Transition State Theory for a reversible elementary reaction that is first order in OH⁻ concentration involving a single critical activated complex. The relationship applies over the experimental ΔG_r range of 0.4–5.5 kJ mol⁻¹ for precipitation and −0.1 to −1.9 kJ mol⁻¹ for dissolution, and the pH_m ≡ −log(mH⁺) range of 6–9.6. The gibbsite precipitation data at 50 °C could also be treated adequately with the same model:R_net,gibbsite=10^-5.86{m_OH^-}{1-exp(ΔG_r/RT)}, over a more limited experimental range of ΔG_r (0.7–3.7 kJ mol⁻¹) and pH_m (8.2–9.7).