Is the Mantel correlogram powerful enough to be useful in ecological analysis? A simulation study

The Mantel correlogram is an elegant way to compute a correlogram for multivariate data. However, recent papers raised concerns about the power of the Mantel test itself. Hence the question: Is the Mantel correlogram powerful enough to be useful? To explore this issue, we compared the performances of the Mantel correlogram to those of other methods, using numerical simulations based on random, normally distributed data. For a single response variable, we compared it to the Moran and Geary correlograms. Type I error rates of the three methods were correct. Power of the Mantel correlogram was nearly as high as that of the univariate methods. For the multivariate case, the test of the multivariate variogram developed in the context of multiscale ordination is in fact a Mantel test, so that the power of the two methods is the same by definition. We devised an alternative permutation test based on the variance, which yielded similar results. Overall, the power of the Mantel test was high, the method successfully detecting spatial correlation at rates similar to the permutation test of the variance statistic in multivariate variograms. We conclude that the Mantel correlogram deserves its place in the ecologist's toolbox.     

Evolution, radiation and chemotaxonomy of Lamellodysidea, a demosponge genus with anti-plasmodial metabolites

Sponges of the family Dysideidae (Dictyoceratida) are renowned for their diversity of secondary metabolites, and its genus Lamellodysidea, particularly Lamellodysidea herbacea, is the most studied taxon biochemically. Despite its importance, the taxonomic status of L. herbacea—whether it is a distinct species or a species complex—has never been assessed. Recent biochemical profiling revealed anti-plasmodial activity of brominated compounds in Lamellodysidea of the Pacific. Here, we present a comparative chemotaxonomic and molecular analysis of selected Dysideidae from the Pacific and the Indian Ocean (New Caledonia, Great Barrier Reef, Fiji, Mayotte, Guam, Palau). We investigated the phylogenetic relationships between the populations and assessed their bioactive (PBDE) compounds in order to unravel the taxonomic status of this commercially important group of sponges and assessed patterns of dispersal and biochemical variation. The molecular phylogeny was based on the internal transcribed ribosomal spacer and compared against a PBDE phylogeny for several specimens. Molecular data revealed a diversity of Indo-Pacific L. herbacea populations, also reflected by different PBDE compound profiles. Molecular and biochemical data also revealed a Lamellodysidea species new to science. Several specimens misidentified as Lamellodysidea were detected based on their position on different clades in the molecular phylogeny and their production of different halogenated compounds (brominated vs. chlorinated). The direct comparison of molecular and biochemical data also provided evidence for the occurrence of a host switch event and support for the theory that abiotic factors, such as sedimentation, affect the chemical constituents produced in L. herbacea.     

Food sources used by sediment meiofauna in an intertidal Zostera noltii seagrass bed: a seasonal stable isotope study

In an intertidal Zostera noltii Hornem seagrass bed, food sources used by sediment meiofauna were determined seasonally by comparing stable isotope signatures (δ¹³C, δ¹⁵N) of sources with those of nematodes and copepods. Proportions of different carbon sources used by consumers were estimated using the SIAR mixing model on δ¹³C values. Contrary to δ¹⁵N values, food source mean δ¹³C values encompassed a large range, from −22.1 ‰ (suspended particulate organic matter) to −10.0 ‰ (Z. noltii roots). δ¹³C values of copepods (from −22.3 to −12.3 ‰) showed that they use many food sources (benthic and phytoplanktonic microalgae, Z. noltii matter). Nematode δ¹³C values ranged from −14.6 to −11.4 ‰, indicating a strong role of microphytobenthos and/or Z. noltii matter as carbon sources. The difference of food source uses between copepods and nematodes is discussed in light of source accessibility and availability.     

Multiscale codependence analysis: an integrated approach to analyze relationships across scales

The spatial and temporal organization of ecological processes and features and the scales at which they occur are central topics to landscape ecology and metapopulation dynamics, and increasingly regarded as a cornerstone paradigm for understanding ecological processes. Hence, there is need for computational approaches which allow the identification of the proper spatial or temporal scales of ecological processes and the explicit integration of that information in models. For that purpose, we propose a new method (multiscale codependence analysis, MCA) to test the statistical significance of the correlations between two variables at particular spatial or temporal scales. Validation of the method (using Monte Carlo simulations) included the study of type I error rate, under five statistical significance thresholds, and of type II error rate and statistical power. The method was found to be valid, in terms of type I error rate, and to have sufficient statistical power to be useful in practice. MCA has assumptions that are met in a wide range of circumstances. When applied to model the river habitat of juvenile Atlantic salmon, MCA revealed that variables describing substrate composition of the river bed were the most influential predictors of parr abundance at 0.4-4.1 km scales whereas mean channel depth was more influential at 200-300 m scales. When properly assessed, the spatial structuring observed in nature may be used purposefully to refine our understanding of natural processes and enhance model representativeness.     

Load Acceleration and Footstep Strategies in Asymmetrical Lifting and Lowering

Accelerated execution effects for lifting and lowering a 12-kg box using two footstep strategies associated with experienced workers were studied. Eight healthy male participants performed a normal and an accelerated execution of a lifting task and a lowering task, using a minimal feet displacement strategy (oblique-step) and a strategy with a step (crossed-step). It was hypothesized that the accelerated executions, as compared to the normal executions, would have a different effect on L5/S1 resultant moment, body posture, and other kinematic variables. A tridimensional dynamic rigid body model was used to compute L5/S1 resultant moments. Results showed that the accelerated condition did not reduce body asymmetry of posture, but it reduced the length of the path of the global center of gravity and the duration of the supporting phase of the box, and it did not significantly affect L5/S1 maximal resultant moments for lifting but increased them for lowering. These results indicate that the net work production for accelerated strategies might be smaller, which may represent an economy of energy. Furthermore, the results showed that the use of an accelerated strategy for lowering should be avoided.     

Aggressions and size-related fecundity of queenless workers in the ant <Emphasis Type="Italic">Cataglyphis cursor</Emphasis>

In social hymenoptera, the reproductive division of labor is often linked to differences in individual body size with the reproductive caste (the queen) being larger than the workers. Likewise, the reproductive potential may vary with size within the worker caste and could affect the evolution of worker size in social insects. Here, we tested the relationship between worker size and reproductive potential in the facultative parthenogenetic ant Cataglyphis cursor. Colonies are headed by a multiply mated queen, but workers can produce gynes (virgin queens) and workers by thelytokous parthenogenesis after the queen's death. We observed the behaviour of workers (n = 357) until the production of gynes (212 h over 3 months) in an orphaned colony (mated queen not present). The size of workers was measured, and their paternal lineage determined using six microsatellite markers, to control for an effect of patriline. Larger workers were more likely to reproduce and lay more eggs indicating that individual level selection could take place. However, paternal lineage had no effect on the reproductive potential and worker size. From the behavioural and genetic data, we also show for the first time in this species, evidence of aggressive interactions among workers and a potential for nepotism to occur in orphaned colonies, as the five gynes produced belonged to a single paternal lineage.     

Water quality assessment by means of HFNI valvometry and high-frequency data modeling

The high-frequency measurements of valve activity in bivalves (e.g., valvometry) over a long period of time and in various environmental conditions allow a very accurate study of their behaviors as well as a global analysis of possible perturbations due to the environment. Valvometry uses the bivalve's ability to close its shell when exposed to a contaminant or other abnormal environmental conditions as an alarm to indicate possible perturbations in the environment. The modeling of such high-frequency serial valvometry data is statistically challenging, and here, a nonparametric approach based on kernel estimation is proposed. This method has the advantage of summarizing complex data into a simple density profile obtained from each animal at every 24-h period to ultimately make inference about time effect and external conditions on this profile. The statistical properties of the estimator are presented. Through an application to a sample of 16 oysters living in the Bay of Arcachon (France), we demonstrate that this method can be used to first estimate the normal biological rhythms of permanently immersed oysters and second to detect perturbations of these rhythms due to changes in their environment. We anticipate that this approach could have an important contribution to the survey of aquatic systems.     

Characterization of municipal solid waste combustion in a grate furnace

The objective of this paper is to evaluate the combustion process of municipal solid waste combustion in a grate furnace both experimentally and numerically by using data of a reference experiment with over-stoichiometric primary air supply. Measurements were carried out inside the combustion chamber of a pilot plant by monitoring temperatures and sampling gaseous combustion products along the bed surface. The data were assessed using elemental and energy balances. Experimental data of the axial temperature profiles of the flue gas, the fuel bed and the grate bars, as well as local gas flows and the flue gas composition measured above the fuel bed along the grate were used to describe the conversion process, including drying and carbon burnout. These data served as input to model the thermo- and fluid dynamic processes of the gas phase above the bed inside the combustion chamber. For this purpose the commercial code FLUENT was employed to carry out the simulations. Thus, the turbulent temperature, flow and species distributions in the combustion chamber of the pilot waste incinerator TAMARA were predicted. The results of the FLUENT modeling showed that under the prevailing conditions the flue gas burnout is almost completed before entering the first flue due to high temperatures, effective mixing and sufficient residence times of the flue gas inside the combustion chamber. This agrees well with the experimental results inside the first flue. On the basis of the above mentioned results, design and parametric studies can be carried out in a more efficient way by saving cost and time.     

Forward selection of explanatory variables

This paper proposes a new way of using forward selection of explanatory variables in regression or canonical redundancy analysis. The classical forward selection method presents two problems: a highly inflated Type I error and an overestimation of the amount of explained variance. Correcting these problems will greatly improve the performance of this very useful method in ecological modeling. To prevent the first problem, we propose a two-step procedure. First, a global test using all explanatory variables is carried out. If, and only if, the global test is significant, one can proceed with forward selection. To prevent overestimation of the explained variance, the forward selection has to be carried out with two stopping criteria: (1) the usual alpha significance level and (2) the adjusted coefficient of multiple determination (Ra(2)) calculated using all explanatory variables. When forward selection identifies a variable that brings one or the other criterion over the fixed threshold, that variable is rejected, and the procedure is stopped. This improved method is validated by simulations involving univariate and multivariate response data. An ecological example is presented using data from the Bryce Canyon National Park, Utah, U.S.A.