Parameterization of the roughness length over the sea in forced and free convection

In the condition of free convection, the Charnock relation is inadequate. In this paper we extend the Charnock relation to include the effect of free convection on the roughness length. As a result, the singularity in the Monin–Obukhov similarity theory can be avoided. This paper shows two approaches to derive the roughness length formula in the forced and free convections. The first approach is based on the mixing length theory and the use of the observational data of the vertical velocity variance. We introduce a new vertical velocity scale based on the vertical velocity variance; this velocity variance is well behaved in the atmospheric boundary layer and easy to obtain from field experiments. The second approach is based on the theoretical framework of Sykes et al. (Q R Met Soc, 119: 409–421). From that framework, we develop a theory to obtain the roughness length formula. The results of these two approaches are in agreement with each other. In the past, a multiplication factor associated with free convection was considered to be a constant. This paper shows that this multiplication factor is, in fact, also dependent on the depth of the mixing height. In previous studies, experimental works were often conducted without taking into account the depth of the mixing height. Not taking into account the mixing height in the estimation of the roughness length in free convection would result in an inaccurate estimate of the roughness length and hence the drag coefficient. An erratum to this article can be found at     

Advection and resulting CO2 exchange uncertainty in a tall forest in central Germany

Potential losses by advection were estimated at Hainich Forest, Thuringia, Germany, where the tower is located at a gentle slope. Three approaches were used: (1) comparing nighttime eddy covariance fluxes to an independent value of total ecosystem respiration by bottom-up modeling of the underlying processes, (2) direct measurements of a horizontal CO2 gradient and horizontal wind speed at 2 m height in order to calculate horizontal advection, and (3) direct measurements of a vertical CO2 gradient and a three-dimensional wind profile in order to calculate vertical advection. In the first approach, nighttime eddy covariance measurements were compared to independent values of total ecosystem respiration by means of bottom-up modeling of the underlying biological processes. Turbulent fluxes and storage term were normalized to the fluxes calculated by the bottom-up model. Below a u(*) threshold of 0.6 m/s the normalized turbulent fluxes decreased with decreasing u(*), but the flux to the storage increased only up to values less than 20% of the modeled flux at low turbulence. Horizontal advection was measured by a horizontal CO2 gradient over a distance of 130 m combined with horizontal wind speed measurements. Horizontal advection occurred at most of the evenings independently of friction velocity above the canopy. Nevertheless, horizontal advection was higher when u(*) was low. The peaks of horizontal advection correlated with changes in temperature. A full mass balance including turbulent fluxes, storage, and horizontal and vertical advection resulted in an increase of spikes and scatter but seemed to generally improve the results from the flux measurements. The comparison of flux data with independent bottom-up modeling results as well as the direct measurements resulted in strong indications that katabatic flows along the hill slope during evening and night reduces the measured apparent ecosystem respiration rate. In addition, anabatic flows may occur during the morning. We conclude that direct measurements of horizontal and vertical advection are highly necessary at sites located even on gentle hill slopes.     

The structure of turbulence near a tall forest edge: the backward-facing step flow analogy revisited

Flow disturbances near tall forest edges are receiving significant attention in diverse disciplines including ecology, forest management, meteorology, and fluid mechanics. Current theories suggest that near a forest edge, when the flow originates from a forest into a large clearing, the flow retains its forest canopy turbulence structure at the exit point. Here, we propose that this framework is not sufficiently general for dense forested edges and suggest that the flow shares several attributes with backward-facing step (BFS) flow. Similar analogies, such as rotor-like circulations, have been proposed by a number of investigators, though the consequences of such circulations on the primary terms in the mean momentum balance at the forest clearing edge have rarely been studied in the field. Using an array of three triaxial sonic anemometers positioned to measure horizontal and vertical gradients of the velocity statistics near a forest edge, we show that the flow structure is more consistent with an intermittent recirculation pattern, rather than a continuous rotor, whose genesis resembles the BFS flow. We also show that the lateral velocity variance, v'2, is the moment that adjusts most slowly with downwind distance as the flow exits from the forest into the clearing. Surprisingly, the longitudinal and vertical velocity variances (u'2 and w'2) at the forest edge were comparable in magnitude to their respective values at the center of a large grass-covered forest clearing, suggesting rapid adjustment at the edge. Discussions on how the forest edge modifies the spectra and co-spectra of momentum fluxes, effective mixing length, and static pressure are also presented.     

Model of turbulence penetration into a suspension layer on a sediment bed and effect on vertical solute transfer

The vertical diffusional mass (solute) transfer through a suspended sediment layer, e.g. at the bottom of a lake, reservoir or estuary, by the propagation of velocity fluctuations from above was investigated. The attenuation of the velocity fluctuations in the suspension layer and the associated effect on solute transfer through the suspension layer was simulated. To represent large eddies traveling downstream in water over a high-concentration suspended sediment layer, a streamwise velocity fluctuation moving in downstream direction was imposed along the upper boundary of the suspension layer. Velocity fluctuations and downstream velocity were normalized by the shearvelocity (U_*) at the top of the suspension layer. Streamwise and vertical velocity components inside the suspension layer, were obtained from the 2-D continuity and the Navier–Stokes equations. The persistence of turbulence with depth—as it penetrates from the overlying water into the suspension layer—was found to depend on its amplitude, its period, and on the apparent viscosity of the suspension. The turbulence was found to propagate efficiently into the suspension layer when its frequency is low, and the apparent viscosity of the suspension is high. Effects on vertical mass transfer were parameterized by penetration depth and effective diffusion coefficient, and related to apparent viscosity of the suspension, Schmidt number and shear velocity on top of the suspension layer. The enhancement of turbulence penetration by viscosity is similar to the flow near an oscillating flat plate (Stokes’ second problem), but is opposite to turbulence penetration into a stationary porous and permeable sediment bed. The information is applicable to water quality modeling mear the sediment/water interface of lakes, river impoundments and estuaries.     

Reef size and isolation determine the temporal stability of coral reef fish populations

Temporal variance in species abundance, a potential driver of extinction, is linked to mean abundance through Taylor's power law, the empirical observation of a linear log-log relationship with a slope between 1 and 2 for most species. Here we test the idea that the slope of Taylor's power law can vary both among species and spatially as a function of habitat area and isolation. We used the world's most extensive database of coral reef fish communities comprising a 15-year series of fish abundances on 43 reefs of Australia's Great Barrier Reef. Greater temporal variances were observed at small and isolated reefs, and lower variances at large and connected ones. The combination of reef area and isolation was associated with an even greater effect on temporal variances, indicating strong empirical support for the idea that populations on small and isolated reefs will succumb more frequently to local extinction via higher temporal variability, resulting in lower resilience at the community level. Based on these relationships, we constructed a regional predictive map of the dynamic fragility of coral reef fish assemblages on the Great Barrier Reef.     

Latitude or biogeographic breaks? Determinants of phenotypic (co)variation in fitness-related traits in Betaeus truncatus along the Chilean coast

Ectothermal organisms distributed along environmental gradients in a wide geographical distribution display extensive phenotypic variation. This is particularly pervasive along latitudinal clines, which are linked to gradual changes in environmental factors. Widespread species may also be distributed among biogeographic breaks, which in contrast to smooth clines, often show abrupt changes in phenotypic traits. In species with widespread latitudinal distribution that also encompass important biogeographical breaks, it is not clear which of those factors prevails on shaping the phenotypic variation or if some traits are particularly more sensitive to one or the other. To evaluate this, we measured 4 fitness-related traits in 6 populations of the intertidal snapping shrimp Betaeus truncatus, as its distribution along Chile expands over 40° in latitude and three major biogeographical provinces. Here, we statistically evaluated the role of both, latitude and biogeographic breaks, on mean population values of fitness-related traits but also on the variances and covariances (i.e., P-matrix) between them. Overall, our results (1) indicate that latitude is more important than breaks in shaping the phenotypic variation of most of these fitness-related traits, (2) show that the differences in the variance–covariance relationship among traits between the extremes of the gradient arises from gradual increases in variance and rather sharp changes in covariance at mid-latitudes and (3) show that at present, it is difficult to unambiguously determine whether natural selection or plasticity is responsible for the observed pattern in means, variances and covariances and only further work might disentangle these possibilities.     

Designing an efficient adaptive cluster sample

The choice of neighborhood definition and critical value in adaptive cluster sampling is critical for designing an efficient survey. In designing an efficient adaptive cluster sample one should aim for a small difference between the initial and final sample size, and a small difference between the within-network and population variances. However, the two aims can be at odds with each other because small differences between initial and final sample size usually means small within-network variance. One way to help in designing an efficient survey is to think in terms of small network sizes since the network size is a function of both critical value and neighborhood definition. One should aim for networks that are small enough to ensure the final sample size is not excessively large compared with the initial sample size but large enough to ensure the within-network variance is a reasonable fraction of the population variance. In this study surveys that had networks that were two to four units in size were the most efficient.     

海水混合和层化对叶绿素a垂直分布的影响

根据长江口海区、近岸浅水区、黄海冷水团海区和水深超过200m的陆架区等四个区域的定点调查获得的数据，分析了温度、盐度与叶绿素a垂直分布的相互关系。结果表明，长江口海区，陆源径流与海水混合不充分，表层营养盐含量较高，表层叶绿素a含量高于中下水层。近岸浅水区的苏北近岸海水垂向没有温跃层和盐跃层，叶绿素a的垂向分布也均匀，东海西部沿岸出现逆温跃层和逆盐跃层，海水垂直混合不充分，叶绿素a含量(6．72mg／m^3)在10m深水域最高。黄海冷水团海区海水的垂直混合不充分，叶绿素a的垂直变化显著，高值区出现在温跃层下方(4．37mg／m^3)。水深超过200m的陆架区，温度的阶梯状结构、营养盐跃层、光照等因素共同导致整个水层叶绿素a含量普遍较低。同时，结合历史资料分析认为，海水的混合、温度和盐度的层化将影响营养盐的浓度和分布，从而影响海水中叶绿素a的垂直分布。     

Thermal behaviour of a small spider (Araneae : Araneidae : Araneinae) on horizontal webs in semi-arid Western Australia

Summary The solar orientation response of a small spider (gen. nov.: Araneae: Araneidae: Araneinae) is examined at a semi-arid location in tropical Western Australia. The spiders rest beneath horizontal dome-shaped webs close to the ground surface on spinifex (Triodia sp.). The spiders are small (c. 5 mm long) and robust, being only twice as long as wide. Spider orientation, elevation, wind velocity, radiation and wet and dry bulb temperatures together with spider and egg-sac temperatures were measured throughout the day. The spiders orientated their long axis to the solar azimuth through much of the day and postured to the solar elevation; in so doing they tracked the apparent movement of the sun throughout the day with considerable precision, far exceeding that of other spiders studied. Stepwise regression indicated that the solar position was the most important environmental factor associated with the posturing and orientation. The variance in the orientation of the population was associated with a suite of environmental variables related to wind velocity and temperature. Spider (abdominal) temperature had a curvilinear relationship with ambient temperature, with a plateau at c. 35° C, whereas egg-sac temperature had a strictly linear relationship with ambient temperature. Stepwise regression indicates that both the spider and egg-sac temperatures are related mainly to net radiation and ambient temperature. The spiders can alter the incident radiation to which they are exposed by c. 2.2 times which represent the limits of the exposed silhouette area between full dorsal sunlight and the posterior (or anterior) aspect. Under hot conditions the spiders posture and orientate such that the long axis of the body is orientated abdomen towards the sun.     

Mixing in the three-dimensional wake of an experimental modelled vehicle

Experimental results on tracer gas diffusion within the near wake of a simplified model car (Ahmed model with a rear slant angle of 25°) are presented. Pollutant emission is simulated using heated air injected through a small pipe at one side of the model base. Fine cold wire thermometry is used to measure instantaneous temperature excess and variance of temperature gradient in the near wake. Measurements of the three mean velocity components were made using a laser Doppler anemometers system. Characteristics of the mean and fluctuating temperature fields, time-averaged flow streamlines and scalar dissipation measurements are presented and discussed. The local mixing time scale is determined from the measured mean dissipation rate of temperature variance. Its value shows that micromixing is not a limiting phenomenon for chemical reactions in the near wake.     

Observations of the cross-wind velocity variance in the stable boundary layer

Nine tower datasets over grassland, brush rangeland, snow covered plain, the ocean, three different pine forests, an aspen forest and an urban site, are used to document the scale-dependence of the cross-wind velocity variance in the stable boundary layer. The turbulence velocity variance scales with the surface momentum flux, as reported in previous studies. Such scaling removes the stability dependence of the variance at a given site, and also removes most of the differences between sites. The scaling is more effective with use of a record-dependent averaging time for defining the turbulent fluctuations. The variable averaging time is the timescale associated with the gap region in the heat flux multiresolution cospectra. On scales larger than turbulence and less than a few hours (mesoscale), variations in the cross-wind velocity variance at a given site are not related to local variables such as the friction velocity. Possible exceptions include suppression of turbulence and mesoscale motions in well-defined drainage flows and enhancement of turbulence and mesoscale motions in stronger winds downstream of a ridge. Larger mesoscale variance is associated with complex terrain and forested sites compared to the more homogeneous sites in flat terrain with short or no vegetation. These differences between sites are related to the absence of a gap region in the velocity spectra at the complex terrain and forested sites. The observed probability distribution functions of the total variance and the mesoscale variance are documented for different averaging times, stability classes and site characteristics.     

Effects of temperature and the presence of benthic predators on the vertical distribution of the ctenophore<Emphasis Type="Italic"> Pleurobrachia pileus</Emphasis>

This study details the investigations of the effect of temperature and the presence of predators on the vertical distribution of the ctenophore Pleurobrachia pileus. Previous studies had suggested that P. pileus moves to the sediment when exposed to colder temperatures as an energy saving strategy. Although we did observe an effect of temperature on the movement of the ciliate combs, the effect of temperature on the average depth of animals in the water column was not significant, but a significant interaction was found between temperature and exposure time. In the presence of benthic predators or of water previously inhabited by these predators P. pileus changed its vertical distribution significantly, and moved away from the sediment, thus reducing the risk of being predated. We conclude that P. pileus can sense the presence of its predators, and changes its behaviour accordingly.Communicated by O. Kinne, Oldendorf/Luhe     

Contribution à l'étude du phototropisme du plancton récolté dans la couche hyponeustonique. II. Anomalocera patersoni—Influence de l'intensité et de la longueur d'onde de la lumière sur les réactions d'orientation et sur l'activité locomotrice

The influence of light intensity and wavelength on orientation and velocity of male and female Anomalocera patersoni is summarized. G test and variance analysis show that light intensity and wavelength can affect the orientation and the velocity of male individuals, but do not significantly affect the polarization coefficient of female individuals. Nevertheless, photopositive reactions are dominant in both sexes, and the parameters studied may be modified by temperature changes.     

Induced triploidy in the soft-shelled clam Mya arenaria: energetic implications

Given that triploid adult bivalves reportedly grow larger and faster than their diploid siblings, such differences should be traceable to variation in energy allocation. In one proposed mechanism, retarded gametogenesis found in triploid adults may allow them more energy for somatic growth. Another hypothesis states that triploids are more heterozygous; increased heterozygosity has been positively correlated with enhanced growth. Juvenile soft-shelled clams, Mya arenaria, were treated with cytochalasin B to induce triploidy and examined with respect to components of a balanced energy equation (C=P+R+E). The variables measured were oxygen uptake (V _{o 2}), filtration rate (FR), dry tissue weight, shell length, shell height and shell inflation. Energy budgets were constructed and diploid and triploid groups compared. Few significant differences were found between diploid and triploid juvenile clams with respect to energy budget components. However, at seven loci assayed electrophoretically the triploid individuals were nearly twice as heterozygous as their diploid siblings. Moreover, tripoloid variances were less than diploid variances for every variable measured. Increased heterozygosity has been correlated with the decreased variance of morphological parameters. This study is believed to be the first to show decreased variance of physiological properties as well as morphological characters. Overall the data clearly indicate that energy allocation in juvenile M. arenaria is not related to ploidy.     

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.     

Including model uncertainty in estimating variances in multiple capture studies

We present bootstrap-based methods which incorporate model uncertainty in estimating variances in multiple capture studies. Each of our three methods has a specific set of properties, and we discuss when each method should be used. Our first method can be used in any multiple capture setting, but it gives an estimate of the variance conditional on the number of observed animals. Our other two methods yield estimates of the unconditional variance; they require good estimates of part or all of the specific probability model, respectively. Smoothed estimated cell probabilities are utilized by the latter method. We contrast the three methods on a real-life data set, and then conduct simulations for a simple setting. Finally, we detail the use of our methodology for specific settings and discuss adaptations for tag-return studies.     

Drivers of euphausiid species abundance and numerical abundance in the Atlantic Ocean

Mid-ocean ridges are common features of the world’s oceans but there is a lack of understanding as to how their presence affects overlying pelagic biota. The Mid-Atlantic Ridge (MAR) is a dominant feature of the Atlantic Ocean. Here, we examined data on euphausiid distribution and abundance arising from several international research programmes and from the continuous plankton recorder. We used a generalized additive model (GAM) framework to explore spatial patterns of variability in euphausiid distribution on, and at either side of, the MAR from 60°N to 55°S in conjunction with variability in a suite of biological, physical and environmental parameters. Euphausiid species abundance peaked in mid-latitudes and was significantly higher on the ridge than in adjacent waters, but the ridge did not influence numerical abundance significantly. Sea surface temperature (SST) was the most important single factor influencing both euphausiid numerical abundance and species abundance. Increases in sea surface height variance, a proxy for mixing, increased the numerical abundance of euphausiids. GAM predictions of variability in species abundance as a function of SST and depth of the mixed layer were consistent with present theories, which suggest that pelagic niche availability is related to the thermal structure of the near surface water: more deeply-mixed water contained higher euphausiid biodiversity. In addition to exposing present distributional patterns, the GAM framework enables responses to potential future and past environmental variability including temperature change to be explored.     

Ontogenetic diel vertical migration of the planktonic copepod Calanus sinicus in the Inland Sea of Japan

The ontogenetic diel vertical migration of the planktonic copepod Calanus sinicus was investigated in the Inland Sea of Japan in November 1988 and March 1989, when the water temperature was weakly stratified in a reversed manner. In both investigations a pronounced ontogenetic difference in vertical distribution was found. Spawning always occurred during nighttime, being confined to the upper 40 m water column in November but to the layer below 35 m in March. The distribution of pre-feeding nauplius stages, NI and NII, was more or less similar to that of the eggs. The first-feeding NIII performed a marked upward migration, and late nauplius stages (NIV to NVI) and early copepodite stages (CI and CII) continuously aggregated in the upper water column where phytoplankton was abundant. CIII to CVI (adult female and male) tended to disperse in the whole water column. In November, however, they avoided the upper 10 m strate during daytime and some individuals migrated upward to the surface during nighttime. In March, CV and CVI aggregated in the layer between 5 and 15 m deep in the daytime and migrated both upward and downward at dusk, resulting in homogeneous distributions during the nighttime.     

Characteristics of flow structure of free-surface flow in a partly obstructed open channel with vegetation patch

Free-surface flows over patchy vegetation are common in aquatic environments. In this study, the hydrodynamics of free-surface flow in a rectangular channel with a bed of rigid vegetation-like cylinders occupying half of the channel bed was investigated and interpreted by means of laboratory experiments and numerical simulations. The channel configurations have low width-to-depth aspect ratio (1.235 and 2.153). Experimental results show that the adjustment length for the flow to be fully developed through the vegetation patch in the present study is shorter than observed for large-aspect-ratio channels in other studies. Outside the lateral edge of the vegetation patch, negative velocity gradient (\(\partial \overline{u}/\partial z < 0\)) and a local velocity maximum are observed in the vertical profile of the longitudinal velocity in the near-bed region, corresponding to the negative Reynolds stress (\(- \overline{{u^{\prime}w^{\prime}}} < 0\)) at the same location. Assuming coherent vortices to be the dominant factor influencing the mean flow field, an improved Spalart–Allmaras turbulence model is developed. The model improvement is based on an enhanced turbulence length scale accounting for coherent vortices due to the effect of the porous vegetation canopy and channel bed. This particular flow characteristic is more profound in the case of high vegetation density due to the stronger momentum exchange of horizontal coherent vortices. Numerical simulations confirmed the local maximum velocity and negative gradient in the velocity profile due to the presence of vegetation and bed friction. This in turn supports the physical interpretation of the flow processes in the partly obstructed channel with vegetation patch. In addition, the vertical profile of the longitudinal velocity can also be explained by the vertical behavior of the horizontal coherent vortices based on a theoretical argument.     

Fish assemblage on natural versus vertical artificial reefs: the rehabilitation perspective

Artificial reefs have been suggested as a potential tool for the restoration of marine habitats. In the present study, the fish assemblage established around the oil jetties of Eilat (northern Red Sea, Israel) was compared to those found in three adjacent natural reef habitats: two in a nature reserve (one shallow and one deep) and a third deep site located near the city. Both species richness and fish abundance were found to be significantly higher around the vertical structures of the jetty's pillars than at all three natural sites, with the lowest values at the site closest to the city. The higher species richness at the jetties may be explained by (1) the vertical relief and high complexity of the jetty which offers a variety of niches for both shallow and deep coral reef species, and (2) by the reduction in available niches at the natural sites as a result of coral destruction due to anthropogenic activity. The pronounced difference in fish abundance is attributed mainly to the high seasonal recruitment at the jetty which was much lower at the natural sites. We therefore suggest that vertical structures are more attractive to fish settlement and recruitment than moderately sloped bottoms such as those found at the fringing reefs of Eilat. High similarity (51 to 56%) was found between fish assemblages at the natural sites while relatively low similarity (27 to 37%) was found between the jetty and the natural reefs. The jetty's complex vertical artificial structures can serve as a model for future construction of artificial reefs designed to restore the fish community in areas where the natural reefs have been damaged. It should be taken into account, however, that these do not necessarily mimic the natural environment but may rather establish a community of their own, which is influenced by the spatial orientation and complexity of the structure. Received: 30 December 1998 / Accepted: 9 December 1999