ALBERT

Description: Author Posting. © Cambridge University Press, 2000. This article is posted here by permission of Cambridge University Press for personal use, not for redistribution. The definitive version was published Journal of the Marine Biological Association of the UK 80 (2000): 827-834, doi:10.1017/S0025315400002800.

Description: The shell and internal anatomy of the montacutid bivalve Mysella verrilli is described for the first time. The species is remarkable in that the oesophagus has developed into a suctorial proboscis. This has been accompanied by the loss of the palps. In addition the gonads have been extended from the dorsal part of the body to form two gill-like extensions to which the reduced inner demibranchs attach along the postero–ventral margin. Mysella verrilli broods its young in the mantle cavity to the late veliger stage before releasing them. It is believed that the species is probably a suctorial ectoparasite on a soft-bodied benthic invertebrate.

Description: Author Posting. © Cambridge University Press, 2006. This article is posted here by permission of Cambridge University Press for personal use, not for redistribution. The definitive version was published in Journal of Fluid Mechanics 564 (2006): 435-454, doi:10.1017/S0022112006001522.

Description: Motivated by work on tilted convection (Sheremet, J. Fluid Mech., vol. 506, 2004, p. 217), a set of experiments is presented here using the same set-up of a tilted tank attached to a rotating centrifuge with a 2.5 m arm. Within the tank small, almost neutrally buoyant, spheres are released, and their trajectories are recorded. Thus the forces acting on a sphere can be analysed in the case of misalignment between the buoyancy force and the axis of rotation. The angles of descent characterizing the trajectory are compared with inviscid linear theory developed by Stewartson (Q. J. Math. Appl. Mech., vol. 6, 1953, p. 141), and the agreement is found to be good. The angles should be independent of the density anomaly of the spheres compared to their environment. Using the descent velocity from non-rotating experiments, the density of the spheres is estimated and used to determine the drag acting on them in the rotating experiments. It is found that the drag is up to 50% larger than expected from Stewartson's theory. The agreement is best, not for infinitesimal, but for small Rossby numbers. The results are consistent with observations recorded by Maxworthy (J. Fluid Mech., vol. 40, 1970, p. 453).

Description: Author Posting. © Cambridge University Press, 2006. This article is posted here by permission of Cambridge University Press for personal use, not for redistribution. The definitive version was published in Journal of Fluid Mechanics 561 (2006):103–112, doi:10.1017/S0022112006000991

Description: A self-consistent formalism to estimate baroclinic energy densities and fluxes resulting from the propagation of internal waves of arbitrary amplitude is derived using the concept of available potential energy. The method can be applied to numerical, laboratory or field data.

Description: MBIWE98 was supported by the US Geological Survey and the Office of Naval Research. A.S. received support from the Office of Naval Research (N00014-05-1-0361), R.B. from the Walter A. and Hope Noyes Smith Chair on Coastal Oceanography and B.B. from the US Geological Survey.

Description: The article addresses old ‘west-east discourses’ and how they continue to develop in the high north, and, not least, in the Norwegian petroleum debate. Adopting a discourse analytical perspective the author shows how environmental safety is used as an argument in favour of Norway producing oil in the Barents Sea at the earliest possible moment. This is only feasible if a connection is made in the public mind between Russia and the environment. These views, it is argued, stem from ideas about Russia that gained currency after the demise of the Soviet Union. While they perhaps have less to do with Russia's petroleum industry and environmental performance today, they nevertheless have a strong impact on how challenges in the high north and Arctic region are perceived. And, perhaps even more importantly, they define freedom of action and available options. In this paper all references to government departments refer to those of Norway.

Description: A study of the morphology of the vortical skeleton behind a flapping NACA0030 wing with a finite aspect ratio of 3, is undertaken. The motivation for this work originates with the proposal that thrust can be efficiently produced by flapping aerofoils. The test condition corresponds to a Strouhal number of 0.35, Reynolds number, based on aerofoil chord, of 600 and an amplitude of flapping, equal to the chord length of the wing. This test condition corresponds to the optimal thrust-producing case in infinite-span flapping wings. This study investigates the effect of wing three-dimensionality on the structure of the wake-flow. This is accomplished here, by quantitatively describing the spatio-temporal variations in the velocity, vorticity and Reynolds stresses for the finite-span-wing case.Preliminary flow visualizations suggest that the presence of wingtip vortices for the three-dimensional-wing case, create a different vortical structure to the two-dimensional-wing case. In the case of a two-dimensional-wing, the flow is characterized by the interaction of leading- and trailing-edge vorticity, resulting in the formation of a clear reverse Kármán vortex street at the selected test condition. In the case of a three-dimensional-wing, the flow exhibits a high degree of complexity and three-dimensionality, particularly in the midspan region. Using phase-averaged particle image velocimetry measurements of the forced oscillatory flow, a quantitative analysis in the plane of symmetry of the flapping aerofoil was undertaken. Using a triple decomposition of the measured velocities, the morphological characteristics of the spanwise vorticity is found to be phase correlated with the aerofoil kinematics. Reynolds stresses in the direction of oscillation are the dominant dissipative mechanism. The mean velocity profiles resemble ajet, indicative of thrust production. Pairs of strong counter-rotating vortices from the leading- and trailing-edge of the aerofoil are shed into the flow at each half-cycle. The large-scale structure of the flow is characterized by constructive merging of spanwise vorticity. The midspan region is populated by cross-sections of interconnected vortex rings.

Description: The effectiveness of international environmental cooperation is a theoretical and methodological challenge that has captured the interest of scholars and students of international relations in the recent years. This article is an evaluation of international environmental cooperation in northwest Russia by applying an approach developed in the Institutional Dimensions of Global Environmental Change (IDGEC) research programme of the International Human Dimension Programme (IHDP). The IDGEC approach emphasises the importance of complex diagnostics of institutional environmental arrangements in terms of performance. The approach relates to the effects of international environmental cooperation in terms of sustainability, efficiency, fairness, and robustness. The article focuses on Russian and Nordic experiences of international environmental cooperation, and the assessment of performance has been done by the participants themselves. According to the results of interviews and questionnaires, the cooperation is considered institutionally effective and robust in general. However, views about the environmental effects vary. Some problems exist that are related to sustainability and fairness, but the most critical issue is the lack of domestic environmental capacity in Russia.

Description: In this study, we performed simulations of turbulent flow over rectangular ribs transversely mounted on one side of a plane in a channel, with the other side being smooth. The separation between ribs is large enough to avoid forming stable vortices in the spacing, which exhibitsk-type, or sand-grain roughness. The Reynolds numberReτof our representative direct numerical simulation case is 460 based on the smooth-wall friction velocity and the channel half-width. The roughness heighthis estimated as 110 wall units based on the rough-wall friction velocity. The velocity profile and kinetic energy budget verify the presence of an equilibrium, logarithmic layer aty≳2h. In the roughness sublayer, however, a significant turbulent energy flux was observed. A high-energy region is formed by the irregular motions just above the roughness. Visualizations of vortical streaks, disrupted in all three directions in the roughness sublayer, indicate that the three-dimensional flow structure of sand-grain roughness is replicated by the two-dimensional roughness, and that this vortical structure is responsible for the high energy production. The difference in turbulence structure between smooth- and rough-wall layers can also be seen in other flow properties, such as anisotropy and turbulence length scales.

Description: A model for simulating the process of growth, collapse and rebound of a cavitation bubble travelling along the flow through a convergent–divergent nozzle producing a cavitating water jet is established. The model is based on the Rayleigh–Plesset bubble dynamics equation using as inputs ambient pressure and velocity profiles calculated with the aid of computational fluid dynamics (CFD) flow modelling. A variable time-step technique is applied to solve the highly nonlinear second-order differential equation. This technique successfully solves the Rayleigh–Plesset equation for wide ranges of pressure variation and bubble original size and saves considerable computing time. Inputs for this model are the pressure and velocity data from CFD calculation. To simulate accurately the process of bubble growth, collapse and rebound, a heat transfer model, which includes the effects of conduction plus radiation, is developed to describe the thermodynamics of the incondensable gas inside the bubble. This heat transfer model matches previously published experimental data well. Assuming that single bubble behaviour also applies to bubble clouds, the calculated distance from the nozzle exit travelled by the bubble to the point where the bubble size becomes invisible is taken to be equal to the bubble cloud length observed. The predictions are compared with experiments carried out in a cavitation cell and show good agreement for different nozzles operating at different pressure conditions.

Description: An experimental study, investigating the mean flow and turbulence in the wind drift layer formed beneath short wind waves was conducted. The degree to which these flows resemble the flows that occur in boundary layers adjacent to solid walls (i.e. wall-layers) was examined. Simultaneous DPIV (digital particle image velocimetry) and infrared imagery were used to investigate these near-surface flows at a fetch of 5.5 m and wind speeds from 4.5 to 11 m s−1. These conditions produced short steep waves with dominant wavelengths from 6 cm to 18 cm. The mean velocity profiles in the wind drift layer were found to be logarithmic and the flow was hydrodynamically smooth at all wind speeds. The rate of dissipation of turbulent kinetic energy was determined to be significantly greater in magnitude than would occur in a comparable wall-layer. Microscale breaking waves were detected using the DPIV data and the characteristics of breaking and non-breaking waves were compared. The percentage of microscale breaking waves increased abruptly from 11% to 80% as the wind speed increased from 4.5 to 7.4 m s−and then gradually increased to 90% as the wind speed increased to 11 m s−. At a depth of 1 mm, the rate of dissipation was 1.7 to 3.2 times greater beneath microscale breaking waves compared to non-breaking waves. In the crest–trough region beneath microscale breaking waves, 40% to 50% of the dissipation was associated with wave breaking. These results demonstrated that the enhanced near-surface turbulence in the wind drift layer was the result of microscale wave breaking. It was determined that the rate of dissipation of turbulent kinetic energy due to wave breaking is a function of depth, friction velocity, wave height and phase speed as proposed by Terrayet al. (1996). Vertical profiles of the rate of dissipation showed that beneath microscale breaking waves there were two distinct layers. Immediately beneath the surface, the dissipation decayed as ζ−0.7and below this in the second layer it decayed as ζ−2. The enhanced turbulence associated with microscale wave breaking was found to extend to a depth of approximately one significant wave height. The only similarity between the flows in these wind drift layers and wall-layers is that in both cases the mean velocity profiles are logarithmic. The fact that microscale breaking waves were responsible for 40%–50% of the near-surface turbulence supports the premise that microscale breaking waves play a significant role in enhancing the transfer of gas and heat across the air–sea interface.

Description: Results are presented from experimental investigations into the motion of a heavy ellipsoid in a horizontal rotating cylinder, which has been completely filled with highly viscous fluid. The motion can be conveniently classified using the ratio between the maximum radius of curvature of the ellipsoid κmaxand the radius of the drumRd. If κmax

Description: Electric conduction from an electrolyte solution into a charge selective solid, such as ion exchange membrane or electrode, becomes unstable when the electrolyte concentration near the interface approaches zero owing to diffusion limitation. The sequence of events leading to instability is as follows: upon the decrease of the interface concentration, the electric double layer at the interface transforms from its common quasi-equilibrium structure to a different, non-equilibrium one. The key feature of this new structure is an extended space charge added to the usual one of the quasi-equilibrium electric double layer. The non-equilibrium electro-osmotic slip related to this extended space charge renders the quiescent conductance unstable. A unified asymptotic picture of the electric double-layer undercurrent, encompassing all regimes from quasi-equilibrium to the extreme non-equilibrium one, is developed and employed for derivation of a universal electro-osmotic slip formula. This formula is used for a linear stability study of quiescent electric conduction, yielding the precise parameter range of instability, compared with that in the full electroconvective formulation. The physical mechanism of instability is traced both kinematically, in terms of non-equilibrium electro-osmotic slip, and dynamically, in terms of forces acting in the electric double layer.

Description: The water-shipping problem is modelled in a two-dimensional framework and studied experimentally and numerically for the case of a fixed barge-shaped structure. The analysis represents the second step of the research discussed in Grecoet al. (J. Fluid Mech., vol. 525, 2005, p. 309). The numerical investigation is performed by using both a boundary element method and a domain-decomposition strategy. The model tests highlight the occurrence of dam-breaking-type water on deck, (a) with and (b) without an initial plunging phase, and (c) an unusual type of water shipping connected with blunt water–deck impacts here called a hammer-fist type event never documented before. Cases (a) and (c) are connected with the most severe events and the related features and green-water loads are discussed in detail. A parametric analysis of water-on-deck phenomena has also been carried out in terms of the local incoming waves and bow flow features. We classify such phenomena in a systematic way to provide a basis for further investigations of water-on-deck events. The severity of (a)-type water-on-deck events is analysed in terms of initial cavity area and water-front velocity along the deck. The former increases as the square power of the modified incoming-wave (front-crest) steepness while the latter scales with its square-root. The two-dimensional investigation gives useful quantitative information in terms of water-front velocity for comparison with three-dimensional water-on-deck experiments on fixed bow models interacting with wave packets.

Description: The irrotational flow past two slender bodies of revolution at angles of yaw, translating in parallel paths in very close proximity, is analysed by extending the classical slender body theory. The flow far away from the two bodies is shown to be a direct problem, which is represented in terms of two line sources along their longitudinal axes, at the strengths of the variation rates of their cross-section areas. The inner flow near the two bodies is reduced to the plane flow problem of the expanding (contracting) and lateral translations of two parallel circular cylinders with different radii, which is then solved analytically using conformal mapping. Consequently, an analytical flow solution has been obtained for two arbitrary slender bodies of revolution at angles of yaw translating in close proximity. The lateral forces and yaw moments acting on the two bodies are obtained in terms of integrals along the body lengths. A comparison is made among the present model for two slender bodies in close proximity, Tuck & Newman's (1974) model for two slender bodies far apart, and VSAERO (AMI)–commercial software based on potential flow theory and the boundary element method (BEM). The attraction force of the present model agrees well with the BEM result, when the clearance,h0, is within 20% of the body length, whereas the attraction force of Tuck & Newman is much smaller than the BEM result whenh0is within 30% of the body length, but approaches the latter whenh0is about half the body length. Numerical simulations are performed for the three typical manoeuvres of two bodies: (i) a body passing a stationary body, (ii) two bodies in a meeting manoeuvre (translating in opposite directions), and (iii) two bodies in a passing manoeuvre (translating in the same direction). The analysis reveals the orders of the lateral forces and yaw moments, as well as their variation trends in terms of the manoeuvre type, velocities, sizes, angles of yaw of the two bodies, and their proximity, etc. These irrotational dynamic features are expected to provide a basic understanding of this problem and will be beneficial to further numerical and experimental studies involving additional physical effects.

Description: The flow associated with a synthetic jet transitioning to turbulence in an otherwise quiescent external flow is examined using time-accurate three-dimensional numerical simulations. The incompressible Navier–Stokes solver uses a second-order accurate scheme for spatial discretization and a second-order semi-implicit fractional step method for time integration. The simulations are designed to model the experiments of C. S. Yaoet al. (Proc. NASA LaRC Workshop, 2004) which have examined, in detail, the external evolution of a transitional synthetic jet in quiescent flow. Although the jet Reynolds and Stokes numbers in the simulations match with the experiment, a number of simplifications have been made in the synthetic jet actuator model adopted in the current simulations. These include a simpler representation of the cavity and slot geometry and diaphragm placement. Despite this, a reasonably good match with the experiments is obtained in the core of the jet and this indicates that for these jets, matching of these key non-dimensional parameters is sufficient to capture the critical features of the external jet flow. The computed results are analysed further to gain insight into the dynamics of the external as well as internal flow. The results indicate that near the jet exit plane, the flow field is dominated by the formation of counter-rotating spanwise vortex pairs that break down owing to the rapid growth of spanwise instabilities and transition to turbulence a short distance from the slot. Detailed analyses of the unsteady characteristics of the flow inside the jet cavity and slot provide insights that to date have not been available from experiments.

Description: A general proof that more energy flows upscale than downscale in two-dimensional turbulence and barotropic quasi-geostrophic (QG) turbulence is given. A proof is also given that in surface QG turbulence, the reverse is true. Though some of these results are known in restricted cases, the proofs given here are pedagogically simpler, require fewer assumptions and apply to both forced and unforced cases.

Description: The local linear stability of forced, stationary long waves produced by topography or potential vorticity (PV) sources is examined using a quasi-geostrophic barotropic model. A multiple scale analysis yields coupled equations for the background stationary wave and low-frequency (LF) disturbance field. Forcing structures for which the LF dynamics are Hamiltonian are shown to yield conservation laws that provide necessary conditions for instability and a constraint on the LF structures that can develop. Explicit knowledge of the forcings that produce the stationary waves is shown to be crucial to predicting a unique LF field. Various topographies or external PV sources can be chosen to produce stationary waves that differ by asymptotically small amounts, yet the LF instabilities that develop can have fundamentally different structures and growth rates. If the stationary wave field is forced solely by topography, LF oscillatory modes always emerge. In contrast, if the stationary wave field is forced solely by PV, two LF structures are possible: oscillatory modes or non-oscillatory envelope modes. The development of the envelope modes within the context of a linear LF theory is novel.An analysis of the complex WKB branch points, which yields an analytical expres-sion for the leading-order eigenfrequency, shows that the streamwise distribution of absolute instability and convective growth is central to understanding and predicting the types of LF structures that develop on the forced stationary wave. The location of the absolute instability region with respect to the stationary wave determines whether oscillatory modes or envelope modes develop. In the absence of absolute instability, eastward propagating wavetrains generated in the far field can amplify via local convective growth in the stationary wave region. If the stationary wave region is streamwise symmetric (asymmetric), the local convective growth results in a local change in wave energy that is transient (permanent).

Description: This paper presents combined theoretical and experimental studies of the two-dimensional piston-like steady-state motions of a fluid in a moonpool formed by two rectangular hulls (e.g. a dual pontoon or catamaran). Vertical harmonic excitation of the partly submerged structure in calm water is assumed. A high-precision analytically oriented linear-potential-flow method, which captures the singular behaviour of the velocity potential at the corner points of the rectangular structure, is developed. The linear steady-state results are compared with new experimental data and show generally satisfactory agreement. The influence of vortex shedding has been evaluated by using the local discrete-vortex method of Graham (1980). It was shown to be small. Thus, the discrepancy between the theory and experiment may be related to the free-surface nonlinearity.

Description: This paper presents a study of the global response of a fluid to impulsive and localized forcing; it has been motivated by the recent laboratory experiments on the locomotion of water-walking insects reported in Hu, Chan & Bush (Nature, vol. 424, 2003, p. 663). These insects create both waves and vortices by their rapid leg strokes and it has been a matter of some debate whether either form of motion predominates in the momentum budget. The main result of this paper is to argue that generically both waves and vortices are significant, and that in linear theory they take up the horizontal momentum with share 1/3 and 2/3, respectively.This generic result, which depends only on the impulsive and localized nature of the forcing, is established using the classical linear impulse theory, with adaptations to weakly compressible flows and flows with a free surface. Additional general comments on experimental techniques for momentum measurement and on the wave emission are given and then the theory is applied in detail to water-walking insects.Owing to its generality, this kind of result and the methods used to derive it should be applicable to a wider range of wave–vortex problems in the biolocomotion of water-walking animals and elsewhere.

Description: The Regenerative Enclosed Life Support Module Simulator (REMS) was designed to simulate the conditions aboard the International Space Station (ISS). This unique terrestrial, encapsulated environment for humans and their associated organisms allowed investigations into the microbial communities within an enclosed habitat system, primarily with respect to diversity, phylogeny and the possible impact on human health. To assess time- and/or condition-dependent changes in microbial diversity within REMS, a total of 27 air samples were collected during three consecutive months. The microbial burden and diversity were elucidated using culture-dependent and culture-independent molecular methods. The results indicate that during controlled conditions the total microbial burden detected by culture-dependent techniques (below a detectable level to 102 cells m−3of air) and intracellular ATP assay was significantly low (102–103 cells m−3of air), but increased during the uncontrolled post-operation phase (∼104 cells m−3of air). Only Gram-positive and α-proteobacteria grew under tested culture conditions, with a predominant occurrence ofMethylobacterium radiotolerans, andSphingomonas yanoikuyae. Direct DNA extraction and 16S rDNA sequencing methodology revealed a broader diversity of microbes present in the REMS air (51 species). Unlike culture-dependent analysis, both Gram-positive and proteobacteria were equally represented, while members of a few proteobaterial groups dominated (Rhodopseudomonas,Sphingomonas,Acidovorax,Ralstonia,Acinetobacter,Pseudomonas, andPsychrobacter). Although the presence of several opportunistic pathogens warrants further investigation, the results demonstrated that routine maintenance such as controlling the humidity, crew’s daily cleaning, and air filtration were effective in reducing the microbial burden in the REMS.

Description: Experimentation and theory are used to study the long-term dynamics of a two-dimensional density current flowing into a two-layer stratified basin. When the initial Richardson number,$hbox{it Ri}_{ ho}^{hbox{scriptsizeit in}}$, characterizing the ratio of the background stratification to the buoyancy flux of the density current, is less than the critical value of$hbox{it Ri}_{ ho}^{*} ,{=}, 21-27$, it is found that the density current penetrates the stratified interface. This result is ostensibly independent of slope for angles between 30° and 90°. If the current does not initially penetrate the interface, then it slowly increases the density of the top layer until the interfacial density difference is reduced sufficiently to drive penetration. The time scale for this to occur,$t_{p} ,{=}, (hbox{it Ri}^{hbox{scriptsizeit in}}_{ ho} - hbox{it Ri}_{ ho}^{*}) L/B^{1/3}$, is explicitly a function of the buoyancy fluxBand the length of the basinL. The initial Richardson number,$hbox{it Ri}^{hbox{scriptsizeit in}}_{ ho}$, is a function of depth, the initial reduced gravity of the interface and a weak function of slope angle. In the absence of initial penetration for very steep slopes of 75° and 90°, we observe that penetrative convection at the interface leads to significant local entrainment. In consequence, the top layer thickens and the interfacial entrainment rate increases as the fifth power of the interfacial Froude number. In contrast, such a process is not observed at comparable interfacial Froude numbers on lower slopes of 30°, 45° and 60°, thereby demonstrating the important role of impact angle on penetrative convection. We attribute the increased interfacial entrainment by the steep density currents as the result of the transition from an undular bore to a turbulent hydraulic jump at the point where the density current intrudes. We discuss the applicability of the observed circulation to the stability of the Arctic halocline where we find$0.56,{lesssim}, t_{p} ,{lesssim},1.2$years for a range of contemporary oceanographic conditions.

Description: Direct numerical simulations (DNS) of stenotic flows under conditions of steady inlet flow were discussed in Part 1 of this study. DNS of pulsatile flow through the 75% stenosed tube (by area) employed for the computations in Part 1 is examined here. Analogous to the steady flow results, DNS predicts a laminar post-stenotic flow field in the case of pulsatile flow through the axisymmetric stenosis model, in contrast to previous experiments, in which intermittent disturbed flow regions and turbulent breakdown were observed in the downstream region. The introduction of a stenosis eccentricity, that was 5% of the main vessel diameter at the throat, resulted in periodic, localized transition to turbulence. Analysis in this study indicates that the early and mid-acceleration phases of the time period cycle were relatively stable, with no turbulent activity in the post-stenotic region. However, towards the end of acceleration, the starting vortex, formed earlier as the fluid accelerated through the stenosis at the beginning of acceleration, started to break up into elongated streamwise structures. These streamwise vortices broke down at peak flow, forming a turbulent spot in the post-stenotic region. In the early part of deceleration there was intense turbulent activity within this spot. Past the mid-deceleration phase, through to minimum flow, the inlet flow lost its momentum and the flow field began to relaminarize. The start of acceleration in the following cycle saw a recurrence of the entire process of a starting structure undergoing turbulent breakdown and subsequent relaminarization of the post-stenotic flow field. Peak wall shear stress (WSS) levels occurred at the stenosis throat, with the rest of the vessel experiencing much lower levels. Turbulent breakdown at peak flow resulted in a sharp amplification of instantaneous WSS magnitudes across the region corresponding to the turbulent spot, accompanied by large axial and circumferential fluctuations, even while ensemble-averaged axial shear stresses remained mostly low and negative. WSS levels dropped rapidly after the mid-deceleration phase, when the relaminarization process took over, and were almost identical to laminar, axisymmetric shear levels through most of the acceleration phase.

Description: We investigate experimentally the force generated by the unsteady vortex formation of low-aspect-ratio normal flat plates with one end free. The objective of this study is to determine the role of the free end, or tip, vortex. Understanding this simple case provides insight into flapping-wing propulsion, which involves the unsteady motion of low-aspect-ratio appendages. As a simple model of a propulsive half-stroke, we consider a rectangular normal flat plate undergoing a translating start-up motion in a towing tank. Digital particle image velocimetry is used to measure multiple perpendicular sections of the flow velocity and vorticity, in order to correlate vortex circulation with the measured plate force. The three-dimensional wake structure is captured using flow visualization. We show that the tip vortex produces a significant maximum in the plate force. Suppressing its formation results in a force minimum. Comparing plates of aspect ratio six and two, the flow is similar in terms of absolute distance from the tip, but evolves faster for aspect ratio two. The plate drag coefficient increases with decreasing aspect ratio.

Description: Using the Boussinesq approximation, a set of depth-integrated wave equations for long-wave propagation over a mud bed is derived. The wave motions above the mud bed are assumed to be irrotational and the mud bed is modelled as a highly viscous fluid. The pressure and velocity are required to be continuous across the water–mud interface. The resulting governing equations are differential–integral equations in terms of the depth-integrated horizontal velocity and the free-surface displacement. The effects of the mud bed appear in the continuity equation in the form of a time integral of weighted divergence of the depth-averaged velocity. Damping rates for periodic waves and solitary waves are calculated. For the solitary wave case, the velocity profiles in the water column and the mud bed at different phases are discussed. The effects of the viscous boundary layer above the mud–water interface are also examined.

Description: During extended deployment at an ocean observatory off the coast of New Jersey, a bottom-mounted five-beam acoustic Doppler current profiler measured large-scale velocity structures that we interpret as Langmuir circulations filling the entire water column. These circulations are the large-eddy structures of wind-wave-driven turbulent flows that occur episodically when a shallow water column experiences prolonged strong wind forcing. Many observational characteristics agree with former descriptions of Langmuir circulations in deep water. The three-dimensional velocity field reveals quasi-organized structures consisting of pairs of surface-intensified counter-rotating vortices, aligned approximately downwind. Maximum downward velocities are stronger than upward velocities, and the downwelling region of each cell, defined as a pair of vortices, is narrower than the upwelling region. Maximum downward vertical velocity occurs at or above mid-depth, and scales approximately with wind speed. The estimated crosswind scale of cells is roughly 3–6 times their vertical scale, set under these conditions by water depth. The long axis of the cells appears to lie at an angle ∼10°–20° to the right of the wind. A major difference from deep-water observations is strong near-bottom intensification of the downwind ‘jets’ found typically centred over downwelling regions. Accessible observational features such as cell morphology and profiles of mean velocities, turbulent velocity variances, and shear stress components are compared with the results of associated large-eddy simulations (reported in Part 2) of shallow water flows driven by surface stress and the Craik–Leibovich vortex forcing generally used to represent generation of Langmuir cells. A particularly sensitive diagnostic for identification of Langmuir circulations as the energy-containing eddies of the turbulent flow is the depth trajectory of invariants of the turbulent stress tensor, plotted in the Lumley ‘triangle’ corresponding to realizable turbulent flows. When Langmuir structures are present in the observations, the Lumley map is distinctly different from that of surface-stress-driven Couette flow, again in agreement with the large-eddy simulations (LES). Unlike the LES, observed velocity fields contain two distinct and significant scales of variability, documented by wavelet analysis of observational records of vertical velocity. Variability with periods of many minutes is that expected from Langmuir cells drifting past the instrument at the slowly time-varying crosswind velocity. Shorter period variability, of the order of 1–2 min, has roughly the observed periodicity of surface wave groups, suggesting a connection with the wave groups themselves and/or the wave breaking associated with them in high wind conditions.

Description: Lagrangian auto- and cross-correlation functions of the rate of strains2, enstrophy ω2, their respective production terms −sijsjkskiand ωiωjsij, and material derivatives, Ds2/Dtand Dω2/Dtare estimated using experimental results obtained through three-dimensional particle tracking velocimetry (three-dimensional-PTV) in homogeneous turbulence atReλ=50. The autocorrelation functions are used to estimate the Lagrangian time scales of different quantities, while the cross-correlation functions are used to clarify some aspects of the interaction mechanisms between vorticity ω and the rate of strain tensorsij, that are responsible for the statistically stationary, in the Eulerian sense, levels of enstrophy and rate of strain in homogeneous turbulent flow. Results show that at the Reynolds number of the experiment these quantities exhibit different time scales, varying from the relatively long time scale of ω2to the relatively shorter time scales ofs2, ωiωjsijand −sijsjkski. Cross-correlation functions suggest that the dynamics of enstrophy and strain, in this flow, is driven by a set of different-time-scale processes that depend on the local magnitudes ofs2and ω2. In particular, there are indications that, in a statistical sense, (i) strain production anticipates enstrophy production in low-strain–low-enstrophy regions (ii) strain production and enstrophy production display high correlation in high-strain–high-enstrophy regions, (iii) vorticity dampening in high-enstrophy regions is associated with weak correlations between −sijsjkskiands2and between −sijsjkskiand Ds2/Dt, in addition to a marked anti-correlation between ωiωjsijand Ds2/Dt. Vorticity dampening in high-enstrophy regions is thus related to the decay ofs2and its production term, −sijsjkski.

Description: Transverse jets arise in many applications, including propulsion, effluent dispersion, oil field flows, and V/STOL aerodynamics. This study seeks a fundamental, mechanistic understanding of the structure and evolution of vorticity in the transverse jet. We develop a high-resolution three-dimensional vortex simulation of the transverse jet at large Reynolds number and consider jet-to-crossflow velocity ratiosrranging from 5 to 10. A new formulation of vorticity-flux boundary conditions accounts for the interaction of channel wall vorticity with the jet flow immediately around the orifice. We demonstrate that the nascent jet shear layer contains not only azimuthal vorticity generated in the jet pipe, but wall-normal and azimuthal perturbations resulting from the jet–crossflow interaction. This formulation also yields analytical expressions for vortex lines in the near field as a function ofr.Transformation of the cylindrical shear layer emanating from the orifice begins with axial elongation of its lee side to form sections of counter-rotating vorticity aligned with the jet trajectory. Periodic roll-up of the shear layer accompanies this deformation, creating complementary vortex arcs on the lee and windward sides of the jet. Counter-rotating vorticity then drives lee-side roll-ups in the windward direction, along the normal to the jet trajectory. Azimuthal vortex arcs of alternating sign thus approach each other on the windward boundary of the jet. Accordingly, initially planar material rings on the shear layer fold completely and assume an interlocking structure that persists for several diameters above the jet exit. Though the near field of the jet is dominated by deformation and periodic roll-up of the shear layer, the resulting counter-rotating vorticity is a pronounced feature of the mean field; in turn, the mean counter-rotation exerts a substantial influence on the deformation of the shear layer. Following the pronounced bending of the trajectory into the crossflow, we observe a sudden breakdown of near-field vortical structures into a dense distribution of smaller scales. Spatial filtering of this region reveals the persistence of counter-rotating streamwise vorticity initiated in the near field.

Description: The decay of a passive scalar in a sinusoidal shear flow translating in the cross-stream direction at a constant speeduis studied in the limit of small diffusivity κ. The decay rate, obtained by solving an eigenvalue problem, is found to tend to a non-zero constant as κ→0 whenuis of order κ1/2. This result, establishing that fast decay is possible in shear flows, is fragile however: because of the existence of pseudomodes, the addition of a small noise leads to decay rates that decrease to zero with κ as κ2/5.

Description: Direct numerical simulation is used to study a round turbulent jet in a laminar crossflow. The ratio of bulk jet velocity to free-stream crossflow velocity is 5.7 and the Reynolds number based on the bulk jet velocity and the jet exit diameter is 5000. The mean velocity and turbulent intensities from the simulations are compared to data from the experiments by Su & Mungal (2004) and good agreement is observed. Additional quantities, not available from experiments, are presented. Turbulent kinetic energy budgets are computed for this flow. Examination of the budgets shows that the near field is far from a state of turbulent equilibrium – especially along the jet edges. Also – in the near field – peak kinetic energy production is observed close to the leading edge, while peak dissipation is observed toward the trailing edge of the jet. The results are used to comment upon the difficulty involved in predicting this flow using RANS computations. There exist regions in this flow where the pressure transport term, neglected by some models and poorly modelled by others, is significant. And past the jet exit, the flow is not close to established canonical flows on which most models appear to be based.

Description: The theory of Wagner from 1932 for the normal symmetric impact of a two-dimensional body of small deadrise angle on a half-space of ideal and incompressible liquid is extended to derive the second-order corrections for the locations of the higher-pressure jet-root regions and for the upward force on the impactor using a systematic matched-asymptotic analysis. The second-order predictions for the upward force on an entering wedge and parabola are compared with numerical and experimental data, respectively, and it is concluded that a significant improvement in the predictive capability of Wagner's theory is afforded by proceeding to second order.

Description: SUMMARYThe purpose of the current work was to study phosphorus (P) metabolism in growing sheep supplemented with different levels of dicalcium phosphate using an extant mathematical model. Twelve male non-castrated Santa Inês sheep, weighing 23 (±2·2) kg, received a basal diet unsupplemented or supplemented with dicalcium phosphate to provide 1·5, 3·0, 4·5 g of P/animal per day (treatments T1 to T4, respectively). After 3 weeks adaptation, 7·4 MBq of32P was injected into the jugular vein of each animal. Samples of blood, faeces and urine were collected every day during a 7-day period and thereafter the animals were sacrificed and samples from liver, kidney, heart, muscle and bone were collected for specific activity and inorganic P determinations. The flows between gut and plasma were similar for each treatment except for T1, which showed the lowest values for both flows (P

Description: SUMMARYMedial claws from the right hind feet were obtained post mortem from four 19–20-month-old beef heifers and from four 28-month-old first-calving dairy heifers 3–4 days postpartum. X-ray diffraction (XRD) studies were undertaken on samples of soft and hard (cornified) integumental tissues of dorsal wall, sole and heel (bulb) for varying sites and planes of exposure. The measurements were interpreted as defining diffraction patterns and intermolecular spacings of cytoskeletal and extracellular fibrous structural proteins. The orientation of these proteins was examined in relation to physical characteristics and function including bearing of body weight by these tissues.Physical measurements taken included impression hardness which showed typically greater values for wall than sole and variable differences between horn of dairy and beef origin and husbandry systems. Claws from dairy heifers had significantly smaller values for toe (dorsal wall) angle, claw height and heel height and thickness of solear horn and heel soft tissue. Although few were studied, the results reflected typical husbandry origins and indicated the susceptibility to the lesion formation well recognized in postpartum dairy cattle.Typical XRD patterns for horn samples showed defined arcs of reflectance on the equatorial axis consistent with findings for the presence of α-helices in fibrils reported to occur in other hard-keratin-containing integumental tissues. However, reflectance on the meridional axis also reported for these other tissues was not detected. A similar defined pattern was obtained for less than 0·10 of samples of internal soft pre-cornified epidermal and attached dermal tissue although the values for intermolecular ‘d’ spacing for these were consistent with those reported for type I collagen. Diffuse reflection patterns were thus evident for the majority of samples of soft tissue epidermis and dermis and also for adipose tissue of the digital cushion.The formation of defined arcs of reflectance allowed the determination of fibril alignment in wall and solear horn. For the orientated samples of dorsal wall horn tissue, the outer layer showed a longitudinal angle of orientation essentially maintained proximal to distal. This pattern was maintained throughout the depth of horn at the proximal site. In contrast, layers in mid-wall and towards the distal edge showed a greater circumferential (horizontal) orientation in sections collected anterior to posterior towards the inner corial, including laminar, tissues. The orientation of fibrils in inner wall horn appears to relate to the direction of load-bearing forces in connecting horn to the distal phalanx. Horizontal alignment of fibrils was observed in the sole. In presenting the long axis of cells to the ground surface this orientation may facilitate erosive forces and contribute to the thinning of cornified sole horn under adverse underfoot conditions.

Description: The flapping coupling between two filaments is studied theoretically and experimentally in this paper. A temporal linear instability analysis is carried out based on a simplified hydrodynamic model. The dispersion relationship between the eigen-frequency ω and wavenumberkis expressed by a quartic equation. Two special cases of flapping coupling, i.e. two identical filaments having the same length and two filaments having different lengths, are studied in detail. In the case of two identical filaments, the theoretical analysis predicts four coupling modes, i.e. the stretched-straight mode, the antisymmetrical in-phase mode, the symmetrical out-of-phase mode and the indefinite mode. The theory also predicts the existence of an eigenfrequency jump during transition between the in-phase and out-of-phase modes, which has been observed in previous experiments and numerical simulations. In the case of two filaments having different lengths, four modes similar to those in the former case are identified theoretically. The distribution of coupling modes for both the cases is shown in two planes. One is a dimensionless plane ofSvs.U, whereSis the density ratio of solid filament to fluid andU2is the ratio of fluid kinetic energy to solid elastic potential energy. The other is a dimensional plane of the half-distance (h) between two filaments vs. the filament length (L). Relevant experiments are carried out in a soap-film tunnel and the stable and unstable modes are observed. Theory and experiment are compared in detail. It should be noted that the model used in our analysis is a very simplified one that can provide intuitional analytical results of the coupling modes as well as their qualitative distributions. The factors neglected in our model, such as vortex shedding, viscous and nonlinear effects, do not allow the model to predict results precisely consistent with the experiments. Moreover, the Strouhal numbers of the flapping filaments are found to be generally around a fixed value in the experiments for both cases, implying that the filaments try to maintain a lower potential energy state.

Description: Numerical simulations on a β-plane are used to further understand the formation of zonal flows from small-scale fluctuations. The dynamics of ‘reduced models’ are computed by restricting the nonlinear term to include a subset of triad interactions in Fourier space. Reduced models of near-resonant triads are considered, as well as the complement set of non-resonant triads. At moderately small values of the Rhines number, near-resonant triad interactions are shown to be responsible for the generation of large-scale zonal flows from small-scale random forcing. Without large-scale drag, both the full system and the reduced model of near resonances produce asymmetry between eastward and westward jets, in favour of stronger westward jets. When large-scale drag is included, the long-time asymmetry is reversed in the full system, with eastward jets that are thinner and stronger than westward jets. Then the reduced model of near resonances exhibits a weaker asymmetry, but there are nevertheless more eastward jets stronger than a threshold value.

Description: Free surface potential flows past disturbances in a channel are considered. Three different types of disturbance are studied: (i) a submerged obstacle on the bottom of a channel; (ii) a pressure distribution on the free surface; and (iii) an obstruction in the free surface (e.g. a sluice gate or a flat plate). Surface tension is neglected, but gravity is included in the dynamic boundary condition. Fully nonlinear solutions are computed by boundary integral equation methods. In addition, weakly nonlinear solutions are derived. New solutions are found when several disturbances are present simultaneously. They are discovered through the weakly nonlinear analysis and confirmed by numerical computations for the fully nonlinear problem.

Description: A numerical description of heterogeneous propellant combustion enables us to examine the spatial and temporal fluctuations in the flow field arising from the heterogeneity. Particular focus is placed on the fluctuations in a zone intermediate between the combustion field (where reaction is important) and the chamber flow domain, for these define boundary conditions for simulations of the turbulent chamber flow. The statistics of the temperature field and the normal velocity field are described, and characteristic length scales and time scales are identified. The length scales are small compared to any relevant length scale of the chamber flow, and so the boundary conditions for this flow at any mesh point are statistically independent of those at any other mesh point. But the temporal correlations at a fixed point are significant, and affect the nature of the chamber flow in a variety of ways. We describe the fluctuations in the head-end pressure that arise because of them, and contrast these results with those calculated using a white-noise assumption.

Description: We consider the flapping stability and response of a thin two-dimensional flag of high extensional rigidity and low bending rigidity. The three relevant non-dimensional parameters governing the problem are the structure-to-fluid mass ratio, μ = ρsh/(ρfL); the Reynolds number,Rey=VL/ν; and the non-dimensional bending rigidity,KB=EI/(ρfV2L3). The soft cloth of a flag is represented by very low bending rigidity and the subsequent dominance of flow-induced tension as the main structural restoring force. We first perform linear analysis to help understand the relevant mechanisms of the problem and guide the computational investigation. To study the nonlinear stability and response, we develop a fluid–structure direct simulation (FSDS) capability, coupling a direct numerical simulation of the Navier–Stokes equations to a solver for thin-membrane dynamics of arbitrarily large motion. With the flow grid fitted to the structural boundary, external forcing to the structure is calculated from the boundary fluid dynamics. Using a systematic series of FSDS runs, we pursue a detailed analysis of the response as a function of mass ratio for the case of very low bending rigidity (KB= 10−4) and relatively high Reynolds number (Rey= 103). We discover three distinct regimes of response as a function of mass ratio μ: (I) a small μ regime of fixed-point stability; (II) an intermediate μ regime of period-one limit-cycle flapping with amplitude increasing with increasing μ; and (III) a large μ regime of chaotic flapping. Parametric stability dependencies predicted by the linear analysis are confirmed by the nonlinear FSDS, and hysteresis in stability is explained with a nonlinear softening spring model. The chaotic flapping response shows up as a breaking of the limit cycle by inclusion of the 3/2 superharmonic. This occurs as the increased flapping amplitude yields a flapping Strouhal number (St= 2Af/V) in the neighbourhood of the natural vortex wake Strouhal number,St≃ 0.2. The limit-cycle von Kármán vortex wake transitions in chaos to a wake with clusters of higher intensity vortices. For the largest mass ratios, strong vortex pairs are distributed away from the wake centreline during intermittent violent snapping events, characterized by rapid changes in tension and dynamic buckling.

Description: The governing equations for the two-point correlations of the turbulent fluctuating velocity in the temporally evolving wake were analysed to determine whether they could have equilibrium similarity solutions. It was found that these equations could have such solutions for a finite-Reynolds-number wake, where the two-point velocity correlations could be written as a product of a time-dependent scale and a function dependent only on similarity variables. It is therefore possible to collapse the two-point measures of all the scales of motions in the temporally evolving wake using a single set of similarity variables. As in an earlier single-point analysis, it was found that the governing equations for the equilibrium similarity solutions could not be reduced to a form that was independent of a growth-rate dependent parameter. Thus, there is not a single ‘universal’ solution that describes the state of the large-scale structures, so that the large-scale structures in the far field may depend on how the flow is generated.The predictions of the similarity analysis were compared to the data from two direct numerical simulations of the temporally evolving wakes examined previously. It was found that the two-point velocity spectra of these temporally evolving wakes collapsed reasonably well over the entire range of scales when they were scaled in the manner deduced from the equilibrium similarity analysis. Thus, actual flows do seem to evolve in a manner consistent with the equilibrium similarity solutions.

Description: This paper examines the consistency of the exact scaling laws for isotropic magnetohydrodynamic (MHD) turbulence in numerical simulations with large magnetic Prandtl numbersPmand withPm= 1. The exact laws are used to elucidate the structure of the magnetic and velocity fields. Despite the linear scaling of certain third-order correlation functions, the situation is not analogous to the case of Kolmogorov turbulence. The magnetic field is adequately described by a model of a stripy (folded) field with direction reversals at the resistive scale. At currently available resolutions, the cascade of kinetic energy is short-circuited by the direct exchange of energy between the forcing-scale motions and the stripy magnetic fields. This non-local interaction is the defining feature of isotropic MHD turbulence.

Description: SUMMARYThe development of the genetic control of nematode resistance in growing lambs is of biological interest, as well as being important in terms of designing practical strategies to breed for increased nematode resistance. The current paper demonstrates the use of random regression techniques for quantifying the development of the heritability of faecal egg count (Fec), the indicator of nematode resistance, in growing lambs and predicted inter-age genetic and phenotypic correlations for Fec. Fec data from 732 lambs, collected at 4-week intervals fromc. 8–24 weeks of age, were analysed using random regression techniques. Random effects fitted in the model included genetic, individual animal environmental, litter and residual random effects. Output (co)variance components were interpolated to weekly time points. Individual variance components showed complex patterns of change over time; however, the estimated heritability increased smoothly with age, from 0·10 to 0·38, and showed more stable time trends than were obtained from univariate analyses of Fec at individual time points. Inter-age correlations decreased as the time interval between measurements increased. Genetic correlations were always positive, with 0·6 of all possible inter-age correlations being greater than 0·80. Phenotypic correlations were lower, and decreased more quickly as the time interval between measurements increased. The results presented confirm biological understanding of the development of immunity to nematode infections in growing lambs. Additionally, they provide a tool to determine optimal sampling ages when assessing lambs' relative resistance to nematode infections.

Description: SUMMARYSince many soils used for growing potatoes in the UK are likely to be close to their plastic limit for cultivation during early spring, there exists the potential for soil compaction to occur during planting which will restrict root penetration. A series of experiments showed that soil compaction delayed emergence, reduced rate of leaf appearance and ground cover expansion, shortened canopy cover duration and restricted light interception, which combined to reduce tuber yield. Rooting density and maximum depth of rooting were reduced, particularly where compaction was shallow. In some soils, irrigation helped alleviate some of the effects of compaction but in others it exacerbated their severity. Using a cone penetrometer, relationships between rate of root penetration and soil resistance (Ω) were established from a number of experiments and replicated blocks in commercial fields and an overall relationship of the formy=16·3–4·08Ω mm/day was produced. Root penetration rates ofc. 20 mm/day were measured in the intensively-cultivated ridge zone but growth rates were halved at a Ω of 1·5 MPa. A survey of 602 commercial fields showed that two thirds of fields had Ωs ⩾3 MPa (where root growth rates would be

Description: The hypothesis of the explosion of a number of planets and moons of our Solar System during its 4.6-billion-year history is in excellent accord with all known observational constraints, even without adjustable parameters orad hochelper hypotheses. Many of its boldest predictions have been fulfilled. In most instances, these predictions were judged highly unlikely by the current standard models. Moreover, in several cases, the entire exploded planet model was at risk of being falsified if the predictions failed. The successful predictions include: (1) satellites of asteroids; (2) satellites of comets; (3) salt water in meteorites; (4) ‘roll marks’ leading to boulders on asteroids; (5) the time and peak rate of the 1999 Leonid meteor storm; (6) explosion signatures for asteroids; (7) the strongly spiked energy parameter for new comets; (8) the distribution of black material on slowly rotating airless bodies; (9) splitting velocities of comets; (10) the asteroid-like nature ofDeep Impacttarget Comet Tempel 1; and (11) the presence of high-formation-temperature minerals in theStardustcomet dust sample return. In physics and astronomy, hypotheses are either falsified if their predictions fail, or proved to be of value if they succeed. By all existing evidence, the exploded planet hypothesis has proved far more useful than the half-dozen or so hypotheses it would replace. Among the many important corollaries are these. (a) Perhaps as many as six former planets of our Solar System have exploded over its 4.6-billion-year history. (b) In particular, Mars is not an original planet, but a former moon of an exploded planet. (c) As a major player in Solar System evolution, the exploded planet scenario must be considered as a likely propagation vehicle for the spread of biogenic organisms. We conclude with a brief mention of three possible planetary explosion mechanisms.

Description: The degree and nature of association between trace metals (Cu, Pb, Zn, Ni, Ag, & Cd) and cyanobacterial mats, phytoplankton and sediments has been assessed in the Lake Vanda region of the Wright Valley, Victoria Land. Trace metal:Fe ratios and SEM imaging confirmed that apparent trace metal enrichment in cyanobacterial mats, relative to the sediment beneath, was due to incorporation of fine (sub-micron) sediment particles in the muciligenous matrix of the mat. In suspended particulate material (SPM) filtered from the oxic water of Lake Vanda and the Onyx River, the degree of trace metal binding to the SPM did not appear to correlate with phytoplankton content. Instead a positive correlation was observed between Fe and trace metal content. The SPM at the top of the lake water column, where only the finest sediment remains suspended, had the highest trace metal concentrations. It is concluded that the trace metal content of cyanobacterial mats and phytoplankton samples is primarily due to incorporation of fine sediment particles of high surface area which therefore enhance trace metal adsorption capacity. This reinforces the existing hypothesis that trace metal solubility in this environment is primarily controlled by abiotic processes.

Description: The statistical properties of fully developed planar turbulent Couette–Poiseuille flow result from the simultaneous imposition of a mean wall shear force together with a mean pressure force. Despite the fact that pure Poiseuille flow and pure Couette flow are the two extremes of Couette–Poiseuille flow, the statistical properties of the latter have proved resistant to scaling approaches that coherently extend traditional wall flow theory. For this reason, Couette–Poiseuille flow constitutes an interesting test case by which to explore the efficacy of alternative theoretical approaches, along with their physical/mathematical ramifications. Within this context, the present effort extends the recently developed scaling framework of Weiet al. (2005a) and associated multiscaling ideas of Fifeet al. (2005a,b) to fully developed planar turbulent Couette–Poiseuille flow. Like Poiseuille flow, and contrary to the structure hypothesized by the traditional inner/outer/overlap-based framework, with increasing distance from the wall, the present flow is shown in some cases to undergo abalance breakingandbalance exchangeprocess as the mean dynamics transition from a layer characterized by a balance between the Reynolds stress gradient and viscous stress gradient, to a layer characterized by a balance between the Reynolds stress gradient (more precisely, the sum of Reynolds and viscous stress gradients) and mean pressure gradient. Multiscale analyses of the mean momentum equation are used to predict (in order of magnitude) the wall-normal positions of the maxima of the Reynolds shear stress, as well as to provide an explicit mesoscaling for the profiles near those positions. The analysis reveals a close relationship between the mean flow structure of Couette–Poiseuille flow and two internal scale hierarchies admitted by the mean flow equations. The averaged profiles of interest have, at essentially each point in the channel, a characteristic length that increases as a well-defined ‘outer region’ is approached from either the bottom or the top of the channel. The continuous deformation of this scaling structure as the relevant parameter varies from the Poiseuille case to the Couette case is studied and clarified.

Description: In this paper, we use numerical simulation and laboratory experimental observation to show that fire whirls can be generated spontaneously through the interaction between a central flame and surrounding organized or randomly distributed flames. The momentum of the air stream entrained by the main flame decreases as it crosses a surrounding flame, so that the main flame rotates if surrounding flames are arranged in such a way as to block the passage of the air stream directed towards the centre of the main flame and to favour flows in a particular circumferential direction. An analysis is performed to study the role of the rotation speed in the flame height. It is found that the flame height initially decreases to a minimum owing to the inflow boundary layer wind reducing the initial vertical velocity of gas for low rotation speed and to entrainment enhancement reducing the rising time, and then it increases owing to the pressure reduction at the centre of the rotating vortex and entrainment suppression extending the rising time.

Description: We investigate aspects of hovering insect flight by finding the optimal wing kinematics which minimize power consumption while still providing enough lift to maintain a time-averaged constant altitude over one flapping period. In particular, we study the flight of three insects whose masses vary by approximately three orders of magnitude: fruitfly (Drosophila melanogaster), bumblebee (Bombus terrestris), and hawkmoth (Manduca sexta). Here, we model an insect wing as a rigid body with three rotational degrees of freedom. The aerodynamic forces are modelled via a quasi-steady model of a thin plate interacting with the surrounding fluid. The advantage of this model, as opposed to the more computationally costly method of direct numerical simulation via computational fluid dynamics, is that it allows us to perform optimization procedures and detailed sensitivity analyses which require many cost function evaluations. The optimal solutions are found via a hybrid optimization algorithm combining aspects of a genetic algorithm and a gradient-based optimizer. We find that the results of this optimization yield kinematics which are qualitatively and quantitatively similar to previously observed data. We also perform sensitivity analyses on parameters of the optimal kinematics to gain insight into the values of the observed optima. Additionally, we find that all of the optimal kinematics found here maintain the same leading edge throughout the stroke, as is the case for nearly all insect wing motions. We show that this type of stroke takes advantage of a passive wing rotation in which aerodynamic forces help to reverse the wing pitch, similar to the turning of a free-falling leaf.

Description: Direct numerical simulations (DNS) of steady and pulsatile flow through 75% (by area reduction) stenosed tubes have been performed, with the motivation of understanding the biofluid dynamics of actual stenosed arteries. The spectral-element method, providing geometric flexibility and high-order spectral accuracy, was employed for the simulations. The steady flow results are examined here while the pulsatile flow analysis is dealt with in Part 2 of this study. At inlet Reynolds numbers of 500 and 1000, DNS predict a laminar flow field downstream of an axisymmetric stenosis and comparison to previous experiments show good agreement in the immediate post-stenotic region. The introduction of a geometric perturbation within the current model, in the form of a stenosis eccentricity that was 5% of the main vessel diameter at the throat, resulted in breaking of the symmetry of the post-stenotic flow field by causing the jet to deflect towards the side of the eccentricity and, at a high enough Reynolds number of 1000, jet breakdown occurred in the downstream region. The flow transitioned to turbulence about five diameters away from the stenosis, with velocity spectra taking on a broadband nature, acquiring a -5/3 slope that is typical of turbulent flows. Transition was accomplished by the breaking up of streamwise, hairpin vortices into a localized turbulent spot, reminiscent of the turbulent puff observed in pipe flow transition, within which r.m.s. velocity and turbulent energy levels were highest. Turbulent fluctuations and energy levels rapidly decayed beyond this region and flow relaminarized. The acceleration of the fluid through the stenosis resulted in wall shear stress (WSS) magnitudes that exceeded upstream levels by more than a factor of 30 but low WSS levels accompanied the flow separation zones that formed immediately downstream of the stenosis. Transition to turbulence in the case of the eccentric stenosis was found to be manifested as large temporal and spatial gradients of shear stress, with significant axial and circumferential variations in instantaneous WSS.

Description: Panspermia, the dissemination of life through space, would require resistance to the conditions found in space, including UV light. All known life forms depend on DNA to store information. In an effort to understand the liabilities of DNA to UV light and modes of DNA protection in terrestrial life forms, we established UV–VUV (125–340 nm) absorption spectra for dry DNA and its polymerized components and mononucleotides, as well as for a selection of potential UV screens ubiquitous in all organisms, including proteins, selected amino acids and amines (polyamines and tyramine). Montmorillonite clay was included as a potential abiotic UV screen. Among the potential screens tested, adenosine triphosphate (ATP) appeared to be particularly attractive, because its UV absorption spectrum was similar to that of DNA. We suggest that the use of ATP in UV protection could have pre-dated its current role in energy transfer. Spectroscopy also showed that UV absorption varied according to nucleotide content, suggesting that base pair usage could be a factor in adaptation to given UV environments and the availability of UV screens.

Description: About 5400 cal yr BP, a large landslide formed a 〉 400-m-tall dam in the upper Marsyandi River, central Nepal. The resulting lacustrine and deltaic deposits stretched 〉 7 km upstream, reaching a thickness of 120 m.14C dating of 7 wood fragments reveals that the aggradation and subsequent incision occurred remarkably quickly (∼ 500 yr). Reconstructed volumes of lacustrine (∼ 0.16 km3) and deltaic (∼ 0.09 km3) deposits indicate a bedload-to-suspended load ratio of 1:2, considerably higher than the ≤ 1:10 that is commonly assumed. At the downstream end of the landslide dam, the river incised a new channel through ≥ 70 m of Greater Himalayan gneiss, requiring a minimum bedrock incision rate of 13 mm/yr over last 5400 yr. The majority of incision presumably occurred over a fraction of this time, suggesting much higher rates. The high bedload ratio from such an energetic mountain river is a particularly significant addition to our knowledge of sediment flux in orogenic environments.

Description: We present both analytical and numerical solutions describing seepage flows in an unsaturated permeable seabed induced by transient long waves. The effects of compressibility of pore water in the seabed due to a small degree of unsaturation are considered in the investigation. To make the problem tractable analytically, we first focus our attention on situations where the horizontal scale of the seepage flow is much larger than the vertical scale. With this simplification the pore-water pressure in the soil column is governed by a one-dimensional diffusion equation with a specified pressure at the water–seabed interface and the no-flux condition at the bottom of the seabed. Analytical solutions for pore-water pressure and velocity are obtained for arbitrary transient waves. Special cases are studied for periodic waves, cnoidal waves, solitary waves and bores. Numerical solutions are also obtained by simultaneously solving the Navier–Stokes equations for water wave motions and the exact two-dimensional diffusion equation for seepage flows in the seabed. The analytical solutions are used to check the accuracy of the numerical methods. On the other hand, numerical solutions extend the applicability of the analytical solutions. The liquefaction potential in a permeable bed as well as the energy dissipation under various wave conditions are then discussed.

Description: The Selborne Group comprises two metamorphic rock units, the muscovite±dolomite bearing Madison Marble and the biotite-muscovite±quartz-calcite Contortion Schist, which contains thick lenses of variably deformed metabasalts and metaconglomerates. Petrological and structural data indicate a polyphase metamorphic evolution including: i) an early stage of upper greenschist regional metamorphism (P = ~0.15–0.3 GPa; T = ~380–450°C), ii) prograde metamorphism during D1up to amphibolite facies peak conditions (P = 0.58–0.8 GPa, T = ~560–645°C), iii) syn-D2unloading-cooling retrograde metamorphism, iv) a post-D2contact metamorphic overprint at variable T between 450 and 550°C and ~0.2 GPa connected to the emplacement of granitic plutons and felsic dyke swarms. Geochronological data constrain the polyphase syn-D1/D2evolution between ~ 510 and 492 Ma. A similar metamorphic path, including a medium P stage but at lower T conditions, is documented in greenschist facies metabasalts within the Byrd Group in the Mount Dick area. The metamorphic pattern and close lithostratigraphic matching between Selborne Group and Byrd Group sharply contrast with the high-grade Horney Formation that is exposed north of the Byrd Glacier and corroborate the hypothesis that the Byrd Glacier discontinuity marks a first-order crustal tectonic boundary crossing the Ross orogen.

Description: Recent experiments have shown intriguing regions of intense luminescence or ‘hotspots’ in the vicinity of triple-point shear layers in propagating gaseous detonation waves. Localized explosions have also been observed to develop in these fronts. These features were observed in higher effective activation energy mixtures, but not in lower effective activation energy mixtures. The increased lead shock oscillation through a cell cycle in higher activation energy mixtures may result in a significantly increased disparity in the induction time on either side of the triple-point shear layer, and thus an enhanced mixing between reacted and non-reacted streams supported by Kelvin–Helmholtz instability. The relation between the shear-layer instability and the mixture effective activation energy is analysed by carrying out a spatial linear stability study for three mixtures with different activation energies and injection conditions that correspond to the experimental conditions. The role of vortical structures associated with Kelvin–Helmholtz instability in the formation of localized ignition is investigated by performing two-dimensional Navier–Stokes simulations with detailed chemical kinetics and transport. In the low activation energy mixture, large-scale vortical structures are observed to occur downstream of the induction distance; these structures do not have a noticeable effect on the reaction. In higher effective activation energy mixtures, a thin transverse ignition front develops near the interface between the two gas streams and results in a combustion structure decoupled from the entrainment region. The decoupling leads to attenuation of the instability growth rate when compared to frozen calculations, and a reduced heat release in the high vorticity region. The analysis indicates the instability plays a modest role in ignition events for high activation energy mixtures. The formation of localized explosions observed in high activation energy systems is instead linked to the impossibility of a one-dimensional reactive combustion wave supported by the injection conditions. In the absence of curvature effects and stream-tube divergence, a system of shock waves is formed which spreads the ignition to the cold gas stream.

Description: A dinoflagellate cyst record with strong Mediterranean/Lusitanian affinities is described from marine deposits of Eemian age (Last Interglacial; Late Pleistocene) at Ristinge Klint, Denmark, revealing new information about the hydrographic evolution of the southwestern Baltic Sea. A revised correlation of the pollen record at Ristinge Klint with that of the annually laminated site at Bispingen in northern Germany provides temporal control. Approximately the first quarter of Eemian time is represented. A marine ingression into a lake took place during theQuercusrise, about 300 years into the interglacial, and is marked by low (

Description: The Tiñu Formation of Oaxaca State is the only fossiliferous lower Palaeozoic unit between the Laurentian platform in northwest Mexico and Gondwanan successions in Andean South America. The Tiñu traditionally has been referred to the Lower Ordovician (Tremadoc) and regarded as having a provincially mixed fauna with Laurentian, Avalonian, and Gondwanan elements. Bio- and lithostratigraphic re-evaluation demonstrates that the Tiñu is a Gondwanan, passive margin succession. It includes a lower, thin (to 16 m), condensed, uppermost Cambrian Yudachica Member (new). The Yudachica nonconformably overlies middle Proterozoic basement as a result of very high late Late Cambrian eustatic levels. The Yudachica changes from storm-dominated, but slightly dysoxic, shelf facies (fossil hash limestone and shale) in the south to an upper slope facies with debris flows 50 km to the north. Three biostratigraphically distinct depositional sequences comprise the Yudachica. The Yudachica has Gondwanan-aspect trilobites with low-diversity conodonts characteristic of unrestricted marine/temperate facies. The upper Tiñu, or Río Salinas Member (new), is a Lower Ordovician (Tremadoc) depositional sequence that records strong early, but not earliest, Tremadoc eustatic rise marked by graptolite- and olenid-bearing dysoxic mudstones. Higher strata shoal upward into shell-hash limestones and proximal tempestite sandstones with upper lower Tremadocian unrestricted marine/temperate conodonts. New taxa includeOrminskia rexroadaeLanding gen. et sp. nov. from theCordylodus andresiZone; this euconodont is related to hyaline coniform genera best known from Ordovician tropical platform successions.Cornuodus?clarkeiLanding sp. nov. resembles the coeval, upper lower Tremadoc tropical speciesScalpellodus longipinnatus(Ji & Barnes).

Description: Mean flow profiles, skin friction, and integral parameters for boundary layers developing naturally over a wide variety of fully aerodynamically rough surfaces are presented and discussed. The momentum thickness Reynolds numberReθextends to values in excess of 47000 and, unlike previous work, a very wide range of the ratio of roughness element height to boundary-layer depth is covered (0.03 0.5). Comparisons are made with some classical formulations based on the assumption of a universal two-parameter form for the mean velocity profile, and also with other recent measurements. It is shown that appropriately re-written versions of the former can be used to collapse all the data, irrespective of the nature of the roughness, unless the surface is very rough, meaning that the typical roughness element height exceeds some 50% of the boundary-layer momentum thickness, corresponding to about$h/delta,{widetilde{〉}},0.2$.

Description: The effects of initial conditions on grid turbulence are investigated for low to moderate Reynolds numbers. Four grid geometries are used to yield variations in initial conditions and a secondary contraction is introduced to improve the isotropy of the turbulence. The hot-wire measurements, believed to be the most detailed to date for this flow, indicate that initial conditions have a persistent impact on the large-scale organization of the flow over the length of the tunnel. The power-law coefficients, determined via an improved method, also depend on the initial conditions. For example, the power-law exponentmis affected by the various levels of large-scale organization and anisotropy generated by the different grids and the shape of the energy spectrum at low wavenumbers. However, the results show that these effects are primarily related to deviations between the turbulence produced in the wind tunnel and true decaying homogenous isotropic turbulence (HIT). Indeed, when isotropy is improved and the intensity of the large-scale periodicity, which is primarily associated with round-rod grids, is decreased, the importance of initial conditions on both the character of the turbulence andmis diminished. However, even in the case where the turbulence is nearly perfectly isotropic,mis not equal to −1, nor does it show an asymptotic trend inxtowards this value, as suggested by recent analysis. Furthermore, the evolution of the second- and third-order velocity structure functions satisfies equilibrium similarity only approximately.

Description: The fragmentation of a laminar undulating liquid sheet flowing in quiescent air is investigated. Combining various observations and measurements we propose a sequential atomization scenario describing the overall sheet–drop transition in this configuration. The undulation results from a controlled primary Kelvin–Helmholtz instability. As the liquid travels through the undulating pattern, it experiences transient accelerations perpendicular to the sheet. These accelerations trigger a secondary instability responsible for the amplification of spanwise thickness modulations of the sheet. This mechanism, called the ‘wavy corridor’, is responsible for the sheet free edge indentations from which liquid ligaments emerge and break, forming drops. The final drop size distribution is of a Gamma type characterized by a unique parameter independent of the operating conditions once drop sizes are normalized by their mean.

Description: High levels of an exotic amino acid (Aib) and enhanced levels of iridium are evident in sediments pre-dating the Chicxulub impact by several tens of millennia. The source is thought to be debris from the break-up of a giant comet or trans-Neptunian body, a large fragment of which was the 10 km sized impactor that caused the famous iridium spike identified with the K/T boundary. In this paper it is argued that the Aib is not extra-terrestrial but the indicator of exotic pathogenic microfungi that flourished through this era. Its abundance implies a significant role for the fungi in the ecology, in species extinction and in driving evolution as the Tertiary period got underway. Microfungi containing the complex of genes that underlie the synthesis of Aib peptides flourished early on in the K/T transition and attacked species – including dinosaurs – that lacked counter immune mechanisms. Species (including mammals) that possessed or developed effective defence mechanisms won through in the early Tertiary-period flowering of new species. The genetic coding for Aib peptide synthesis might have evolved by natural selection. However, the coincidence in the boundary record between Aib peptides and the rise of iridium is indicative of the Aib blueprint arriving from space, in some carrier-organism or in microfungi themselves.

Description: Pleistocene Lake Bonneville of western Utah contains a variety of spits associated with shorelines and other features that formed between 21,000 and 12,00014C yr BP. Field studies in the low-lying mountain ranges of the central portion of Lake Bonneville identified 17 spits of various types. The spits are connected to small mountain ranges and islands, vary in size from 0.02 to 0.5 km2, and are composed of coarse-grained, well-rounded, poorly-sorted sedimentary material. Sixteen of the 17 spits have a northeasterly to southwesterly orientation implying that winds were from the northwest to northeast, approximately 180° out of phase with modern winds in the eastern Great Basin. Lake Bonneville spit orientation is best explained as the result of persistent northerly winds caused by the high atmospheric pressure cell of the continental ice sheet and passage of low pressure extratropical storms south of the lake. Similar, strong persistent winds are a common feature of modern continental ice sheets and passing low pressure systems. If so, the North American jet stream tracked south of Lake Bonneville as recently as 12,00014C yr BP, well past the height of the last glacial maximum.

Description: Accelerator mass spectrometry (AMS) radiocarbon measurements of organic samples require combustion to obtain CO2 for graphitization. Furthermore, determination of δ13C values is required in order to correct the 14C age due to carbon isotope fractionation effects. δ13C analysis is commonly carried out by stable isotope mass spectrometry because most applications demand high-precision δ13C values in addition to the requirements of 14C dating. A simplifying step is therefore to combine the combustion for stable isotope analysis with cryogenic trapping of CO2 for AMS graphite targets. Presented here is a simple CO2 trapping device based on a modified Gilson 220XL sampling (manifold) robot coupled to the inlet manifold system of a GV Instruments IsoPrime stable isotope mass spectrometer. The system is capable of batch combustion and analysis of up to 40 samples and is under full computer control by the mass spectrometer software. All trapping parameters such as flush time prior to trapping and total trap time are adjustable through the standard software user interface. A low 14C activity of background materials and high precision and accuracy of stable isotope analysis of carbon and nitrogen are demonstrated.

Description: We present an overview of accelerator mass spectrometry (AMS) radiocarbon sample preparation and measurements, describing the technical upgrades that now allow us to routinely obtain 0.2–0.3% precision for 1-mg carbon samples. A precision of ∼1% on samples with 100 μg of carbon can also be achieved. We have also developed graphitization techniques and AMS procedures for ultra-small samples (down to 0.002 mg of carbon). Detailed time series are presented for large and small aliquots of standards such as NIST OX-I and OX-II; FIRI-C and -D; IAEA-C6, -C7 and -C8; and 14C-free samples.

Description: The effect of nitrous oxide (N2O) upon the graphitization of small (∼40 μg of carbon) CO2 samples at the ANSTO and University of California, Irvine, radiocarbon laboratories was investigated. Both laboratories produce graphite samples by reduction of CO2 over a heated iron catalyst in the presence of an excess of H2. Although there are significant differences between the methods employed at each laboratory, it was found that N2O has no effect upon the reaction at levels of up to 9.3% by volume of CO2. Further, it was systematically determined that more effective water vapor trapping resulted in faster reaction rates. Using larger amounts of the Fe catalyst generally resulted in higher yields or reaction rates (but not both). The effects of changing the type of Fe catalyst on the final yield and reaction rate were less clear.

Description: Marine upwelling along coastal Peru can be intense and variable, making radiocarbon dating marine and coastal systems complex. Historical and proxy records of upwelling along coastal Peru are few, and long-lived species such as corals do not grow in the cold coastal waters. Mollusk shell carbonate, however, can record both the magnitude of the local marine reservoir correction, ΔR, and of seasonal oscillations in the ventilation age of coastal waters. If large, these seasonal oscillations would complicate radiocarbon dating of marine organisms. To examine this possibility, we sampled for δ13C, δ18O, and 14C content a set of pre-bomb Argopecten purpuratus shells collected from coastal Peru during 1908 and 1926. Intrashell variations of up to 216 14C yr were noted, but these were not consistently correlated with seasonal changes in δ18O or δ13C. Only an 11 yr difference was observed in the weighted average ΔR of Callao Bay shells collected during normal (1908) and El Niño (1926) years. Despite the intrashell 14C variation noted, weighted average ΔR values from all 3 sample sites and from normal and El Niño years all overlap at 1 σ. We report ΔR values of 183 ± 18 and 194 ± 23 yr from Callao Bay (12°4′S), 165 ± 24 yr from Salaverry (8°14′S), and 189 ± 23 yr from Sechura Bay (5°45′S).

Description: Changes in soil organic carbon stocks were simulated with the Rothamsted carbon (RothC) model. We evaluated the calculation of a major input variable, the amount of inert organic matter (IOM), using measurable data. Three different approaches for quantifying IOM were applied to soils with mainly recent organic matter and with carbon contribution from fossil fuels: 1) IOM estimation via total soil organic carbon (SOC); 2) through bulk soil radiocarbon and a mass balance; and 3) by quantifying the portion of black carbon via a specific marker. The results were highly variable in the soil containing lignite-derived carbon and ranged from 8% to 52% inert carbon of total SOC, while nearly similar amounts of 5% to 8% were determined in the soil with mainly recent organic matter. We simulated carbon dynamics in both soils using the 3 approaches for quantifying IOM in combination with carbon inputs derived from measured crop yields. In the soil with recent organic matter, all approaches gave a nearly similar good agreement between measured and modeled data, while in the soil with a fossil carbon admixture, only the 14C approach was successful in matching the measured data. Although 14C was useful for initializing RothC, care should be taken when interpreting SOC dynamics in soils containing carbon from fossil fuels, since these reflect the contribution from both natural and anthropogenic carbon sources.

Description: In Siberia, the accumulation of radiocarbon dates from archaeological sites since the 1960s makes it possible to compile a general Paleolithic 14C database, which contains about 440 entries as of late 2005. With these data, we can reveal the main chronological patterns of Paleolithic complexes, with a focus on the late Middle Paleolithic (Mousterian) and Upper Paleolithic. The 14C dates for late Middle Paleolithic industries in Siberia are quite “young,” up to about 30,000–28,500 BP and perhaps ∼27,000 BP. The emergence of the Upper Paleolithic in Siberia took place relatively early compared with Eastern Europe. At about 43,000–35,000 BP, blade-dominated industries existed in the Altai Mountains and Lake Baikal region, and numerous adornments are known from several sites of that age. The late Upper Paleolithic complexes with microblade technology from the Altai Mountains are 14C dated to about 35,000–28,000 BP, and represent the earliest unequivocal evidence of microblade manufacture in northern Eurasia. The end of the Paleolithic in Siberia is related to the appearance of pottery, which indicates the beginning of the Neolithic period. In northern Transbaikal, the earliest pottery complexes are dated to about 12,000–11,000 BP and in the Russian Far East even to ∼13,000 BP, while in most of Siberia they date to approximately 8000–6000 BP. The most important features of the Siberian Paleolithic chronology are: a) the long persistence of Middle Paleolithic complexes, until about 30,000–27,000 BP; b) very early Middle to Upper Paleolithic transition, ∼43,000 BP, closely connected with the emergence of art and symbolic behavior in the earliest Upper Paleolithic at this time; c) the very early origin of microblade complexes, at least at about 35,000 BP; and d) a gradual Paleolithic–Neolithic transition, beginning in the Russian Far East at ∼13,000 BP and in Transbaikal about 12,000–11,000 BP, in most of Siberia at about 8000–6000 BP, and even later in some northern regions.

Description: The Tavan Tolgoi (Five Holy Hills) site, located in Ongon sum, Sukhbaatar aimag, in southeastern Mongolia, consists of about 20 burials. During the preliminary 2004 excavations conducted by the Department of Anthropology and Archaeology, National University of Mongolia, 7 graves were unearthed. In grave 1 (2004), the skeleton of a woman 40 yr old, wearing golden rings with the inscription of a Siberian falcon, was found together with other ornamental artifacts. In grave 2 (2004), a man with a gold-gilded saddle and a horse were buried. Adornments strongly indicate that these burials date to the Great Mongol Empire period and may relate to the Golden Horde lineage of Genghis Khan. Initial accelerator mass spectrometry (AMS) radiocarbon dating of wood from a coffin at burial 2004–6 (Table 1) gave an age of 860 ± 60 BP, and the age of a human bone sample from burial 2004–1 was determined as 890 ± 40 BP. Subsequent excavations yielded 13 samples for 14C dating, and 7 of them have been dated thus far. The calibrated dates were in the range of AD 1130–1250, which is in agreement with Genghis Khan's life span. Artifacts strongly suggest that these burials belong to nobility or members of the royal family. Given that such burials are hard to find, the Tavan Tolgoi site is expected to yield important archaeological and historical information. In this paper, the historical importance of the artifacts recovered is discussed in light of 14C dating and the results of additional scientific analyses.

Description: We present a broad and detailed study of the background components of a liquid scintillation (LS) detector, using a simple laboratory-built system, ICELS. It was specifically designed for radiocarbon dating and is compact and easily transportable (total weight 35 kg). Its flexible LS detector unit has a dome-shaped vial with 3 mL benzene to which 45 mg butyl-PBD is added. The vial sits on the top of a vertical 28-mm-diameter phototube. The gamma radiation, to which the benzene is exposed under varying conditions, was measured by replacing the vial with a 38-mm-diameter Nal crystal. The pulse-height spectra of the 14C LS background and the Nal gamma background were measured in a surface laboratory and in a deep underground counting room with: 1) a lead shield of varying thickness; 2) lead of normal and low 210Pb concentration; 3) phototubes of 2 different types; and 4) varying benzene volume. The beta emission from the face of the tubes was measured with a low-level Geiger counter.

Description: Previous work has suggested that seasonal and interannual upwelling of deep, cold, radiocarbon-depleted waters from the South Atlantic has caused variations in the reservoir effect (R) through time along the southern coast of Brazil. This work aims to examine the possible upwelling influence on the paleoreservoir age of Brazilian surficial coastal waters based on paired terrestrial/marine samples obtained from archaeological remains. On the Brazilian coast, there are hundreds of shell middens built up by an ancient culture that lived between 6500 to 1500 yr ago, but there are few shell middens located on open-coast sites with a known upwelling influence. Three archaeological sites located in a large headland in Arraial do Cabo and Una de Cabo Frio on the southeastern coast of Brazil, with open-ocean conditions and a well-known strong and large upwelling of the Malvinas/Falkland current, were chosen for this study. The 14C age differences between carbonized seed and marine samples varied from 281 ± 44 to 1083 ± 51 14C yr. There are also significant age differences between carbonized seed samples (977 14C yr) and marine samples (200 and 228 14C yr) from the same archaeological layer that cannot be explained by a reservoir effect or an old-wood effect for charcoal. Therefore, the present data from the southeastern Brazilian coast are inconclusive for identifying an upwelling effect on R. To do so, it would be necessary to more precisely define the present–pre-bomb R in upwelling regions, and to analyze paired marine/terrestrial samples that are contemporaneous beyond doubt.

Description: A prominent variant of the compensation hypothesis rests on the premise that increased trade exposure heightens domestic economic volatility, prompting demands for compensation via generous systems of transfers and services. Economic theory suggests that because the expansion of international trade entails integration into larger, deeper, more stable markets, and may entail risk diversification, it may actually promote rather than reduce stability. By the same token, however, economic theory also suggests that smaller economies should experience greater levels of volatility than larger economies, and thus also greater levels of insecurity. The evidence presented here suggests that the level of domestic economic volatility in the developed economies, during the latter half of the twentieth century, may indeed have been driven by the size and depth of markets. And critically, for these countries international trade integration may have eased rather than accentuated domestic economic volatility.

Description: In this paper we analyze the distribution of pay and changing trends of inequality in Argentina and Brazil, illuminating the specific winners and losers, by region and by economic activity (sector). In both countries we find that inequality rose in the neoliberal period, but that it declined following the severe crises of neoliberal policy, in 1993 in Brazil and in late 2001 in Argentina. This period of post-neoliberalism is characterized in both countries by a decline in the economic weight of the financial sector and a recovery of the position of the civil service. In both countries, the rise in inequality leading to the crisis produced an increase in the relative position of the major metropolitan centers; this positional advantage also declined modestly in the post-crisis recovery period.

Description: The introduction to the special issue develops a systematic and theoretically grounded framework for assessing business power in global governance. It is shown that power is said to have shifted from the world of states to the world of business. However, in order to evaluate such a claim first a differentiation of power in its instrumental, structural, and discursive facets is necessary. It is furthermore explained that the strength of such a three-dimensional assessment is that it combines different levels of analysis and considers actor-specific and structural dimensions and their material and ideational sources. Following a short introduction to the more empirical articles is provided summarizing their commonalities and differences.

Description: Research on conflict in resource-rich countries suggests that resource extraction companies contribute to tension but not development. In recent times, public-private partnerships (PPPs) have flourished, in which set up regulation is not against business but in joint cooperation with corporate actors. Yet PPPs are criticized for serving business self-interest and increasing business power rather than the common good. The paper takes the Kimberley Process and the diamond industry as an example to examine the multi-faceted nature of business power when this PPP was negotiated. The core of the argument is that realist-informed perspectives about business power in PPPs and constructivist accounts emphasizing socialization and social learning processes only tell one part of the story. While demonstrating that the diamond industry acted as a both a socializing and socialized agent, the analysis of the different facets of power shows that structural and discursive power were crucial elements in making socialization happen in the first place.

Description: Jiangsu and Zhejiang are of two of China's most prosperous and dynamic provinces. This paper first presents a factual account of two empirical phenomena: 1) FDI has played a more substantial role in the economic development of Jiangsu than in Zhejiang, and 2) ownership biases against domestic private firms in Jiangsu were more substantial than in Zhejiang. The paper hypothesizes that there is a connection between these two empirical phenomena. Specifically, ownership biases against domestic private firms increase preferences for FDI because FDI provides a measure of relative property rights security. Thus a biased domestic private firm has an incentive to move its assets and/or future growth opportunities to the foreign sector. The paper uses two private-sector surveys—one conducted in 1993 and the other in 2002—to provide an empirical test of this hypothesis. Our analysis shows, controlling for a variety of firm-level attributes and industry and regional characteristics, those private firms which perceive ownership biases to be more severe are more likely to form joint ventures with foreign firms.

Description: The concluding article to the special issue critically reflects on arguments and analysis presented in the preceding articles. It argues that globalisation, new forms of private authority and the increased power of transnational business have not generally weakened the state, but rather advanced a business-oriented transformation of statehood. To understand this transformation the article first provides a very short overview of the state-globalisation debate. Subsequently, it deals more explicitly with the state theoretical debate. In particular, it brings together neo-Marxist and post-Weberian conceptualisations in order to address both the social nature of the state and the particular forms and processes by which it is interactively embedded in the economy and society. After an outline of the transformation of statehood and the strategic options for non-state actors, the article concludes with some critical remarks on the future of democratic politics.

Description: Environmental policies of providers of international finance – namely the World Bank, export credit agencies, and Equator Principles banks – provide interesting cases within which to study the power of business not as only an input to the political process or as a constraint on politics, but also as a conduit for both state and non-state actors.This paper shows how targeting financial actors has allowed NGOs to transform their rather weak discursive power base into instrumental power over business actors in other sectors. NGOs have channeled their power through states, consumers, and financial institutions; this has allowed them to augment discursive power over their targets with additional indirect, yet more immediate, forms of structural and instrumental power. As a consequence of both direct and indirect NGO pressure, financial institutions have adopted environmental policies. This article posits a theoretical explanation of the underspecified power relationships in NGO strategies that allow NGOs to exploit weak links in commodity chains for their campaigns.This paper argues that financial institutions wield considerable structural power through their ability to control access to finance. It is particularly this power base which has made them prime targets for NGOs campaigning for the greening of infrastructure development projects. As a consequence of NGO pressure, financial institutions have adopted environmental policies which in turn have provided the World Bank and Equator banks with additional sources of discursive power.

Description: By analyzing the interaction between a business firm and multiple government institutions (including a regulatory agency, an executive and a bicameral legislature), we develop predictions about how firms target their political strategies at different branches of government when seeking more favorable public policies. The core of our argument is that firms will target their resources at the institution that is ‘pivotal’ in the policy-making process. We develop a simple framework, drawing on the political science literature, which identifies pivotal institutions in different types of political environments. We find empirical support for our thesis in an analysis of how U.S. accounting firms shifted their political campaign contributions between the House and Senate in response to the threat of new regulations governing auditor independence during the 1990s.

Description: Recent work suggests that the most fruitful approach to accounting for variations in interest system diversity of any type lies in understanding variations in interest system density (Lowery, Gray and Fellowes 2005). We build on this insight by examining the sources of variation in the substantive diversity of health interests in the American states, focusing on how the densities of several sub-guilds of health interest organizations vary in their responses to changes in the sizes of the constituencies that give rise to them and variations in the policy and political energy supporting their mobilization. We discuss the concept of interest system diversity in the first section of the paper, highlighting its multiple meanings and the limits of prior research. This is followed by a close empirical examination of 14 sub-guilds of state health interest organizations. We conclude by discussing the inherent difficulties of understanding interest system diversity.

Description: We have determined the radiocarbon ages for 40-yr-interval tree rings in 2 fossil trees of the Towada Hachinohe buried forest, northeastern Honshu Island, Japan. The 14C ages range from 13.0 to 13.3 kyr BP (about 15.5 cal kyr BP). The weighted average of the 14C age of the outermost 5 rings is 13,133 ± 33 BP, which can be calibrated to 15,363–15,679 cal BP by using the IntCal04 standard curve (Reimer et al. 2004). The estimated δ14C values range between 265 and 300% and show approximately sinusoidal fluctuation of an indicated ∼200-yr cycle, perhaps reflecting contemporary solar activity change. Comparison between the tree 14C profile and the Cariaco Basin 14C record provides further information on the accurate date of the Towada Hachinohe buried forest and the eruption that produced it. 14C analysis of tree rings from the buried forest may contribute to the construction of a better 14C calibration curve and the elucidation of solar activity change during the last glacial period, as well as possible global and regional impacts of the huge eruption from Towada Volcano.

Description: Two sediment sequences from Big Kyzykul Lake and the Shushenskoe paleolake in the Minusinsk depression, Southern Siberia, were studied by pollen, microfossil, and geochemical analyses, as well as radiocarbon dating. The records indicate the persistence of an arid period between ∼11.7–7.6 cal kyr BP, increased effective moisture since ∼7.6 cal kyr BP, 2 humid impulses at ∼5.1 and 2.8 cal kyr BP separated by a dry interval, and the return to generally drier conditions after ∼1.5 cal kyr BP. This is contrary to the findings noted for the Eurasian temperate zone, but agrees with proxy data reported for arid and semi-arid zones of Central Asia. Reconstructed changes in climate and environment are in good agreement with archaeological data. Almost no evidence of the Mesolithic-Neolithic cultures has been reported for the depression, which is consistent with a dry early and mid-Holocene. Effective moisture started to rise from ∼7.6 cal kyr BP, followed by the beginning of human occupation at ∼6 cal kyr BP. Two maxima of humidity are recorded in the late Holocene, corresponding to the arrival of trees in the depression. No gap was to be found from the Early Bronze to the Iron ages cultures at this time, with the exception of a dry interval at ∼3.6–3.3 cal kyr BP, when the Minusinsk depression was sparsely occupied. The data obtained suggest a close relationship between climate change and cultural dynamics in the steppe zone of Southern Siberia.

Description: High-precision radiocarbon age calibration for different terrestrial samples allows us to establish accurate boundaries for many climatic time series. At the same time, the fluctuations of 14C content reflect solar variability. A bispectrum analysis of long-term series of the 14C content deduced from decadal measurements in tree rings demonstrates the existence of amplitude modulation, with a period of main modulation of ∼2400 yr. In 14C time series for the last 11 kyr, major oscillations are distinguished at 8.5–7.8, 5.4–4.7, 2.6–2.2, and 1.1–0.4 cal kyr BP with ∼2400-yr periodicity. High amplitudes in cosmogenic isotope content with a periodicity of about 2400 yr appear synchronous to cooling events documented in Greenland ice cores, to the timing of worldwide Holocene glacier expansion, and to the periods of lake-level changes. This paper focuses on revealing solar forcing on the Earth's climate and about the nature, significance, and impact of sharp Holocene climate variability on human societies and civilizations.

Description: A basin-scale repeat hydrography of the WOCE-P06 line along approximately 32°S has revealed decadal changes of bomb-produced radiocarbon in the subtropical South Pacific Ocean between 1992 and 2003. Surface δ14C decreased by about 30%. A δ14C decrease was also found in the upper thermocline from the surface to 400 m water depth. In the lower thermocline, from 400 m to 1100 m depth, however, δ14C increased, with a maximum increase of 25% at 700 m depth. This contrast between the upper and lower thermoclines resulted in an overall 10% increase of the specific column inventory of bomb 14C from 1992 to 2003. The global ocean inventory of bomb 14C was assessed based on the specific inventory increase in the subtropical South Pacific Ocean. The meridional distribution of bomb 14C in the early 1990s suggests that the bomb 14C increases observed along 32°S in 2003 were primarily caused by mixing along isopycnals.

Description: Rivers and streams in many parts of the world contain high concentrations of dissolved carbon dioxide, which is lost to the atmosphere by evasion (outgassing). Recent methodological advances now enable the carbon isotopic composition of this evaded CO2 to be measured directly, with early results from peatland streams showing it to be depleted in 13C relative to the atmosphere. The first direct measurements of the radiocarbon age of evaded CO2 for a stream draining a peatland site in the United Kingdom covered a large range, ranging from modern to 〉1400 yr BP. We investigated whether a proportion of the carbon fixed by plants growing adjacent to the stream was derived from evasion; this would have implications for the cycling of carbon in such streams, and lead to riparian plants having older, rather than contemporary, 14C ages.13C analysis of riparian plants at the site suggested that up to 20% of the carbon they fixed was derived from evasion. Although the 14C content of the same samples suggested a similar proportion of plant carbon was derived from evaded CO2, the uncertainty in these estimates was greater. Fixation of evaded CO2 was greatest in plants growing within 1 m of the stream. Although the results for this site suggest a small amount of CO2 lost by evasion from the stream surface is incorporated in plant material, it may be a more important factor at other sites characterized by higher evasion rates, lower wind speeds, and different stream bank morphology.

Description: This paper examines 2 potential sources of the radiocarbon offset between human and terrestrial mammal (horse) bones recovered from Norse (∼AD 870–1000) pagan graves in Mývatnssveit, north Iceland. These are the marine and freshwater 14C reservoir effects that may be incorporated into human bones from dietary sources. The size of the marine 14C reservoir effect (MRE) during the Norse period was investigated by measurement of multiple paired samples (terrestrial mammal and marine mollusk shell) at 2 archaeological sites in Mývatnssveit and 1 site on the north Icelandic coast. These produced 3 new δR values for the north coast of Iceland, indicating a δR of 106 ± 10 14C yr at AD 868–985, and of 144 ± 28 14C yr at AD 1280–1400. These values are statistically comparable and give an overall weighted mean δR of 111 ± 10 14C yr.The freshwater reservoir effect was similarly quantified using freshwater fish bones from a site in Mývatnssveit. These show an offset of between 1285 and 1830 14C yr, where the fish are depleted in 14C relative to the terrestrial mammals. This is attributed to the input of geothermally derived CO2 into the groundwater and subsequently into Lake Mývatn. We conclude the following: i) some of the Norse inhabitants of Mývatnssveit incorporated non-terrestrial resources into their diet that may be identified from the stable isotope composition of their bone collagen; ii) the MRE off the north Icelandic coast during the Norse period fits a spatial gradient of wider North Atlantic MRE values with increasing values to the northwest; and iii) it is important to consider the effect that geothermal activity could have on the 14C activity of samples influenced by groundwater at Icelandic archaeological sites.

Description: Radiocarbon concentrations in the northernmost region of the Japan Sea were observed during the summer of 2002. The averaged surface δ14C (above 100 m depth) was 52 ± 8%, which is significantly higher compared with the values of the Pacific Ocean and Okhotsk Sea. The δ14C in the deep water decreased with density, and the minimum value was −70%. By analyzing 14C and other hydrographic data, we found that i) the Tsushima Warm Current Water reaches to the surface layer in the southern Tatarskiy Strait; ii) deep convection did not occur in the northernmost region, at least not after the winter of 2001–2002; and iii) the bottom water that was previously formed in this region may step down southward along the bottom slope and mix with the Japan Sea Bottom Water. Furthermore, a new water mass characterized by high salinity (〉34.09 psu) was found in the subsurface layer in the area north of 46°N.

Description: Approximately 10% of Hungary is covered by dust sequences of the Quaternary period. Samples have been taken from more than 50 outcrops and boreholes during fieldwork in the past 20 yr. Some 81 bulk samples taken from 27 profiles of the Hungarian loess regions have been analyzed for radiocarbon. Based on the 14C results, loess layers that accumulated between 33,000 and 12,000 BP were selected for further investigation. The sedimentation rates of the 27 loess profiles suggest periods of slow and quick dust accumulation in the Carpathian Basin during the Upper Weichselian period. It seems to us that some soil development and intense weathering periods interrupted the loess development during the Upper Weichselian in Hungary. According to the 14C dates, the estimated average rate of sedimentation was 0.31 mm/yr in both the northern and southern parts of the Carpathian Basin between 33,000 and 12,000 BP.

Description: Lacustrine carbonate deposition in Hungary has been traditionally interpreted as the outcome of the dry, hot climate prevailing between 7500 and 5000 14C yr BP (hereafter BP) (∼6400 and 3800 BC), triggering the partial desiccation of minor ponds and lakes. A comparative analysis of 5 14C results from the site of Csólyospálos, central Hungary, with those of other Hungarian lacustrine carbonates yielded stunning new results. According to these new dates, carbonate deposition must have initiated much earlier, possibly around 10,000–11,000 BP (9500–11,000 BC) in the Carpathian Basin. Furthermore, the formation of lacustrine carbonates must have come to an end at very different times in different parts of the basin, contrasting previous views on the uniform and synchronous cessation of lacustrine carbonate formation in Hungary. According to the newest results presented here, carbonate deposition in the southern and southeastern parts of the basin ceased around 6000 BP (∼4900 BC). Meanwhile, in the central parts, deposition continued as long as the terminal Bronze Age (∼1300 BC).

Description: High-resolution radiocarbon calibration for the last 14,000 cal yr has been developed in large part using European oaks and pines. Recent subfossil wood collections from the Great Lakes region provide an opportunity to measure 14C activity in decadal series of rings in North America prior to the White Mountains bristlecone record. We developed decadal 14C series from wood at the classic Two Creeks site (∼11,850 BP) in east-central Wisconsin, the Liverpool East site (∼10,250 BP) in northwestern Indiana, and the Gribben Basin site (∼10,000 BP) in the Upper Peninsula of Michigan. Initial AMS dates on holocellulose produced younger-than-expected ages for most Two Creeks subsamples and for a few samples from the other sites, prompting a systematic comparison of chemical pretreatment using 2 samples from each site, and employing holocellulose, AAA-treated holocellulose, alpha-cellulose, and AAA-treated whole wood. The testing could not definitively reveal the source of error in the original analyses, but the “best” original ages together with new AAA-treated holocellulose and α-cellulose ages were visually fitted to the IntCal04 calibration curve at ages of 13,760–13,530 cal BP for the Two Creeks wood, 12,100–12,020 cal BP for Liverpool East, and 11,300–11,170 cal BP for Gribben Basin. The Liverpool East age falls squarely within the Younger Dryas (YD) period, whereas the Gribben Basin age appears to postdate the YD by ∼300 yr, although high scatter in the decadal Gribben Basin results could accommodate an older age nearer the end of the YD.

Description: Oxbow lake deposits of the Neris River at the Valakupiai site in Vilnius (Lithuania) have been studied by different methods including radiocarbon dating. A timescale was attained for the development of the oxbow lake and climatic events recorded in the sediments. 14C dates obtained for 24 samples cover the range 990–6500 BP (AD 580 to 5600 BC). Medieval human activity was found in the upper part of the sediments. Mollusk fauna found in the basal part of the terrace indicate contact between people living in the Baltic and the Black Sea basins. Mean rates were calculated for erosion of the river and for accumulation during the formation of the first terrace.

Description: At 2 locations of the southern Spanish coast, we found indications for tsunamis induced by submarine earthquakes. Charcoal, which we sampled in 2 stratified, assumed tsunamigenic sediment (tsunamites) layers at the shore outside the ruins of the Roman city of Baelo Claudia, close to the Strait of Gibraltar (province of Cádiz), and radiocarbon dated, surprisingly turned out to be of identical age, about 465 BP. In the laguna of the Cabo de Gata (province of Almería), we found possible remnants of tsunamites in drilling cores deposited above organic sediments, 14C dated as 680 ± 30 BP.

Description: Time series are presented of radiocarbon and 13C contents in atmospheric carbon dioxide over eastern Europe (southern Poland), covering the periods 1983–1994 and 2000–2004. The carbon isotope composition was measured in biweekly composite samples of atmospheric CO2, collected about 20 m above the local ground level. The data for 2 observational sites are presented: i) city of Kraków (50°04′N, 19°55′E; 220 m asl; for 1983–1994 and 2000–2004); and ii) Kasprowy Wierch, Tatra Mountains (49°14′N, 19°56′E; 1989 m asl; for 2000–2004). The latter site is considered a regional reference station, relatively free of anthropogenic influences. During the period 1983–1994, observations in the Kraków area revealed a gradual decrease of 14C content with a broad minimum around 1991 and a small increase by about 10% in the subsequent years. δ13C also changes with time, showing a decreasing trend from approximately −9.6% in 1983, with a slope of −0.02%/yr. The observed trends for both isotopes coincide well with a substantial reduction of coal consumption in Poland and partial replacement of coal by natural gas, especially in urban regions. After 2000, the δ13C slightly increases, reaching a mean value of −10% in 2004, while Δ14C is below the reference level by ∼3.5%. Observations at Kasprowy Wierch (regional reference station) also reflect a diminishing input of fossil carbon into the regional atmosphere. The fossil component in atmospheric CO2, calculated with the aid of 14C data available for the 2 study periods, shows a reduction of anthropogenic input by a factor of 2, which is confirmed by annual statistics of coal consumption.

Description: New records of δ13C and Δ14C values in annual rings of pine and oak from different sites around the world were obtained with a time resolution of 1 yr. The results obtained for Europe (Poland), east Asia (Japan), and South America (Peru) are presented in this paper. The δ13C and radiocarbon concentration of α-cellulose from annual tree rings of pine and of the latewood of oak were measured by both accelerator mass spectrometry (AMS) and liquid scintillation spectrometry (LSC). The values of 14S, which represent decreasing 14C concentrations caused by the emission of CO2 from fossil fuel use (Suess effect; Suess 1955), were calculated for each site. Low average 14S (about −0.4 to 0.8%) values for clean areas and high values (about 3.4–3.6%) for industrial and/or urbanized areas were noted. Records of the δ13C values obtained for pine and oak from Poland were used to reconstruct climate changes during the last 400 yr. The results clearly indicate the climate cooling during the periods of the Maunder minimum (1645–1715) and the Dalton minimum (1790–1820). The anti-correlation between the δ13C and Δ14C records during those 2 periods is clear if the 14C record is shifted toward older ages by 24 yr.

Description: Ajdovska Jama (The Pagan's Cave) in southeast Slovenia lies within the catchment of the River Sava, a major tributary of the Danube. The site is well known for its Neolithic burials and has been excavated to a high standard on various occasions since 1884. The human remains at the site occurred as distinct clusters of mainly disarticulated bones belonging to at least 31 individuals. Hitherto, dating of the burials has been based on the associated archaeological finds, including a few low-precision radiometric radiocarbon measurements on charred plant material. In the present study, bones from 15 individuals were subsampled for accelerator mass spectrometry (AMS) and stable isotope analyses. These comprised adults and children from 3 of the clusters. The results of the study indicate that the burials all belong to a relatively short time interval, while the stable isotope data indicate a mixed diet based on C3 plant and animal food sources. These interpretations differ somewhat from those of previous researchers. The AMS 14C and stable isotope analyses form part of a wider investigation of dietary and demographic change from the Mesolithic to the Iron Age in the Danube Basin.

Description: Sohr Damb/Nal, the type site of the Nal complex, is located in Balochistan, Pakistan. After 1 season of excavation by H Hargreaves in 1924, which made the polychrome Nal pottery widely known, no further work took place until the Joint German-Pakistani Archaeological Mission to Kalat resumed excavations in 2001. So far, 4 seasons of excavations have been undertaken, which have revealed 4 periods of occupation, dated from about 3800 to 2000 BC. The well-stratified assemblages provide new insights into cultural processes and developments, and enhance the comparative frameworks through typological series and a comprehensive set of radiocarbon dates. This information is essential for assessing cross-cultural relations and the date of urbanization. In this paper, the 14C dates from Sohr Damb/Nal are presented and their cultural context is discussed. Period III has several links to sites in Pakistan, Afghanistan, and Iran, such as Miri Qalat IIIb–c, Mehrgarh VI–VII/Nausharo I, Quetta III, Mundigak IV, and Shahr-e Sokhta II–III. Period IV represents a Kulli-Harappan occupation, which is dated to the second half of the 3rd millennium BC.

Description: The famous single-arch stone bridge over the Neretva River in Mostar, Bosnia and Herzegovina, was built in 1557–1566 by the order of Ottoman Sultan Süleyman the Magnificent. During the reconstruction of the Old Bridge, which was destroyed during the war in 1993, remnants of 2 older wooden bridges were found at the same location. Six wood and 2 charcoal samples were dated using the radiocarbon method. Wood samples with visible tree rings were taken for dendrochronological dating. The results point to several periods of construction and repairs to the bridge and the towers at each end, spanning from the 12th to 18th century. Calibrated 14C and dendrochronological ages were in good agreement.