Nd isotopic composition of Paleoproterozoic volcanic and granitoid rocks of Vila Riozinho: implications for the crustal evolution of the Tapajós gold province, Amazon craton

Nd data from the Paleoproterozoic magmatic rocks of Vila Riozinho and Jamanxim (Tapajós gold province) indicate that original magmas were not produced exclusively by the remelting of Archean sialic crust and point to dominant Paleoproterozoic sources. ε_Nd(T) values preclude derivation from mantle sources for the ca. 2.0 Ga Vila Riozinho volcanics and older São Jorge granite. They may represent a subduction-related magmatic arc with magmas modified by interaction with crust or a post- to late-orogenic remelting of an older Paleoproterozoic juvenile arc with minimal contribution from the Archean crust. The origin of the 1.88 Ga Parauari, Maloquinha, Iriri, and Moraes Almeida igneous associations and the Jamanxim rhyolites has been attributed to large-scale taphrogenesis that marked the breakup of a large Paleoproterozoic continent. Derivation of the original magmas from the remelting of crustal sources older than ca. 1.9 Ga is consistent with geochemical and Nd isotopic data. Archean remnants probably occur between the Paleoproterozoic terrains of the Ventuari-Tapajós province. Archean terrains of the Amazon craton extend from the Xingu to the Itaituba region but have not been identified in the southern Guyana shield. Thus, data reveal that the boundaries between the central Amazon and Ventuari-Tapajós provinces need better definition and more detailed field and geochronological work.     

Rainfall interception by an isolated evergreen oak tree in a Mediterranean savannah

Redistribution of ground‐level rainfall and interception loss by an isolated Quercus ilex tree were measured over 2 years in a Mediterranean oak savannah. Stemflow, meteorological variables and sap flow were also monitored. Rainfall at ground level was measured by a set of rain‐gauges located in a radial layout centred on the tree trunk and extending beyond the crown limits. Interception loss was computed as the difference between the volume of rainwater that would reach the ground in the absence of the tree and the volume of water that actually fell on the ground sampling area (stemflow included). This procedure provided correct interception loss estimates, irrespective of rainfall inclination. Results have shown a clear non‐random spatial distribution of ground‐level rainfall, with rainwater concentrations upwind beneath the crown and rain‐shadows downwind. Interception loss amounted to 22% of gross rainfall, per unit of crown‐projected area. Stand interception loss, per unit of ground area, was only 8% of gross rainfall and 28% of tree evapotranspiration. These values reflect the low crown cover fraction of the stand (0·39) and the specific features of the Mediterranean rainfall regime (predominantly with few large storms). Nevertheless, it still is an important component of the water balance of these Mediterranean ecosystems. Copyright © 2005 John Wiley & Sons, Ltd.     

A cosmological hydrodynamic code based on the piecewise parabolic method

We present a hydrodynamical code for cosmological simulations which uses the piecewise parabolic method (PPM) to follow the dynamics of the gas component and an N -body particle–mesh algorithm for the evolution of the collisionless component. The gravitational interaction between the two components is regulated by the Poisson equation which is solved by a standard fast Fourier transform (FFT) procedure. In order to simulate cosmological flows we have introduced several modifications to the original PPM scheme which we describe in detail. Various tests of the code are presented including adiabatic expansion, single and multiple pancake formation and three-dimensional cosmological simulations with initial conditions based on the cold dark matter scenario.     

A parallelizable method for two‐phase flows in naturally‐fractured reservoirs

A parallelizable, semi‐implicit numerical method is proposed for the study of naturally‐fractured reservoir systems. It has proved to be computationally efficient in producing accurate numerical solutions for the dual‐porosity model for immiscible, two‐phase flow in such reservoirs. The method combines hybridized mixed finite elements, a new version of the modified method of characteristics, a sophisticated operator‐splitting procedure for separating the pressure calculation in the fractures from that of the saturation, another operator splitting to handle the interaction of the matrix blocks and the fractures, and domain decomposition iterative procedures for both the pressure and the saturation. It permits moderately long time steps for the pressure and the saturation in the fractures and matrix blocks by using short, inexpensive microsteps to treat the transport portion of the saturation equation in the fractures. This paper is devoted to the formulation of the method and a discussion of numerical results for five‐spot and vertical cross‐section examples.     

Geochemistry of granitic rocks from the Hercynian Sardinia-Corsica batholith: Implication for magma genesis

The Sardinia-Corsica batholith was structured in the late stage of the Hercynian orogenesis. The granitoids intrude mainly metamorphic complexes grading from zeolite up to amphibolite facies. The batholith is heterogeneous consisting of complexes with different affinity, chemical composition, age and degree of deformation. The present paper reports major-and trace-element data for selected samples coming only from Sardinian outcrops.

The rocks range from gabbro-diorite to tonalite, monzogranite and leucogranite. The two latter lithologies are the most abundant, gabbrodiorites and tonalites occurring in minor amounts and mainly in northern-central Sardinia. Over 75% of the granitoids contain microgranular enclaves of magmatic origin. The age of the rocks falls in the interval between 307 and 281 Ma. Sr isotope initial ratios are high, ranging between 0.7083 and 0.7107.

REE, Rb, Sr, Ba, Zr, Th, Ta, Hf, Co and Sc abundances were determined on selected samples. All elements follow three types of trends vs. CaO, which is used as differentiation index. Two trends show positive and negative correlations while the third one shows a bell-shaped pattern. LREE have different degrees of enrichment (La = 20−120× ch) and HREE show variable fractionation with prevailing (Tb/Yb)_n<1. The two peraluminous samples have very different geochemical characteristics.

From the geochemical point of view all the rocks coming from the Sardinian segment of the batholith display a typical calc-alkaline chemical character showing the imprint of both “normal and mature” continental arc geodynamic environments.

Geochemical trends suggest some petrogenetic constraints. The complete sequence of differentiation can be neither the product of crystal/liquid fractionation processes starting from a single basic parent magma nor the product of an AFC process. On the contrary, a two-stage model can be proposed. In the first stage a mafic melt of subcrustal origin interacted with monzogranitic magmas derived from 25–35% degree of melting of a crustal biotite amphibolitic source. Such a mixing process acted together with a crystal/liquid fractionation process to give tonalites and granodiorites. In the second stage lesser degrees of melting of the same crustal source could give the late-stage leucogranitic masses.

A possible scenario, able to take into account field and geochemical data, can be suggested for the genesis of this suite and we propose it as a working model for future investigations.     

Water and sediment movement around a coastal headland: Portland Bill,southern UK

Numerical simulations of tidal flow and sand transport around a coastal headland (Portland Bill, southern UK) were undertaken to investigate patterns of sand transport during the development of tidally induced transient eddies. Results obtained from a 2-D finite-element hydrodynamic model (TELEMAC-2D) were combined with a sediment transport model (SEDTRANS), to simulate the sand transport processes around the headland. Simulation of the tidal flow around Portland Bill has shown the formation and evolution of tidally induced transient eddies, around the headland. During the evolution of these transient eddies, no current-induced bedload (transport) eddy is formed for either side of the headland. Net bedload sand transport direction, around a coastal headland, is the result of instantaneous gradients in bedload transport rates, during flood and ebb flows, rather than the average (residual) flow. Thus, the use of residual (water) circulation to describe patterns of sediment movement as bedload is not an appropriatedapproach. In the case study presented here, the distinct characteristics of the coastal and seabed morphology around the Isle of Portland (i.e. headland shape and the bathymetry) indicate that these parameters can be influencing tidal (flow) and sediment dispersion around the headland. Such an interpretation has broader implications and applications to headland-associated sandbanks elsewhere.Responsible Editor: Hans Burchard     

Morphology and internal structure of sand shoals and sandbanks off the Dorset coast, English Channel

The morphology and internal structure of sand shoals and sandbanks around a coastal headland (Portland Bill, southern UK) are described on the basis of sidescan sonar and high-resolution seismic data sets. Morphological and architectural evidence, combined with the spatial distribution and nature of the bedrock surface, indicates that the evolution of these deposits, especially the sandbanks, may not only be hydrodynamically controlled but also morphologically controlled. The internal structure of the sand shoals reveals a simple pattern of clinoforms dipping in the same direction as their steeper profile. In contrast, the internal structure of the sandbanks reveals that their lower seismic units consist of a sedimentary core interpreted as the remains of a lowstand deposit or an early stage of the bank development under different prevailing hydrodynamic conditions. The present morphology of the sandbanks is represented by upper seismic units characterized by large-scale sandwave foresets dipping at 6–8°. The nature and morphology of the bedrock surface may also influence the development of sandbanks in some cases. For example, to the west of Portland Bill, the morphology of the bedrock surface (basal reflector) restricts the development of the Portland Bank. This information may explain the asymmetric evolution of sandbanks around Portland Bill.     

Enviromental genotoxicity evaluation: Bayesian approach for a mixture statistical model

 The data analyzed in this paper are part of the results described in Bueno et al. (2000). Three cytogenetics endpoints were analyzed in three populations of a species of wild rodent – Akodon montensis – living in an industrial, an agricultural, and a preservation area at the Itajaí Valley, State of Santa Catarina, Brazil. The polychromatic/normochromatic ratio, the mitotic index, and the frequency of micronucleated polychromatic erythrocites were used in an attempt to establish a genotoxic profile of each area. It was assumed that the three populations were in the same conditions with respect to the influence of confounding factors such as animal age, health, nutrition status, presence of pathogens, and intra- and inter-populational genetic variability. Therefore, any differences found in the endpoints analyzed could be attributed to the external agents present in each area. The statistical models used in this paper are mixtures of negative-binomials and Poisson variables. The Poisson variables are used as approximations of binomials for rare events. The mixing distributions are beta densities. The statistical analyzes are under the bayesian perspective, as opposed to the frequentist ones often considered in the literature, as for instance in Bueno et al. (2000).     

The Buffon's needle problem and the design of a geophysical survey

A logical procedure for designing a geophysical survey when sampling an area with a regular grid can be summarized as follows: model the expected anomaly, estimate the expected noise level, estimate the area of the anomaly above the noise level, choose the spacing, in both the x - and y -directions, of the measurement grid. This last step can be approached according to two main strategies: either when applying the sampling theorem to the shortest dimension of the anomaly or when using a coarser grid, leaving a more complete definition of the anomaly to a later fitting. When following this second option, it can be constructive to estimate the probability of intercepting a given anomaly with a specific segment of profile and a given profile spacing. This latter procedure is analysed by considering a rectangle approximating the plane projection of the anomaly shape and taking into account various ratios between the grid spacing and the rectangle sides. The formulae for estimating the probability of intersecting a given anomaly with a given segment of a given profile spacing are calculated. To demonstrate the accuracy of the results, a Monte-Carlo simulation on a synthetic magnetic map was performed, obtaining, for different ratios between the sides, segment length and profile interval, an agreement better than 0.1% with the analytical formulae.