Velocities of compressional and shear waves in limestones

Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional- and shear-wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse-echo method. The measurements were made both in vacuum-dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin-impregnated polished thin sections, X-ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional- and shear-wave velocities (V-p and V-s, respectively) decrease with increasing porosity and that V-p decreases approximately twice as fast as V-s. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore-structure-dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot-Gassmann velocity values are greater than the measured velocity values due to the rock-fluid interaction. This is not accounted for in the Biot-Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time-average relationships overestimated the measured velocities even more than the Biot model.

Fate of prions in soil: Detergent extraction of PrP from soils

The transmissible spongiform encephalopathies (TSEs) are caused by infectious agents whose structures have not been fully characterized but include abnormal forms of the host protein PrP, designated PrPSc, which are deposited in infected tissues. The transmission routes of scrapie and chronic wasting disease (CWD) seem to include environmental spread in their epidemiology, yet the fate of TSE agents in the environment is poorly understood. There are concerns that, for example, buried carcasses may remain a potential reservoir of infectivity for many years. Experimental determination of the environmental fate requires methods for assessing binding/elution of TSE infectivity, or its surrogate marker PrPSc, to and from materials with which it might interact. We report a method using Sarkosyl for the extraction of murine PrPSc, and its application to soils containing recombinant ovine PrP (recPrP). Elution properties suggest that PrP binds strongly to one or more soil components. Elution from a clay soil also required proteinase K digestion, suggesting that in the clay soil binding occurs via the N-terminal of PrP to a component that is absent from the sandy soils tested.

Fate of prions in soil: Longevity and migration of recPrP in soil columns

Research into transmissible spongiform encephalopathy (TSE) diseases has become a high priority worldwide in recent years yet remarkably little is known about the behaviour of TSE infectivity in the environment. The resilience and stability of prion proteins could lead to soils becoming a potential reservoir of TSE infectivity as a result of contamination from activities such as infected carcass burial or the dispersion of effluents from slaughter houses, or by contamination of pastures by infected animals, (e.g. scrapie in sheep). Knowledge of the fate of prion proteins in soils, and associated physico-chemical conditions which favour migration, can be used to help prevent re-infection of animals through grazing, to protect watercourses and develop good management practices. In two consecutive experiments of 9 and 6 months, the migration of recombinant ovine PrP (recPrP) in soil columns was followed under contrasting levels of microbial activity (normal versus reduced), under varying regimes of soil water content and redox potential, and in two different soil types (loamy sand and clay loam). At each analysis time (1, 3, 6 or 9 months), in both soil types, full-length recPrP was detected in the original contaminated layer, indicating the resilience and stability of recPrP under varied soil conditions, even in the presence of active soil microbial populations. Evidence of protein migration was found in every soil column at the earliest analysis time (1 or 3 months), but was restricted to a maximum distance of 1 cm, indicative of limited initial mobility in soils followed by strong adsorption over the following days to weeks. The survival of recPrP in the soil over a period of at least 9 months was demonstrated. In this study, recPrP was used as an indicator for potential TSE infectivity, although infectivity tests should be carried out before conclusions can be drawn regarding the infection risk posed by prions in soil. However, it has been demonstrated that soil is likely to act as a significant barrier to the dispersion of contaminated material at storage or burial sites. (c) 2007 Elsevier Ltd. All rights reserved.

Seasonal deposition of Holocene banded sediments in the Severn Estuary Levels (southwest Britain): palynological and sedimentological evidence

Banded sediments outcrop widely in the intertidal zone of the Severn Estuary and have been suggested, on the basis of textural analysis, to have formed in response to seasonal variations in sea temperature and windiness (Holocene, 14 (2004) 536). Here palynological and sedimentological analyses of banded sediments of mid-Holocene date from Gold Cliff, on the Welsh side of the Severn Estuary, are combined to test and further develop the hypothesis of seasonal deposition. Pollen percentage and concentration data are presented from a short sequence of bands to establish whether textural variations in the bands coincide with variations in pollen content reflecting seasonal flowering patterns. It is shown that fine-grained band parts contain higher total pollen concentrations, and a higher proportion of pollen from late spring- to summer-flowering plants, than coarse-grained band parts. Pollen in the coarser deposits appears primarily to reflect deposition from the buffering `reservoir' of suspended pollen in the estuarine water-body and from rivers, when there is little pollen in the air in winter, while the finer sediments contain pollen deposited from the atmosphere during the flowering season, superimposed on these `background' sources. The potential of such deposits for refining chronologies and identifying seasonality of coastal processes is noted, and the results of charcoal particle analysis of the bands presented as an example of how they have the potential to shed light on seasonal and annual patterns of human activity. (C) 2004 Elsevier Ltd. All rights reserved.

A cross-calibration of chlorine isotopic measurements and suitability of seawater as the international reference material

A collection of 24 seawaters from various worldwide locations and differing depth was culled to measure their chlorine isotopic composition (delta(37)Cl). These samples cover all the oceans and large seas: Atlantic, Pacific, Indian and Antarctic oceans, Mediterranean and Red seas. This collection includes nine seawaters from three depth profiles down to 4560 mbsl. The standard deviation (2sigma) of the delta(37)Cl of this collection is +/-0.08 parts per thousand, which is in fact as large as our precision of measurement ( +/- 0.10 parts per thousand). Thus, within error, oceanic waters seem to be an homogeneous reservoir. According to our results, any seawater could be representative of Standard Mean Ocean Chloride (SMOC) and could be used as a reference standard. An extended international cross-calibration over a large range of delta(37)Cl has been completed. For this purpose, geological fluid samples of various chemical compositions and a manufactured CH3Cl gas sample, with delta(37)Cl from about -6 parts per thousand to +6 parts per thousand have been compared. Data were collected by gas source isotope ratio mass spectrometry (IRMS) at the Paris, Reading and Utrecht laboratories and by thermal ionization mass spectrometry (TIMS) at the Leeds laboratory. Comparison of IRMS values over the range -5.3 parts per thousand to +1.4 parts per thousand plots on the Y=X line, showing a very good agreement between the three laboratories. On 11 samples, the trend line between Paris and Reading Universities is: delta(37)Cl(Reading)= (1.007 +/- 0.009)delta(37)Cl(Paris) - (0.040 +/- 0.025), with a correlation coefficient: R-2 = 0.999. TIMS values from Leeds University have been compared to IRMS values from Paris University over the range -3.0 parts per thousand to +6.0 parts per thousand. On six samples, the agreement between these two laboratories, using different techniques is good: delta(37)Cl(Leeds)=(1.052 +/- 0.038)delta(37)Cl(Paris) + (0.058 +/- 0.099), with a correlation coefficient: R-2 = 0.995. The present study completes a previous cross-calibration between the Leeds and Reading laboratories to compare TIMS and IRMS results (Anal. Chem. 72 (2000) 2261). Both studies allow a comparison of IRMS and TIMS techniques between delta(37)Cl values from -4.4 parts per thousand to +6.0 parts per thousand and show a good agreement: delta(37)Cl(TIMS)=(1.039 +/- 0.023)delta(37)Cl(IRMS)+(0.059 +/- 0.056), with a correlation coefficient: R-2 = 0.996. Our study shows that, for fluid samples, if chlorine isotopic compositions are near 0 parts per thousand, their measurements either by IRMS or TIMS will give comparable results within less than +/- 0.10 parts per thousand, while for delta(37)Cl values as far as 10 parts per thousand (either positive or negative) from SMOC, both techniques will agree within less than +/- 0.30 parts per thousand. (C) 2004 Elsevier B.V. All rights reserved.

The Upper Valanginian (Early Cretaceous) positive carbon-isotope event recorded in terrestrial plants

Our understanding of the ancient ocean-atmosphere system has focused on oceanic proxies. However, the study of terrestrial proxies is equally necessary to constrain our understanding of ancient climates and linkages between the terrestrial and oceanic carbon reservoirs. We have analyzed carbon-isotope ratios from fossil plant material through the Valanginian and Lower Hauterivian from a shallow-marine, ammonite-constrained succession in the Crimean Peninsula of the southern Ukraine in order to determine if the Upper Valanginian positive carbon-isotope excursion is expressed in the atmosphere. delta(13)C(plant) values fluctuate around -23% to -22% for the Valanginian-Hauterivian, except during the Upper Valanginian where delta(13)C(plant) values record a positive excursion to similar to-18%. Based upon ammonite biostratigraphy from Crimea, and in conjunction with a composite Tethyan marine delta(13)C(carb) curve, several conclusions can be drawn: (1) the delta(13)C(plant) record indicates that the atmospheric carbon reservoir was affected; (2) the defined ammonite correlations between Europe and Crimea are synchronous; and (3) a change in photosynthetic carbon-isotope fractionation, caused by a decrease in atmospheric PCO2, occurred during the Upper Valanginian Positive delta(13)C excursion. Our new data, combined with other paleoenvironmental and paleoclimatic information, indicate that the Upper Valanginian was a cool period (icehouse) and highlights that the Cretaceous period was interrupted by periods of cooling and was not an equable climate as previously thought. (C) 2005 Elsevier B.V. All rights reserved.

Marine04 marine radiocarbon age calibration, 0-26 cal kyr BP

New radiocarbon calibration curves, IntCal04 and Marine04, have been constructed and internationally ratified to replace the terrestrial and marine components of IntCal98. The new calibration data sets extend an additional 2000 yr, from 0-26 cal kyr BP (Before Present, 0 cal. BP = AD 1950), and provide much higher resolution, greater precision, and more detailed structure than IntCal98. For the Marine04 curve, dendrochronologically-dated tree-ring samples, converted with a box diffusion model to marine mixed-layer ages, cover the period from 0-10.5 call kyr BR Beyond 10.5 cal kyr BP, high-resolution marine data become available from foraminifera in varved sediments and U/Th-dated corals. The marine records are corrected with site-specific C-14 reservoir age information to provide a single global marine mixed-layer calibration from 10.5-26.0 cal kyr BR A substantial enhancement relative to IntCal98 is the introduction of a random walk model, which takes into account the uncertainty in both the calendar age and the C-14 age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine04) are discussed here. The tree-ring data sets, sources of uncertainty, and regional offsets are presented in detail in a companion paper by Reimer et al. (this issue).

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.

Sonic to ultrasonic Q of sandstones and limestones: Laboratory measurements at in situ pressures

Laboratory measurements of the attenuation and velocity dispersion of compressional and shear waves at appropriate frequencies, pressures, and temperatures can aid interpretation of seismic and well-log surveys as well as indicate absorption mechanisms in rocks. Construction and calibration of resonant-bar equipment was used to measure velocities and attenuations of standing shear and extensional waves in copper-jacketed right cylinders of rocks (30 cm in length, 2.54 cm in diameter) in the sonic frequency range and at differential pressures up to 65 MPa. We also measured ultrasonic velocities and attenuations of compressional and shear waves in 50-mm-diameter samples of the rocks at identical pressures. Extensional-mode velocities determined from the resonant bar are systematically too low, yielding unreliable Poisson's ratios. Poisson's ratios determined from the ultrasonic data are frequency corrected and used to calculate the sonic-frequency compressional-wave velocities and attenuations from the shear- and extensional-mode data. We calculate the bulk-modulus loss. The accuracies of attenuation data (expressed as 1000/Q, where Q is the quality factor) are +/- 1 for compressional and shear waves at ultrasonic frequency, +/- 1 for shear waves, and +/- 3 for compressional waves at sonic frequency. Example sonic-frequency data show that the energy absorption in a limestone is small (Q(P) greater than 200 and stress independent) and is primarily due to poroelasticity, whereas that in the two sandstones is variable in magnitude (Q(P) ranges from less than 50 to greater than 300, at reservoir pressures) and arises from a combination of poroelasticity and viscoelasticity. A graph of compressional-wave attenuation versus compressional-wave velocity at reservoir pressures differentiates high-permeability (> 100 mD, 9.87 X 10(-14) m(2)) brine-saturated sandstones from low-permeability (< 100 mD, 9.87 X 10 (14) m(2)) sandstones and shales.

Chlorophyll-a in the rivers of eastern England

Chlorophyll-a concentration variations are described for two major river basins in England, the Humber and the Thames and related to catchment characteristics and nutrient concentrations across a range of rural, agricultural and urban/industrial settings. For all the rivers there are strong seasonal variations, with concentrations peaking in the spring and summer time when biological activity is at its highest. However, there are large variations in the magnitude of the seasonal effects across the rivers. For the spring-summer low-flow periods, average concentrations of chlorophyll-a correlate with soluble reactive phosphor-us (SRP). Chlorophyll-a is also correlated with particulate nitrogen (PN), organic carbon (POC) and suspended sediments. However, the strongest relationships are with catchment area and flow, where two straight line relationships are observed. The results indicate the importance of residence times for determining planktonic growth within the rivers. This is also indicated by the lack of chlorophyll-a response to lowering of SRP concentrations in several of the rivers in the area due to phosphorus stripping of effluents at major sewage treatment works. A key control on chlorophyll-a concentration may be the input of canal and reservoir waters during the growing period: this too relates to issues of residence times. However, there may well be a complex series of factors influencing residence time across the catchments due to features such as inhomogeneous flow within the catchments, a fractal distribution of stream channels that leads to a distribution of residence times and differences in planktonic inoculation sources. Industrial pollution on the Aire and Calder seems to have affected the relationship of chlorophyll-a with PN and POC. The results are discussed in relation to the Water Framework Directive. (c) 2006 Elsevier B.V. All rights reserved.

IntCal04 terrestrial radiocarbon age calibration, 0-26 cal kyr BP

A new calibration curve for the conversion of radiocarbon ages to calibrated (cal) ages has been constructed and internationally ratified to replace ImCal98, which extended from 0-24 cal kyr BP (Before Present, 0 cal BP = AD 1950). The new calibration data set for terrestrial samples extends from 0-26 cal kyr BP, but with much higher resolution beyond 11.4 cal kyr BP than ImCal98. Dendrochronologically-dated tree-ring samples cover the period from 0-12.4 cal kyr BP. Beyond the end of the tree rings, data from marine records (corals and foraminifera) are converted to the atmospheric equivalent with a site-specific marine reservoir correction to provide terrestrial calibration from 12.4-26.0 cal kyr BP. A substantial enhancement relative to ImCal98 is the introduction of a coherent statistical approach based on a random walk model, which takes into account the uncertainty in both the calendar age and the C-14 age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The tree-ring data sets, sources of uncertainty, and regional offsets are discussed here. The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine 04) are discussed in brief, but details are presented in Hughen et al. (this issue a). We do not make a recommendation for calibration beyond 26 cal kyr BP at this time; however, potential calibration data sets are compared in another paper (van der Plicht et al., this issue).

Shifts in rhizosphere microbial communities and enzyme activity of Poa alpina across an alpine chronosequence

This study quantifies the influence of Poa alpina on the soil microbial community in primary succession of alpine ecosystems, and whether these effects are controlled by the successional stage. Four successional sites representative of four stages of grassland development (initial, 4 years (non-vegetated); pioneer, 20 years; transition, 75 years; mature, 9500 years old) on the Rotmoos glacier foreland, Austria, were sampled. The size, composition and activity of the microbial community in the rhizosphere and bulk soil were characterized using the chloroform-fumigation extraction procedure, phospholipid fatty acid (PLFA) analysis and measurements of the enzymes beta-glucosidase, beta-xylosidase, N-acetyl-beta-glucosaminidase, leucine aminopeptidase, acid phosphatase and sulfatase. The interplay between the host plant and the successional stage was quantified using principal component (PCA) and multidimensional scaling analyses. Correlation analyses were applied to evaluate the relationship between soil factors (C-org, N-t, C/N ratio, pH, ammonium, phosphorus, potassium) and microbial properties in the bulk soil. In the pioneer stage microbial colonization of the rhizosphere of P. alpina was dependent on the reservoir of microbial species in the bulk soil. As a consequence, the rhizosphere and bulk soil were similar in microbial biomass (ninhydrin-reactive nitrogen (NHR-N)), community composition (PLFA), and enzyme activity. In the transition and mature grassland stage, more benign soil conditions stimulated microbial growth (NHR-N, total amount of PLFA, bacterial PLFA, Gram-positive bacteria, Gram-negative bacteria), and microbial diversity (Shannon index H) in the rhizosphere either directly or indirectly through enhanced carbon allocation. In the same period, the rhizosphere microflora shifted from a G(-) to a more G(+), and from a fungal to a more bacteria-dominated community. Rhizosphere beta-xylosidase, N-acetyl-beta-glucosaminidase, and sulfatase activity peaked in the mature grassland soil, whereas rhizosphere leucine aminopeptidase, beta-glucosidase, and phosphatase activity were highest in the transition stage, probably because of enhanced carbon and nutrient allocation into the rhizosphere due to better growth conditions. Soil organic matter appeared to be the most important driver of microbial colonization in the bulk soil. The decrease in soil pH and soil C/N ratio mediated the shifts in the soil microbial community composition (bacPLFA, bacPLFA/fungPLFA, G(-), G(+)/G(-)). The activities of beta-glucosidase, beta-xylosidase and phosphatase were related to soil ammonium and phosphorus, indicating that higher decomposition rates enhanced the nutrient availability in the bulk soil. We conclude that the major determinants of the microllora vary along the successional gradient: in the pioneer stage the rhizosphere microflora was primarily determined by the harsh soil environment; under more favourable environmental conditions, however, the host plant selected for a specific microbial community that was related to the dynamic interplay between soil properties and carbon supply. (C) 2004 Elsevier Ltd. All rights reserved.

Lava production at Soufriere Hills Volcano, Montserrat: 1995-2009

[1] We estimate that about 1 km3 of andesitic lava has been produced at Soufrière Hills Volcano, Montserrat from 1995 to 2009. There were three major episodes of extrusion, each lasting about 2 to 3.5 years and producing about 280 to 340 M m3 of lava, and one minor episode. Our estimates account for the dense rock equivalent volumetric contributions from the core and talus components of the lava dome, pyroclastic flow deposits and air-fall deposits. By 2005 at least two thirds of the erupted mass has already entered the sea. The average lava flux across the major extrusion episodes has been 3–5 m3s−1, with short-period (10–15 days) pulses up to 10–20 m3s−1. The first and third episodes of extrusion show similar flux histories suggesting similar behaviour of the system ten years apart. Waning flux towards the end of each episode may be caused by declining overpressure in the magma reservoir.

The unsteady flow of a weakly compressible fluid in a thin porous layer. I: Two-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a two-dimensional reservoir in an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting or extracting fluid. Numerical solution of this problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l. This is a situation which occurs frequently in the application to oil reservoir recovery. Under the assumption that epsilon=h/l<<1, we show that the pressure field varies only in the horizontal direction away from the wells (the outer region). We construct two-term asymptotic expansions in epsilon in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive analytical expressions for all significant process quantities. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the reservoir, epsilon, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighborhood of wells and away from wells.

A Bayesian approach for multi-model downscaling: Seasonal forecasting of regional rainfall and river flows in South America

This study addresses three issues: spatial downscaling, calibration, and combination of seasonal predictions produced by different coupled ocean-atmosphere climate models. It examines the feasibility Of using a Bayesian procedure for producing combined, well-calibrated downscaled seasonal rainfall forecasts for two regions in South America and river flow forecasts for the Parana river in the south of Brazil and the Tocantins river in the north of Brazil. These forecasts are important for national electricity generation management and planning. A Bayesian procedure, referred to here as forecast assimilation, is used to combine and calibrate the rainfall predictions produced by three climate models. Forecast assimilation is able to improve the skill of 3-month lead November-December-January multi-model rainfall predictions over the two South American regions. Improvements are noted in forecast seasonal mean values and uncertainty estimates. River flow forecasts are less skilful than rainfall forecasts. This is partially because natural river flow is a derived quantity that is sensitive to hydrological as well as meteorological processes, and to human intervention in the form of reservoir management.