Magmatic modification of the uppermost mantle beneath the Basin and Range to Colorado Plateau Transition Zone; Evidence from xenoliths, Wikieup, Arizona

Upper mantle xenoliths from Wikieup, AZ, provide abundant evidence for magmatic modification of the uppermost mantle beneath the Transition Zone between the Colorado Plateau and the southern Basin and Range province. Upper mantle lithologies in this xenolith suite are represented by spinel peridotite, wehrlite, plagioclase peridotite, and Al-augite group pyroxenites. Isotopic data for these xenoliths yield relatively uniform values and suggest a common petrogenesis. Al-augite-bearing gabbro and pyroxenite xenoliths from this locality are interpreted to have formed by crystal fractionation processes from parent alkali basalts similar to the Wikieup host basalt. Mineral and whole rock compositions show consistent trends of increasing incompatible element contents (Fe, Al, Ca, Na, K, LIL, and LREE), and decreasing compatible element contents (Mg, Cr, Ni) from spinel peridotite to wehrlite to plagioclase peridotite to the host basalt composition. These compositional trends are interpreted as resulting from varying degrees of magma-mantle wall rock interaction as ascending mafic magmas infiltrated upper mantle peridotite. Small degrees of melt infiltration resulted in slightly modified spinel peridotite compositions while moderate degrees metasomatized spinel peridotite to wehrlite, and the highest degrees metasomatized it to plagioclase peridotite. Whole rock compositions and clinopyroxene, plagioclase, and whole rock isotopic data suggest that the infiltrating magmas were the same as those from which the gabbros and pyroxenites crystallized, and that they were alkalic in composition, similar to the Wikieup host alkali olivine basalts. Relatively uniform ¹⁴³Nd/¹⁴⁴Nd for the mineral separates and whole rocks in spite of the significantly wide range in their ¹⁴⁷Sm/¹⁴⁴Nd (0.71–0.23 in clinopyroxene) suggests that the Wikieup xenoliths including gabbro, pyroxenite, peridotite, wehrlite, and plagioclase peridotite, are all relatively young rocks formed or metasomatized by a relatively recent magmatic episode. Received: 21 May 1996 / Accepted: 23 December 1996     

Estimation of Sill Matrices in the Linear Model of Coregionalization

When using least squares to fit the linear model of coregionalization to multivariate geostatistical data, the sill matrices for the different regions must be estimated, subject to the constraint that they be non-negative definite. In 1992, Goulard and Voltz proposed and empirically examined an iterative algorithm for doing this. Although no proof was given for its convergence or for the uniqueness of the solution to the problem, the algorithm has subsequently been extensively and successfully used. In this paper, we prove that the minimization problem, in fact, has a unique solution and that the algorithm is guaranteed to converge to it from any starting point. We also discuss the effect of the starting point on the speed of convergence.     

Regional-scale back-analysis using TRIGRS: an approach to advance landslide hazard modeling and prediction in sparse data regions

Landslides in Kerala, India, have been shown to be preceded not only by critical rainfall over a short period but also a much longer period of elevated pore pressure. Such rainfall-triggered landslides are difficult to monitor due to a lack of adequate data on the locations of failures and precipitation. Here, a method is presented using Transient Rainfall Infiltration and Grid-based Regional Slope stability (TRIGRS) as a tool to model the relationship between critical rainfall and antecedent pore pressure as they relate to slope stability, which can be useful for hazard assessment in sparse data regions. This is demonstrated by parameterizing the model with a combination of regional data sources, remote sensing, and temporal back-analysis based on two known failure events (June 2004 and July 2007). Ranges of possible geotechnical and hydraulic parameters were obtained from various local and regional sources, and soil thickness was modeled as a function of slope angle. Rainfall was estimated using satellite microwave radiometry data. For back-analysis, combinations of cohesion, friction angle, and water table depth were then tested in TRIGRS using trial and error until the predicted and observed failure times coincided for the two failure events. While the spatial prediction accuracy of the model is low and multiple solution sets are expected to exist, the results confirm that information regarding the critical pre-failure conditions and stability changes over time can be derived despite data-poor circumstances. Future studies can be undertaken extending this method to characterize many parameter combinations and incorporate more failure cases to develop probabilistic early-warning thresholds.     

Sea level forcing in the Mediterranean Sea between 1960 and 2000

Sea level trends and inter-annual variability in the Mediterranean Sea for the period 1960–2000 is explored by comparing observations from tide gauges with sea level hindcasts from a barotropic 2D circulation model, and two full primitive equation 3D ocean circulation models, a regional one and the Mediterranean component of a global one,. In the 2D model, 50% of the sea level variance was found to result from the wind and atmospheric pressure forcing. In the 3D models, 20% of the sea level variance was explained by the steric effects. The sea level residuals at the tide gauges locations, calculated by subtraction of the 2D model output from the sea level observations are significantly correlated (r = 0.4) with the steric signals from the 3D models. After the removal of the atmospheric and the steric contributions the tide-gauge sea level records indicate a period where sea level was stable (1960–1975) and a period where sea level was rising (1975–2000) with rates in the range 1.1–1.8 mm/yr. A part of the residual trend can be explained by the contribution of local land movements (0.3 mm/yr) while its major part indicates a global signal, probably mass addition, appearing after 1975.     

A Comparison of Hyporheic Transport at a Cross‐Vane Structure and Natural Riffle

While restoring hyporheic flowpaths has been cited as a benefit to stream restoration structures, little documentation exists confirming that constructed restoration structures induce comparable hyporheic exchange to natural stream features. This study compares a stream restoration structure (cross‐vane) to a natural feature (riffle) concurrently in the same stream reach using time‐lapsed electrical resistivity (ER) tomography. Using this hydrogeophysical approach, we were able to quantify hyporheic extent and transport beneath the cross‐vane structure and the riffle. We interpret from the geophysical data that the cross‐vane and the natural riffle induced spatially and temporally unique hyporheic extent and transport, and the cross‐vane created both spatially larger and temporally longer hyporheic flowpaths than the natural riffle. Tracer from the 4.67‐h injection was detected along flowpaths for 4.6 h at the cross‐vane and 4.2 h at the riffle. The spatial extent of the hyporheic zone at the cross‐vane was 12% larger than that at the riffle. We compare ER results of this study to vertical fluxes calculated from temperature profiles and conclude significant differences in the interpretation of hyporheic transport from these different field techniques. Results of this study demonstrate a high degree of heterogeneity in transport metrics at both the cross‐vane and the riffle and differences between the hyporheic flowpath networks at the two different features. Our results suggest that restoration structures may be capable of creating sufficient exchange flux and timescales of transport to achieve the same ecological functions as natural features, but engineering of the physical and biogeochemical environment may be necessary to realize these benefits.     

Application of geochemical and stable isotopic tracers to investigate groundwater salinity in the Ochi-Narkwa Basin,Ghana

Rainwater, groundwater and soil-water samples were analysed to assess groundwater geochemistry and the origin of salinity in the Ochi-Narkwa basin of the Central Region of Ghana. The samples were measured for major ions and stable isotopes (δ¹⁸O, δ²H and δ¹³C). The Cl^? content in rainwater decreased with distance from the coast. The major hydrochemical facies were Na-Cl for the shallow groundwaters and Ca-Mg-HCO₃, Na-Cl and Ca-Mg-Cl-SO₄ for the deep groundwaters. Groundwater salinization is caused largely by halite dissolution and to a minor extent by silicate weathering and seawater intrusion. Stable isotope composition of the groundwaters followed a slope of 3.44, suggesting a mixing line. Chloride profiles in the soil zone revealed the existence of salt crusts, which support halite dissolution in the study area. A conceptual flow model developed to explain the mechanism of salinization showed principal groundwater flow in the NW–SE direction.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR K. Heal     

High-resolution mapping of large gas emitting mud volcanoes on the Egyptian continental margin (Nile Deep Sea Fan) by AUV surveys

Two highly active mud volcanoes located in 990–1,265 m water depths were mapped on the northern Egyptian continental slope during the BIONIL expedition of R/V Meteor in October 2006. High-resolution swath bathymetry and backscatter imagery were acquired with an autonomous underwater vehicle (AUV)-mounted multibeam echosounder, operating at a frequency of 200 kHz. Data allowed for the construction of ~1 m pixel bathymetry and backscatter maps. The newly produced maps provide details of the seabed morphology and texture, and insights into the formation of the two mud volcanoes. They also contain key indicators on the distribution of seepage and its tectonic control. The acquisition of high-resolution seafloor bathymetry and acoustic imagery maps with an AUV-mounted multibeam echosounder fills the gap in spatial scale between conventional multibeam data collected from a surface vessel and in situ video observations made from a manned submersible or a remotely operating vehicle.     

Hydrogen and carbon isotopes of petroleum and related organic matter

$\text{D}\text{H}$ and ${}^{13}\text{C}{}^{12}\text{C}$ ratios were measured for 114 petroleum samples and for several samples of related organic matter. δD of crude oil ranges from ?85 to ?181‰, except for one distillate (?250‰) from the Kenai gas field; δ¹³C of crude oil ranges from ?23.3 to ?32.5‰, Variation in δD and δ¹³C values of compound-grouped fractions of a crude oil is small, 3 and 1.1%., respectively, and the difference in δD and δ¹³C between oil and coeval wax is slight. Gas fractions are 53–70 and 22.6–23.2‰ depleted in D and ¹³C, respectively, relative to the coexisting oil fractions.The δD and δ¹³C values of the crude oils appear to be largely determined by the isotopic compositions of their organic precursors. The contribution of terrestrial organic debris to the organic precursors of most marine crude oils may be significant.     

Diffusion of a multi-component plasma

Diffusion in a weakly ionized plasma composed of negative ions as well as positive ions is examined using appropriate linearized fluid equations. For initial, electrically neutral, density perturbations of the form exp(ikx) the diffusion process is characterized by electrical non-neutrality and by three stages or time scales. For equal positive and negative ion diffusion coefficients these stages are in general (1) equilibration of the electron gas so that pressure gradient and electric forces are balanced (2) ambipolar-like diffusion of all three species, and concluding with (3) free ion diffusion. The details of the process are governed by the product kλ_e (wave number times electron Debye length) and the ambient ratio of negative ion to electron number density. Numerical and analytic results for separate positive and negative diffusion coefficients show added complexity which is briefly described. These results or the more complete numerical solutions find application to the lower D region of the ionosphere.