Lead, zinc, and antimony contamination of the Rio Chilco-Rio Tupiza drainage system, Southern Bolivia

An intense, but localized rainfall event in February 2003, led to the severe erosion and failure of a tailings disposal impoundment at the Abarόa Antimony Mine in southern Bolivia. The failure released approximately 5,500 m³ of contaminated tailings into the Rio Chilco-Rio Tupiza drainage system. The impacts of the event on sediment quality are examined and compared to contamination resulting from historic mining operations in the headwaters of the basin. Of primary concern are contaminated floodplain soils located along downstream reaches of the Rio Tupiza which were found to contain lead (Pb), zinc (Zn), and antimony (Sb) concentrations that locally exceed Canadian, German, and Dutch guidelines for agricultural use. Spatial patterns in sediment-borne trace metal concentrations, combined with Pb isotopic data, indicate that Pb, Zn, and Sb are derived from three tributary basins draining the Abarόa, Chilcobija, and Tatasi-Portugalete mining districts. Downstream of each tributary, geographical patterns in trace metal concentrations reflect local geomorphic changes throughout the drainage system. Trace metal concentrations within the Rio Chilco decrease rapidly downstream as a result of dilution by uncontaminated sediments and storage of metal enriched particles (e.g., sulfide minerals) in the channel bed as a result of ongoing aggradation. Storage in the floodplains is limited. These processes significantly reduced the dispersal and, thus, the relative environmental affects of tailings eroded from the Abarόa Mine during the 2003 flood. In contrast, storage of Pb, Zn, and Sb in floodplains along the Rio Tupiza is significant, the majority of which is derived from historic mining operations, particularly mining within the Tatasi-Portugalete district.     

Comparison of ERS-SAR and spot data for sugar beet crop cover assessment

The present work was aimed to compare the abilities of radar and optical satellite data to estimate crop canopy cover, which is a key component of productivity estimates. Three ERS-1 SAR images were obtained of East Anglia (UK) in 1995 and one ERS-2 SAR image in 1996. The images covered a study area around the IACR Brooms Barn Sugar Beet Research Institute. Field data comprising radiometric and biophysical measurements of the crop canopy were collected in two fields from June 22 to August 3, 1995 to coincide with ERS-1 SAR overpass dates. In 1996, field data were collected in two fields from June 11 to July 29 on a weekly basis. A previously calibrated version of the water cloud model was inverted to estimate Leaf Area Index (LAI) from ERS-1 and ERS-2 SAR backscatter and soil moisture samples. Canopy cover was estimated from the radar-estimated LAI using a standard exponential relationship that has a well-established coefficient for sugar beet. Radio-metrically and atmospherically corrected data from three SPOT images in 1995 and one SPOT image in 1996 were used to calculate the Optimised Soil Adjusted Vegetation Index (OSAVI), from which crop canopy cover was estimated using a relationship determined previously by canopy modelling. The crop cover values estimated by satellite were in good agreement with those measured on ground with the Parkinson radiometer. Radar data may be able to provide useful estimates of canopy cover for crop production modelling, especially in the case of loss of optical data due to cloud.     

The Port Mouton Shear Zone: intersection of a regional fault with a crystallizing granitoid pluton

The peraluminous tonalite–monzogranite Port Mouton Pluton is a petrological, geochemical, structural, and geochronological anomaly among the many Late Devonian granitoid intrusions of the Meguma Lithotectonic Zone of southern Nova Scotia. The most remarkable structural feature of this pluton is a 4-km-wide zone of strongly foliated (040/subvertical) monzogranites culminating in a narrow (10–30 m), straight, zone of compositionally banded rocks that extends for at least 3 km along strike. The banded monzogranites consist of alternating melanocratic and leucocratic compositions that are complementary to the overall composition of that part of the pluton, suggesting an origin by mineral–melt and mineral–mineral sorting. Biotite and feldspar are strongly foliated in the plane of the compositional bands. These compositional variations and foliations originated by a process of segregation flow during shearing of the main magma with a crystallinity of 55–75%. Subsequent minor brittle fracturing of feldspars, twinning of microcline, development of blocky sub-grains in quartz, and kinking of micas demonstrate overprinting by a high-temperature deformation straddling the monzogranite solidus. Small folds and late sigmoidal dykes indicate dextral movement on the shear zone. This Port Mouton Shear Zone (PMSZ) is approximately co-linear with the only outcrops of Late Devonian mafic intrusions in the area, two of which are syn-plutonic with well-developed mingling textures in the marginal tonalite of the Port Mouton Pluton. Also closely co-linear with the mafic intrusions are a granitoid dyke that extends well beyond the outer contact of the Port Mouton Pluton, a swarm of large aligned angular xenolithic slabs, a zone of thin wispy schlieren banding, a large Be-bearing pegmatite, and a breccia pipe with abundant garnetiferous metapelitic xenoliths. In various ways, the shear zone may control all of these features. The Port Mouton Shear Zone is parallel to many other NE-trending faults and shear zones in the northern Appalachians, probably related to the docking of the Meguma Zone along the Cobequid–Chedabucto Fault system.     

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.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.