A model for oxygen and sulfur isotope fractionation in sulfate during bacterial sulfate reduction processes

We present a model of bacterial sulfate reduction that includes equations describing the fractionation relationship between the sulfur and the oxygen isotope composition of residual sulfate (δ³⁴S_{SO4_residual}, δ¹⁸O_{SO4_residual}) and the amount of residual sulfate. The model is based exclusively on oxygen isotope exchange between cell-internal sulfur compounds and ambient water as the dominating mechanism controlling oxygen isotope fractionation processes. We show that our model explains δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} patterns observed from natural environments and from laboratory experiments, whereas other models, favoring kinetic isotope fractionation processes as dominant process, fail to explain many (but not all) observed δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} patterns. Moreover, we show that a “typical” δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} slope does not exist. We postulate that measurements of δ³⁴S_{SO4_residual} and δ¹⁸O_{SO4_residual} can be used as a tool to determine cell-specific sulfate reduction rates, oxygen isotope exchange rates, and equilibrium oxygen isotope exchange factors. Data from culture experiments are used to determine the range of sulfur isotope fractionation factors in which a simplified set of equations can be used. Numerical examples demonstrate the application of the equations. We postulate that, during denitrification, the oxygen isotope effects in residual nitrate are also the result of oxygen isotope exchange with ambient water. Consequently, the equations for the relationship between δ³⁴S_{SO4_residual}, δ¹⁸O_{SO4_residual}, and the amount of residual sulfate could be modified and used to calculate the fractionation-relationship between δ¹⁵N_{NO3_residual}, δ¹⁸O_{NO3_residual}, and the amount of residual nitrate during denitrification.     

Longitudinal and depth variation of bacterioplankton productivity and related factors in a temperate estuary

Bacterioplankton productivity (BP) spatial variation was investigated in relation to potential resources, including primary productivity and dissolved organic matter, in the micro-tidal Neuse River–Pamlico Sound estuarine system, North Carolina, USA. Estuarine BP was predicted to correlate with the trophic gradient, decreasing along the salinity gradient in parallel with the decrease in organic matter and primary productivity. This prediction was tested over four years at spatial scales ranging from kilometers to meters along the riverine axis and with depth. The general pattern of BP across the salinity gradient was unimodal and matched the phytoplankton patterns in peak location and variability. Peak locations varied with discharge, especially in 2003 when above average discharge moved peaks downstream. Spatial coherence of BP with other variables was much less at short time scales. The effect of temperature, nutrients, and phytoplankton on BP varied by location, especially fresh versus brackish stations, although only temperature explained more than 20% of the BP variation. Depth variation of BP was as great as longitudinal variation and bottom samples were often higher than surface. BP was strongly correlated with particulate organic carbon at the pycnocline and bottom, highlighting the importance of particulate matter as a resource. Station-averaged BP and phytoplankton data corresponded well with two published meta-analyses, although the offset of the freshwater station suggested longitudinal differences in community composition or resource availability.     

Climate feedbacks determined using radiative kernels in a multi-thousand member ensemble of AOGCMs

The use of radiative kernels to diagnose climate feedbacks is a recent development that may be applied to existing climate change simulations. We apply the radiative kernel technique to transient simulations from a multi-thousand member perturbed physics ensemble of coupled atmosphere-ocean general circulation models, comparing distributions of model feedbacks with those taken from the CMIP-3 multi GCM ensemble. Although the range of clear sky longwave feedbacks in the perturbed physics ensemble is similar to that seen in the multi-GCM ensemble, the kernel technique underestimates the net clear-sky feedbacks (or the radiative forcing) in some perturbed models with significantly altered humidity distributions. In addition, the compensating relationship between global mean atmospheric lapse rate feedback and water vapor feedback is found to hold in the perturbed physics ensemble, but large differences in relative humidity distributions in the ensemble prevent the compensation from holding at a regional scale. Both ensembles show a similar range of response of global mean net cloud feedback, but the mean of the perturbed physics ensemble is shifted towards more positive values such that none of the perturbed models exhibit a net negative cloud feedback. The perturbed physics ensemble contains fewer models with strong negative shortwave cloud feedbacks and has stronger compensating positive longwave feedbacks. A principal component analysis used to identify dominant modes of feedback variation reveals that the perturbed physics ensemble produces very different modes of climate response to the multi-model ensemble, suggesting that one may not be used as an analog for the other in estimates of uncertainty in future response. Whereas in the multi-model ensemble, the first order variation in cloud feedbacks shows compensation between longwave and shortwave components, in the perturbed physics ensemble the shortwave feedbacks are uncompensated, possibly explaining the larger range of climate sensitivities observed in the perturbed simulations. Regression analysis suggests that the parameters governing cloud formation, convection strength and ice fall speed are the most significant in altering climate feedbacks. Perturbations of oceanic and sulfur cycle parameters have relatively little effect on the atmospheric feedbacks diagnosed by the kernel technique.     

Quantifying internal variability in a regional climate model: a case study for Southern Africa

This study aims at presenting various methodologies to separate the reproducible and irreproducible components of seasonal and intraseasonal climate variability simulated by a regional climate model over Southern Africa (south of 15°S), during an austral summer rainy season representative of the climatology. To that end, a 30-member ensemble simulation is performed using WRF laterally forced by the ERA40 reanalyses. Retained metrics include the analysis of weather regimes, signal-to-noise ratio, inter-member standard deviation and coefficient of variation. At the seasonal timescale, simulated rainfall amounts generally show a strong reproducibility, except in the subtropics and over the southern part of the Mozambique Channel. There, the number of rainy days is roughly similar in all members, while their average intensity varies extensively. At the intraseasonal timescale, the chronology of weather regimes, derived from the 500?hPa geopotential height, is highly reproducible. Rainfall variability is much less reproducible, especially in the central parts of the domain and near its outflow boundaries. Analysis of a South African regional index nonetheless indicates that both wet and dry spells tend to be accurately simulated and occur in phase in most members, demonstrating that they are embedded in large-scale variability patterns. Internal variability is lastly related to the lateral forcings along the domain boundaries. An objective classification of inflow/outflow mass fluxes allows identification of the recurrent synoptic configurations that favor strong or weak regional reproducibility. The main uncertainties concern the basic morphological features of rain-bearing systems (i.e., their spatial extension, location and propagation speed). Consequences for tropical-temperate interactions are then discussed.     

Distinguishing volcanic debris avalanche deposits from their reworked products: the Perrier sequence (French Massif Central)

Debris avalanches associated with volcanic sector collapse are usually high-volume high-mobility phenomena. Debris avalanche deposit remobilisation by cohesive debris flows and landslides is common, so they can share textural characteristics such as hummocks and jigsaw cracks. Distinguishing original deposits from reworked products is critical for geological understanding and hazard assessment because of their different origin, frequency and environmental impact. We present a methodology based on field evidence to differentiate such epiclastic breccias. Basal contact mapping constrained by accurate altitude and location data allows the reconstruction of deposit stratigraphy and geometry. Lithological analysis helps to distinguish the different units. Incorporation structures, kinematic indicators and component mingling textures are used to characterise erosion and transport mechanisms. We apply this method to the enigmatic sequence at Perrier (French Massif Central), where four units (U1–U4) have been interpreted either as debris flow or debris avalanche deposits. The sequence results from activity on the Monts Dore Volcano about 2 Ma ago. The epiclastic units are matrix supported with an almost flat top. U2 and U3 have clear debris flow deposit affinities such as rounded clasts and intact blocks (no jigsaw cracks). U1 and U4 have jigsaw cracked blocks with matrix injection and stretched sediment blocks. U1 lacks large blocks (>10 m wide) and has a homogenous matrix with an upward increase of trapped air vesicle content and size. This unit is interpreted as a cohesive debris flow deposit spawned from a debris avalanche upstream. In contrast, U4 has large mega-blocks (up to 40 m wide), sharp contacts between mixed facies zones with different colours and numerous jigsaw fit blocks (open jigsaw cracks filled by monogenic intra-clast matrix). Mega-blocks are concentrated near the deposit base and are spatially associated with major substratum erosion. This deposit has a debris avalanche distal facies with local debris flow affinities due to partial water saturation. We also identify two landslide deposits (L1 and L2) resulting from recent reworking that has produced a similar facies to U1 and U4. These are distinguishable from the original deposits, as they contain blocks of mixed U1/U4 facies, a distinctly less consolidated and more porous matrix and a fresh hummocky topography. This work shows how to differentiate epiclastic deposits with similar characteristics, but different origins. In doing so, we improve understanding of present and past instability of the Monts Dore and identify present landslide hazards at Perrier.     

Climate during the last glacial maximum in the Wasatch and southern Uinta Mountains inferred from glacier modeling

Recent improvements in understanding glacial extents and chronologies in the Wasatch and Uinta Mountains and other mountain ranges in the western U.S. call for a more detailed approach to using glacier reconstructions to infer paleoclimates than commonly applied AAR-ELA-ÄT methods. A coupled 2-D mass balance and ice-flow numerical modeling approach developed by [Plummer, M.A., Phillips, F.M., 2003. A 2-D numerical model of snow/ice energy balance and ice flow for paleoclimatic interpretation of glacial geomorphic features. Quaternary Science Reviews 22, 1389–1406] allows exploration of the combined effects of temperature, precipitation, shortwave radiation and many secondary parameters on past ice extents in alpine settings. We apply this approach to the Little Cottonwood Canyon in the Wasatch Mountains and the Lake Fork and Yellowstone Canyons in the south-central Uinta Mountains. Results of modeling experiments indicate that the Little Cottonwood glacier required more precipitation during the local Last Glacial Maximum (LGM) than glaciers in the Uinta Mountains, assuming lapse rates were similar to modern. Model results suggest that if temperatures in the Wasatch Mountains and Uinta Mountains were  6 °C to 7 °C colder than modern, corresponding precipitation changes were  3 to 2× modern in Little Cottonwood Canyon and  2 to 1× modern in Lake Fork and Yellowstone Canyons. Greater amounts of precipitation in the Little Cottonwood Canyon likely reflect moisture derived from the surface of Lake Bonneville, and the lake may have also affected the mass balance of glaciers in the Uinta Mountains.     

Submarine channel levee shape and sediment waves from physical experiments

Submarine channel levees aggrade through repeated overspill events from the channel axis. The shape of the levees may therefore reflect some characteristic(s) of the overspilling flow. It has been noted that basin floor levees typically have a relatively low-relief and taper exponentially to their termination; in contrast slope channel levees may be much steeper close to the channel. A simple physical experiment was performed where a surge-like sediment-laden current flowed through a curved channel. Significant overspill occurred and generated a deposit flanking the channel on either side. The experiment was repeated 25 times to build up low-relief channel-levees. It was found that in proximal areas, levees were steep and characterised by power-law decays, a transitional zone of logarithmically thinning levee was found a little further down-channel, followed by exponential decays in medial to distal areas. The style of levee decay is a function of spatial variation in overbank sedimentation rates. Where flows rapidly lose momentum and deposit across the grain-size spectrum, i.e., in proximal areas, levees tend to be steep; farther down the channel, the steep levee slope gives way to a more gradually tapering deposit. In more distal parts of the channel, deposition is directly related to sediment settling velocity (rather than the suspended load exceeding flow transport capacity as is the case in proximal areas), the deposit reflects this with relatively simple exponential thickness decays. Additionally, small-scale sediment waves developed under lee wave conditions on the inner-bend overbank. The waves initially migrated slightly towards the channel, but as the style of overspill evolved due to intra-channel deposition, flows moved out of the lee wave window and sedimentation became out of phase with the wavelength of the features and the topography was healed.     

3D seismic analysis of the gas hydrate system and the fluid migration paths in part of the Niger Delta Basin,Nigeria

Comprehensive qualitative and semi-quantitative seismic analysis was carried out on 3-dimensional seismic data acquired in the deepwater compressional and shale diapiric zone of the Niger Delta Basin using an advanced seismic imaging tool. The main aim of this work is to obtain an understanding of the forming mechanism of the gas hydrate system, and the fluid migration paths associated with this part of the basin. The results showed the presence of pockmarks on the seafloor and bottom simulating reflectors (BSRs) in the field, indicating the active fluid flux and existence of gas hydrate system in the area. In the area of approximately 195 km² occupying nearly 24% of the entire study field, three major zones with continuous or discontinuous BSRs of 3 to 7 km in length which are in the northeastern, southern and eastern part of the field respectively were delineated. The BSR is interpreted to be the transition between the free gas zone and the gas hydrate zone. The geologic structures including faults (strike-slip and normal faults), chimneys and diapirs were deduced to be the main conduits for gas migration. It is concluded that the biogenic gases generated in the basin were possibly transported via faults and chimneys by advection processes and subsequently accumulated under low temperature and high pressure conditions in the free gas zone below the BSR forming gas hydrate. A plausible explanation for the presence of the ubiquitous pockmarks of different diameters and sizes in the area is the transportation of the excessive gas to the seafloor through these mapped geologic structures.     

Homogeneous nucleation of amorphous solid water particles in the upper mesosphere

Condensed water particles are known to exist in the high latitude upper mesosphere during the summer months. However, the mechanism or mechanisms through which they nucleate remains uncertain. It is postulated here that particles of amorphous solid water (ASW, condensed water with a non-crystalline structure) may nucleate homogeneously in the summer mesosphere. Using classical nucleation theory and a one-dimensional model, it is shown that more than 10⁵ cm^?3 amorphous solid water particles can nucleate homogeneously under mesopause conditions. Furthermore, it is shown that homogeneous nucleation competes with heterogeneous nucleation on meteoric smoke particles when the cooling rate is >0.5 K/h. The homogeneous nucleation of amorphous solid water could provide an explanation for the high density of ice particles (many thousands per cm³) thought to be required for electron depletions in the upper mesosphere. A parameterisation for homogeneous nucleation is presented which can be used in other mesospheric cloud models.