Substitution Random Fields with Gaussian and Gamma Distributions: Theory and Application to a Pollution Data Set

This paper presents random field models with Gaussian or gamma univariate distributions and isofactorial bivariate distributions, constructed by composing two independent random fields: a directing function with stationary Gaussian increments and a stationary coding process with bivariate Gaussian or gamma distributions. Two variations are proposed, by considering a multivariate directing function and a coding process with a separable covariance, or by including drift components in the directing function. Iterative algorithms based on the Gibbs sampler allow one to condition the realizations of the substitution random fields to a set of data, while the inference of the model parameters relies on simple tools such as indicator variograms and variograms of different orders. A case study in polluted soil management is presented, for which a gamma model is used to quantify the risk that pollutant concentrations over remediation units exceed a given toxicity level. Unlike the multivariate Gaussian model, the proposed gamma model accounts for an asymmetry in the spatial correlation of the indicator functions around the median and for a spatial clustering of high pollutant concentrations.     

New Insights in Shallow Gas Generation from Very High Resolution Seismic and Bathymetric Surveys in the Marennes-Oléron Bay, France

Very high resolution seismic profiles, ground-truthed by vibrocores, have revealed the occurrence of kilometre-scale acoustic turbidity in the Marennes-Oléron Bay, France. Such acoustic turbidity is commonly interpreted as gas-charged sediments. Comparison between accurate historical bathymetric data and the present day bathymetry has shown high sediment accretion zones in the study area (locally up to 8 m since 1824). The superimposition of seismic and bathymetric datasets displays a striking correlation between the high sedimentation rate area and the boundaries of the acoustic turbidity, i.e. gas-charged sediments. The key role of sedimentation rates in shallow gas generation in the study area is pointed out. It is also concluded that shallow gas is probably generated at short (decadal to secular) time scales.     

Spatial and temporal zoning of hydrothermal alteration and mineralization in the Sossego iron oxide–copper–gold deposit,Carajás Mineral Province,Brazil: paragenesis and stable isotope constraints

The Sossego iron oxide–copper–gold deposit (245 Mt @ 1.1% Cu, 0.28 g/t Au) in the Carajás Mineral Province of Brazil consists of two major groups of orebodies (Pista–Sequeirinho–Baiano and Sossego–Curral) with distinct alteration assemblages that are separated from each other by a major high angle fault. The deposit is located along a regional WNW–ESE-striking shear zone that defines the contact between metavolcano–sedimentary units of the ∼2.76 Ga Itacaiúnas Supergroup and tonalitic to trondhjemitic gneisses and migmatites of the ∼2.8 Ga Xingu Complex. The deposit is hosted by granite, granophyric granite, gabbro, and felsic metavolcanic rocks. The Pista–Sequeirinho–Baiano orebodies have undergone regional sodic (albite–hematite) alteration and later sodic–calcic (actinolite-rich) alteration associated with the formation of massive magnetite–(apatite) bodies. Both these alteration assemblages display ductile to ductile–brittle fabrics. They are cut by spatially restricted zones of potassic (biotite and potassium feldspar) alteration that grades outward to chlorite-rich assemblages. The Sossego–Curral orebodies contain weakly developed early albitic alteration and very poorly developed subsequent calcic–sodic alteration. These orebodies contain well-developed potassic alteration assemblages that were formed during brittle deformation that resulted in the formation of breccia bodies. Breccia matrix commonly displays coarse mineral infill suggestive of growth into open space. Sulfides in both groups of deposits were precipitated first with potassic alteration and more importantly with a later assemblage of calcite–quartz–epidote–chlorite. In the Sequeirinho orebodies, sulfides range from undeformed to deformed; sulfides in the Sossego–Curral orebodies are undeformed. Very late, weakly mineralized hydrolytic alteration is present in the Sossego/Currral orebodies. The sulfide assemblage is dominated by chalcopyrite with subsidiary siegenite, and millerite. Pyrrhotite and pyrite are minor constituents of ore in the Sequerinho orebodies while pyrite is relatively abundant in the Sossego–Curral bodies. Oxygen isotope partitioning between mineral pairs constrains temperatures in the deposit spatially and through time. In the Sequeirinho orebody, the early sodic–calcic alteration stage was characterized by temperatures exceeding 500°C and values for the alteration fluid of 6.9 ± 0.9‰. Temperature declines outward and upward from the zone of most intense alteration. Paragenetically later copper–gold mineralization displays markedly lower temperatures (<300°C) and was characterized by the introduction of ¹⁸O-depleted hydrothermal fluids −1.8 ± 3.4‰. The calculated δD_H2O and values suggest that the fluids that formed the early calcic–sodic alteration assemblage were of formational/metamorphic or magmatic origin. The decrease of values through time may reflect influx of surficially derived waters during later alteration and mineralization events. Influx of such fluids could be related to episodic fluid overpressure, resulting in dilution and cooling of the metalliferous fluid, causing deposition of metals transported as metal chloride complexes.     

New geologic, fluid inclusion and stable isotope studies on the controversial Igarapé Bahia Cu–Au deposit, Carajás Province, Brazil

The Igarapé Bahia Cu–Au deposit in the Carajás Province, Brazil, is hosted by steeply dipping metavolcano-sedimentary rocks of the Igarapé Bahia Group. This group consists of a low greenschist grade unit of the Archean (∼2,750 Ma) Itacaiúnas Supergroup, in which other important Cu–Au and iron ore deposits of the Carajás region are also hosted. The orebody at Igarapé Bahia is a fragmental rock unit situated between chloritized basalt, with associated hyaloclastite, banded iron formation (BIF), and chert in the footwall and mainly coarse- to fine-grained turbidites in the hanging wall. The fragmental rock unit is a nearly concordant, 2 km long and 30–250 m thick orebody made up of heterolithic, usually matrix-supported rocks composed mainly of coarse basalt, BIF, and chert clasts derived from the footwall unit. Mineralization is confined to the fine-grained matrix and comprises disseminated to massive chalcopyrite accompanied by magnetite, gold, U- and light rare earth element (LREE)-minerals, and minor other sulfides like bornite, molybdenite, cobaltite, digenite, and pyrite. Gangue minerals include siderite, chlorite, amphibole, tourmaline, quartz, stilpnomelane, epidote, and apatite. A less important mineralization style at Igarapé Bahia is represented by late quartz–chalcopyrite–calcite veins that crosscut all rocks in the deposit area. Fluid inclusions trapped in a quartz cavity in the ore unit indicate that saline aqueous fluids (5 to 45 wt% NaCl + CaCl₂ equiv), together with carbonic (CO₂ ± CH₄) and low-salinity aqueous carbonic (6 wt% NaCl equiv) fluids, were involved in the mineralization process. Carbonates from the fragmental layer have δ¹³C values from −6.7 to −13.4 per mil that indicate their origin from organic and possibly also from magmatic carbon. The δ³⁴S values for chalcopyrite range from −1.1 to 5.6 per mil with an outlier at −10.8 per mil, implying that most sulfur is magmatic or leached from magmatic rocks, whereas a limited contribution of reduced and oxydized sulfur is also evident. Oxygen isotopic ratios in magnetite, quartz, and siderite yield calculated temperatures of ∼400°C and δ¹⁸O-enriched compositions (5 to 16.5 per mil) for the ore-forming fluids that suggest a magmatic input and/or an interaction with ¹⁸O-rich, probably sedimentary rocks. The late veins of the Igarapé Bahia deposit area were formed from saline aqueous fluids (2 to 60 wt% NaCl + CaCl₂ equiv) with δ¹⁸O_fluid compositions around 0 per mil that indicate contribution from meteoric fluids. With respect to geological features, Igarapé Bahia bears similarity with syngenetic, volcanic-hosted massive sulfide (VHMS)-type deposits, as indicated by the volcano-sedimentary geological context, stratabound character, and association with submarine volcanic flows, hyaloclastite, and exhalative beds such as BIF and chert. On the other hand, the highly saline ore fluids and the mineral assemblage, dominated by magnetite and chalcopyrite, with associated gold, U- and LREE-minerals and scarce pyrite, indicate that Igarapé Bahia belongs to the Fe oxide Cu–Au (IOCG) group of deposits. The available geochronologic data used to attest syngenetic or epigenetic origins for the mineralization are either imprecise or may not represent the main mineralization episode but a later, superimposed event. The C, S, and O isotopic results obtained in this study do not clearly discriminate between fluid sources. However, recent B isotope data obtained on tourmaline from the matrix of the fragmental rock ore unit (Xavier, Wiedenbeck, Dreher, Rhede, Monteiro, Araújo, Chemical and boron isotopic composition of tourmaline from Archean and Paleoproterozoic Cu–Au deposits in the Carajás Mineral Province, 1° Simpósio Brasileiro de Metalogenia, Gramado, Brazil, extended abstracts, CD-ROM, 2005) provide strong evidence of the involvement of a marine evaporitic source in the hydrothermal system of Igarapé Bahia. Evaporite-derived fluids may explain the high salinities and the low reduced sulfur mineral paragenesis observed in the deposit. Evaporite-derived fluids also exclude a significant participation of magmatic or mantle-derived fluids, reinforcing the role of nonmagmatic brines in the genesis of Igarapé Bahia. Considering this aspect and the geological features, the possibility that the deposit was generated by a hydrothermal submarine system whose elevated salinity was acquired by leaching of ancient evaporite beds should be evaluated.     

Together,but separate: neighborhood-scale patterns and correlates of spatial segregation between male and female same-sex couples in Melbourne and Sydney

Considerable research has examined the relationships between queer communities and the city, particularly within a qualitative, North American context. This paper provides the first quantitative treatment of the spatial patterning of male and female same-sex couples in Australia’s two largest urban centers, examining differences between the two groups and their potential drivers. Male and female same-sex couples are shown to concentrate in distinct regions of Melbourne and Sydney that are segregated from each other, but nonetheless geographically close: together, but separate. Spatial regression analyses provided an incomplete explanation of these patterns. A number of factors were identified, endogenous and exogenous to the same-sex couple community, which provide a potential explanatory framework for explaining the aggregation of these groups. Very little emerges from this analysis which explains segregation for these groups. The ultimate etiology of segregation between male and female same-sex couples in Melbourne and Sydney remains elusive.     

The single-cell transport problem for two-phase flow with polymer

Polymer injection is a widespread strategy in enhanced oil recovery. Polymer increases the water viscosity and creates a more favorable mobility ratio between the injected water and the displaced oil. The computational cost of simulating polymer injection can be significantly reduced if one splits the governing system of two-phase equations into a pressure equation and a set of saturation/component equations and use a Gauss–Seidel algorithm with optimal cell ordering to solve the nonlinear systems arising from an implicit discretization of the saturation/component equations. This approach relies on a robust single-cell solver that computes the saturation and polymer concentration of a cell, given the total flux and the saturation and polymer concentration of the neighboring cells. In this paper, we consider a relatively comprehensive polymer model used in an industry-standard simulator, and show that, in the case of a discretization using a two-point flux approximation, the single-cell problem always admits a solution that is also unique.     

Lateral variations of attenuation coefficients, group and phase velocities of rayleigh waves in Europe

Rayleigh wave group and phase velocities and attenuation coefficients are investigated in the 15–60-s period range for earthquakes occurring in the European area and registered at European stations. Observed group and phase velocities show clear differences in the period range 20–40 s, whereas attenuation coefficients are differentiated from 45 s to larger periods. From the shape of dispersion and attenuation curves, two broad zones have been isolated and designated as Western Europe and Southeastern Europe. Although some scatter is present, specially for attenuation, the data suggest that in general the shear velocity and Q_β values are higher for Southeastern Europe than for Western Europe. For Southeastern Europe there is a channel of low velocity and low Q_β, apparently correlated, located at the lower crust: in both regions there is a channel of low velocity and low Q_β values located at the uppermost mantle. In each case there also exists a channel of high Q_β values located at the lower crust but at slightly different depths.     

Geoid,Pangea and convection

We show that the present geoid has a simple low-order configuration with an axis of symmetry in the equatorial plane. We show further that it is a “tennis-ball” pattern, with an equatorial high belt and a polar low one, which is clearly controlled by the rotation of the Earth. Finally, we show that the outline of Pangea between at least 200 Ma and 125 Ma ago lay along a great circle passing through the paleo-poles of rotation. Thus, it also had an axis of symmetry in the equatorial plane. This hemispheric super-continent configuration ended in Middle Cretaceous time during a major geologic catastrophe which was accompanied by high rates of spreading, hotspot outbreaks and high sea-level stands. We interpret this evidence in terms of separate steady state lower mantle convection, responsible for the present geoid, weakly coupled to the upper mantle one. This weak coupling leads to the hemispheric continent configuration which ends when excessive heating of the upper mantle due to the insulating continental cap leads to continent dispersal. The complete cycle, from one supercontinent to the next, might be of the order of 400 Ma.