The inertial range “outer scale” and optical turbulence

The crucial parameter used to calculate turbulence effects upon light waves propagating through the atmosphere is known as the structure constant, . As Tatarski has shown, this parameter depends upon the “outer scale” of the inertial sub-range of the turbulence. Recently there have been successful predictions of astronomical “seeing” conditions at Mauna Kea Astronomical Observatory which have increased interest in this subject and in the use of the so-called “Dewan Optical Turbulence Model”. In the case of the Air Force, there has been a longstanding need for such optical turbulence prediction, especially in the stratosphere. In the past researchers have used a relation due to Tatarski, (which plays a prominent role in this model) in order to deduce values of the “outer scale” from measurements. When doing this, they have been surprised to find values very much smaller than expected. The goal of the paper is to explain this unexpected result. As we will show, this result can be explained by two factors: (a) the average turbulent layer thicknesses are smaller than originally believed, and, more importantly, (b) only a minor fraction of the stratosphere is turbulent. In order to arrive at this conclusion, we used the high-resolution (10 m) wind profiles that were originally used to formulate the previously mentioned optical turbulence model. The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

A statistical model to forecast the profile of the index structure constant $${C_N^2}$$

A statistical study was performed over 145 profiles of meteorological balloons, equipped with microthermal sensors, from the ground to the midstratosphere. This study put into evidence the lognormal distribution of the fluctuations of the structure constant of temperature, the fluctuations of the buoyancy force and the vertical shear of the horizontal wind speed. We show that these quantities, computed over a large scale (100 m), are correlated up to the midstratosphere. A model is adjusted to estimate the optical turbulence strength from the macroscopic meteorological parameters. The model performances, to estimate and forecast the strength and the altitude of the optical turbulent layers, are quantified and compared with other already defined models. This model was shown to have the best performances and put into evidence a new relation to describe the turbulence on a large scale.     

Comments on “turbulent velocity profile in fully-developed open channel flows”

In a recent paper published in this journal, Bonakdari et al. (Environ Fluid Mech 8:1–17, 2008) presented a new formulation of the vertical velocity profile in the central portion of steady fully developed turbulent open-channel flows which is based on an analysis of the Navier–Stokes equations. The predicted mean streamwise velocity profiles presented by the authors represent dip phenomenon (maximum velocity below the free surface). The discusser would like to point out some contradictions and an error in the main demonstration. This error appears in the right-hand side term of the main Eq. 22 (Environ Fluid Mech 8:1–17, 2008) in the parameter α. Our corrected demonstration shows that α is equal to and not to as obtained by the authors (Environ Fluid Mech 8:1–17, 2008). The values of α used by the authors are under-estimated by 1, this difference has a significant effect on velocities and therefore will involve a difference between measured and predicted velocity profiles. Finally, the assumption of at the free surface is in contradiction with predicted velocity profiles.     

A performance comparison between nonlinear similarity functions in bulk parameterization for very stable conditions

Three sets of nonlinear similarity functions for strong stability are selected to compare their performance in bulk parameterization. To uncover their advantages and disadvantages, theoretical and measurement analyses are made with four profile metrics and the Deacon number technique. Main disadvantages include the negligence of the different transfer efficiency between momentum and heat, the flux cutoff due to the upper limit in gradient Richardson number (Ri) and the ignorance of limited stability range where the dimensionless gradient functions ( and ) approach constants. Accordingly, three suggestions are made for future improvement. First, the functions for wind velocity and potential temperature should have the same function form, but with different coefficients. Second, and should approach constants only within a certain stability range. Third, the limit value in Ri should be avoided to widen their applicability in flux modeling. Furthermore, quantitative comparisons in transfer coefficients for moment and sensible heat (C _D and C _H) are made among the similarity functions in the bulk Richardson number (Ri _B) range 0 < Ri _B < 1. Generally, significant discrepancy is found, which may approach a factor of two and three at large Ri _B in C _D and C _H, respectively. Finally, a new recommendation is made to one of the three sets, mainly because of its ability to predict C _D and C _H that decrease rather slowly in very stable conditions.     

Trophic variation with length in two ommastrephid squids, <Emphasis Type="Italic">Ommastrephes bartramii</Emphasis> and <Emphasis Type="Italic">Sthenoteuthis oualaniensis</Emphasis>

From 1998 to 2001 a total of 200 Ommastrephes bartramii (27 paralarvae) and 170 Sthenoteuthis oualaniensis (14 paralarvae) were sampled from the Central North Pacific. One group of non-paralarval O. bartramii (n = 30) was sampled from farther northwest in 1996. The δ¹⁵N of mantle muscle of non-paralarval O. bartramii ( = 12.4‰) was significantly greater than that of non-parlarval S. oualaniensis ( = 8.1‰) (P < 0.001). The δ¹⁵N of whole paralarvae of O. bartramii ( = 6.4‰) was not significantly different than parlarvalae of S. oualaniensis ( = 6.1‰) (P = 0.528). There was no significant difference between the mantle muscle δ¹⁵N values of male (n = 95, = 13.3‰) and female (n = 18, = 12.9 ‰) O. bartramii greater than 300 mm mantle length (ML) (P = 0.15). There was also no significant difference between the mantle muscle δ¹⁵N values of male (n = 15, = 7.2‰) and female (n = 26, = 7.3 ‰) S. oualaniensis in the same size range (P = 0.41). Overall there was a distinct logistic increase in δ¹⁵N with mantle length for O. bartramii, whereas S. oualaniensis showed an exponential increase in δ¹⁵N with mantle length that was stronger within individual years than with all samples combined. In general, adult O. bartramii are more than a trophic level above S. oualaniensis (4.3‰, 1.3 TLs). Because of the nature of the sampling protocol, this study could not separate spatial and temporal effects on the δ¹⁵N signals from each squid species. This study demonstrates the ability of stable isotope analyses to differentiate trophic levels between squid species as well as track trophic changes across size ranges from paralarvae to adults. Additional research is needed to validate these trophic changes across size within individuals.     

Hypoxia increases the behavioural activity of schooling herring: a response to physiological stress or respiratory distress?

Atlantic herring, Clupea harengus, increase their swimming speed during low O₂ (hypoxia) and it has been hypothesised that the behavioural response is modulated by the degree of “respiratory distress” (i.e. a rise in anaerobic metabolism and severe physiological stress). To test directly whether a deviation in physiological homeostasis is associated with any change in behavioural activity, we exposed C. harengus in a school to a progressive stepwise decline in water oxygen pressure and measured fish swimming speed and valid indicators of primary and secondary stress (i.e. blood cortisol, lactate, glucose and osmolality). Herring in hypoxia increased their swimming speed by 11–39% but only when was <8.5 kPa and in an unsteady (i.e. declining) state. In parallel with the shift in behaviour, plasma cortisol also exhibited an increase with plasma osmolality was subject to a transient rise at 8.5 kPa and plasma glucose was generally reduced at However, without any rise in anaerobically derived lactate levels, there was no evidence of respiratory distress at any set We show that a shift in physiological homeostasis is indeed linked with an increase in the swimming speed of herring but the physiological response reflects a hypoxia-induced shift in metabolic fuel-use rather than respiratory distress per se. The significance of this behavioural–physiological reaction is discussed in terms of behavioural-energetic trade-offs, schooling dynamics and the hypoxia tolerance of herring.     

Oxygen consumption of the semi-terrestrial crab <Emphasis Type="Italic">Pachygrapsus marmoratus</Emphasis> in relation to body mass and temperature: an information theory approach

Pachygrapsus marmoratus is a semi-terrestrial crab and the most common grapsid crab in the intertidal belt of rocky shores throughout the Mediterranean Sea, Black Sea and northeastern Atlantic. In this study, the combined effects of temperature (T), body mass (M), and sex (S) on the routine oxygen consumption rate (R) in P. marmoratus were quantified. The blotted wet body mass of the specimens ranged between 43 mg and 18.0 g, and five test temperatures were used between 13.5 and 28.0°C. Six candidate models that reflected different assumptions regarding the dependence of R on S and T were compared. Model selection was based on Kullback–Leibler’s information theory and Akaike’s information criterion (AIC). The model had the highest support by the data (E is the activation energy, B = 8.618 × 10⁻⁵ eV K⁻¹ is Boltzmann’s constant, T _a is the absolute temperature in Kelvin, and b the allometric scaling exponent); for P. marmoratus it was found that No sex dependence of R was supported by the data. Following a multi-model inference (MMI) approach, the mean (± SE) allometric exponent was 0.750 (± 0.013) having a 95% (bootstrap) confidence interval of 0.726–0.774. Thus, it was established that P. marmoratus follows Kleiber’s 3/4 law, as seems to be generally true for intertidal crabs. The allometric exponent was independent of temperature as has also been reported for many other marine invertebrates (at normal temperatures). Q ₁₀ values were relatively low, indicating wide thermal tolerance of the species. Model selection based on information theory is recommended for respiration studies, as an effective method in finding a parsimonious approximating model. MMI by model averaging, based on Akaike weights, is an effective way to make robust parameter estimations and deal with model selection uncertainty.     

OHP–APT 2002 gravity wave campaign: waves,turbulence and forecasts

In July 2002, a multi-sensor campaign was conducted in southern France to investigate the hypothesized connection between gravity waves and optical turbulence. A generalized scidar (GS) was mounted on the 1.93 m diameter telescope at l’Observatoire de Haute-Provence (OHP). The GS provides continuous profiles of optical turbulence with 300 m vertical resolution from telescope altitude up to 25 km. Thermosondes, which provide in situ measurement of optical turbulence by measuring temperature variance, were launched at the OHP site and at a site approximately 20 km west-northwest of OHP. Gravity wave activity was deduced from temperature and wind velocity measured by radiosondes, which are part of the thermosonde system. In this paper, gravity waves were analyzed using techniques for simple two-dimensional mountain waves, with only fair results. Mesoscale models were run at moderately high resolution for the period. The forecasts were analyzed for wave activity, and a post processor model was used to diagnose the optical turbulence. Mountain waves were evident in the forecast, but quantitative comparison showed the forecast to be inadequate in predicting wave strength. The forecast optical turbulence was in fair agreement with measurements with notable exceptions. This paper is declared a work of the U.S. Government and is not subject to copyright protection in the United States. The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged. Presented at the Turbulence and Waves in Stably Stratified Atmospheric Shear Flows: Measurements, High Resolution Simulations and Numerical Prediction Challenges Conference, 13–15 September 2004, Lighthill Institute of Mathematical Sciences, University College London, London, UK.     

Energy dissipation estimates in oscillating grid setup: LDV and PIV measurements

The dissipation of turbulent kinetic energy has been increasingly used as a scaling parameter to integrate microbiological accrual and metabolic rates with fluid-flow motion in natural and engineered aquatic ecosystems. The estimation of turbulent kinetic energy under field conditions and the generation of energy dissipation rates under controlled laboratory conditions with microbiological organisms are necessities required to integrate environmental/ecological laboratory protocols with a moving fluid in the environment. Turbulent fluid-flow conditions were generated in an oscillating grid setup, and turbulence variables were quantified using laser-Doppler velocimetry (LDV) and particle image velocimetry (PIV) measuring techniques. The rate of dissipation of the turbulent kinetic energy in the setup ranged from 10⁻⁹ to 10⁻⁴ m²/s³ and was similar to the levels of energy dissipation commonly reported in engineered and natural aquatic ecosystems. Time-averaged velocities were close to zero with the root-mean-square velocity ratios about 1, indicating nearly isotropic fluid-flow conditions in the setup. The velocity spectra, obtained by stationary LDV measurements for the vertical and horizontal velocity components across the setup revealed the existence of inertial subrange with the frequency power scaling law of “ω ^−5/3.” The estimated Eulerian frequency spectrum followed the theoretical functional relation and confirmed the applicability of inertial dissipation method for the estimation of turbulent kinetic energy dissipation rates. PIV was used for a direct estimation of dissipation by evaluating spatially distributed velocity gradients. The direct dissipation estimate in conjunction with the estimated Eulerian frequency spectrum provided evaluation of a “universal” constant, α, commonly used for the estimation of an energy dissipation rate over the inertial subrange of the Eulerian spectrum. The results demonstrated a range of values, rather than a universal constant, of α with a lognormal probability distribution for vertical and horizontal velocity components. In order to encompass a 0.955 probability range under the lognormal distribution the universal constant, α, should be in the range 2.91 ≥ α _u ≥ 0.43 and 4.44 ≥ α _w ≥ 0.42 for horizontal and vertical velocity components, respectively.     

Effect of tidal currents,seston, and bottom sediments on growth of Mercenaria mercenaria: results of a field experiment

Tidal currents, seston, and sediments separately influence growth of the hard clam, Mercenaria mercenaria, but it is uncertain how these factors may interact. A 3×3 factorial field experiment, carried out in Great Sound, a coastal lagoon in Southern New Jersey, USA, between May and September 1986, determined the relative effects of three sediment types and three site-specific seston/tidal current regimes on the individual growth of M. mercenaria. Analysis of variance of the change in shell length after 15 wk (differences in initial and final lengths) demonstrated a significant difference (P=0.0064) in growth among sites, but no significant differences (P=0.1331) for growth in different sediments, although trends were evident. Effects of sites were independent of sediment type (P=0.2621). Shell growth rates differed by 10.7% between the slowest and fastest sites, but only differed by 5.7% between sediment types, with fastest growth in sand and slowest in mud. Tidal current speeds and four measures of seston (chlorophyll a, particulate inorganic and organic matter, PIM and POM, and energy content) were measured >20 times in near-bottom waters at each site. Horizontal fluxes of POM exhibited higher correlation coefficients with growth rates, than did seston concentrations or current speeds alone. We attribute significant site differences to differences in horizontal seston fluxes fluxes among sites. We suggest that horizontal seston fluxes may be a major factor affecting individual growth of suspension-feeding bivalves.No. 216 of the Jackson Estuarine Laboratory Contribution Series     

Photosynthetic characteristics of nanoplankton and picoplankton from the surface mixed layer

Nanoplankton and picoplankton primary production has been studied at two oceanic stations in the Porcupine Sea-bight and at one shelf station in the Celtic Sea. At both sites, low wind conditions in June and July 1985 resulted in greatly reduced vertical turbulent mixing and a secondary, temporary thermocline developed in what is usually a well-mixed surface layer; as a result, there was physical separation of the phytoplankton within two zones of the surface mixed layer. The photosynthetic characteristics of three size fractions (>5 m, <5 to >1 m and <1 to >0.2 m) of phytoplankton populations from the two zones have been measured. Phytoplankton was more abundant at the oceanic stations and chlorophyll a values were between 1.3 and 2.2 mg chlorophyll a m^-3, compared with 0.3 to 0.6 mg chlorophyll a m^-3 at the shelf station; at both stations, numbers of cyanobacteria were slightly higher in the lower zone of the surface mixed layer. There was no effect of the temporary thermocline on the vertical profiles of primary production and most phtosynthesis occurred in the surface 10 m. Photosynthetic parameters of the three size fractions of phytoplankton have been determined; there was considerable day-to-day variation in the measured photosynthetic parameters. Assimilation number (P _m ^B ) of all >5 m phytoplankton was lower for the deeper than for the surface populations, but there was little change in initial slope (a ^B ). The small oceanic nanoplankton (<5 to >1 m) showed changes similar to the >5 m phytoplankton, but the same size fraction from the shelf station showed changes that were more like those shown by the picoplankton (<1 m) viz, little change in P _m ^B but an increase in a ^B with depth. Values of a ^B were generally greater for the picoplankton fraction than for the larger phytoplankton, but values of adaptation parameter (I _k )(=P _m ^B /) were not always less. There was little evidence to support the hypothesis that these populations of picoplankton were significantly more adapted to low light conditions than the larger phytoplankton cells. When photosynthetic parameters of the picoplankton were normalised to cell number (P _m ^C /a ^C ) rather than chlorophyll a, P _m ^C was comparable to other published data for picoplankton, but a ^C was much lower. The maximum doubling time of the picoplankton at saturating irradiance is calculated to be ca. 8.5 h for the oceanic population and ca. 6.2 h for the shelf population.     

Spatial and temporal distribution of atmospheric mercury species over the Adriatic Sea

Field measurements of atmospheric mercury and related species were carried out during an intensive cruise campaign performed over the Adriatic sea from October 26th to November 12th, 2004 on board the R/V Urania. Hg⁰ ranged between 0.8 and 3.3 ng/m³ with an average of 1.6 ± 0.4 ng/m³ over the entire period. Hg(II) concentrations ranged from 0.1 to 62.8 pg/m³ with an average of 6.7 ± 11.7 pg/m³ whereas Hg-p levels were in a range of 0.04 and 51 pg/m³ with an average of 4.5 ± 8 pg/m³. Higher Hg⁰ and Hg-p concentrations were observed in the Gulf of Venice and Gulf of Trieste due primarily to air masses transported from the mainland reflecting the contribution from anthropogenic sources. In contrast, higher Hg(II) concentrations observed during the first period of the cruise campaign were likely due to the occurrence of photo-oxidants production which are the main players of the gas phase oxidation of to Hg(II)_(g). These findings have been confirmed by the backward trajectories analysis of air masses crossing the studied area during selected days.     

Distribution and phytoavailability of antimony at an antimony mining and smelting area,Hunan, China

An investigation of the distribution, fractionation and phytoavailability of antimony (Sb) and other heavy metals in soil sampled at various locations in the vicinity of a Sb mine revealed elevated levels of Sb, most certainly due to the mining activities. The concentration of Sb in the soil samples was 100.6–5045 mg kg⁻¹; in comparison, the maximum permissible concentration for Sb in soil in The Netherlands is 3.5 mg kg⁻¹, and the maximum permissible concentration of pollutant Sb in receiving soils recommended by the World Health Organization is 36 mg kg⁻¹. The soil sampled near the Sb mine areas had also contained high concentrations of As and Hg. Root and leaf samples from plants growing in the Sb mine area contained high concentrations of Sb, with the concentration of Sb in the leaves of radish positively correlating with Sb concentrations in soil. The distribution of Sb in the soil showed the following order: strongly bound to the crystalline matrix > adsorbed on Fe/Mn hydrous oxides, complexed to organic/sulfides, bound to carbonates > weakly bound and soluble. Solvents showed varying levels of effectiveness in extracting Sb (based on concentration) from the soil, with , in decreasing order. The concentration of easily phytoavailable Sb was high and varied from 2.5 to 13.2 mg kg⁻¹, the percentage of moderately phytoavailable Sb ranged from 1.62 to 8.26%, and the not phytoavailable fraction represented 88.2–97.9% of total Sb in soils.     

The response of Atlantic cod, <Emphasis Type="Italic">Gadus morhua</Emphasis>, to progressive hypoxia: fish swimming speed and physiological stress

Atlantic cod, Gadus morhua, were exposed to a progressive stepwise decline in water oxygen pressure Fish swimming speed and indicators of primary and secondary stress (e.g. blood cortisol and lactate) were measured to assess whether a severe shift in physiological homeostasis (i.e. stress) preceded any change in behaviour or vice versa. Swimming speed increased by 18% when was reduced rapidly from 19.9 kPa to 13.2 kPa and was interpreted as an initial avoidance response. However, swimming speed was reduced by 21% at a moderate level of steady (8.4 kPa) and continued to drop by 41% under progressively deep hypoxia (4.3 kPa). Elevations in plasma cortisol and blood lactate indicated major physiological stress but only at 4.3 kPa, which corresponds to the critical oxygen tension of this species. We propose that the drop in speed during hypoxia aids to offset major stress and is adaptive for the survival of cod in extensive areas of low oxygen.     

Ecophysiology of the benthic amphipod Monoporeia affinis in an open-sea area of the northern Baltic Sea: seasonal variations in oxygen consumption and ammonia excretion

Seasonal variations in the oxygen uptake rate ( ) and ammonia excretion rate ( ) of the benthic deposit-feeding amphipod Monoporeia affinis vere studied in an open-sea population of the northern Baltic Sea during three consecutive years (1991 to 1993). In the early summer an increase in the weightspecific by 22% (standard 1-mg individual) was registered; an even more striking change was noted in the , which showed an increase in the early summer of ca. six times the extremely low winter/spring valnes. The metabolic changes observed could not be attributed to temperature, which remained at steadystate throughout the year, nor could seasonal variations in the biochemical composition of the amphipods explain the results. Therefor, it is evident that quality-food abundance following the sedimentation of the spring phytoplankton bloom induces the observed early-summer elevations in and . In turn, the marked decrease in the towards autumn reflects the gradual depletion of quality-food in the benthic environment. Atomic O: N ratios indicated that at the sampling location, M. affinis uses lipid almost exclusively for metabolic energy throughout the year, except in early summer when the use of lipid and protein was observed to be almost equal. Results also indicated a time-lag between sedimentation of the phytoplankton bloom and metabolic response in the amphipod population.     

The effect of temperature-salinity combinations on survival and growth of juvenile Patiriella pseudoexigua (Echinodermata: Asteroidea)

Adult Patiriella pseudoexigua were collected in October 1989 from Wanlitung, Taiwan and then induced to spawn in the laboratory. Post-metamorphosed juvenile P. pseudoexigua were reared on a diet of benthic algae Navicula sp. at 25°C and salinity (34). Six weeks after metamorphosis, juvenile P. pseudoexigua at ca. 400 m in radius were reared on a diet of benthic algae Navicula sp. at different combinations of temperatures (20, 25, 30°C) and salinities (26, 30, 34) for 40 d. Both temperature and salinity had a significant effect on juvenile survival and growth. Juveniles survived best (>90%) at 25°C and 34 and grew best (to ca. 750 m in radius) at 30°C and 34. Variation in juvenile size was small immediately after metamorphosis and increased with time.     

Spatial patterns and temporal succession in soft-bottom macroinvertebrate assemblages surrounding an ocean outfall on the southern San Diego shelf: relation to anthropogenic and natural events

Spatial patterns and temporal succession in soft-bottom macroinvertebrate assemblages are described for an area of the southern San Diego shelf which is influenced by both anthropogenic and natural events. The study transect covered 19 km along the 60 m-depth contour and spanned the terminus of the Point Loma wastewater outfall (32°40N; 117°17W). Ordination and classification analyses revealed that eight major benthic assemblages occurred in the study area from 1986–1990. These analyses further indicated that: (1) the impact of wastewater discharge on the benthos was asymmetric about the outfall, and (2) primary and secondary impact zones existed within the outfall footprint. The primary impact extended only 0.5 km south, but usually up to 1 km north, of the point of discharge, reflecting the prevailing northerly flow of currents along the San Diego coast. Stations within the primary impact zone usually supported similar faunal assemblages at any point in time. A lessened, or secondary, impact was apparent at stations 1 to 2 km south and 2 km north of the outfall. The two stations 9 km to the north appeared to be outside the outfall footprint. Assemblages in the three zones did not differ greatly in species composition, but did differ with respect to organismal density, relative abundances (dominance hierarchy) of component species, and temporal persistence. Biostimulation of the benthos, manifest as increased species richness and organismal abundance, was apparent in the primary impact zone from 1986–1988, although dominance was relatively low. Dominance increased at these sites during 1989 and 1990 due to large increases in populations of a few small species, such as the polychaete Myriochele sp. M and the ostracod Euphilomedes carcharodonta. Biomass was slightly elevated in the primary impact zone during 1986–1987; however, this trend was not sustained. Communities within the outfall footprint were much less persistent over time than those beyond it. The observed shifts in benthic assemblages within the outfall footprint corresponded to natural oceanographic cycles (i.e., the transition from El Niño to La Niña conditions) and anthropogenic changes (i.e., changes in wastewater output). These observations suggest that the Point Loma outfall has a localized effect on the stability of infaunal communities. Possible destabilization mechanisms are explored. Finally, comparisons with pre-discharge data taken during 1956–1957 revealed that about half of the species present prior to construction of the outfall were still common in 1990, after 27 yr of wastewater discharge, although large increases in species diversity and abundance seem to have occurred at all stations.     

Abundance and biomass of macrobenthos in the vicinity of Carson Submarine Canyon,northwest Atlantic Ocean

The abundance, biomass and mean weight of macrofauna in the Carson Canyon region (Lat. 45°30N, Long. 48°40W) of the Grand Banks were estimated from 40 quantitative infaunal samples taken at 15 stations in June 1980. In contrast to other areas, there were no significant differences in these values between samples taken in and outside of the canyon. Abundance ( = 716 individuals m^-2), biomass ( = 10.9 gm^-2), and mean weight ( = 0.03 g individual^-1) displayed no consistent relationship with depth, but stations grouped according to sediment type had differing dominant taxa. Biomass and abundance at the edge of the Grand Banks are unexpectedly low. Larger animals, including fishes, that forage there must depend more on pelagic than on benthic food sources.     

Nuclear and mitochondrial DNA markers indicate unidirectional gene flow of Indo-Pacific to Atlantic bigeye tuna (<Emphasis Type="Italic">Thunnus obesus</Emphasis>) populations,and their admixture off southern Africa

A sharp genetic break separates Atlantic from Indo-Pacific bigeye tuna (Thunnus obesus) populations, as the frequencies of two major mitochondrial (mt) DNA types ( and ) found in this species are different across the tip of southern Africa. The level of nucleotide divergence between mtDNA types and is of the same order as that between reproductively isolated taxa. To further investigate the genetic structure of bigeye tuna over its distribution range and in the contact zone off southern Africa, bigeye tuna samples collected between 1992 and 2001 (including samples from a previous mtDNA survey) were characterized for four nuclear DNA loci and for mtDNA. Nuclear markers did not support the hypothesis that and mitochondria characterize sibling species. Significant allele-frequency differences at one intronic locus (GH2) and one microsatellite locus (µ208) were found between Atlantic and Indo-Pacific samples, although the level of nuclear genetic differentiation (Weir and Cockerhams =0.025 to 0.042) was much lower than in mtDNA ( =0.664 to 0.807). Probabilistic Bayesian assignment of individuals to a population confirmed that southern African bigeye tuna samples represent a simple mixture of individuals from Atlantic and Indian stocks that do not interbreed, with a higher contribution from Indian Ocean individuals (about 2/3 vs. 1/3).Communicated by M.S. Johnson, Crawley