Distribution patterns of chlorophyll a in spring and autumn in association with hydrological features in the southern Yellow Sea and northern East China Sea

Two field studies were conducted to measure pigments in the Southern Yellow Sea (SYS) and the northern East China Sea (NECS) in April (spring) and September (autumn) to evaluate the distribution pattern of phytoplankton stock (Chl a concentration) and the impact of hydrological features such as water mass, mixing and tidal front on these patterns. The results indicated that the Chl a concentration was 2.43 +/- 2.64 (Mean +/- SD) mg m(-3) in April (range, 0.35 to 17.02 mg m(-3)) and 1.75 +/- 3.10 mg m(-3) in September (from 0.07 to 36.54 mg m(-3)) in 2003. Additionally, four areas with higher Chl a concentrations were observed in the surface water in April, while two were observed in September, and these areas were located within or near the point at which different water masses converged (temperature front area). The distribution pattern of Chl a was generally consistent between onshore and offshore stations at different depths in April and September. Specifically, higher Chl a concentrations were observed along the coastal line in September, which consisted of a mixing area and a tidal front area, although the distributional pattern of Chl a concentrations varied along transects in April. The maximum Chl a concentration at each station was observed in the surface and subsurface layer (0-10 m) for onshore stations and the thermocline layer (10-30 m) for offshore stations in September, while the greatest concentrations were generally observed in surface and subsurface water (0-10 m) in April. The formation of the Chl a distributional pattern in the SYS and NECS and its relationship with possible influencing factors is also discussed. Although physical forces had a close relationship with Chl a distribution, more data are required to clearly and comprehensively elucidate the spatial pattern dynamics of Chl a in the SYS and NECS.

喜马拉雅西构造结区岩石圈演化三维有限元动力学数值模拟

The past three decades have seen numerous attempts to numerically model stress and strain patterns in the lithosphere of the Earth on both global and regional scales. This efforts have been indispensable in identifying the features we need to include in our endeavour to develop better models of our planet’s lithosphere and they have also raised our awareness for the many unresolved issue in the deep geodynamical issues that need to be addressed in the future. Nonetheless, in most models, the lithosphere is treated as a single layer with depth-averaged properties, and as the same distribution in the stress and strain fields, and as deforming under plane strain. All these above make a great hander for its reality and degree of recognition. As the beginning in this paper, some principal numerical models and results on the evolution of Tibetan plateau are reviewed and analyzed. Then, the geological and geophysical expedition on the Western Himalayan Syntaxis is briefly reviewed. Furthermore, we analysis the feature in deep geophysical field studies in this area and adjacent regions. Because, for most continents, stress models driven by plate boundary forces have successfully reproduced the main characteristics of the stress and strain field, we present a set of three-dimensional models of lithosphere system for a simplified geometry of the Western Himalayan Syntaxis area and its adjacent regions, where we try to match the first-order characteristics of the stress and strain fields of lithosphere since 10 Ma, and deformation and geodynamical evolution process in former 2Ma. Of course, the kinematic boundary conditions of the stress models driven by plate boundary forces were applied. The rheology plays a significant role in the lithospheric tectonics, which lead to different rheological parameters were used in different works although the have the same constitutive equations in models. So, in this paper we do not aim to produce all characteristics of the Western Himalayan Syntaxis areas’ stress and strain fields by the choices of various parameters, but rather the dynamic response between various rheological parameters and stress and strain fields. We have chosen to concentrate on the importance of rheology and lateral strength variations for lithospheric stress and strain patterns and use our findings to build a model of the Western Himalayan Syntaxis areas. In doing so, we want to go beyond purely elastic models or purely viscoelastic models. Compared the results of the crust viscosity in the Western Himalayan Syntaxis areas, we believed that, when various viscoelastic models are adopted, the selection of the coefficient of viscosity in the Western Syntaxis area has important influence on the its uplifts and evolutions. A wider uplift ranges and gently elevation was observed at the same time when a lower viscosity was used in our models, and vice versa. Data of stress magnitudes are not available, but it is clear that the stress levels must be at or below the failure threshold of rock under compression. Under these criteria, the calculation results show that the viscosity in the Western Syntaxis area should be smaller than 1023Pa.s When elastic model is adopted in relatively rigid Tarim basin, obvious changes are induced to the stress and strain fields of the whole Western Syntaxis area. We found that rigid block of lithosphere reduced stress levels within its interior and that, at the edges of such regions, stress orientation can change. Furthermore there is no evidence that such rigid regions act as stress barriers in that they shield areas in opposite sides of the structure from the influence of one another. In our models, the upper crustal material of the Western Syntaxis area does not turns to move westward. Whereas, because of the stress and strain fields have been decoupling at the interior of the lithosphere, we can get the results that the deep material must not move westward.

Exploring alternatives for milk quality improvement and more efficient dairy production in a smallholder farming context – case study: Mantaro Valley (Peru)

Consumer demand is revolutionizing the way products are being produced, distributed and marketed. In relation to the dairy sector in developing countries, aspects of milk quality are receiving more attention from both society and the government. However, milk quality management needs to be better addressed in dairy production systems to guarantee the access of stakeholders, mainly small-holders, into dairy markets. The present study is focused on an analysis of the interaction of the upstream part of the dairy supply chain (farmers and dairies) in the Mantaro Valley (Peruvian central Andes), in order to understand possible constraints both stakeholders face implementing milk quality controls and practices; and evaluate “ex-ante” how different strategies suggested to improve milk quality could affect farmers and processors’ profits. The analysis is based on three complementary field studies conducted between 2012 and 2013. Our work has shown that the presence of a dual supply chain combining both formal and informal markets has a direct impact on dairy production at the technical and organizational levels, affecting small formal dairy processors’ possibilities to implement contracts, including agreements on milk quality standards. The analysis of milk quality management from farms to dairy plants highlighted the poor hygiene in the study area, even when average values of milk composition were usually high. Some husbandry practices evaluated at farm level demonstrated cost effectiveness and a big impact on hygienic quality; however, regular application of these practices was limited, since small-scale farmers do not receive a bonus for producing hygienic milk. On the basis of these two results, we co-designed with formal small-scale dairy processors a simulation tool to show prospective scenarios, in which they could select their best product portfolio but also design milk payment systems to reward farmers’ with high milk quality performances. This type of approach allowed dairy processors to realize the importance of including milk quality management in their collection and manufacturing processes, especially in a context of high competition for milk supply. We concluded that the improvement of milk quality in a smallholder farming context requires a more coordinated effort among stakeholders. Successful implementation of strategies will depend on the willingness of small-scale dairy processors to reward farmers producing high milk quality; but also on the support from the State to provide incentives to the stakeholders in the formal sector.

Biodeposition and Biogeochemical Processes in Shallow, Mesohaline Sediments of Chesapeake Bay

Stocks of the eastern oyster, Crassostrea virginica, have been declining in Chesapeake Bay since the late 19th century, and current strategies involve restoring culture of Crassostrea virginica on-bottom and in devices suspended within the water column. Sub-tidal suspension culture of Crassostrea virginica in Chesapeake Bay occurs mostly in sheltered inlets and tidal creeks and, thereby, has the potential to influence shallow water biogeochemical processes. To assess the influence of Crassostrea virginica biodeposits and benthic microalgae on sediment nitrogen and phosphorus exchange, field studies with Crassostrea virginica held in aquaculture floats and laboratory experiments were conducted. Enhanced organic nitrogen deposition from Crassostrea virginica biodeposits led to gradual increases in surface sediment nitrogen and pore water ammonium concentrations; however, modifications to pore water concentrations were not always expressed at the sediment-water interface. Benthic microalgae often modulated the influence of biodeposits on sediment nitrogen exchange but, as observed in laboratory experiments, the supply of nitrogen from Crassostrea virginica biodeposits may exceed their biological demand. Organic carbon from biodeposits had varying influences on aerobic respiration but consistently stimulated anaerobic metabolism. Shifts in net phosphorus exchange were driven by this anaerobic remineralization and concentrations of iron and manganese oxy(hydr)oxides, with transitions in fluxes coinciding with changes in benthic photosynthesis and oxidation of surface sediments. Manganese and iron oxy(hydr)oxides from biodeposits supported incorporation of added phosphorus and prevented exchange at the sediment-water interface in the absence of iron-sulfide mineral formation. Differences in the response of shallow water sediments to Crassostrea virginica biodeposits were due to the quality and quantity of biodeposits supplied, as well as the spatial and temporal variability within these sediments. Initial conditions and corresponding reference sediments illustrated the potential for sediment biogeochemistry and nutrient exchange from tidal creek sediments to vary spatially and temporally on relatively small scales. Factors influencing variability within tidal creek sediments were related to shifts in riverine freshwater inputs, macroalgal blooms, nutrient concentrations in overlying waters, and bioirrigation from the clam, Macoma balthica.

Global Rates of Free Hydrogen (H2) Production by Serpentinization and other Abiogenic Processes within Young Ocean Crust

The main conclusion of this dissertation is that global H2 production within young ocean crust (<10 Mya) is higher than currently recognized, in part because current estimates of H2 production accompanying the serpentinization of peridotite may be too low (Chapter 2) and in part because a number of abiogenic H2-producing processes have heretofore gone unquantified (Chapter 3). The importance of free H2 to a range of geochemical processes makes the quantitative understanding of H2 production advanced in this dissertation pertinent to an array of open research questions across the geosciences (e.g. the origin and evolution of life and the oxidation of the Earth’s atmosphere and oceans).

The first component of this dissertation (Chapter 2) examines H2 produced within young ocean crust [e.g. near the mid-ocean ridge (MOR)] by serpentinization. In the presence of water, olivine-rich rocks (peridotites) undergo serpentinization (hydration) at temperatures of up to ~500°C but only produce H2 at temperatures up to ~350°C. A simple analytical model is presented that mechanistically ties the process to seafloor spreading and explicitly accounts for the importance of temperature in H2 formation. The model suggests that H2 production increases with the rate of seafloor spreading and the net thickness of serpentinized peridotite (S-P) in a column of lithosphere. The model is applied globally to the MOR using conservative estimates for the net thickness of lithospheric S-P, our least certain model input. Despite the large uncertainties surrounding the amount of serpentinized peridotite within oceanic crust, conservative model parameters suggest a magnitude of H2 production (~1012 moles H2/y) that is larger than the most widely cited previous estimates (~1011 although previous estimates range from 1010-1012 moles H2/y). Certain model relationships are also consistent with what has been established through field studies, for example that the highest H2 fluxes (moles H2/km2 seafloor) are produced near slower-spreading ridges (<20 mm/y). Other modeled relationships are new and represent testable predictions. Principal among these is that about half of the H2 produced globally is produced off-axis beneath faster-spreading seafloor (>20 mm/y), a region where only one measurement of H2 has been made thus far and is ripe for future investigation.

In the second part of this dissertation (Chapter 3), I construct the first budget for free H2 in young ocean crust that quantifies and compares all currently recognized H2 sources and H2 sinks. First global estimates of budget components are proposed in instances where previous estimate(s) could not be located provided that the literature on that specific budget component was not too sparse to do so. Results suggest that the nine known H2 sources, listed in order of quantitative importance, are: Crystallization (6x1012 moles H2/y or 61% of total H2 production), serpentinization (2x1012 moles H2/y or 21%), magmatic degassing (7x1011 moles H2/y or 7%), lava-seawater interaction (5x1011 moles H2/y or 5%), low-temperature alteration of basalt (5x1011 moles H2/y or 5%), high-temperature alteration of basalt (3x1010 moles H2/y or <1%), catalysis (3x108 moles H2/y or <<1%), radiolysis (2x108 moles H2/y or <<1%), and pyrite formation (3x106 moles H2/y or <<1%). Next we consider two well-known H2 sinks, H2 lost to the ocean and H2 occluded within rock minerals, and our analysis suggests that both are of similar size (both are 6x1011 moles H2/y). Budgeting results suggest a large difference between H2 sources (total production = 1x1013 moles H2/y) and H2 sinks (total losses = 1x1011 moles H2/y). Assuming this large difference represents H2 consumed by microbes (total consumption = 9x1011 moles H2/y), we explore rates of primary production by the chemosynthetic, sub-seafloor biosphere. Although the numbers presented require further examination and future modifications, the analysis suggests that the sub-seafloor H2 budget is similar to the sub-seafloor CH4 budget in the sense that globally significant quantities of both of these reduced gases are produced beneath the seafloor but never escape the seafloor due to microbial consumption.

The third and final component of this dissertation (Chapter 4) explores the self-organization of barchan sand dune fields. In nature, barchan dunes typically exist as members of larger dune fields that display striking, enigmatic structures that cannot be readily explained by examining the dynamics at the scale of single dunes, or by appealing to patterns in external forcing. To explore the possibility that observed structures emerge spontaneously as a collective result of many dunes interacting with each other, we built a numerical model that treats barchans as discrete entities that interact with one another according to simplified rules derived from theoretical and numerical work, and from field observations: Dunes exchange sand through the fluxes that leak from the downwind side of each dune and are captured on their upstream sides; when dunes become sufficiently large, small dunes are born on their downwind sides (“calving”); and when dunes collide directly enough, they merge. Results show that these relatively simple interactions provide potential explanations for a range of field-scale phenomena including isolated patches of dunes and heterogeneous arrangements of similarly sized dunes in denser fields. The results also suggest that (1) dune field characteristics depend on the sand flux fed into the upwind boundary, although (2) moving downwind, the system approaches a common attracting state in which the memory of the upwind conditions vanishes. This work supports the hypothesis that calving exerts a first order control on field-scale phenomena; it prevents individual dunes from growing without bound, as single-dune analyses suggest, and allows the formation of roughly realistic, persistent dune field patterns.

Oceanographic responses to climate in the northwest Atlantic

Situated in an oceanographic transition zone, the Gulf of Maine/Western Scotian Shelf (GOM/WSS) region of the Northwest Atlantic is especially susceptible to changes in the climate system. Recent studies have shown that a coupled slope water system (CSWS) operates in the Northwest Atlantic and responds in a similar manner to climatic forcing over a broad range of time scales. These studies further suggest that it may be possible to associate different modes of the CSWS with the different phases of the North Atlantic Oscillation (NAO). Results from recent GLOBEC field studies in the Northwest Atlantic provide strong evidence linking physical responses of the CSWS to basin-scale forcing associated with the NAO. By placing these results in the context of time-series data collected from the GOM/WSS over the past half century, we show that we show that: (i) the region’s shelf ecosystems respond both physically and biologically to modal shifts in the CSWS; (ii) the CSWS mediates the effects on these ecosystems of basin-scale climatic forcing associated with the NAO and (iii) certain planktonic species can be good indicators of the CSWS’s modal state on inter-annual to interdecadal time scales.

Some Ecological Consequences Of The Physiological And Biochemical Effects Of Petroleum Compounds On Marine Mollusks

The problems of relating the results of experiments in the laboratory to events in nature are twofold: to equate the response to a single variable (hydrocarbons) with the natural variability in the biological material in a multivariate environment, and to consider whether the response established experimentally has any relevance to the animal's chances of survival and reproduction (i.e. its fitness) in the natural population. Recent studies of the effects of petroleum hydrocarbons on marine invertebrates are reviewed, with an emphasis on the physiological and cytochemical responses by bivalve molluscs. The dose-response relations that emerge suggest the intensity of the 'signal' that must be detected in nature if the chronic, sublethal effects of petroleum pollution are to be measured. The natural variability in these physiological and cytochemical processes are then reviewed and the main causes of variability in natural populations, both endogenous and exogenous, discussed. These results indicate the extent of the `noise' above which the signal from possible pollution effects must be detected. The results from recent field studies on the common mussel, Mytilus edulis, are discussed. The results are as complex as expected, but it proves possible to reduce the variance in the measured responses so that pollution effects, including those due to hydrocarbons, can be detected. The ecological consequences of the observed effects of petroleum hydrocarbons are then discussed in terms of reproductive effort and reproductive value. Considerable variation between populations exists here also and this can be used to help in the interpretation of the extent of the impact of the environment on the ecology of the population. The result is to place the findings of the laboratory experiments in an ecological context of natural variability and of the physiological costs of adaptation.

Strong preference for decapod prey by the western rock lobster Panulirus cygnus

Western rock lobsters, Panulirus cygnus are an abundant benthic consumer distributed along the temperate west coast of Australia and constitute the largest single species fishery in Australia. As a dominant consumer, it is important to understand their predator-prey interactions as they can potentially exert strong trophic effects, and may influence ecosystem function as seen in other spiny lobster species. While previous field studies have focused on the diet composition of P. cygnus, this study investigated their preference for various benthic invertebrate prey to better understand the likely predator-prey interactions of P. cygnus. Prey preferences of small sub-legal juvenile lobsters, as well as medium and large legal-sized mature lobsters were investigated using laboratory feeding trials to identify size-associated differences in lobster prey preference. Handling time and diet quality were investigated to estimate energetic cost and gain from consuming different prey which may explain prey choice by lobsters. It was found that large lobsters preferred crabs and mussels while medium and small lobsters preferred crabs over mussels, gastropods, and sea urchins. This suggests that strong predator-prey interactions between P. cygnus and crabs may occur in the wild.

Indications of a negative impact of herring on recruitment of Norway pout

The Norway pout (Trisopterus esmarkii) stock in the North Sea has experienced poor recruitment recently. Herring (Clupea harengus) has been suggested to be a major predator on fish larvae in the North Sea. We investigated possible interactions between herring and Norway pout using a simple statistical analysis and a modified stock - recruit relationship. There was a significant negative relationship (linear regression, r = 20.44, p < 0.05) between total herring biomass and recruitment of Norway pout. The spawning stock of Norway pout is typically dominated by 2-year-olds, and there was a strong negative relationship (linear regression, r = 20.79, p < 0.01) between herring biomass and Norway pout spawning-stock biomass (SSB) 2 years later. A Beverton-Holt model fitted to stock recruit data of Norway pout produced a rather poor correlation (r(2) = 0.04). However, when only the Norway pout SSB not overlapping with herring is considered, the fit between the model and the stock - recruit data improves (r(2) = 0.31). The analyses indicate a negative impact by herring on recruitment of Norway pout, the most plausible cause for this being herring predation on Norway pout larvae, but field studies are needed to verify such predation.

Satellite-detected fluorescence reveals global physiology of ocean phytoplankton

Phytoplankton photosynthesis links global ocean biology and climate-driven fluctuations in the physical environment. These interactions are largely expressed through changes in phytoplankton physiology, but physiological status has proven extremely challenging to characterize globally. Phytoplankton fluorescence does provide a rich source of physiological information long exploited in laboratory and field studies, and is now observed from space. Here we evaluate the physiological underpinnings of global variations in satellite-based phytoplankton chlorophyll fluorescence. The three dominant factors influencing fluorescence distributions are chlorophyll concentration, pigment packaging effects on light absorption, and light-dependent energy-quenching processes. After accounting for these three factors, resultant global distributions of quenching-corrected fluorescence quantum yields reveal a striking consistency with anticipated patterns of iron availability. High fluorescence quantum yields are typically found in low iron waters, while low quantum yields dominate regions where other environmental factors are most limiting to phytoplankton growth. Specific properties of photosynthetic membranes are discussed that provide a mechanistic view linking iron stress to satellite-detected fluorescence. Our results present satellite-based fluorescence as a valuable tool for evaluating nutrient stress predictions in ocean ecosystem models and give the first synoptic observational evidence that iron plays an important role in seasonal phytoplankton dynamics of the Indian Ocean. Satellite fluorescence may also provide a path for monitoring climate-phytoplankton physiology interactions and improving descriptions of phytoplankton light use efficiencies in ocean productivity models.

Bridging the gap between marine biogeochemical and fisheries sciences; configuring the zooplankton link

Exploring climate and anthropogenic impacts on marine ecosystems requires an understanding of how trophic components interact. However, integrative end-to-end ecosystem studies (experimental and/or modelling) are rare. Experimental investigations often concentrate on a particular group or individual species within a trophic level, while tropho-dynamic field studies typically employ either a bottom-up approach concentrating on the phytoplankton community or a top-down approach concentrating on the fish community. Likewise the emphasis within modelling studies is usually placed upon phytoplankton-dominated biogeochemistry or on aspects of fisheries regulation. In consequence the roles of zooplankton communities (protists and metazoans) linking phytoplankton and fish communities are typically under-represented if not (especially in fisheries models) ignored. Where represented in ecosystem models, zooplankton are usually incorporated in an extremely simplistic fashion, using empirical descriptions merging various interacting physiological functions governing zooplankton growth and development, and thence ignoring physiological feedback mechanisms. Here we demonstrate, within a modelled plankton food-web system, how trophic dynamics are sensitive to small changes in parameter values describing zooplankton vital rates and thus the importance of using appropriate zooplankton descriptors. Through a comprehensive review, we reveal the mismatch between empirical understanding and modelling activities identifying important issues that warrant further experimental and modelling investigation. These include: food selectivity, kinetics of prey consumption and interactions with assimilation and growth, form of voided material, mortality rates at different age-stages relative to prior nutrient history. In particular there is a need for dynamic data series in which predator and prey of known nutrient history are studied interacting under varied pH and temperature regimes.

Long-term annual primary production in the Ulleung Basin as a biological hot spot in the East/Japan Sea

Although the Ulleung Basin is an important biological hot spot in East/Japan Sea (hereafter the East Sea), very limited knowledge for seasonal and annual variations in the primary productivity exists. In this study, a recent decadal trend of primary production in the Ulleung Basin was analyzed based on MODIS-derived monthly primary production for a better annual production budget. Based on the MODIS-derived primary production, the mean daily primary productivity was 766.8 mg C m-2 d-1 (SD=+/- 196.7 mg C m-2 d-1) and the annual primary productivity was 280.2 g C m-2 yr-1 (SD=+/- 14.9 g C m-2 yr-1) in the Ulleung Basin during the study period. The monthly contributions of primary production were not largely variable among different months, and a relatively small interannual production variability was also observed in the Ulleung Basin, which indicates that the Ulleung Basin is a sustaining biologically productive region called as hot spot in the East Sea. However, a significant recent decline in the annual primary production was observed in the Ulleung Basin after 2006. Although no strong possibilities were found in this study, the current warming sea surface temperature and a negative phase PDO index were suggested for the recent declining primary production. For a better understanding of subsequent effects on marine ecosystems, more intensive interdisciplinary field studies will be required in the Ulleung Basin.

Coastal upwelling and downwelling forcing of circulation in a semi-enclosed bay: Ria de Vigo

Semi-enclosed bays in upwelling regions are exposed to forcing related to winds, currents and buoyancy over the shelf. The influence of this external forcing is moderated by factors such as connectivity to the open ocean, shelter by surrounding topography, dimensions of the bay, and freshwater outflows. Such bays, preferred locations for ports, mariculture, marine industry, recreational activities and coastal settlement, present a range of characteristics, understanding of which is necessary to their rational management. Observations in such a semi-enclosed bay, the Ria de Vigo in Spain, are used to characterize the influence of upwelling and downwelling pulses on its circulation. In this location, near the northern limit of the Iberian upwelling system, upwelling events dominate during a short summer season and downwelling events the rest of the year. The ria response to the external forcing is central to nutrient supply and resultant plankton productivity that supports its high level of cultured mussel production. Intensive field studies in September 2006 and June 2007 captured a downwelling event and an upwelling event, respectively. Data from eight current profiler moorings and boat-based MiniBat/ADCP surveys provided an unprecedented quasi-synoptic view of the distribution of water masses and circulation patterns in any ria. In the outer ria, circulation was dominated by the introduction of wind-driven alongshore flow from the external continental shelf through the ria entrances and its interaction with the topography. In the middle ria, circulation was primarily related to the upwelling/downwelling cycle, with a cool, salty and dense lower layer penetrating to the inner ria during upwelling over the shelf. A warmer, lower salinity and less dense surface layer of coastal waters flowed inward during downwelling. Without external forcing, the inner ria responded primarily to tides and buoyancy changes related to land runoff. Under both upwelling and downwelling conditions, the flushing of the ria involved shelf responses to wind pulses. Their persistence for a few days was sufficient to allow waters from the continental shelf to penetrate the innermost ria. Longer term observations supported by numerical modeling are required to confirm the generality of such flushing events in the ria and determine their typical frequency, while comparative studies should explore how these scenarios fit into the range of conditions experienced in other semi-enclosed bays.

Biomonitoring of air pollutants by using lichens (Evernia prunastri) in areas between Kenitra and Mohammedia cities in Morocco.

In this study, Evernia prunastri, a lichen growing in its natural habitat in Morocco was analysed for the concentration of five heavy metals (Fe, Pb, Zn, Cu and Cr) from eleven sites between Kenitra and Mohammedia cities. The control site was Dar Essalam, an isolated area with low traffic density and dense vegetation. In the investigated areas, the concentration of heavy metals was correlated with vehicular traffic, industrial activity and urbanization. The total metal concentration was highest in Sidi Yahya, followed by Mohammedia and Bouznika. The coefficient of variation was higher for Pb and lower for Cu, Zn and Fe. The concentrations of most heavy metals in the thalli differed significantly between sites (p<0.01). Principal component analysis (PCA) revealed a significant correlation between heavy metal accumulation and atmospheric purity index. This study demonstrated also that the factors most strongly affecting the lichen flora were traffic density, the petroleum industry and paper factories in these areas. Overall, these results suggest that the index of atmospheric purity and assessment of heavy metals in lichen thalli are good indicators of the air quality at the studied sites.

Sexual cannibalism and mate choice

Sexual cannibalism, where a female kills and consumes a courting male, represents an extreme form of sexual conflict and has been proposed as a mechanism of mate choice. We evaluate the evidence for mate choice through premating sexual cannibalism via mate rejection, other indirect mechanisms of mate 'choice' and choice in postmating sexual cannibalism. We highlight a paucity of investigations, particularly of field studies, and note gaps in our knowledge. There is empirical support for the size-dependent sexual cannibalism predicted by mate choice through premating sexual cannibalism. This may represent mate choice operating on absolute male size but it could be a by-product of female foraging behaviour and greater vulnerability of relatively smaller males. Thus, indirect mate choice is as plausible an explanation of size-dependent sexual cannibalism as is direct mate choice based on discrimination of male traits. Direct female choice, mediated through premating sexual cannibalism, has yet to be demonstrated. We suggest a framework for distinguishing direct and indirect choice and note an absence of information on which to test it. There is evidence for sequential mate choice in postmating sexual cannibalism, but the nature or basis of the female's discriminatory behaviour remains unclear. Costs and long-term fitness benefits of the putative mate choice have been largely ignored. Reversed sexual cannibalism, in which the male eats the female, presumably occurs when the gain from food is high and potential gain from mating low and probably has little to do with mate choice. (c) 2005 The Association for the Study of Animal Behaviour Published by Elsevier Ltd. All rights reserved.