Guide and examples for users of the ‘SPACEKAP’ style file

Guide and examples for users of the SPACEKAP style fileBasic instructions     

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.     

Anaerobic methane oxidation and a deep H2S sink generate isotopically heavy sulfides in Black Sea sediments

The main terminal processes of organic matter mineralization in anoxic Black Sea sediments underlying the sulfidic water column are sulfate reduction in the upper 2-4 m and methanogenesis below the sulfate zone. The modern marine deposits comprise a ca. 1-m-deep layer of coccolith ooze and underlying sapropel, below which sea water ions penetrate deep down into the limnic Pleistocene deposits from >9000 years BP. Sulfate reduction rates have a subsurface maximum at the SO₄²⁻-CH₄ transition where H₂S reaches maximum concentration. Because of an excess of reactive iron in the deep limnic deposits, most of the methane-derived H₂S is drawn downward to a sulfidization front where it reacts with Fe(III) and with Fe²⁺ diffusing up from below. The H₂S-Fe²⁺ transition is marked by a black band of amorphous iron sulfide above which distinct horizons of greigite and pyrite formation occur. The pore water gradients respond dynamically to environmental changes in the Black Sea with relatively short time constants of ca. 500 yr for SO₄²⁻ and 10 yr for H₂S, whereas the FeS in the black band has taken ca. 3000 yr to accumulate. The dual diffusion interfaces of SO₄²⁻-CH₄ and H₂S-Fe²⁺ cause the trapping of isotopically heavy iron sulfide with δ³⁴S = +15 to +33‰ at the sulfidization front. A diffusion model for sulfur isotopes shows that the SO₄²⁻ diffusing downward into the SO₄²⁻-CH₄ transition has an isotopic composition of +19‰, close to the +23‰ of H₂S diffusing upward. These isotopic compositions are, however, very different from the porewater SO₄²⁻ (+43‰) and H₂S (−15‰) at the same depth. The model explains how methane-driven sulfate reduction combined with a deep H₂S sink leads to isotopically heavy pyrite in a sediment open to diffusion. These results have general implications for the marine sulfur cycle and for the interpretation of sulfur isotopic data in modern sediments and in sedimentary rocks throughout earth’s history.     

Phytoplankton production and nutrient distributions in a subtropical estuary: Importance of freshwater flow

The relationships between phytoplankton productivity, nutrient distributions, and freshwater flow were examined in a seasonal study conducted in Escambia Bay, Florida, USA, located in the northeastern Gulf of Mexico. Five sites oriented along the salinity gradient were sampled 24 times over the 28-mo period from 1999 to 2001. Water column profiles of temperature and salinity were measured along with surface chlorophyll and surface inorganic nutrient concentrations. Primary productivity was measured at 2 sites on 11 dates, and estimated for the remaining dates and sites using an empirical regression model relating phytoplankton net production to the product of chlorophyll, euphotic zone depth, and daily solar insolation. Freshwater flow into the system varied markedly over the study period with record low flow during 2000, a flood event in March 2001, and subsequent resumption of normal flow. Flushing times ranged from 1 d during the flood to 20 d during the drought. Freshwater input strongly affected surface salinity distributions, nutrient flux, chlorophyll, and primary productivity. The flood caused high turbidity and rapid flushing, severely reducing phytoplankton production and biomass accumulation. Following the flood, phytoplankton biomass and productivity sharply increased. Analysis of nutrient distributions suggested Escambia Bay phytoplankton alternated between phosphorus limitation during normal flow and nitrogen limitation during low flow periods. This study found that Escambia Bay is a moderately productive estuary, with an average annual integrated phytoplankton production rate of 290 g C m⁻² yr⁻¹.     

The Kaidun meteorite: Mineralogy of an unusual CM1 lithology

Abstract Kaidun is a breccia of disparate enstatite and carbonaceous chondrite clasts that continues to provide real surprises. Many Kaidun clasts have been intensely altered by aqueous fluids, as evidenced by the widespread occurrence of ferromagnesian phyllosilicates and by the presence of carbonate- and phyllo-silicate-filled veins. In this report, we describe an unusual CM lithology containing many mineralogical features not previously reported from any meteorite, including pyrrhotite, with exclusive needlelike morphologies and thick mantles of phyllosilicate, and complex aggregates of phyllosilicate, melanite garnet, crosscut by pentlandite veins. The latter features appear to be due in large part to extensive hydrothermal alteration at temperatures on the order of 450 °C, which is significantly higher than that attained during secondary processing from other known CM material.     

Structure of the crust and upper mantle in Canada

The seismic probing of the crust and upper mantle in Canada started in 1938 and since then has involved many government and university groups using a wide variety of techniques. These have included simple profiling with both wide and narrow station spacing, areal time-term surveys, detailed deep reflection experiments, very long-range refraction studies and the analysis of surface wave dispersion between stations of the Canadian Standard Network.

A review of the published interpretation leads to the general conclusion that:

1. (1) P_n-velocities vary from a value possibly as low as 7.7 km/sec under Vancouver Island to 8.6 km/sec and higher in the extreme eastern part of the shield and some parts of the Atlantic coast.

2. (2) Large areas of Canada have a crustal thickness of 30–40 km, with Vancouver Island, the southwestern Prairies, the Lake Superior basin and parts of the eastern shield of Quebec being thicker. No continental area in Canada is known to have a crust thinner than 29 km.

3. (3) The Riel discontinuity — a deep intra-crustal reflector and sometime refractor, is widely reported in the Prairies and Manitoba. It is not seen to the north in the vicinity of Great Slave Lake, nor in the Hudson Bay, Lake Superior and Maritime regions, nor in the interior of British Columbia. It may be present in some areas of the eastern shield.

4. (4) As experiments have become more detailed, crustal structures of greater complexity have been revealed. The concept that crustal structure becomes simpler with increasing depth is apparently unfounded.

Long-range refraction studies suggest that the Gutenberg P-wave low-velocity channel is poorly developed under the Canadian Shield. The analysis of the dispersion of surface waves, however, suggests that the channel is better developed for S-waves, and is present throughout the country. The lid of the channel is deepest under the central shield and shallowest under the Cordillera.     

Petrophysics and mineral exploration: a workflow for data analysis and a new interpretation framework

As mineral exploration seeks deeper targets, there will be a greater reliance on geophysical data and a better understanding of the geological meaning of the responses will be required, and this must be achieved with less geological control from drilling. Also, exploring based on the mineral system concept requires particular understanding of geophysical responses associated with altered rocks. Where petrophysical datasets of adequate sample size and measurement quality are available, physical properties show complex variations, reflecting the combined effects of various geological processes. Large datasets, analysed as populations, are required to understand the variations. We recommend the display of petrophysical data as frequency histograms because the nature of the data distribution is easily seen with this form of display. A petrophysical dataset commonly contains a combination of overlapping sub-populations, influenced by different geological factors. To understand the geological controls on physical properties in hard rock environments, it is necessary to analyse the petrophysical data not only in terms of the properties of different rock types. It is also necessary to consider the effects of processes such as alteration, weathering, metamorphism and strain, and variables such as porosity and stratigraphy. To address this complexity requires that much more supporting geological information be acquired than in current practice. The widespread availability of field portable instruments means quantitative geochemical and mineralogical data can now be readily acquired, making it unnecessary to rely primarily on categorical rock classification schemes. The petrophysical data can be combined with geochemical, petrological and mineralogical data to derive explanations for observed physical property variations based not only on rigorous rock classification methods, but also in combination with quantitative estimates of alteration and weathering. To understand how geological processes will affect different physical properties, it is useful to define three end-member forms of behaviour. Bulk behaviour depends on the physical properties of the dominant mineral components. Density and, to a lesser extent, seismic velocity show such behaviour. Grain and texture behaviour occur when minor components of the rock are the dominate controls on its physical properties. Grain size and shape control grain properties, and for texture properties the relative positions of these grains are also important. Magnetic and electrical properties behave in this fashion. Thinking in terms of how geological processes change the key characteristics of the major and minor mineralogical components allows the resulting changes in physical properties to be understood and anticipated.