Simulating property exchange in estuarine ecosystem models at ecologically appropriate scales

A computational scheme has been developed and tested to simulate property exchange by advection and dispersion in estuaries at time and space scales that are well suited to ecological and management simulations, but are coarse relative to the demands of physical hydrodynamic models. An implementation of the Regional Ocean Model System (ROMS) for the Providence River and Narragansett Bay (RI, USA) was used to determine property exchanges between the spatial elements of an ecological box model. The basis for the method is the statistical tabulation of numerical dye experiments done with the full ROMS physical model. The ROMS model domain was subdivided into fifteen coarse boxes, each with two vertical layers, defining 30 elements that were used for the box model simulations. Dye concentrations were set to arbitrary initial concentrations for all ROMS grids in the large elements, and the ROMS model was run for 24 h. The final distribution of the dye among the elements was used as a tracer for property exchange over that day and was used to develop an exchange matrix. Box model predictions of salinity over 77 days in each element compared favorably with ROMS simulated salinity averaged over the same spatial elements, although the disparity was greater in areas where large river inflows caused strong gradients in ROMS within elements assumed to be homogeneous in the box model. The 77-day simulation included periods of high and low river flow. Despite the large size of the spatial elements, dispersion artifacts were small, much less than the modeled daily exchanges. While others have taken a similar approach, we found a number of theoretical and practical considerations deserved careful attention for this approach to perform satisfactorily. Whereas the full ROMS model takes 9 days on a powerful computing cluster to compute the physics simulation for 77 days, the box model simulates physics and biology for the same interval in 5 s on a personal computer, and a full year in under 1 min. The exchange matrix mixing model is a fast, cost effective, and convenient way to simulate daily variation of complex estuarine physics in ecological modeling at appropriate scales of space and time.     

Transmitter attachment and release methods for short-term shark and stingray tracking on coral reefs

This study details a simple and cost-effective means of attaching acoustic transmitters to coral reef sharks and stingrays, which potentially allows for retrieval and reuse on completion of tracks. Between 2008 and 2011, galvanised timed releases were trialled in both static field tests and on blacktip reef sharks Carcharhinus melanopterus, cowtail Pastinachus atrus, and porcupine Urogymnus asperrimus rays in Coral Bay (?23°08′41″, 113°45′53″), Western Australia. The timed releases remained attached to animals for the duration required for tracking and in four out of five deployments transmitters were recovered after release from the animals. The use of modified Rototags for sharks and stainless steel darts for stingrays allowed rapid and effective attachment to animals, with limited impact on their welfare in the short term. External attachment for short-term tracking of coral reef-associated elasmobranchs should be considered as a complementary option to internal placement of transmitters in animals either by surgery or by ingestion.     

Transatlantic migration by post-breeding puffins: a strategy to exploit a temporarily abundant food resource?

The distribution of Atlantic puffins (Fratercula arctica) from Skellig Michael, south-west Ireland, was investigated using geolocation loggers between the 2010 and 2011 breeding seasons. All tracked birds travelled rapidly west into the North Atlantic at the end of the breeding season in August, with the majority undertaking transatlantic trips from Ireland to the Newfoundland-Labrador shelf. The furthest distance from the colony reached by each bird was not influenced by body mass or sex and was achieved in approximately 20 days. By October, all birds had moved back to the mid Atlantic where they remained resident until returning to the breeding colony. The most parsimonious explanation for the rapid, directed long-distance migration is that birds exploit the seasonally high abundance of prey [e.g., fish species such as capelin (Mallotus villosus) and sandlance (Ammodytes spp.)] off the Canadian coast, which is also utilised by large populations of North American seabirds at this time. Once the availability of this short-term prey resource has diminished, the tracked puffins moved back towards the north-east Atlantic. A relationship between relative abundance of puffins and zooplankton was found in all winter months, but after correcting for spatial autocorrelation, was only significant in November and January. Nevertheless, these results suggest a potential switch in diet from mainly fish during the breeding and early post-breeding periods to zooplankton over the remaining winter period. This study suggests that puffins from south-west Ireland have a long-distance migration strategy that is rare in breeding puffins from the UK and identifies a key non-breeding destination for puffins from Ireland. This has implications for the susceptibility of different breeding populations to the effects of possible climatic or oceanographic change.     

Tolerance of Hyas araneus zoea I larvae to elevated seawater PCO2 despite elevated metabolic costs

Early life stages of marine crustaceans respond sensitively to elevated seawater PCO₂. However, the underlying physiological mechanisms have not been studied well. We therefore investigated the effects of elevated seawater PCO₂ on oxygen consumption, dry weight, elemental composition, median developmental time (MDT) and mortality in zoea I larvae of the spider crab Hyas araneus (Svalbard 79°N/11°E; collection, May 2009; hatch, December 2009). At the time of moulting, oxygen consumption rate had reached a steady state level under control conditions. In contrast, elevated seawater PCO₂ caused the metabolic rate to rise continuously leading to a maximum 1.5-fold increase beyond control level a few days before moulting into the second stage (zoea II), followed by a pronounced decrease. Dry weight of larvae reared under high CO₂ conditions was lower than in control larvae at the beginning of the moult cycle, yet this difference had disappeared at the time of moulting. MDT of zoea I varied between 45 ± 1 days under control conditions and 42 ± 2 days under the highest seawater CO₂ concentration. The present study indicates that larval development under elevated seawater PCO₂ levels results in higher metabolic costs during premoulting events in zoea I. However, H. araneus zoea I larvae seem to be able to compensate for higher metabolic costs as larval MDT and survival was not affected by elevated PCO₂ levels.     

Female preferences for spectral call properties in the western genetic lineage of Cope's gray treefrog (Hyla chrysoscelis)

Female frogs discriminate among potential mates based on individual variation in male advertisement calls. While considerable data have accumulated allowing comparisons of female preference functions among species, we still lack fundamental knowledge about how and why the shapes of preference functions for particular call properties vary among populations within all but a few species. Here, we report results from a study aimed at describing female preference functions for spectral call properties in Cope's gray treefrog (Hyla chrysoscelis). Widespread throughout the eastern half of North America, Cope's gray treefrog is the diploid member of the cryptic diploid–tetraploid Hyla versicolor species complex, and its populations are divided into two distinct genetic lineages (eastern and western). In this study of a western lineage population, we recorded and analyzed the spectral properties of 1,000 advertisement calls from 50 males and conducted two-stimulus phonotaxis experiments to estimate a population-level preference function. Females preferred calls with average frequencies over calls with frequencies that were 2 or 3 semitones (1.4 or 2.1 standard deviations, respectively) lower than the population mean. We observed no behavioral discrimination between calls with average and higher-than-average frequencies. Preferences discriminating against low-frequency calls were weak and were abolished by attenuating the preferred average call by 3?dB. We discuss these results in light of previous studies of eastern lineage populations, geographic variation in female preference functions, and the potential adaptive value of discriminating against calls with low frequencies.     

Multiple predators indirectly alter community assembly across ecological boundaries

Models of habitat selection often assume that organisms choose habitats based on their intrinsic quality, regardless of the position of these habitats relative to low-quality habitats in the landscape. We created a habitat matrix in which high-quality (predator-free) aquatic habitat patches were positioned adjacent to (predator-associated) or isolated from (control) patches with single or two species of caged predators. After 16 days of colonization, larval insect abundance was reduced by 50% on average in both the predator and predator-associated treatments relative to isolated controls. Effects were largely similar among predator treatments despite variation in number of predator species, predator biomass, and whether predators were native or nonnative. Importantly, the strength of effects did not depend on whether predators were physically present. These results demonstrate that predator cues can cascade with equal strength across ecological boundaries, indirectly altering community assembly via habitat selection in intrinsically high-quality habitats.     

Habitat selection predicts genetic relatedness in an alpine ungulate

Landscape heterogeneity plays an integral role in shaping ecological and evolutionary processes. Despite links between the two disciplines, ecologists and population geneticists have taken different approaches to evaluating habitat selection, animal movement, and gene flow across the landscape. Ecologists commonly use statistical models such as resource selection functions (RSFs) to identify habitat features disproportionately selected by animals, whereas population genetic approaches model genetic differentiation according to the distribution of habitat variables. We combined ecological and genetic approaches by using RSFs to predict genetic relatedness across a heterogeneous landscape. We constructed sex- and season-specific resistance surfaces based on RSFs estimated using data from 102 GPS (global positioning system) radio-collared mountain goats (Oreamnos americanus) in southeast Alaska, USA. Based on mountain goat ecology, we hypothesized that summer and male surfaces would be the best predictors of relatedness. All individuals were genotyped at 22 microsatellite loci, which we used to estimate genetic relatedness. Summer resistance surfaces derived from RSFs were the best predictors of genetic relatedness, and winter models the poorest. Mountain goats generally selected for areas close to escape terrain and with a high heat load (a metric related to vegetative productivity and snow depth), while avoiding valleys. Male- and female-specific surfaces were similar, except for winter, for which male habitat selection better predicted genetic relatedness. The null models of isolation-by-distance and barrier only outperformed the winter models. This study merges high-resolution individual locations through GPS telemetry and genetic data, that can be used to validate and parameterize landscape genetics models, and further elucidates the relationship between landscape heterogeneity and genetic differentiation.     

Being young in a changing world: how temperature and salinity changes interactively modify the performance of larval stages of the barnacle <Emphasis Type="Italic">Amphibalanus improvisus</Emphasis>

The fate of key species, such as the barnacle Amphibalanus improvisus, in the course of global change is of particular interest since any change in their abundance and/or performance may entail community-wide effects. In the fluctuating Western Baltic, species typically experience a broad range of environmental conditions, which may preselect them to better cope with climate change. In this study, we examined the sensitivity of two crucial ontogenetic phases (naupliar, cypris) of the barnacle toward a range of temperature (12, 20, and 28°C) and salinity (5, 15, and 30 psu) combinations. Under all salinity treatments, nauplii developed faster at intermediate and high temperatures. Cyprid metamorphosis success, in contrast, was interactively impacted by temperature and salinity. Survival of nauplii decreased with increasing salinity under all temperature treatments. Highest settlement rates occurred at the intermediate temperature and salinity combination, i.e., 20°C and 15 psu. Settlement success of “naive” cyprids, i.e., when nauplii were raised in the absence of stress (20°C/15 psu), was less impacted by stressful temperature/salinity combinations than that of cyprids with a stress history. Here, settlement success was highest at 30 psu particularly at low and high temperatures. Surprisingly, larval survival was not highest under the conditions typical for the Kiel Fjord at the season of peak settlement (20°C/15 psu). The proportion of nauplii that ultimately transformed to attached juveniles was, however, highest under these “home” conditions. Overall, only particularly stressful combinations of temperature and salinity substantially reduced larval performance and development. Given more time for adaptation, the relatively smooth climate shifts predicted will probably not dramatically affect this species.     

A 200-year perspective on alternative stable state theory and lake management from a biomanipulated shallow lake

Multiple stressors to a shallow lake ecosystem have the ability to control the relative stability of alternative states (clear, macrophyte-dominated or turbid, algal-dominated). As a consequence, the use of remedial biomanipulations to induce trophic cascades and shift a turbid lake to a clear state is often only a temporary solution. Here we show the instability of short-term manipulations in the shallow Lake Christina (Minnesota, USA) is governed by the long-term state following a regime shift in the lake. During the modern, managed period of the lake, three top-down manipulations (fish kills) were undertaken inducing temporary (5-10 years) unstable clear-water states. Paleoecological remains of diatoms, along with proxies of primary production (total chlorophyll a and total organic carbon accumulation rate) and trophic state (total P) from sediment records clearly show a single regime shift in the lake during the early 1950s; following this shift, the functioning of the lake ecosystem is dominated by a persistent turbid state. We find that multiple stressors contributed to the regime shift. First, the lake began to eutrophy (from agricultural land use and/or increased waterfowl populations), leading to a dramatic increase in primary production. Soon after, the construction of a dam in 1936 effectively doubled the depth of the lake, compounded by increases in regional humidity; this resulted in an increase in planktivorous and benthivorous fish reducing phytoplankton grazers. These factors further conspired to increase the stability of a turbid regime during the modern managed period, such that switches to a clear-water state were inherently unstable and the lake consistently returned to a turbid state. We conclude that while top-down manipulations have had measurable impacts on the lake state, they have not been effective in providing a return to an ecosystem similar to the stable historical period. Our work offers an example of a well-studied ecosystem forced by multiple stressors into a new long-term managed period, where manipulated clear-water states are temporary, managed features.     

Randomized graph sampling

Randomized graph sampling (RGS) is an approach for sampling populations associated with or describable as graphs, when the structure of the graph is known and the parameter of interest is the total weight of the graph. RGS is related to, but distinct from, other graph-based approaches such as snowball and network sampling. Graph elements are clustered into walks that reflect the structure of the graph, as well as operational constraints on sampling. The basic estimator in RGS can be constructed as a Horvitz-Thompson estimator. I prove it to be design-unbiased, and also show design-unbiasedness of an estimator of the sample variance when walks are sampled with replacement. Covariates can be employed for variance reduction either through improved assignment of selection probabilities to walks in the design step, or through the use of alternative estimators during analysis. The approach is illustrated with a trail maintenance example, which demonstrates that complicated approaches to assignment of selection probabilities can be counterproductive. I describe conditions under which RGS may be efficient in practice, and suggest possible applications.