Paleomagnetic estimate of the emplacement temperature of the long-runout Nevado de Colima volcanic debris avalanche deposit, Mexico

Stoopes and Sheridan have mapped a volcanic debris avalanche of Nevado de Colima which has an exceptionally long runout (120 km) and low fall-height to length ratio (H/L = 0.04). We present paleomagnetic results from this volcanic debris avalanche deposit which provide evidence that this avalanche was emplaced at elevated temperatures. The majority of samples, collected from lithic clasts in the volcanic debris avalanche deposit, exhibit two-component remanent magnetizations with a low-temperature component (25–350°C) which is well grouped about the geomagnetic field direction at Colima and a high-temperature component (350–580°C) which is randomly oriented. Although the temperature of the deposit most likely varied with distance from the volcanic source and the thickness of the deposit, our results suggest an emplacement temperature of approximately 350°C at intermediate distances (18–26 km) from the source. In order for the rock clasts (20–40 cm diameter) to be heated to these temperatures, the avalanche was most likely the results of a magmatic, Bezymianny-type eruption. The mixing of hot, juvenile gases with the clasts provides an explanation for the high degree of fluidization of this material, as evidenced by the long runout of this avalanche deposit.     

Predation by jellyfish on large and emergent zooplankton: Implications for benthic–pelagic coupling

Stable carbon isotopes were used to determine the contribution of emergent demersal zooplankton to the diet of the scyphozoan jellyfish Catostylus mosaicus at Smiths Lake, New South Wales, Australia. A preliminary study in 2004 indicated that there was no difference in the δ¹³C of ectodermal tissue and mesoglea of the medusae. In 2005, medusae and zooplankton present during the day and night were sampled and isotopic signatures were modelled using IsoSource. Modelling indicated that: (1) mollusc veligers and copepods sampled during the day contributed <13% of the carbon to the jellyfish; (2) copepods sampled at night contributed up to 25%; and (3) the large, emergent decapod Lucifer sp. contributed 88–94%. We hypothesised, therefore, that medusae derive most of their carbon from emergent species of zooplankton. In 2006, sampling done in 2005 was repeated three times over a period of 4 weeks to measure short-term temporal variation in isotopic signatures of medusae and zooplankton, and emergent demersal zooplankton was specifically sampled using emergence traps. Short-term temporal variation in isotopic signatures was observed for some taxa, however, actual variations were small (<1.5‰) and the values of medusae and zooplankton remained consistent relative to each other. IsoSource modelling revealed that mysid shrimp and emergent copepods together contributed 79–100% of the carbon to the jellyfish, and that the maximum possible contribution of daytime copepods and molluscs was only 22%. Jellyfish apparently derive most of their carbon from emergent zooplankton and by capturing small numbers of relatively large taxa, such as Lucifer sp. or mysid shrimp. Small but abundantly captured zooplankton (such as mollusc veligers) contribute only minor amounts of carbon. Jellyfish have a major role in the transfer of carbon between benthic and pelagic food webs in coastal systems.     

A record of the upper Cochiti polarity transition from the Southern Hemisphere (FIJI)

We sampled the upper Cochiti polarity transition recorded in the Suva Marl in Fiji (18° S, 178° E). The Suva Marl accumulated at an average rate of 83.2 m/m.y. and provides a unique opportunity to address the question of whether present-day overprints bias the transition records. Fiji has undergone a 30° counterclockwise rotation since the deposition of the Suva Marl and hence a present day overprint onto the transitions would bias the VGPs toward the east of the site. Replicate transition records yield VGPs over both Asia and the Americas. Progressive demagnetization reveals a normal polarity overprint which was acquired prior to the rotation of Fiji, and therefore does not produce an easterly bias to the transition data even after correction for the rotation. We collected an oriented block sample across a portion of the transitional interval in the upper Cochiti reversal. Subsampling of this block into 1.0 cm thick wafers cut parallel to bedding provides considerably greater detail during the transition. These detailed data suggest that the oscillatory movement of the VGPs in these sections may be a result of the averaging caused by standard, detailed minicoring and not by the sediments or the remanence acquisition process.     

Precambrian plate tectonics: The midcontinent gravity high

Examination of the shape of the midcontinent gravity high of central North America has led to the hypothesis that the Keweenawan rift system that caused it is the result of plate tectonic interaction. A numerical test has been carried out on the width and postulated transform fault offsets of the gravity high. The exactness of fit to a plate tectonic geometry implies that the continental lithosphere behaved as rigid plates during the Late Precambrian, about 1.1 by ago. This exactness of fit also suggests that the total amount of separation on the Keweenawan rifts is equal to the width of the gravity high. Gravity modelling studies bear out the plausibility of a major amount of rifting, up to 90 km under central Lake Superior. The midcontinent gravity high may represent an intermediate stage of continental rifting, since similar gravity highs and strong associated magnetic anomalies are found on the modern rifted margins of the Atlantic Ocean.     

The effects of wildfires on the magnetic properties of soils in the Everglades

We present results of a rock‐magnetic study of soils that were affected by wildfires that burned portions of the Everglades in the Spring of 2008. Soils at sites that were extensively burned exhibit a pronounced surface magnetic enhancement effect with magnetizations of surface samples up to 16 times greater than that observed at depth (>7 cm) at these sites. The increase in magnetization results from an increased abundance of a low‐coercivity phase (maghemite) that occurs at the expense of the abundance of a high‐coercivity phase (goethite). These results indicate that fire‐induced heating caused goethite in the surface soils to convert into a more magnetic, low‐coercivity phase, such as maghemite. Goethite is an excellent adsorber of phosphorus, and therefore we hypothesize that the destruction of goethite as a result of burning may have important implications for phosphorus cycling in the Everglades ecosystem. Copyright © 2010 John Wiley & Sons, Ltd.     

An interpolation method taking into account inequality constraints: I. Methodology

Several kinds of data can provide information about a variable measured on a one- or two-dimensional space; at some points, the value is known to be equal to a certain number. At other points, the only information may be that the variable is greater or smaller than a given value. The theory of splines provides interpolating functions that can take into account both equality and inequality data. These interpolating functions are presented. The parallel between splines and kriging is reviewed, using the formalism of dual kriging. Coefficients of dual kriging can be obtained directly by minimizing a quadratic form. By adding some inequality constraints to this minimization, an interpolating function may be calculated which takes into account inequality data and is more general than a spline. The method is illustrated by some simple one-dimensional examples.Work performed at Sohio Petroleum Company     

Water resource impacts of climate change in southwestern Bulgaria

This study assesses the regional impact of climate change on runoff in a mountainous region of southwestern Bulgaria. A GIS-based distributed hydrologic model and two climate change scenarios – HadCM2 and CCC – were employed for years around 2025 and 2085. Results from both scenarios demonstrate the basin's sensitivity of runoff to climate change, which produce significant spatial and temporal changes in the basin's water yield with maximum runoff shift into early spring and further decreases in summer runoff. There could be no reduction in mean annual runoff, except under the Hadley scenario by 2085. Changes in the magnitude of mean monthly and peak flow are associated with early snowmelt and a reduction in snow cover in spring. Increases in spatial variability of runoff reflect the basin's complex physiographic characteristics. The increase of spatial and temporal variability in runoff points to different strategies for future water resource management. This revised version was published online in July 2006 with corrections to the Cover Date.     

Empirical likelihood estimation of the spatial quantile regression

The spatial quantile regression model is a useful and flexible model for analysis of empirical problems with spatial dimension. This paper introduces an alternative estimator for this model. The properties of the proposed estimator are discussed in a comparative perspective with regard to the other available estimators. Simulation evidence on the small sample properties of the proposed estimator is provided. The proposed estimator is feasible and preferable when the model contains multiple spatial weighting matrices. Furthermore, a version of the proposed estimator based on the exponentially tilted empirical likelihood could be beneficial if model misspecification is suspect.     

Increasing illumination and sensitivity of reverse‐time migration with internal multiples

Reverse‐time migration is a two‐way time‐domain finite‐frequency technique that accurately handles the propagation of complex scattered waves and produces a band‐limited representation of the subsurface structure that is conventionally assumed to be linear in the contrasts in model parameters. Because of this underlying linear single‐scattering assumption, most implementations of this method do not satisfy the energy conservation principle and do not optimally use illumination and model sensitivity of multiply scattered waves. Migrating multiply scattered waves requires preserving the non‐linear relation between the image and perturbation of model parameters. I modify the extrapolation of source and receiver wavefields to more accurately handle multiply scattered waves. I extend the concept of the imaging condition in order to map into the subsurface structurally coherent seismic events that correspond to the interaction of both singly and multiply scattered waves. This results in an imaging process referred to here as non‐linear reverse‐time migration. It includes a strategy that analyses separated contributions of singly and multiply scattered waves to a final non‐linear image. The goal is to provide a tool suitable for seismic interpretation and potentially migration velocity analysis that benefits from increased illumination and sensitivity from multiply scattered seismic waves. It is noteworthy that this method can migrate internal multiples, a clear advantage for imaging challenging complex subsurface features, e.g., in salt and basalt environments. The results of synthetic seismic imaging experiments, including a subsalt imaging example, illustrate the technique.