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.

Public values for badgers, bovine TB reduction and management strategies

This paper applies an attribute-based stated choice experiment approach to estimate the value that society places on changes to the size of the badger population in England and Wales. The study was undertaken in the context of a rising incidence of bovine tuberculosis (bTB) in cattle and the government's review of current bTB control policy. This review includes consideration of culling badgers to reduce bTB in cattle, since badgers are thought to be an important wildlife reservoir for the disease. The design of the CE involved four attributes (size of badger population, cattle slaughtered due to bTB, badger management strategy and household tax) at four levels with eight choice sets of two alternatives presented to respondents. Telephone interviews were undertaken with over 400 respondents, which elicited their attitudes and preferences concerning badgers, bTB in cattle and badger management strategies. The study estimated a willingness to pay of 0.10 pound per household per year per 100,000 badgers and 1.52 pound per household per year per 10,000 cattle slaughtered due to bTB which aggregated to 22 per badger and 3298 pound per bTB slaughtered animal for all households in England and Wales. Management strategy toward badgers had a very high valuation, highlighting the emotive issue of badger culling for respondents and the importance of government policy towards badgers.

The relationship between multi-species catch and effort: among fishery comparisons

The control of fishing mortality via fishing effort remains fundamental to most fisheries management strategies even at the local community or co-management level. Decisions to support such strategies require knowledge of the underlying response of the catch to changes in effort. Even under adaptive management strategies, imprecise knowledge of the response is likely to help accelerate the adaptive learning process. Data and institutional capacity requirements to employ multi-species biomass dynamics and age-structured models invariably render their use impractical particularly in less developed regions of the world. Surplus production models fitted to catch and effort data aggregated across all species offer viable alternatives. The current paper seeks models of this type that best describe the multi-species catch–effort responses in floodplain-rivers, lakes and reservoirs and reef-based fisheries based upon among fishery comparisons, building on earlier work. Three alternative surplus production models were fitted to estimates of catch per unit area (CPUA) and fisher density for 258 fisheries in Africa, Asia and South America. In all cases examined, the best or equal best fitting model was the Fox type, explaining up to 90% of the variation in CPUA. For lake and reservoir fisheries in Africa and Asia, the Schaefer and an asymptotic model fitted equally well. The Fox model estimates of fisher density (fishers km−2) at maximum yield (iMY) for floodplain-rivers, African lakes and reservoirs and reef-based fisheries are 13.7 (95% CI [11.8, 16.4]); 27.8 (95% CI [17.5, 66.7]) and 643 (95% CI [459,1075]), respectively and compare well with earlier estimates. Corresponding estimates of maximum yield are also given. The significantly higher value of iMY for reef-based fisheries compared to estimates for rivers and lakes reflects the use of a different measure of fisher density based upon human population size estimates. The models predict that maximum yield is achieved at a higher fishing intensity in Asian lakes compared to those in Africa. This may reflect the common practice in Asia of stocking lakes to augment natural recruitment. Because of the equilibrium assumptions underlying the models, all the estimates of maximum yield and corresponding levels of effort should be treated with caution.

Comparative micromorphology of nectariferous and nectarless labellar spurs in selected clades of subtribe Orchidinae (Orchidaceae)

Floral nectar spurs are widely considered to influence pollinator behaviour in orchids. Spurs of 21 orchid species selected from within four molecularly circumscribed clades of subtribe Orchidinae (based on Platanthera s.l., Gymnadenia-Dactylorhiza s.l., Anacamptis s.l., Orchis s.s.) were examined under light and scanning electron microscopes in order to estimate correlations between nectar production (categorized as absent, trace, reservoir), interior epidermal papillae (categorized as absent, short, medium, long) and epidermal cell striations (categorized as apparently absent, weak, moderate, strong). Closely related congeneric species scored similarly, but more divergent species showed less evidence of phylogenetic constraints. Nectar secretion was negatively correlated with striations and positively correlated with papillae, which were especially frequent and large in species producing substantial reservoirs of nectar. We speculate that the primary function of the papillae is conserving energy through nectar resorption and explain the presence of large papillae in a minority of deceit-pollinated species by arguing that the papillae improve pollination because they are a tactile expectation of pollinating insects. In contrast, the prominence of striations may be a 'spandrel', simply reflecting the thickness of the overlying cuticle. Developmentally, the spur is an invagination of the labellum; it is primarily vascularized by a single 'U'-shaped primary strand, with smaller strands present in some species. Several suggestions are made for developing further, more targeted research programmes. (C) 2009 The Linnean Society of London, Botanical Journal of the Linnean Society, 2009, 160, 369-387.

Corrosion by sulfate reducing bacteria that utilize nitrate

16S rRNA gene sequencing was used to identify a sulfate-reducing bacterium (SRB) from a Danish North Sea oilfield water injection system. This species was cultivated, purified and subsequently identified as being >97% similar to Desulfovibrio gracilis. Like some other Desulfovibrio species this SRB, strain OP102, could reduce nitrate as an electron acceptor and produce ammonia in the absence of sulfate. In addition, in the presence of sulfate, when nitrate was dosed at 100 mg/l it was again reduced by the bacterium, with some ammonium production. Therefore, this mechanism could be important in oilfield systems where nitrate is applied to prevent sulfide generation by SRB which leads to reservoir souring. In static tests the influence of this Desulfovibrio on corrosion was assessed using carbon steel coupons, in the presence of sulfate and in the presence of sulfate with 100 mg/l nitrate. Corrosion rates were less than 1.5 mpy when coupons were incubated in the same water, with sulfate and with nitrate. Furthermore, the occurrence of pitting corrosion was fairly low under all circumstances.

Solutions of a two-dimensional dam-break problem

Solutions of a two-dimensional dam break problem are presented for two tailwater/reservoir height ratios. The numerical scheme used is an extension of one previously given by the author [J. Hyd. Res. 26(3), 293–306 (1988)], and is based on numerical characteristic decomposition. Thus approximate solutions are obtained via linearised problems, and the method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids non-physical, spurious oscillations.

Last Glacial Maximum CO2 and δ13C successfully reconciled

During the Last Glacial Maximum (LGM, ∼21,000 years ago) the cold climate was strongly tied to low atmospheric CO2 concentration (∼190 ppm). Although it is generally assumed that this low CO2 was due to an expansion of the oceanic carbon reservoir, simulating the glacial level has remained a challenge especially with the additional δ13C constraint. Indeed the LGM carbon cycle was also characterized by a modern-like δ13C in the atmosphere and a higher surface to deep Atlantic δ13C gradient indicating probable changes in the thermohaline circulation. Here we show with a model of intermediate complexity, that adding three oceanic mechanisms: brine induced stratification, stratification-dependant diffusion and iron fertilization to the standard glacial simulation (which includes sea level drop, temperature change, carbonate compensation and terrestrial carbon release) decreases CO2 down to the glacial value of ∼190 ppm and simultaneously matches glacial atmospheric and oceanic δ13C inferred from proxy data. LGM CO2 and δ13C can at last be successfully reconciled.

Stress field control of eruption dynamics at a rift volcano: Nyamuragira, D.R.Congo

Using the record of 30 flank eruptions over the last 110 years at Nyamuragira, we have tested the relationship between the eruption dynamics and the local stress field. There are two groups of eruptions based on their duration (< 80days >) that are also clustered in space and time. We find that the eruptions fed by dykes parallel to the East African Rift Valley have longer durations (and larger volumes) than those eruptions fed by dykes with other orientations. This is compatible with a model for compressible magma transported through an elastic-walled dyke in a differential stress field from an over-pressured reservoir (Woods et al., 2006). The observed pattern of eruptive fissures is consistent with a local stress field modified by a northwest-trending, right lateral slip fault that is part of the northern transfer zone of the Kivu Basin rift segment. We have also re-tested with new data the stochastic eruption models for Nyamuragira of Burt et al. (1994). The time-predictable, pressure-threshold model remains the best fit and is consistent with the typically observed declining rate of sulphur dioxide emission during the first few days of eruption with lava emission from a depressurising, closed, crustal reservoir. The 2.4-fold increase in long-term eruption rate that occurred after 1977 is confirmed in the new analysis. Since that change, the record has been dominated by short-duration eruptions fed by dykes perpendicular to the Rift. We suggest that the intrusion of a major dyke during the 1977 volcano-tectonic event at neighbouring Nyiragongo volcano inhibited subsequent dyke formation on the southern flanks of Nyamuragira and this may also have resulted in more dykes reaching the surface elsewhere. Thus that sudden change in output was a result of a changed stress field that forced more of the deep magma supply to the surface. Another volcano-tectonic event in 2002 may also have changed the magma output rate at Nyamuragira.

Direct Computation of Stochastic Flow in Reservoirs with Uncertain Parameters

A direct method is presented for determining the uncertainty in reservoir pressure, flow, and net present value (NPV) using the time-dependent, one phase, two- or three-dimensional equations of flow through a porous medium. The uncertainty in the solution is modelled as a probability distribution function and is computed from given statistical data for input parameters such as permeability. The method generates an expansion for the mean of the pressure about a deterministic solution to the system equations using a perturbation to the mean of the input parameters. Hierarchical equations that define approximations to the mean solution at each point and to the field covariance of the pressure are developed and solved numerically. The procedure is then used to find the statistics of the flow and the risked value of the field, defined by the NPV, for a given development scenario. This method involves only one (albeit complicated) solution of the equations and contrasts with the more usual Monte-Carlo approach where many such solutions are required. The procedure is applied easily to other physical systems modelled by linear or nonlinear partial differential equations with uncertain data.