Evolution of droplets in subsea oil and gas blowouts: Development and validation of the numerical model VDROP-J

The droplet size distribution of dispersed phase (oil and/or gas) in submerged buoyant jets was addressed in this work using a numerical model, VDROP-J. A brief literature review on jets and plumes allows the development of average equations for the change of jet velocity, dilution, and mixing energy as function of distance from the orifice. The model VDROP-J was then calibrated to jets emanating from orifices ranging in diameter, D, from 0.5 mm to 0.12 m, and in cross-section average jet velocity at the orifice ranging from 1.5 m/s to 27 m/s. The d₅₀/D obtained from the model (where d₅₀ is the volume median diameter of droplets) correlated very well with data, with an R² = 0.99. Finally, the VDROP-J model was used to predict the droplet size distribution from Deepwater Horizon blowouts. The droplet size distribution from the blowout is of great importance to the fate and transport of the spilled oil in marine environment.     

Focal depths and mechanisms of Tohoku-Oki aftershocks from teleseismic P wave modeling

Aftershocks of the 2011 Tohoku-Oki great earthquake have a wide range of focal depths and fault plane mechanisms. We constrain the focal depths and focal mechanisms of 69 aftershocks with M _w > 5.4 by modeling the waveforms of teleseismic P and its trailing near-surface reflections pP and sP. We find that the “thrust events” are within 10 km from the plate interface. The dip angles of these thrust events increase with depth from ~5° to ~25°. The “non-thrust events” vary from 60 km above to 40 km below the plate interface. Normal and strike-slip events within the overriding plate point to redistribution of stress following the primary great earthquake; however, due to the spatially variable stress change in the Tohoku-Oki earthquake, an understanding of how the mainshock affected the stresses that led to the aftershocks requires accurate knowledge of the aftershock location.     

Ambient noise as the new source for urban engineering seismology and earthquake engineering: a case study from Beijing metropolitan area

In highly populated urban centers, traditional seismic survey sources can no longer be properly applied due to restrictions in modern civilian life styles. The ambient vibration noise, including both microseisms and microtremor, though are generally weak but available anywhere and anytime, can be an ideal supplementary source for conducting seismic surveys for engineering seismology and earthquake engineering. This is fundamentally supported by advanced digital signal processing techniques for effectively extracting the useful information out from the noise. Thus, it can be essentially regarded as a passive seismic method. In this paper we first make a brief survey of the ambient vibration noise, followed by a quick summary of digital signal processing for passive seismic surveys. Then the applications of ambient noise in engineering seismology and earthquake engineering for urban settings are illustrated with examples from Beijing metropolitan area. For engineering seismology the example is the assessment of site effect in a large area via microtremor observations. For earthquake engineering the example is for structural characterization of a typical reinforced concrete high-rise building using background vibration noise. By showing these examples we argue that the ambient noise can be treated as a new source that is economical, practical, and particularly valuable to engineering seismology and earthquake engineering projects for seismic hazard mitigation in urban areas.     

Simulating the diurnal cycle of rainfall in global climate models: resolution versus parameterization

The effects of horizontal resolution and the treatment of convection on simulation of the diurnal cycle of precipitation during boreal summer are analyzed in several innovative weather and climate model integrations. The simulations include: season-long integrations of the Non-hydrostatic Icosahedral Atmospheric Model (NICAM) with explicit clouds and convection; year-long integrations of the operational Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts at three resolutions (125, 39 and 16 km); seasonal simulations of the same model at 10 km resolution; and seasonal simulations of the National Center for Atmospheric Research (NCAR) low-resolution climate model with and without an embedded two-dimensional cloud-resolving model in each grid box. NICAM with explicit convection simulates best the phase of the diurnal cycle, as well as many regional features such as rainfall triggered by advancing sea breezes or high topography. However, NICAM greatly overestimates mean rainfall and the magnitude of the diurnal cycle. Introduction of an embedded cloud model within the NCAR model significantly improves global statistics of the seasonal mean and diurnal cycle of rainfall, as well as many regional features. However, errors often remain larger than for the other higher-resolution models. Increasing resolution alone has little impact on the timing of daily rainfall in IFS with parameterized convection, yet the amplitude of the diurnal cycle does improve along with the representation of mean rainfall. Variations during the day in atmospheric prognostic fields appear quite similar among models, suggesting that the distinctive treatments of model physics account for the differences in representing the diurnal cycle of precipitation.     

Advances in geochemistry during the last four decades: A personal perspective

This is the author’s speech at the meeting in Cologne (2007) to celebrate the 40th anniversary of the International Association of Geochemistry and Cosmochemistry, which the author served as the President in 2000 to 2004. The paper narrates the author’s personal involvement in important scientific programs during the last 4 decades, including implementation of isotope techniques, oil-and-gas research, diamond research, deep-sea drilling, space research, molecular biology and the origin of life.     

A coquinoid tsunamite from the Pliocene of Salento (SE Italy)

A Pliocene composite shell concentration from a shallow-marine carbonate succession located in the small, slowly subsiding Novoli graben (Salento peninsula, Puglia, S-Italy) is interpreted as a tsunamite. This interpretation is based on taphonomic, palaeoecological and sedimentologic features. The deposit appears to drape the previous substrate, with thinning on the topographic highs and thickening in the lows, and consists of up to six sub-layers showing repeated grading (normal and reverse), suggesting that deposition occurred during up to six surges. The shell assemblage results from a chaotic mixing of taxa spanning a large spectrum of environments and taphonomic features. Although a fraction of broken shells is present, most shells are integral, suggesting a transport mechanism preventing strong impacts. On the other hand, the high comminution of skeletal remains occurring in the matrix of the deposit is thought to reflect the composition of sea-floor sediment before the inferred high-energy event.The predominance of displaced, but still articulated valves indicates that bivalves were mostly transported while still alive, after having been lifted from the sea bed. Exhumation, displacement, and large-scale transport of large-sized, deep-infaunal organisms (such as Panopea) implies that across relatively large areas the upper few decimetres of the middle–inner shelf floor were eroded by highly turbulent flows of exceptional intensity. A transport as turbulent suspension, very rapid deposition and immediate and permanent burial are strongly suggested by the general integrity of the shells indicating limited effects of physical damage and lack of biological alteration. The variable faunal characteristics of the unit lying below the inferred tsunamite suggest that deposition of the latter took place in settings ranging from a somewhat restricted, relatively low-energy bay, to a more open inner-shelf environment.     

Weakly bound water structure,bond valence saturation and water dynamics at the goethite (100) surface/aqueous interface: ab initio dynamical simulations

Background

Many important geochemical and biogeochemical reactions occur in the mineral/formation water interface of the highly abundant mineral, goethite [α-Fe(OOH)]. Ab initio molecular dynamics (AIMD) simulations of the goethite α-FeOOH (100) surface and the structure, water bond formation and dynamics of water molecules in the mineral/aqueous interface are presented. Several exchange correlation functionals were employed (PBE96, PBE96 + Grimme, and PBE0) in the simulations of a (3 × 2) goethite surface with 65 absorbed water molecules in a 3D-periodic supercell (a = 30 Å, FeOOH slab ~12 Å thick, solvation layer ~18 Å thick).

Results

The lowest energy goethite (100) surface termination model was determined to have an exposed surface Fe³⁺ that was loosely capped by a water molecule and a shared hydroxide with a neighboring surface Fe³⁺. The water molecules capping surface Fe³⁺ ions were found to be loosely bound at all DFT levels with and without Grimme corrections, indicative that each surface Fe³⁺ was coordinated with only five neighbors. These long bonds were supported by bond valence theory calculations, which showed that the bond valence of the surface Fe³⁺ was saturated and surface has a neutral charge. The polarization of the water layer adjacent to the surface was found to be small and affected only the nearest water. Analysis by density difference plots and localized Boys orbitals identified three types of water molecules: those loosely bound to the surface Fe³⁺, those hydrogen bonded to the surface hydroxyl, and bulk water with tetrahedral coordination. Boys orbital analysis showed that the spin down lone pair orbital of the weakly absorbed water interact more strongly with the spin up Fe³⁺ ion. These weakly bound surface water molecules were found to rapidly exchange with the second water layer (~0.025 exchanges/ps) using a dissociative mechanism.

Conclusions

Water molecules adjacent to the surface were found to only weakly interact with the surface and as a result were readily able to exchange with the bulk water. To account for the large surface Fe–OH₂ distances in the DFT calculations it was proposed that the surface Fe³⁺ atoms, which already have their bond valence fully satisfied with only five neighbors, are under-coordinated with respect to the bulk coordination.

Open image in new window

Graphical abstract All first principle calculations, at all practically achievable levels, for the goethite 100 aqueous interface support a long bond and weak interaction between the exposed surface Fe³⁺ and water molecules capping the surface. This result is supported by bond valence theory calculations and is indicative that each surface Fe³⁺ is coordinated with only 5 neighbors.

     

Cambrian–early Ordovician volcanism across the South Armorican and Occitan domains of the Variscan Belt in France: Continental break-up and rifting of the northern Gondwana margin

The Cambrian–lower Ordovician volcanic units of the South Armorican and Occitan domains are analysed in a tectonostratigraphic survey of the French Variscan Belt. The South Armorican lavas consist of continental tholeiites in middle Cambrian–Furongian sequences related to continental break-up. A significant volcanic activity occurred in the Tremadocian, dominated by crustal melted rhyolitic lavas and initial rifting tholeiites. The Occitan lavas are distributed into five volcanic phases: (1) basal Cambrian rhyolites, (2) upper lower Cambrian Mg-rich tholeiites close to N-MORBs but crustal contaminated, (3) upper lower–middle Cambrian continental tholeiites, (4) Tremadocian rhyolites, and (5) upper lower Ordovician initial rift tholeiites. A rifting event linked to asthenosphere upwelling took place in the late early Cambrian but did not evolve. It renewed in the Tremadocian with abundant crustal melting due to underplating of mixed asthenospheric and lithospheric magmas. This main tectono-magmatic continental rift is termed the “Tremadocian Tectonic Belt” underlined by a chain of rhyolitic volcanoes from Occitan and South Armorican domains to Central Iberia. It evolved with the setting of syn-rift coarse siliciclastic deposits overlain by post-rift deep water shales in a suite of sedimentary basins that forecasted the South Armorican–Medio-European Ocean as a part of the Palaeotethys Ocean.     

Local-scale characteristics of high-marsh communities next to developed and undeveloped shorelines in an ocean-dominated estuary, Murrells Inlet, SC

Anthropogenic alteration of terrestrial shorelines can have pronounced effects on marine environments at the upland-marsh boundary. Possible terrestrial development effects on several physical and biological variables of high-marsh habitats were examined along developed and undeveloped shorelines in an ocean-dominated, southeastern US estuary. Analyses of sediment characteristics suggested development of the upland boundary affected physical conditions within the high-marsh. For example, pore water salinities were greater along undeveloped shorelines during a non-drought period even after rain events. Significant floral and faunal differences also existed between shoreline treatments. Black needle rush stems were significantly taller and marsh periwinkle densities significantly greater, but eastern coffee bean snail densities were significantly reduced along developed shorelines. Benthic infaunal community abundance and composition also were significantly different between shoreline treatments with sand fly larvae, human pest precursors, either only present or present in greater densities along developed shorelines. Sediment respirometry experiments indicated significant differences in heterotrophic and autotrophic processes occurring between shoreline treatments. Greater sediment surface temperatures along developed shorelines provided one possible mechanism driving high-marsh responses to boundary alteration. The history and extent of shoreline development along with a tendency in ocean-dominated southeastern marshes to resist change likely influenced current ecological conditions within our high-marsh study areas. A greater understanding of the driving mechanisms producing localized effects on salt marshes and recognizing regional differences in marsh resistance to change will facilitate predictions of shoreline development consequences and help in proposing effective management strategies for coastal boundaries.