15,000-yr pollen record of vegetation change in the high altitude tropical Andes at Laguna Verde Alta, Venezuela

Pollen analysis of sediments from a high-altitude (4215 m), Neotropical (9°N) Andean lake was conducted in order to reconstruct local and regional vegetation dynamics since deglaciation. Although deglaciation commenced 15,500 cal yr B.P., the area around the Laguna Verde Alta (LVA) remained a periglacial desert, practically unvegetated, until about 11,000 cal yr B.P. At this time, a lycopod assemblage bearing no modern analog colonized the superpáramo. Although this community persisted until 6000 cal yr B.P., it began to decline somewhat earlier, in synchrony with cooling following the Holocene thermal maximum of the Northern Hemisphere. At this time, the pioneer assemblage was replaced by a low-diversity superpáramo community that became established 9000 cal yr B.P. This replacement coincides with regional declines in temperature and/or available moisture. Modern, more diverse superpáramo assemblages were not established until 4600 cal yr B.P., and were accompanied by a dramatic decline in Alnus, probably the result of factors associated with climate, humans, or both. Pollen influx from upper Andean forests is remarkably higher than expected during the Late Glacial and early to middle Holocene, especially between 14,000 and 12,600 cal yr B.P., when unparalleled high values are recorded. We propose that intensification of upslope orographic winds transported lower elevation forest pollen to the superpáramo, causing the apparent increase in tree pollen at high altitude. The association between increased forest pollen and summer insolation at this time suggests a causal link; however, further work is needed to clarify this relationship.     

Long-term dynamical evolution of dusty ejecta from Deimos

We re-assess expected properties of the presumed dust belt of Mars formed by impact ejecta from Deimos. Previous studies have shown that dynamics of Deimos particles are dominated by two perturbing forces: radiation pressure (RP) and Mars’ oblateness (J2). At the same time, they have demonstrated that lifetimes of particles, especially of grains about ten of micrometers in size, may reach more than 10⁴ years. On such timescales, the Poynting-Robertson drag (PR) becomes important. Here we provide a study of the dynamics under the combined action of all three perturbing forces. We show that a PR decay of the semimajor axes leads to an adiabatic decrease of amplitudes and periods of oscillations in orbital inclinations predicted in the framework of the underlying RP+J2 problem. Furthermore, we show that smallest of the long-lived Deimos grains (radius≈5- may reach a chaotic regime, resulting in unpredictable and abrupt changes of their dynamics. The particles just above that size (≈10-) should be the most abundant in the Deimos torus. Our dynamical analysis, combined with a more accurate study of the particle lifetimes, provides corrections to earlier predictions about the dimensions and geometry of the Deimos torus. In addition to a population, appreciably inclined and shifted towards the Sun, the torus should contain a more contracted, less asymmetric, and less tilted component between the orbits of Phobos and Deimos.     

Optimal and sustainable extraction of groundwater in coastal aquifers

Four examples are investigated for the optimal and sustainable extraction of groundwater from a coastal aquifer under the threat of seawater intrusion. The objectives and constraints of these management scenarios include maximizing the total volume of water pumped, maximizing the profit of selling water, minimizing the operational and water treatment costs, minimizing the salt concentration of the pumped water, and controlling the drawdown limits. The physical model is based on the density-dependent advective-dispersive solute transport model. Genetic algorithm is used as the optimization tool. The models are tested on a hypothetical confined aquifer with four pumping wells located at various depths. These solutions establish the feasibility of simulating various management scenarios under complex three-dimensional flow and transport processes in coastal aquifers for the optimal and sustainable use of groundwater.     

Magnetic Twist and Writhe of δ Active Regions

We have selected 104 active regions with a δ magnetic configuration from 1996 to 2002 to study how important a role the kink instability plays in such active regions. In this study, we employ the systematic tilt angle of each active region as a proxy for the writhe of a flux tube and the force-free parameter, α_best, as a proxy for the magnetic field twist in the flux tube. It is found that 65–67% of the active regions have the same sign of twist and writhe. About 34% (32%) of the active regions violate (follow) the Hale-Nicholson and Joy's Laws (HNJL) but follow (violate) the hemispheric helicity rule (HHR). Sixty-one (61) of the 104 active regions studied each produced more than five large flares. Active regions violating HNJL, but following HHR, have a much stronger tendency to produce X-class flares and/or strong proton events. Comparing with previous studies for active regions with well-defined (simpler) bipolar magnetic configuration, it is found that the numbers following both HNJL and HHR are significantly lower in the δ-configuration case, while numbers violating one of the laws and the rule significantly increase with the increase of the magnetic complexity of the active regions. These results support the prediction for the presence of a kink instability, that the twist and writhe of the magnetic fields exhibit the same sign for δ active regions (Linton et al., Astrophys. J. 507, 40, 1998, Astrophys. J. 522, 1205, 1999; Fan et al., Astrophys. J. 521, 460, 1999). Finally, we analyze possible origins of the twist and writhe of the magnetic fields for the active regions studied.     

Acceleration of energetic particles of the outer regions of planetary magnetospheres: Inferences from laboratory and space experiments

High energy particles, with energies above those attainable by adiabatic or steady-state electric field acceleration, have been observed in and around the outer regions of planetary magnetospheres. Acceleration by large amplitude sporadic cross-tail electric fields over an order of magnitude greater than steady-state convection fields is proposed as a source of these particles. It is suggested that such explosive electric fields will occur intermittently in the vicinity of the tail neutral line in the expansive phases of substorms. We use laboratory Double Inverse Pinch Device (DIPD) and satellite evidence to estimate this electric potential for substorms at Earth; values of 500 kV to 2 MV are calculated, in agreement with particle observations. It is further suggested that these particles, which have been accelerated in the night side magnetosphere, drift to the dayside on closed field lines, and under certain interplanetary conditions can escape to regions upstream of the bow shock.     

Amoeboid olivine aggregates and related objects in carbonaceous chondrites: records of nebular and asteroid processes

Amoeboid olivine aggregates (AOAs) are the most common type of refractory inclusions in CM, CR, CH, CV, CO, and ungrouped carbonaceous chondrites Acfer 094 and Adelaide; only one AOA was found in the CB_b chondrite Hammadah al Hamra 237 and none were observed in the CB_a chondrites Bencubbin, Gujba, and Weatherford. In primitive (unaltered and unmetamorphosed) carbonaceous chondrites, AOAs consist of forsterite (Fa_<2), Fe, Ni-metal (5-12 wt% Ni), and Ca, Al-rich inclusions (CAIs) composed of Al-diopside, spinel, anorthite, and very rare melilite. Melilite is typically replaced by a fine-grained mixture of spinel, Al-diopside, and ±anorthite; spinel is replaced by anorthite. About 10% of AOAs contain low-Ca pyroxene replacing forsterite. Forsterite and spinel are always ¹⁶O-rich (δ^17,18O∼−40‰ to −50‰), whereas melilite, anorthite, and diopside could be either similarly ¹⁶O-rich or ¹⁶O-depleted to varying degrees; the latter is common in AOAs from altered and metamorphosed carbonaceous chondrites such as some CVs and COs. Low-Ca pyroxene is either ¹⁶O-rich (δ^17,18O∼−40‰) or ¹⁶O-poor (δ^17,18O∼0‰). Most AOAs in CV chondrites have unfractionated (∼2-10×CI) rare-earth element patterns. AOAs have similar textures, mineralogy and oxygen isotopic compositions to those of forsterite-rich accretionary rims surrounding different types of CAIs (compact and fluffy Type A, Type B, and fine-grained, spinel-rich) in CV and CR chondrites. AOAs in primitive carbonaceous chondrites show no evidence for alteration and thermal metamorphism. Secondary minerals in AOAs from CR, CM, and CO, and CV chondrites are similar to those in chondrules, CAIs, and matrices of their host meteorites and include phyllosilicates, magnetite, carbonates, nepheline, sodalite, grossular, wollastonite, hedenbergite, andradite, and ferrous olivine.Our observations and a thermodynamic analysis suggest that AOAs and forsterite-rich accretionary rims formed in ¹⁶O-rich gaseous reservoirs, probably in the CAI-forming region(s), as aggregates of solar nebular condensates originally composed of forsterite, Fe, Ni-metal, and CAIs. Some of the CAIs were melted prior to aggregation into AOAs and experienced formation of Wark-Lovering rims. Before and possibly after the aggregation, melilite and spinel in CAIs reacted with SiO and Mg of the solar nebula gas enriched in ¹⁶O to form Al-diopside and anorthite. Forsterite in some AOAs reacted with ¹⁶O-enriched SiO gas to form low-Ca pyroxene. Some other AOAs were either reheated in ¹⁶O-poor gaseous reservoirs or coated by ¹⁶O-depleted pyroxene-rich dust and melted to varying degrees, possibly during chondrule formation. The most extensively melted AOAs experienced oxygen isotope exchange with ¹⁶O-poor nebular gas and may have been transformed into magnesian (Type I) chondrules. Secondary mineralization and at least some of the oxygen isotope exchange in AOAs from altered and metamorphosed chondrites must have resulted from alteration in the presence of aqueous solutions after aggregation and lithification of the chondrite parent asteroids.     

Noril’sk–Talnakh Cu–Ni–PGE deposits: a revised tectonic model

The Norilsk mining district is located at the northwest margin of the Tunguska basin, in the centre of the 3,000×4,000 km Siberian continental flood basalt (CFB) province. This CFB province was formed at the Permo-Triassic boundary from a superplume that ascended into the geometric centre of the Laurasian continent, which was surrounded by subducting slabs of oceanic crust. We suggest that these slabs could have reached the core–mantle boundary, and they may have controlled the geometric focus of the superplume. The resulting voluminous magma intruded and erupted in continental rifts and related extensive flood basalt events over a 2–4 Ma period. Cu–Ni–PGE sulfide mineralization is found in olivine-bearing differentiated mafic intrusions beneath the flood basalts at the northwestern margin of the Siberian craton and also in the Taimyr Peninsula, some 300 km east of a triple junction of continental rifts, now buried beneath the Mesozoic–Cenozoic sedimentary basin of western Siberia. The Norilsk-I and Talnakh-Oktyabrsky deposits occur in the Norilsk–Kharaelakh trough of the Tunguska CFB basin. The Cu–Ni–PGE-bearing mineralized intrusions are 2–3 km-wide and 20 km-long differentiated chonoliths. Previous studies suggested that parts of the magma remained in intermediate-level crustal chambers where sulfide saturation and accumulation took place before emplacement. The 5–7-km-thick Neoproterozoic to Palaeozoic country rocks, containing sedimentary Cu mineralization and evaporites, may have contributed additional metal and sulfur to this magma. Classic tectonomagmatic models for these deposits proposed that subvertical crustal faults, such as the northeast-trending Norilsk–Kharaelakh fault, were major trough-parallel conduits providing access for magmas to the final chambers. However, geological maps of the Norilsk region show that the Norilsk–Kharaelakh fault offsets the mineralization, which was deformed into folds and offset by related reverse faults, indicating compressional deformation after mineralization in the Late Triassic to Early Jurassic. In addition, most of the intrusions are sills, not dykes as should be expected if the vertical faults were major conduits. A revised tectonic model for the Norilsk region takes into account the fold structure and sill morphology of the dominant intrusions, indicating a lateral rather than vertical emplacement direction for the magma into final chambers. Taking into account the fold structure of the country rocks, the present distribution of the differentiated intrusions hosting the Norilsk-I and Talnakh–Oktyabrsky deposits may represent the remnants of a single, >60 km long, deformed and eroded palm-shaped cluster of mineralized intrusions, which are perceived as separate intrusions at the present erosional level. The original direction of sill emplacement may have been controlled by a northeast-trending paleo-rise, which we suggest is present at the southeastern border of the Norilsk–Kharaelakh trough based on analysis of the unconformity at the base of the CFB. The mineralized intrusions extend along this rise, which we interpret as a structure that formed above the extensionally tilted block in the metamorphic basement. Geophysical data indicate the presence of an intermediate magma chamber that could be linked with the Talnakh intrusion. In turn, this T-shaped flat chamber may link with the Yenisei–Khatanga rift along the northwest-trending Pyasina transform fault, which may have served as the principal magma conduit to the intermediate chamber. It then produced the differentiated mineralized intrusions that melted through the evaporites with in situ precipitation of massive, disseminated, and copper sulfide ore. The Norilsk–Kharaelakh crustal fault may not relate to mineralization and possibly formed in response to late Mesozoic spreading in the Arctic Ocean.Editorial handling: P. Lightfoot