Significance of quaternary sea-level fluctuations for delta construction along the Brazilian coast

Many important areas of Quaternary deltaic sedimentation along the Brazilian coast are practically unknown in the geologic literature, especially outside Brazil. Our studies show that these areas, previously considered as highly destructive wave-dominated deltas of Holocene age, were affected by a wave-dominated phase during the Pleistocene, succeeded by highly constructive, intralagoonal stage in the Holocene, in turn followed by wave-dominated deltaic sedimentation that continues until today. The geologic evolution of these coastal plains is exemplified here by the Doce River mouth area, State of Espírito Santo, where relative sea-level fluctuations during the Quaternary played an important role in the construction of the deltaic complex.     

Cosmic ray ionization of the Jovian atmosphere

An approximate form of the Boltzmann equation has been used to obtain local ionization rates due to the absorption of galactic cosmic rays in the Jovian atmosphere. It is shown that the muon flux component of the cosmic ray-induced cascade may be especially importannt in ionizing the atmosphere at levels where the total number density exceeds 10¹⁹ cm^?3 (well below the ionospheric layers produced by solar euv). A model containing both positive and negative ion reactions has been employed to compute equilibrium electron and ion number densities. Peak electron number densities on the order of 10³ cm^?3 may be expected even at relatively low magnetic latitudes. The dominant positive ions are NH₄⁺ and C_nH_m⁺ cluster ions, with n ? 2; it is suggested that the absorption of galactic cosmic ray energy at such relatively high pressures in the Jovian atmosphere $\text{(M ? 10}{}^{18}\text{to}\text{10}{}^{20}\text{cm}{}^{\mathrm{?3}}\text{)}$ and the subsequent chemical reactions may be instrumental in the local formation of complex hydrocarbons.     

Laboratory simulation of infrared astrophysical features

Laboratory infrared emission and absorption spectra have been taken of terrestrial silicates, meteorites and lunar soils in the form of micrometer and sub-micrometer grains. The emission spectra were taken in a way that imitates telescopic observations. The purpose was to see which materials best simulate the 10 m astrophysical feature.The emission spectra of dunite, fayalite and Allende give a good fit to the 10 m broadband emission feature of comets Bennett and Kohoutek. A study of the effect of grain size on the presence of the 10 m emission feature of dunite shows that for particles larger than 37 m no feature is seen.The emission spectrum of the Murray meteorite, a Type 2 carbonaceous chondrite, is quite similar to the intermediate resolution spectrum of comet Kohoutek in the 10 m region. Hydrous silicates or amorphous magnesium silicates in combination with high-temperature condensates, such as olivine or anorthite, would yield spectra that match the intermediate resolution spectrum of comet Kohoutek in the 10 m region. Glassy olivine and glassy anorthite in approximately equal proportions would also give a spectrum that is a good fit to the cometary 10 m feature.Invited contribution to the Proceedings of a Workshop onThermodynamics and Kinetics of Dust Formation in the Space Medium held at the Lunar and Planetary Institute, Houston, 6 8 September, 1978.     

Modelling weathering processes at the catchment scale: The WITCH numerical model

A numerical model of chemical weathering in soil horizons and underlying bedrock (WITCH) has been coupled to a numerical model of water and carbon cycles in forest ecosystems (ASPECTS) to simulate the concentration of major species within the soil horizons and the stream of the Strengbach granitic watershed, located in the Vosges Mountains (France). For the first time, simulations of solute concentrations in soil layers and in the catchment river have been performed on a seasonal basis. The model is able to reproduce the concentrations of most major species within the soil horizons, as well as catching the first-order seasonal fluctuations of aqueous calcium, magnesium and silica concentrations. However, the WITCH model underestimates concentrations of Mg²⁺ and silica at the spring of the catchment stream, and significantly underestimates Ca²⁺ concentration. The deficit in calculated calcium can be compensated for by dissolution of trace apatite disseminated in the bedrock. However, the resulting increased Ca²⁺ release yields important smectite precipitation in the deepest model layer (in contact with the bedrock) and subsequent removal of large amount of silica and magnesium from solution. In contrast, the model accurately accounts for the concentrations of major species (Ca, Mg and silica) measured in the catchment stream when precipitation of clay minerals is not allowed. The model underestimation of Mg²⁺ and H₄SiO₄ concentrations when precipitation of well crystallized smectites is allowed strongly suggests that precipitation of well crystallized clay minerals is overestimated and that more soluble poorly crystallized and amorphous materials may be forming. In agreement with observations on other watersheds draining granitic rocks, this study indicates that highly soluble trace calcic phases control the aqueous calcium budget in the Strengbach watershed.     

Optimization of nonconventional wells under uncertainty using statistical proxies

The determination of the optimal type and placement of a nonconventional well in a heterogeneous reservoir represents a challenging optimization problem. This determination is significantly more complicated if uncertainty in the reservoir geology is included in the optimization. In this study, a genetic algorithm is applied to optimize the deployment of nonconventional wells. Geological uncertainty is accounted for by optimizing over multiple reservoir models (realizations) subject to a prescribed risk attitude. To reduce the excessive computational requirements of the base method, a new statistical proxy (which provides fast estimates of the objective function) based on cluster analysis is introduced into the optimization process. This proxy provides an estimate of the cumulative distribution function (CDF) of the scenario performance, which enables the quantification of proxy uncertainty. Knowledge of the proxy-based performance estimate in conjunction with the proxy CDF enables the systematic selection of the most appropriate scenarios for full simulation. Application of the overall method for the optimization of monobore and dual-lateral well placement demonstrates the performance of the hybrid optimization procedure. Specifically, it is shown that by simulating only 10% or 20% of the scenarios (as determined by application of the proxy), optimization results very close to those achieved by simulating all cases are obtained.     

The Removal of Dissolved Metals by Hydroxysulphate Precipitates during Oxidation and Neutralization of Acid Mine Waters, Iberian Pyrite Belt

This study examines the removal of dissolved metals during the oxidation and neutralization of five acid mine drainage (AMD) waters from La Zarza, Lomero, Esperanza, Corta Atalaya and Poderosa mines (Iberian Pyrite Belt, Huelva, Spain). These waters were selected to cover the spectrum of pH (2.2–3.5) and chemical composition (e.g., 319–2,103 mg/L Fe; 2.85–33.3 g/L SO₄⁼) of the IPB mine waters. The experiments were conducted in the laboratory to simulate the geochemical evolution previously recognized in the field. This evolution includes two stages: (1) oxidation of dissolved Fe(II) followed by hydrolysis and precipitation of Fe(III), and (2) progressive pH increase during mixing with fresh waters. Fe(III) precipitates at pH < 3.5 (stages 1 and 2) in the form of schwertmannite, whereas Al precipitates during stage 2 at pH 5.0 in the form of several hydroxysulphates of variable composition (hydrobasaluminite, basaluminite, aluminite). During these stages, trace elements are totally or partially sorbed and/or coprecipitated at different rates depending basically on pH, as well as on the activity of the SO₄⁼ anion (which determines the speciation of metals). The general trend for the metals which are chiefly present as aqueous free cations (Pb²⁺, Zn²⁺, Cu²⁺, Cd²⁺, Mn²⁺, Co²⁺, Ni²⁺) is a progressive sorption at increasing pH. On the other hand, As and V (mainly present as anionic species) are completely scavenged during the oxidation stage at pH < 3.5. In waters with high activities (> 10⁻¹) of the SO ₄⁼ ion, some elements like Al, Zn, Cd, Pb and U can also form anionic bisulphate complexes and be significantly sorbed at pH < 5. The removal rates at pH 7.0 range from around 100% for As, V, Cu and U, and 60–80% for Pb, to less than 20% for Zn, Co, Ni and Mn. These processes of metal removal represent a significant mechanism of natural attenuation in the IPB.     

GGM02 – An improved Earth gravity field model from GRACE

A new generation of Earth gravity field models called GGM02 are derived using approximately 14 months of data spanning from April 2002 to December 2003 from the Gravity Recovery And Climate Experiment (GRACE). Relative to the preceding generation, GGM01, there have been improvements to the data products, the gravity estimation methods and the background models. Based on the calibrated covariances, GGM02 (both the GRACE-only model GGM02S and the combination model GGM02C) represents an improvement greater than a factor of two over the previous GGM01 models. Error estimates indicate a cumulative error less than 1 cm geoid height to spherical harmonic degree 70, which can be said to have met the GRACE minimum mission goals. Electronic Supplementary Material Supplementary material is available in the online version of this article at     

An Early Cretaceous extinct spreading center in the northern Natal valley

We have identified an extinct E–W spreading center in the northern Natal valley on the basis of magnetic anomalies which was active from chron M11 (133 Ma) to 125.3 Ma, just before chron M2 (124 Ma) in the Early Cretaceous. Seafloor spreading in the northern Natal valley accounts for approximately 170 km of north–south motion between the Mozambique Ridge and Africa. This extension resolves the predicted overlap of the continental (central and southern) Mozambique Ridge and Antarctica in the chron M2 to M11 reconstructions from Mesozoic finite rotation parameters for Africa and Antarctica. In addition, the magnetic data reveal that the Mozambique Ridge was an independent microplate from at least 133 to 125 Ma. The northern Natal valley extinct spreading center connects to the spreading center separating the Mozambique Basin and the Riiser-Larsen Sea to the east. It follows that the northern Mozambique Ridge was either formed after the emplacement of the surrounding oceanic crust or it is the product of a very robust spreading center. To the west the extinct spreading center connects to the spreading center separating the southern Natal valley and Georgia Basin via a transform fault. Prior to chron M11, there is still a problem with the overlap of Mozambique Ridge if it is assumed to be fixed with respect to either the African or Antarctic plates. Some of the overlap can be accounted for by Jurassic deformation of the Mozambique Ridge, Mozambique Basin, and Dronning Maud land. It appears though that the Mozambique Ridge was an independent microplate from the breakup of Gondwana, 160 Ma, until it became part of the African plate, 125 Ma.     

Numerical Calculation of Equivalent Cell Permeability Tensors for General Quadrilateral Control Volumes

A new method for upscaling fine scale permeability fields to general quadrilateral-shaped coarse cells is presented. The procedure, referred to as the conforming scale up method, applies a triangle-based finite element technique, capable of accurately resolving both the coarse cell geometry and the subgrid heterogeneity, to the solution of the local fine scale problem. An appropriate averaging of this solution provides the equivalent permeability tensor for the coarse scale quadrilateral cell. The general level of accuracy of the technique is demonstrated through application to a number of flow problems. The real strength of the conforming scale up method is demonstrated when the method is applied in conjunction with a flow-based gridding technique. In this case, the approach is shown to provide results that are significantly more accurate than those obtained using standard techniques.