Spatial dynamics of the knob-tailed gecko Nephrurus stellatus in a fragmented agricultural landscape

In fragmented landscapes, a species' dispersal ability and response to habitat condition are key determinants of persistence. To understand the relative importance of dispersal and condition for survival of Nephrurus stellatus (Gekkonidae) in southern Australia, we surveyed 92 woodland remnants three times. This gecko favours early post-fire succession conditions so may be at risk of extinction in the long-unburnt agricultural landscape. Using N-mixture models, we compared the influence of four measures of isolation, patch area and two habitat variables on the abundance and occurrence of N. stellatus, while taking into account detection probability. Patch occupancy was high, despite the long-term absence of fire from most remnants. Distance to the nearest occupied site was the most informative measure of patch isolation, exhibiting a negative relationship with occupancy. Distance to a nearby conservation park had little influence, suggesting that mainland-island metapopulation dynamics are not important. Abundance and occurrence were positively related to %-cover of spinifex (Triodia), indicating that niche-related factors may also contribute to spatial dynamics. Patterns of patch occupancy imply that N. stellatus has a sequence of spatial dynamics across an isolation gradient, with patchy populations and source-sink dynamics when patches are within 300 m, metapopulations at intermediate isolation, and declining populations when patches are separated by >1-2 km. Considering the conservation needs of the community, habitat condition and connectivity may need to be improved before fire can be reintroduced to the landscape. We speculate that fire may interact with habitat degradation and isolation, increasing the risk of local extinctions.

Simple rules to contain an invasive species with a complex life cycle and high dispersal capacity

Designing practical rules for controlling invasive species is a challenging task for managers, particularly when species are long-lived, have complex life cycles and high dispersal capacities. Previous findings derived from plant matrix population analyses suggest that effective control of long-lived invaders may be achieved by focusing on killing adult plants. However, the cost-effectiveness of managing different life stages has not been evaluated. We illustrate the benefits of integrating matrix population models with decision theory to undertake this evaluation, using empirical data from the largest infestation of mesquite (Leguminosae: Prosopis spp) within Australia. We include in our model the mesquite life cycle, different dispersal rates and control actions that target individuals at different life stages with varying costs, depending on the intensity of control effort. We then use stochastic dynamic programming to derive cost-effective control strategies that minimize the cost of controlling the core infestation locally below a density threshold and the future cost of control arising from infestation of adjacent areas via seed dispersal. Through sensitivity analysis, we show that four robust management rules guide the allocation of resources between mesquite life stages for this infestation: (i) When there is no seed dispersal, no action is required until density of adults exceeds the control threshold and then only control of adults is needed; (ii) when there is seed dispersal, control strategy is dependent on knowledge of the density of adults and large juveniles (LJ) and broad categories of dispersal rates only; (iii) if density of adults is higher than density of LJ, controlling adults is most cost-effective; (iv) alternatively, if density of LJ is equal or higher than density of adults, management efforts should be spread between adults, large and to a lesser extent small juveniles, but never saplings. Synthesis and applications.In this study, we show that simple rules can be found for managing invasive plants with complex life cycles and high dispersal rates when population models are combined with decision theory. In the case of our mesquite population, focussing effort on controlling adults is not always the most cost-effective way to meet our management objective.

No evidence for change in oviposition behaviour of Helicoverpa armigera (Hubner) (Lepidoptera : Noctuidae) after widespread adoption of transgenic insecticidal cotton

Cotton growing landscapes in Australia have been dominated by dual-toxin transgenic Bt varieties since 2004. The cotton crop has thus effectively become a sink for the main target pest, Helicoverpa armigera. Theory predicts that there should be strong selection on female moths to avoid laying on such plants. We assessed oviposition, collected from two cotton-growing regions, by female moths when given a choice of tobacco, cotton and cabbage. Earlier work in the 1980s and 1990s on populations from the same geographic locations indicated these hosts were on average ranked as high, mid and low preference plants, respectively, and that host rankings had a heritable component. In the present study, we found no change in the relative ranking of hosts by females, with most eggs being laid on tobacco, then cotton and least on cabbage. As in earlier work, some females laid most eggs on cotton and aspects of oviposition behaviour had a heritable component. Certainly, cotton is not avoided as a host, and the implications of these finding for managing resistance to Bt cotton are discussed.

Distribution and eradication of an exotic tephritid fruit fly in Australia: Relevance of invasion theory

Data from the eradication of the incursion of Bactrocera papayae Drew and Hancock (Dipt.: Tephritidae) in Australia (1995-1998) are used to assess the significance of various aspects of invasion theory, including the influence of towns on establishment, influence of propagule pressure on the pattern of establishment, and the existence of source-sink dynamics. Because there were no sentinel traps in place, considerable spread had occurred before the eradication campaign started. The distribution of fly density around the epicentre in the town of Cairns and a transect along the main traffic routes to the north and south fitted a Cauchy model with a tail having the same slope as a power model with an exponent of -2.4 extending to 160 km. The Cauchy model indicated that 50% of the flies on the transect would have occurred within 3.2 km of the epicentre, 90% within 13.2 km, and 99% within 60 km. The two major satellites at Mareeba (35 km from the epicentre in Cairns) and Mossman (65 km) were not used for the transect data and had respectively 15 and 30 times the density predicted by the model. The proportion of traps that caught flies (a measure of site occupancy) fell with distance from the epicentre. B. papayae was trapped consistently on only three of the 16 rainforest transects that were surveyed and these were relatively close to urban areas where eradication efforts were intense. Despite there being no eradication effort in the rainforest, the trends to extinction were similar to those in adjacent areas. The strategy of initially concentrating eradication efforts on the core and major satellites while maintaining a quarantine barrier at the airport and the boundaries of the infested area appears to be the key to the containment and rapid eradication of the incursion.

Inter-seasonal population dynamics and pest status of Bemisia tabaci (Gennadius) biotype B in an Australian cropping system

Bemisia tabaci, biotype B, commonly known as the silverleaf whitefly (SLW) is an alien species that invaded Australia in the mid-90s. This paper reports on the invasion ecology of SLW and the factors that are likely to have contributed to the first outbreak of this major pest in an Australian cotton cropping system, population dynamics of SLW within whitefly-susceptible crop (cotton and cucurbit) and non-crop vegetation (sowthistle, Sonchus spp.) components of the cropping system were investigated over four consecutive growing seasons (September-June) 2001/02-2004/05 in the Emerald Irrigation Area (EIA) of Queensland, Australia. Based on fixed geo-referenced sampling sites, variation in spatial and temporal abundance of SLW within each system component was quantified to provide baseline data for the development of ecologically sustainable pest management strategies. Parasitism of large (3rd and 4th instars) SLW nymphs by native aphelinid wasps was quantified to determine the potential for natural control of SLW populations. Following the initial outbreak in 2001/02, SLW abundance declined and stabilised over the next three seasons. The population dynamics of SLW is characterised by inter-seasonal population cycling between the non-crop (weed) and cotton components of the EIA cropping system. Cotton was the largest sink for and source of SLW during the study period. Over-wintering populations dispersed from weed host plant sources to cotton in spring followed by a reverse dispersal in late summer and autumn to broad-leaved crops and weeds. A basic spatial source-sink analysis showed that SLW adult and nymph densities were higher in cotton fields that were closer to over-wintering weed sources throughout spring than in fields that were further away. Cucurbit fields were not significant sources of SLW and did not appear to contribute significantly to the regional population dynamics of the pest. Substantial parasitism of nymphal stages throughout the study period indicates that native parasitoid species and other natural enemies are important sources of SLW mortality in Australian cotton production systems. Weather conditions and use of broad-spectrum insecticides for pest control are implicated in the initial outbreak and on-going pest status of SLW in the region.

Dynamics of Finnish starlings in 1951-2005 : from monitoring to population modelling

Farmland bird species have been declining in Europe. Many declines have coincided with general intensification of farming practices. In Finland, replacement of mixed farming, including rotational pastures, with specialized cultivation has been one of the most drastic changes from the 1960s to the 1990s. This kind of habitat deterioration limits the persistence of populations, as has been previously indicated from local populations. Integrated population monitoring, which gathers species-specific information of population size and demography, can be used to assess the response of a population to environment changes also at a large spatial scale. I targeted my analysis at the Finnish starling (Sturnus vulgaris). Starlings are common breeders in farmland habitats, but severe declines of local populations have been reported from Finland in the 1970s and 1980s and later from other parts of Europe. Habitat deterioration (replacement of pasture and grassland habitats with specialized cultivation areas) limits reproductive success of the species. I analysed regional population data in order to exemplify the importance of agricultural change to bird population dynamics. I used nestling ringing and nest-card data from 1951 to 2005 in order to quantify population trends and per capita reproductive success within several geographical regions (south/north and west/east aspects). I used matrix modelling, acknowledging age-specific survival and fecundity parameters and density-dependence, to model population dynamics. Finnish starlings declined by 80% from the end of the 1960s up to the end of the 1980s. The observed patterns and the model indicated that the population decline was due to the decline of the carrying capacity of farmland habitats. The decline was most severe in north Finland where populations largely become extinct. However, habitat deterioration was most severe in the southern breeding areas. The deteriorations in habitat quality decreased reproduction, which finally caused the decline. I suggest that poorly-productive northern populations have been partly maintained by immigration from the highly-productive southern populations. As the southern populations declined, ceasing emigration caused the population extinction in north. This phenomenon was explained with source sink population dynamics, which I structured and verified on the basis of a spatially explicit simulation model. I found that southern Finnish starling population exhibits ten-year cyclic regularity, a phenomenon that can be explained with delayed density-dependence in reproduction.

Turbulent vertical fluxes and air quality measured in urban air in Helsinki

There is a growing need to understand the exchange processes of momentum, heat and mass between an urban surface and the atmosphere as they affect our quality of life. Understanding the source/sink strengths as well as the mixing mechanisms of air pollutants is particularly important due to their effects on human health and climate. This work aims to improve our understanding of these surface-atmosphere interactions based on the analysis of measurements carried out in Helsinki, Finland. The vertical exchange of momentum, heat, carbon dioxide (CO2) and aerosol particle number was measured with the eddy covariance technique at the urban measurement station SMEAR III, where the concentrations of ultrafine, accumulation mode and coarse particle numbers, nitrogen oxides (NOx), carbon monoxide (CO), ozone (O3) and sulphur dioxide (SO2) were also measured. These measurements were carried out over varying measurement periods between 2004 and 2008. In addition, black carbon mass concentration was measured at the Helsinki Metropolitan Area Council site during three campaigns in 1996-2005. Thus, the analyzed dataset covered far, the most comprehensive long-term measurements of turbulent fluxes reported in the literature from urban areas. Moreover, simultaneously measured urban air pollution concentrations and turbulent fluxes were examined for the first time. The complex measurement surrounding enabled us to study the effect of different urban covers on the exchange processes from a single point of measurement. The sensible and latent heat fluxes closely followed the intensity of solar radiation, and the sensible heat flux always exceeded the latent heat flux due to anthropogenic heat emissions and the conversion of solar radiation to direct heat in urban structures. This urban heat island effect was most evident during winter nights. The effect of land use cover was seen as increased sensible heat fluxes in more built-up areas than in areas with high vegetation cover. Both aerosol particle and CO2 exchanges were largely affected by road traffic, and the highest diurnal fluxes reached 109 m-2 s-1 and 20 µmol m-2 s-1, respectively, in the direction of the road. Local road traffic had the greatest effect on ultrafine particle concentrations, whereas meteorological variables were more important for accumulation mode and coarse particle concentrations. The measurement surroundings of the SMEAR III station served as a source for both particles and CO2, except in summer, when the vegetation uptake of CO2 exceeded the anthropogenic sources in the vegetation sector in daytime, and we observed a downward median flux of 8 µmol m-2 s-1. This work improved our understanding of the interactions between an urban surface and the atmosphere in a city located at high latitudes in a semi-continental climate. The results can be utilised in urban planning, as the fraction of vegetation cover and vehicular activity were found to be the major environmental drivers affecting most of the exchange processes. However, in order to understand these exchange and mixing processes on a city scale, more measurements above various urban surfaces accompanied by numerical modelling are required.

Canopy processes, fluxes and microclimate in a pine forest

Interaction between forests and the atmosphere occurs by radiative and turbulent transport. The fluxes of energy and mass between surface and the atmosphere directly influence the properties of the lower atmosphere and in longer time scales the global climate. Boreal forest ecosystems are central in the global climate system, and its responses to human activities, because they are significant sources and sinks of greenhouse gases and of aerosol particles. The aim of the present work was to improve our understanding on the existing interplay between biologically active canopy, microenvironment and turbulent flow and quantify. In specific, the aim was to quantify the contribution of different canopy layers to whole forest fluxes. For this purpose, long-term micrometeorological and ecological measurements made in a Scots pine (Pinus sylvestris) forest at SMEAR II research station in Southern Finland were used. The properties of turbulent flow are strongly modified by the interaction between the canopy elements: momentum is efficiently absorbed in the upper layers of the canopy, mean wind speed and turbulence intensities decrease rapidly towards the forest floor and power spectra is modulated by spectral short-cut . In the relative open forest, diabatic stability above the canopy explained much of the changes in velocity statistics within the canopy except in strongly stable stratification. Large eddies, ranging from tens to hundred meters in size, were responsible for the major fraction of turbulent transport between a forest and the atmosphere. Because of this, the eddy-covariance (EC) method proved to be successful for measuring energy and mass exchange inside a forest canopy with exception of strongly stable conditions. Vertical variations of within canopy microclimate, light attenuation in particular, affect strongly the assimilation and transpiration rates. According to model simulations, assimilation rate decreases with height more rapidly than stomatal conductance (gs) and transpiration and, consequently, the vertical source-sink distributions for carbon dioxide (CO2) and water vapor (H2O) diverge. Upscaling from a shoot scale to canopy scale was found to be sensitive to chosen stomatal control description. The upscaled canopy level CO2 fluxes can vary as much as 15 % and H2O fluxes 30 % even if the gs models are calibrated against same leaf-level dataset. A pine forest has distinct overstory and understory layers, which both contribute significantly to canopy scale fluxes. The forest floor vegetation and soil accounted between 18 and 25 % of evapotranspiration and between 10 and 20 % of sensible heat exchange. Forest floor was also an important deposition surface for aerosol particles; between 10 and 35 % of dry deposition of particles within size range 10 30 nm occurred there. Because of the northern latitudes, seasonal cycle of climatic factors strongly influence the surface fluxes. Besides the seasonal constraints, partitioning of available energy to sensible and latent heat depends, through stomatal control, on the physiological state of the vegetation. In spring, available energy is consumed mainly as sensible heat and latent heat flux peaked about two months later, in July August. On the other hand, annual evapotranspiration remains rather stable over range of environmental conditions and thus any increase of accumulated radiation affects primarily the sensible heat exchange. Finally, autumn temperature had strong effect on ecosystem respiration but its influence on photosynthetic CO2 uptake was restricted by low radiation levels. Therefore, the projected autumn warming in the coming decades will presumably reduce the positive effects of earlier spring recovery in terms of carbon uptake potential of boreal forests.

Obliquely incident water waves against a vertical cliff

Utilizing the commutativity property of the Cartesian coordinate differential operators arising in the boundary conditions associated with the propagation of surface water waves against a vertical cliff, under the assumptions of linearized theory, the problem of obliquely incident surface waves is considered for solution. The case of normal incidence, handled by previous workers follow as a particular limiting case of the present problem, which exhibits a source/sink type behavior of the velocity potential at the shore-line. An independent method of attack is also presented to handle the case of normal incidence.

Multi-hop Cooperative Wireless Networks: Diversity Multiplexing Tradeoff and Optimal Code Design

We consider single-source single-sink (ss-ss) multi-hop relay networks, with slow-fading links and single-antenna half-duplex relay nodes. While two-hop cooperative relay networks have been studied in great detail in terms of the diversity-multiplexing tradeoff (DMT), few results are available for more general networks. In this paper, we identify two families of networks that are multi-hop generalizations of the two-hop network: K-Parallel-Path (KPP)networks and layered networks.KPP networks, can be viewed as the union of K node-disjoint parallel relaying paths, each of length greater than one. KPP networks are then generalized to KPP(I) networks, which permit interference between paths and to KPP(D) networks, which possess a direct link from source to sink. We characterize the DMT of these families of networks completely for K > 3. Layered networks are networks comprising of layers of relays with edges existing only between adjacent layers, with more than one relay in each layer. We prove that a linear DMT between the maximum diversity dmax and the maximum multiplexing gain of 1 is achievable for single-antenna fully-connected layered networks. This is shown to be equal to the optimal DMT if the number of relaying layers is less than 4.For multiple-antenna KPP and layered networks, we provide an achievable DMT, which is significantly better than known lower bounds for half duplex networks.For arbitrary multi-terminal wireless networks with multiple source-sink pairs, the maximum achievable diversity is shown to be equal to the min-cut between the corresponding source and the sink, irrespective of whether the network has half-duplex or full-duplex relays. For arbitrary ss-ss single-antenna directed acyclic networks with full-duplex relays, we prove that a linear tradeoff between maximum diversity and maximum multiplexing gain is achievable.Along the way, we derive the optimal DMT of a generalized parallel channel and derive lower bounds for the DMT of triangular channel matrices, which are useful in DMT computation of various protocols. We also give alternative and often simpler proofs of several existing results and show that codes achieving full diversity on a MIMO Rayleigh fading channel achieve full diversity on arbitrary fading channels. All protocols in this paper are explicit and use only amplify-and-forward (AF) relaying. We also construct codes with short block-lengths based on cyclic division algebras that achieve the optimal DMT for all the proposed schemes.Two key implications of the results in the paper are that the half-duplex constraint does not entail any rate loss for a large class of cooperative networks and that simple AF protocols are often sufficient to attain the optimal DMT

Removal of unwanted fluid

This work is concerned with the removal of unwanted fluid through the source-sink pair. The source consists of fluid issuing out of a nozzle in the form of a jet and the sink is a pipe that is kept some distance from the source pipe. Of concern is the percentage of source fluid sucked through the sink. The experiments have been carried in a large glass water tank. The source nozzle diameter is 6Â mm and the sink pipe diameter is either 10 or 20Â mm. The horizontal and vertical separations and angles between these source and sink pipes are adjustable. The flow was visualized using KMnO 4 dye, planer laser induced fluorescence and particle streak photographs. To obtain the effectiveness (that is percentage of source fluid entering the sink pipe), titration method is used. The velocity profiles with and without the sink were obtained using particle image velocimetry. The sink flow rate to obtain a certain effectiveness increase dramatically with lateral separation. The sink diameter and the angle between source and the sink axes don't influence effectiveness as much as the lateral separation.

Removal of unwanted fluid

This work is concerned with the removal of unwanted fluid through the source-sink pair. The source consists of fluid issuing out of a nozzle in the form of a jet and the sink is a pipe that is kept some distance from the source pipe. Of concern is the percentage of source fluid sucked through the sink. The experiments have been carried in a large glass water tank. The source nozzle diameter is 6 mm and the sink pipe diameter is either 10 or 20 mm. The horizontal and vertical separations and angles between these source and sink pipes are adjustable. The flow was visualized using KMnO4 dye, planer laser induced fluorescence and particle streak photographs. To obtain the effectiveness (that is percentage of source fluid entering the sink pipe), titration method is used. The velocity profiles with and without the sink were obtained using particle image velocimetry. The sink flow rate to obtain a certain effectiveness increase dramatically with lateral separation. The sink diameter and the angle between source and the sink axes don't influence effectiveness as much as the lateral separation.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

Caribbean connectivity: implications for Marine Protected Area Management. Proceedings of a Special Symposium, 9-11 November 2006, 59th Annual Meeting of the Gulf and Caribbean Fisheries Institute, Belize City, Belize

Executive Summary: Tropical marine ecosystems in the Caribbean region are inextricably linked through the movement of pollutants, nutrients, diseases, and other stressors, which threaten to further degrade coral reef communities. The magnitude of change that is occurring within the region is considerable, and solutions will require investigating pros and cons of networks of marine protected areas (MPAs), cooperation of neighboring countries, improved understanding of how external stressors degrade local marine resources, and ameliorating those stressors. Connectivity can be broadly defined as the exchange of materials (e.g., nutrients and pollutants), organisms, and genes and can be divided into: 1) genetic or evolutionary connectivity that concerns the exchange of organisms and genes, 2) demographic connectivity, which is the exchange of individuals among local groups, and 3) oceanographic connectivity, which includes flow of materials and circulation patterns and variability that underpin much of all these exchanges. Presently, we understand little about connectivity at specific locations beyond model outputs, and yet we must manage MPAs with connectivity in mind. A key to successful MPA management is how to most effectively work with scientists to acquire the information managers need. Oceanography connectivity is poorly understood, and even less is known about the shape of the dispersal curve for most species. Dispersal kernels differ for various systems, species, and life histories and are likely highly variable in space and time. Furthermore, the implications of different dispersal kernels on population dynamics and management of species is unknown. However, small dispersal kernels are the norm - not the exception. Linking patterns of dispersal to management options is difficult given the present state of knowledge. The behavioral component of larval dispersal has a major impact on where larvae settle. Individual larval behavior and life history details are required to produce meaningful simulations of population connectivity. Biological inputs are critical determinants of dispersal outcomes beyond what can be gleaned from models of passive dispersal. There is considerable temporal and spatial variation to connectivity patterns. New models are increasingly being developed, but these must be validated to understand upstream-downstream neighborhoods, dispersal corridors, stepping stones, and source/sink dynamics. At present, models are mainly useful for providing generalities and generating hypotheses. Low-technology approaches such as drifter vials and oceanographic drogues are useful, affordable options for understanding local connectivity. The “silver bullet” approach to MPA design may not be possible for several reasons. Genetic connectivity studies reveal divergent population genetic structures despite similar larval life histories. Historical stochasticity in reproduction and/or recruitment likely has important, longlasting consequences on present day genetic structure. (PDF has 200 pages.)

An examination of spatial and temporal genetic variation in walleye pollock (Theragra chalcogramma) using allozyme, mitochondrial DNA, and microsatellite data

We used allozyme, microsatellite, and mitochondrial DNA (mtDNA) data to test for spatial and interannual genetic diversity in wall-eye pollock (Theragra chalcogramma) from six spawning aggregations representing three geographic regions: Gulf of Alaska, eastern Bering Sea, and eastern Kamchatka. Interpopulation genetic diversity was evident primarily from the mtDNA and two allozyme loci (SOD-2*, MPI*). Permutation tests ˆindicated that FST values for most allozyme and microsatellite loci were not significantly greater than zero. The microsatellite results suggested that high locus polymorphism may not be a reliable indicator of power for detecting population differentiation in walleye pollock. The fact that mtDNA revealed population structure and most nuclear loci did not suggests that the effective size of most walleye pollock populations is large (genetic drift is weak) and migration is a relatively strong homogenizing force. The allozymes and mtDNA provided mostly concordant estimates of patterns of spatial genetic variation. These data showed significant genetic variation between North American and Asian populations. In addition, two spawning aggregations in the Gulf of Alaska, in Prince William Sound, and off Middleton Island, appeared genetically distinct from walleye pollock spawning in the Shelikof Strait and may merit management as a distinct stock. Finally, we found evidence of interannual genetic variation in two of three North American spawning aggregations, similar in magnitude to the spatial variation among North American walleye pol-lock. We suggest that interannual genetic variation in walleye pollock may be indicative of one or more of the following factors: highly variable reproductive success, adult philopatry, source-sink metapopulation structure, and intraannual variation (days) in spawning timing among genetically distinct but spatially identical spawning aggregates.