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.     

Thermal evolution and physics of melt extraction on the ureilite parent body

We develop a physical model of the thermal history of the ureilite parent body (UPB) that numerically tracks the history of its heating, hydration, dehydration, partial melting and smelting as a function of its formation time and the initial values of its composition, formation temperature and water ice content. Petrologic and chemical data from the main group (non-polymict) ureilite meteorites, which sample the interior of the UPB between depths corresponding to pressures in the range 3-10 MPa, are used to constrain the model. We find that to achieve the ∼30% melting inferred for ureilites from all sampled depths, the UPB must have had a radius between ∼80 and ∼130 km and must have accreted about 0.55 Ma after CAI formation. Melting began in the body at ∼1 Ma after CAI, and the time at which 30% melting was reached varied with depth in the asteroid but was always between ∼4.5 and ∼5.8 Ma after CAI. The total rate at which melt was produced in the UPB varied from more than 100 m³ s⁻¹ in the very early stages of melting at ∼1 Ma after CAI to ∼5 m³ s⁻¹ between 2 and 3 Ma after CAI, decreasing to extremely small values as the end of melting was approached beyond ∼5 Ma. Although the initial period of high melt production occupied only a short time around 1 Ma after CAI, it corresponded to ∼half (16%) of total silicate melting, and all strictly basaltic (i.e. plagioclase-saturated) melts must have been produced during this period.A very efficient melt transport network, consisting of a hierarchy of veins and larger pathways (dikes), developed quickly at the start of melting, ensuring rapid (timescales of months) transport of any single parcel of melt to shallow levels, thus ensuring that chemical interaction between melts and the rocks through which they subsequently passed was negligible. Volatile (mainly carbon monoxide) production due to smelting began at the start of silicate melting in the shallowest parts of the UPB and at later times at greater depths. Except at the very start and very end of melting, the volatile content of the melts produced was always high - generally between 15 and 35 mass % - and most of the melt produced was erupted at the surface of the UPB with speeds well in excess of the escape velocity and was lost into space. However, we show that 30% melting at the 3 MPa pressure level was only possible if ∼15% of the total melt produced in the asteroid was retained as a small number (∼5) of very extensive, sill-like intrusions centered at a depth of ∼7 km below the surface, near the base of the ∼8 km thick outer crust of the asteroid that was maintained at temperatures below the basalt solidus by conductive heat loss to the surface. The horizontal extents of these sills occupied about 75% of the surface area of the UPB, and the sills acted as buffers between the steady supply of melt from depth and the intermittent explosive eruption of the melt into space. We infer that samples from these intrusions are preserved as the rare feldspathic (loosely basaltic) clasts in polymict ureilites, and show that the cooling histories of the sills are consistent with these clasts reaching isotopic closure at ∼5 Ma after CAI, as given by ²⁶Al-²⁶Mg, ⁵³Mn-⁵³Cr and Pb-Pb age dates.     

The internal structures and densities of asteroids

Abstract— Four asteroidal bodies (the Martian satellites Phobos and Deimos and the main-belt asteroids 243 Ida and 253 Mathilde) have now been the subjects of sufficiently close encounters by spacecraft that the masses and sizes and, hence, the densities of these bodies can be estimated to ～10%. All of these asteroids are significantly less dense than most members of the classes of meteorites identified as being compositionally most nearly similar to them on the basis of spectral characteristics. We show that two processes can act, independently or in concert, during the evolutionary histories of asteroids to produce a low bulk density. One of these processes is the result of one or more impact events and can affect any asteroid type, whereas the other can occur only for certain types of small asteroids that have undergone aqueous alteration.     

Towards a climate event stratigraphy for New Zealand over the past 30 000 years (NZ‐INTIMATE project)

It is widely recognised that the acquisition of high‐resolution palaeoclimate records from southern mid‐latitude sites is essential for establishing a coherent picture of inter‐hemispheric climate change and for better understanding of the role of Antarctic climate dynamics in the global climate system. New Zealand is considered to be a sensitive monitor of climate change because it is one of a few sizeable landmasses in the Southern Hemisphere westerly circulation zone, a critical transition zone between subtropical and Antarctic influences. New Zealand has mountainous axial ranges that amplify the climate signals and, consequently, the environmental gradients are highly sensitive to subtle changes in atmospheric and oceanic conditions. Since 1995, INTIMATE has, through a series of international workshops, sought ways to improve procedures for establishing the precise ages of climate events, and to correlate them with high precision, for the last 30 000 calendar years. The NZ‐INTIMATE project commenced in late 2003, and has involved virtually the entire New Zealand palaeoclimate community. Its aim is to develop an event stratigraphy for the New Zealand region over the past 30 000 years, and to reconcile these events against the established climatostratigraphy of the last glacial cycle which has largely been developed from Northern Hemisphere records (e.g. Last Glacial Maximum (LGM), Termination I, Younger Dryas). An initial outcome of NZ‐INTIMATE has been the identification of a series of well‐dated, high‐resolution onshore and offshore proxy records from a variety of latitudes and elevations on a common calendar timescale from 30 000 cal. yr BP to the present day. High‐resolution records for the last glacial coldest period (LGCP) (including the LGM sensu stricto) and last glacial–interglacial transition (LGIT) from Auckland maars, Kaipo and Otamangakau wetlands on eastern and central North Island, marine core MD97‐2121 east of southern North Island, speleothems on northwest South Island, Okarito wetland on southwestern South Island, are presented. Discontinuous (fragmentary) records comprising compilations of glacial sequences, fluvial sequences, loess accumulation, and aeolian quartz accumulation in an andesitic terrain are described. Comparisons with ice‐core records from Antarctica (EPICA Dome C) and Greenland (GISP2) are discussed. A major advantage immediately evident from these records apart from the speleothem record, is that they are linked precisely by one or more tephra layers. Based on these New Zealand terrestrial and marine records, a reasonably coherent, regionally applicable, sequence of climatically linked stratigraphic events over the past 30 000 cal. yr is emerging. Three major climate events are recognised: (1) LGCP beginning at ca. 28 000 cal. yr BP, ending at Termination I, ca. 18 000 cal. yr BP, and including a warmer and more variable phase between ca. 27 000 and 21 000 cal. yr BP, (2) LGIT between ca. 18 000 and 11 600 cal. yr BP, including a Lateglacial warm period from ca. 14 800 to 13 500 cal. yr BP and a Lateglacial climate reversal between ca. 13 500 and 11 600 cal. yr BP, and (3) Holocene interglacial conditions, with two phases of greatest warmth between ca. 11 600 and 10 800 cal. yr BP and from ca. 6 800 to 6 500 cal. yr BP. Some key boundaries coincide with volcanic tephras. Copyright © 2007 John Wiley & Sons, Ltd.     

Seston transport and deposition in Pelorus sound,south Island,New Zealand

Transport of seston (suspended sediment) in Pelorus Sound is controlled by tides and freshwater inflow. During high freshwater inflow, a moderately stratified estuarine circulation may be superimposed on the tidal circulation, but the latter dominates and transports seston seawards and landwards with the ebb and flood phases respectively. With extreme freshwater inflow, the estuarine circulation gains impetus and most seston is rapidly transported seaward in the low saline surface layer.

Irrespective of circulation there is a persistent trend in seston concentrations. Highest values occur at the sound's head because of the influence of nearby Pelorus and Kaituna Rivers and because of resuspension of bottom sediment by strong tidal currents. Seston concentrations wane along the sound until near the entrance, where values increase as a result of greater production of biogenic seston and because additional seston is brought in from Cook Strait with the flood tide. This trend parallels variability in the thicknesses of muddy bottom sediments. Muds are thick at the head where an extensive delta extends from the river mouths; muds gradually thin seaward and then thicken markedly in the vicinity of the sound entrance.

Seston weight and composition patterns and 3.5 kHz seismic profiles indicate Pelorus Sound acts as a double‐ended sediment trap. The upper reaches receive and retain river‐derived seston, whereas the sound entrance traps seston derived from Cook Strait. This situation appears to hold for both high and extremely high influxes of sediment.     

Eruptive history of the Colima volcanic complex (Mexico)

The evolution of the Colima volcanic complex can be divided into successive periods characterized by different dynamic and magmatic processes: emission of andesitic to dacitic lava flows, acid-ash and pumice-flow deposits, fallback nuées ardentes leading to pyroclastic flows with heterogeneous magma, plinian air-fall deposits, scoriae cones of alkaline and calc-alkaline nature. Four caldera-forming events, resulting either from major ignimbrite outbursts or Mount St. Helens-type eruptions, separate the main stages of development of the complex from the building of an ancient shield volcano (25 × 30 km wide) up to two summit cones, Nevado and Fuego.The oldest caldera, C1 (7–8 km wide), related to the pouring out of dacitic ash flows, marks the transition between two periods of activity in the primitive edifice called Nevado I: the first one, which is at least 0.6 m.y. old, was mainly andesitic and effusive, whereas the second one was characterized by extrusion of domes and related pyroclastic products. A small summit caldera, C2 (3–3.5 km wide), ended the evolution of Nevado I.Two modern volcanoes then began to grow. The building of the Nevado II started about 200,000 y. ago. It settled into the C2 caldera and partially overflowed it. The other volcano, here called Paleofuego, was progressively built on the southern side of the former Nevado I. Some of its flows are 50,000 y. old, but the age of its first outbursts is not known. However, it is younger than Nevado II. These two modern volcanoes had similar evolutions. Each of them was affected by a huge Mount St. Helens-type (or Bezymianny-type) event, 10,000 y. ago for the Paleofuego, and hardly older for the Nevado II. The landslides were responsible for two horseshoe-shaped avalanche calderas, C3 (Nevado) and C4 (Paleofuego), each 4–5 km wide, opening towards the east and the south. In both cases, the activity following these events was highly explosive and produced thick air-fall deposits around the summit craters.The Nevado III, formed by thick andesitic flows, is located close to the southwestern rim of the C3 caldera. It was a small and short-lived cone. Volcan de Fuego, located at the center of the C4 caldera, is nearly 1500 m high. Its activity is characterized by an alternation of long stages of growth by flows and short destructive episodes related to violent outbursts producing pyroclastic flows with heterogeneous magma and plinian air falls.The evolution of the primitive volcano followed a similar pattern leading to formation of C1 and then C2. The analogy between the evolutions of the two modern volcanoes (Nevado II–III; Paleofuego-Fuego) is described. Their vicinity and their contemporaneous growth pose the problem of the existence of a single reservoir, or two independent magmatic chambers, after the evolution of a common structure represented by the primitive volcano.     

Scaling of the Vertical Spreading of a Plume of a Passive Tracer in a Rough-Wall Neutral Boundary Layer

The influence of surface roughness on the dispersion of a passive scalar in a rough wall turbulent boundary layer has been studied using wind-tunnel experiments. The surface roughness was varied using different sizes of roughness elements, and different spacings between the elements. Vertical profiles of average concentration were measured at different distances downwind of the source, and the vertical spread of the plume was computed by fitting a double Gaussian profile to the data. An estimate of the integral length scale is derived from the turbulence characteristics of the boundary layer and is then used to scale the measured values of plume spread. This scaling reduces the variability in the data, confirming the validity of the model for the Lagrangian integral time scale, but does not remove it entirely. The scaled plume spreading shows significant differences from predictions of theoretical models both in the near and in the far field. In the region immediately downwind of the source this is due to the influence of the wake of the injector for which we have developed a simple model. In the far field we explain that the differences are mainly due to the absence of large-scale motions. Finally, further downwind of the source the scaled values of plume spread fall into two distinct groups. It is suggested that the difference between the two groups may be related to the lack of dynamical similarity between the boundary-layer flows for varying surface roughness or to biased estimates of the plume spread.     

Formation of an eroded lava channel within an Elysium Planitia impact crater: Distinguishing between a mechanical and thermal origin

A lava channel identified on the wall of an Elysium Planitia impact crater is investigated to identify the dominant erosion mechanism, mechanical vs. thermal, acting during channel formation. Observations of channel morphology are used to supplement analytical models of lava channel formation in order to calculate the duration of channel formation, the velocity of the lava flowing through the channel, and the erosion rate in each erosion regime considered. Results demonstrate that the channel observed in the Elysium Planitia impact crater formed primarily due to mechanical erosion. In a more general sense, results of this study suggest that lava channels can form primarily due to thermal erosion in the presence of more gradual slopes and more consolidated substrates whereas lava channels can form primarily due to mechanical erosion in the presence of more energetic flows on steeper slopes and more poorly consolidated substrates. Therefore, both erosion regimes must be considered when analyzing origins of eroded lava channels that cut through strata of different strengths.