Catastrophic landslides,glacier behaviour and moraine formation – A view from an active plate margin

The influence of large bedrock landslides (“rock avalanches”) on the behaviour of glaciers is incompletely recognised. Here we present an example from an active tectonic margin in South Island, New Zealand where large earthquakes leave a significant imprint on glacial records. We demonstrate that terminal moraines on the western side of the Southern Alps record both ‘ordinary’ (i.e. climate-driven) and landslide-initiated glacial advances. Following consideration of the processes involved in rock avalanche-initiated moraine construction we suggest ways of determining the nature of the advance that built the terminal moraine. The implications of these observations are important in breaking the conventional linkage of individual terminal moraines with climate forcing.     

Establishing Quaternary as a formal internationa Period/System

Despite being the most widely used unit infield mapping and having the greatest number of active researchers, the interval known as Quaternary is unique among the chronostratigraphic subdivisions of the Geological Time Scale （GTS） in having the most controversial definition and rank. After more than 100 years of debate, the base of the Quaternary is now widely recognized at -2.6 Ma, marking a dramatic and so-far irreversible shift to the ice-age-dominated world of oscillating glacial advances over the northern continents.     

Metal-residence sites in lavas and tuffs from Volcán Popocatépetl, Mexico: implications for metal mobility in the environment

 Volcan Popocatépetl is a Quaternary stratovolcano located 60 km southeast of Mexico City. The summit crater is the site of recent ash eruptions, excess degassing, and dacite dome growth. The modern cone comprises mainly pyroclastic flow deposits, airfall tephras, debris flows, and reworked deposits of andesitic composition; it is flanked by more mafic monogenetic vents. In least-degassed fallout tuffs and mafic scoria, transition metals are concentrated in phases formed before eruption, during eruption, and after eruption. Preeruptive minerals occur in both lavas and tephra, and include oxides and sulfides in glass and phenocrysts. The magmatic oxides consist of magnetite, ilmenite, and chromite; the sulfides consist of both (Fe,Ni)_1-xS (MSS) and Cu–Fe sulfide (ISS). Syn- and posteruptive phases occur in vesicles in both lavas and tephra, and on surfaces of ash and along fractures. The mineral assemblages in lavas include Cu–Fe sulfide and Fe–Ti oxide in vesicles, and Fe sulfide and Cu–Fe sulfide in segregation vesicles. Assemblages in vesicles in scoria include Fe–Ti oxide and rare Fe–Cu–Sn sulfide. Vesicle fillings of Fe–Ti oxide, Ni-rich chromite, Fe sulfide, Cu sulfide, and barite are common to two pumice samples. The most coarse-grained of the vesicle fillings are Cu–Fe sulfide and Cu sulfide, which are as large as 50 μ in diameter. The youngest Plinian pumice also contains Zn(Fe) sulfide, as well as rare Ag–Cu sulfide, Ag–Fe sulfide, Ag bromide, Ag chloride, and Au–Cu telluride. The assemblage is similar to those typically observed in high-sulfidation epithermal mineralization. The fine-grained nature and abundance of syn- and/or posteruptive phases in porous rocks makes metals susceptible to mobilization by percolating fluids. The abundance of metal compounds in vesicles indicates that volatile exsolution prior to and/or during eruption played an important role in releasing metals to the atmosphere. Received: March 1997 · Accepted: 27 May 1997     

Damped Oscillations of Coronal Loops

A mechanism of damped oscillations of a coronal loop is investigated. The loop is treated as a thin toroidal flux rope with two stationary photospheric footpoints, carrying both toroidal and poloidal currents. The forces and the flux-rope dynamics are described within the framework of ideal magnetohydrodynamics (MHD). The main features of the theory are the following: i) Oscillatory motions are determined by the Lorentz force that acts on curved current-carrying plasma structures and ii) damping is caused by drag that provides the momentum coupling between the flux rope and the ambient coronal plasma. The oscillation is restricted to the vertical plane of the flux rope. The initial equilibrium flux rope is set into oscillation by a pulse of upflow of the ambient plasma. The theory is applied to two events of oscillating loops observed by the Transition Region and Coronal Explorer (TRACE). It is shown that the Lorentz force and drag with a reasonable value of the coupling coefficient (c _d ) and without anomalous dissipation are able to accurately account for the observed damped oscillations. The analysis shows that the variations in the observed intensity can be explained by the minor radial expansion and contraction. For the two events, the values of the drag coefficient consistent with the observed damping times are in the range c _d ≈2 – 5, with specific values being dependent on parameters such as the loop density, ambient magnetic field, and the loop geometry. This range is consistent with a previous MHD simulation study and with values used to reproduce the observed trajectories of coronal mass ejections (CMEs).     

Correspondence analysis of water quality data: Implications for fauna deaths at Stillwater Lakes, Nevada

Over the last 2 years, more than seven million fish and 16,000 birds have died from drought-induced disease and suspected poisoning at Stillwater Lakes, Nevada. Correspondence analysis of water samples taken in June and September 1986, from lakes and reservoirs near Stillwater Lakes and four background lakes/reservoirs (not associated with fauna disease and death), is attempted to investigate the correlation among sites associated with fauna deaths. This analysis identifies high correlation between some Stillwater Lakes sites and boron. A correlation between boron and arsenic is also identified. Leaching and transport of toxins to the Stillwater Lakes system is one hypothesis to explain the wildlife deaths. Other hypotheses include the direct introduction of toxins by groundwater associated with geothermal activity. A natural cause for the wildlife deaths is also suggested, a cause not associated with the introduction of toxins.     

Concentration–discharge relationships reflect chemostatic characteristics of US catchments

Concentration–discharge relationships have been widely used as clues to the hydrochemical processes that control runoff chemistry. Here we examine concentration–discharge relationships for solutes produced primarily by mineral weathering in 59 geochemically diverse US catchments. We show that these catchments exhibit nearly chemostatic behaviour; their stream concentrations of weathering products such as Ca, Mg, Na, and Si typically vary by factors of only 3 to 20 while discharge varies by several orders of magnitude. Similar patterns are observed at the inter‐annual time scale. This behaviour implies that solute concentrations in stream water are not determined by simple dilution of a fixed solute flux by a variable flux of water, and that rates of solute production and/or mobilization must be nearly proportional to water fluxes, both on storm and inter‐annual timescales. We compared these catchments' concentration–discharge relationships to the predictions of several simple hydrological and geochemical models. Most of these models can be forced to approximately fit the observed concentration–discharge relationships, but often only by assuming unrealistic or internally inconsistent parameter values. We propose a new model that also fits the data and may be more robust. We suggest possible tests of the new model for future studies. The relative stability of concentration under widely varying discharge may help make aquatic environments habitable. It also implies that fluxes of weathering solutes in streams, and thus fluxes of alkalinity to the oceans, are determined primarily by water fluxes. Thus, hydrology may be a major driver of the ocean‐alkalinity feedback regulating climate change. Copyright © 2009 John Wiley & Sons, Ltd.     

Biogeochemical impact of a model western iron source in the Pacific Equatorial Undercurrent

Trace element distributions in the source waters of the Pacific Equatorial Undercurrent (EUC) show the existence of elevated total acid-soluble iron concentrations. This region has been suggested to contribute enough bioavailable iron to regulate interannual and interglacial variability in biological productivity downstream in the high-nitrate low-chlorophyll upwelling zone of the eastern equatorial Pacific. We investigated the advection and first-order biogeochemical impact of an imposed, data-based iron maximum in the western Pacific EUC using an ecosystem model forced by a global dynamical model. We imposed two source profiles of iron constrained by total acid-soluble iron measurements. Though the data for total acid-soluble iron included both dissolved and acid-soluble particulate iron species, we treated all of the total acid-soluble iron as if it was dissolved and bioavailable. A deeper (270 m) source was centered in the density horizon of the observed iron maximum and a shallower (180 m) source was located in the core of our model's EUC, where a dissolved iron maximum has been frequently postulated. These source runs were compared with a control run that contained no specific source of iron associated with the EUC. In the source runs elevated iron concentrations were simulated in the EUC across its entire zonal span, evident as a subsurface plume of dissolved iron slightly below the core of the EUC. In the control run there was no iron maximum associated with the EUC. Upwelling of iron-replete water in the central and eastern equatorial Pacific increased integrated primary productivity in the Wyrtki box (180°W:90°W, 5°S:5°N, 0:200 m) by 41% and 66% for the deeper and shallower iron perturbation, respectively. The source runs increased the realism of the zonal extent of HNLC conditions and the meridional distributions of biological productivity, relative to the control run. However, in the source simulations surface chlorophyll concentrations were too high by a factor of two and maximum surface nitrate concentrations were too low, relative to climatologies. The relative abundance of diatoms roughly doubled upon the input of additional iron, exceeding field observations. Though biogeochemical data are limited and we did not adjust parameters to optimize the model fits to observations, these results suggest that acid-soluble particulate iron supplied to the EUC in the western equatorial Pacific is unlikely to be entirely bioavailable.     

Simulation of the Magnetothermal Instability

In many magnetized, dilute astrophysical plasmas, thermal conduction occurs almost exclusively parallel to magnetic field lines. In this case, the usual stability criterion for convective stability, the Schwarzschild criterion, which depends on entropy gradients, is modified. In the magnetized long mean free path regime, instability occurs for small wavenumbers when (∂ P/∂z) (∂ ln T/∂ z) > 0, which we refer to as the Balbus criterion. We refer to the convective-type instability that results as the magnetothermal instability (MTI). We use the equations of MHD with anisotropic electron heat conduction to numerically simulate the linear growth and nonlinear saturation of the MTI in plane-parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured from the simulations are in excellent agreement with the weak field dispersion relation. The addition of isotropic conduction, e.g. radiation, or strong magnetic fields can damp the growth of the MTI and affect the nonlinear regime. The instability saturates when the atmosphere becomes isothermal as the source of free energy is exhausted. By maintaining a fixed temperature difference between the top and bottom boundaries of the simulation domain, sustained convective turbulence can be driven. MTI-stable layers introduced by isotropic conduction are used to prevent the formation of unresolved, thermal boundary layers. We find that the largest component of the time-averaged heat flux is due to advective motions as opposed to the actual thermal conduction itself. Finally, we explore the implications of this instability for a variety of astrophysical systems, such as neutron stars, the hot intracluster medium of galaxy clusters, and the structure of radiatively inefficient accretion flows. J. M. Stone: Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ 08544