Spatial and temporal hydrodynamic and water quality modeling analysis of a large reservoir on the South Carolina (USA) coastal plain

Two-dimensional, 31-segment, 61-channel hydrodynamic and water quality models of Lake Marion (surface area 330.7 km²; volume 1548.3×10⁶ m³) were developed using the WASP5 modeling system. Field data from 1985 to 1990 were used to parameterize the models. Phytoplankton kinetic rates and constants were obtained from a related in situ study; others from modeling literature. The hydrodynamic model was calibrated to estimates of daily lake volume; the water quality model was calibrated for ammonia, nitrate, ortho-phosphate, dissolved oxygen, chlorophyll-a, biochemical oxygen demand, organic nitrogen, and organic phosphorus. Water quality calibration suggested the model characterized phytoplankton and nutrient dynamics quite well. The model was validated (Kolmogorov–Smirnov two-sample goodness-of-fit test at P<0.05) by reparameterizing the nutrient loading functions using an independent set of field data. The models identified several factors that may contribute to the spatial variability previously reported from other research in the reservoir, despite the superficial absence of complex structure. Sensitivity analysis of the phytoplankton kinetic rates suggest that study site-specific estimates were important for obtaining model fit to field data. Sediment sources of ammonia (10–60 mg m⁻² day⁻¹) and phosphate (1–6 mg m⁻² day⁻¹) were important to achieve model calibration, especially during periods of high temperatures and low dissolved oxygen. This sediment flux accounted for 78% (nitrogen) and 50% (phosphorus) of the annual load. Spatial and temporal variability in the lake, reflected in the calibrated and validated models, suggest that ecological factors that influence phytoplankton productivity and nutrient dynamics are different in various parts of the lake. The WASP5 model as implemented here does not fully accommodate the ecological variability in Lake Marion due to model constraints on the specification of rate constants. This level of spatial detail may not be appropriate for an operational reservoir model, but as a research tool the models are both versatile and useful.     

Bayesian hierarchical mixture models for otolith microchemistry analysis

The field of fisheries research commonly uses classical statistical classification methods to estimate the proportion of fish that return to natal spawning grounds to spawn. With the advent of otolith microchemical analysis, researchers are able to extract information from fish ear stones (otoliths) about the chemical composition of water in which fish have spent distinct periods of their lives. Here we present a method of analysis set in the Bayesian statistical paradigm which enables explicit incorporation of habitat information into the analysis. The ecological system is seen as arising from a mixture of disparate fish populations and information from the biological relationships inherent in otolith formation is exploited through the hierarchical model structure. We present the model and motivation, demonstrate the validity of the model through simulation studies, and conclude with an analysis of a data set from Lake Erie.     

A hydrodynamic and water quality modeling study of spatial and temporal patterns of phytoplankton growth in a stratified lake with buoyant incoming flow

A three-dimensional hydrodynamic and water quality model was applied to Lake Paldang, a lake in South Korea that is stratified by incoming flows. The spatial and temporal patterns of phytoplankton growth in this lake were determined from the model. The model was calibrated and verified using data measured under different hydrological conditions. The model results were in reasonable agreement with the field measurements, in both the calibration and verification phases. The distributions of water quality and residence time in the lake and phytoplankton response to changes in nutrient loads were examined with the model, and the influence of the hydrodynamics on phytoplankton response was analyzed. The simulation results indicated that Lake Paldang is an essentially phosphorus-limited system, but that phytoplankton growth is limited by low water temperature and short residence time during the winter and the summer monsoon period, respectively. The results of sensitivity analyses also suggested that the hydrodynamics within the lake may have an indirect influence on phytoplankton responses to changes in the limiting nutrient loads, and that reducing phosphorus loading from Kyoungan Stream should be a high priority policy for controlling algal blooms during the pre- and post-monsoon periods. From this study, it was concluded that the three-dimensional water quality model incorporating hydrodynamic processes could successfully simulate phytoplankton response to changes in nutrient loads and that it could become a useful tool for identifying the essential factors determining phytoplankton growth and for developing the best management policy for algal blooms in Lake Paldang.     

Study of the nutrient and plankton dynamics in Lake Tanganyika using a reduced-gravity model

An eco-hydrodynamic (ECOH) model is proposed for Lake Tanganyika to study the plankton productivity. The hydrodynamic sub-model solves the non-linear, reduced-gravity equations in which wind is the dominant forcing. The ecological sub-model for the epilimnion comprises nutrients, primary production, phytoplankton biomass and zooplankton biomass. In the absence of significant terrestrial input of nutrients, the nutrient loss is compensated for by seasonal, wind-driven, turbulent entrainment of nutrient-rich hypolimnion water into the epilimnion, which gives rise to high plankton productivity twice in the year, during the transition between two seasons. Model simulations predict well the seasonal contrasts of the measured physical and ecological parameters. Numerical tests indicate that the half saturation constant for grazing by zooplankton and the fish predation rate on zooplankton affect the zooplankton biomass measurably more than that of phytoplankton biomass. This work has implications for the application of this model to predict the climatological biological productivity of Lake Tanganyika.     

Structural changes in lake functioning induced from nutrient loading and climate variability

Climate variability is increasingly recognized as an important regulatory factor, capable of influencing the structural properties of aquatic ecosystems. Lakes appear to be particularly sensitive to the ecological impacts of climate variability, and several long time series have shown a close coupling between climate, lake thermal properties and individual organism physiology, population abundance, community structure, and food web dynamics. Thus, understanding the complex interplay among meteorological forcing, hydrological variability, and ecosystem functioning is essential for improving the credibility of model-based water resources/fisheries management. Our objective herein is to examine the relative importance of the ecological mechanisms underlying plankton seasonal variability in Lake Washington, Washington State (USA), over a 35-year period (1964–1998). Our analysis is founded upon an intermediate complexity plankton model that is used to reproduce the limiting nutrient (phosphate)–phytoplankton–zooplankton–detritus (particulate phosphorus) dynamics in the lake. Model parameterization is based on a Bayesian calibration scheme that offers insights into the degree of information the data contain about model inputs and allows obtaining predictions along with uncertainty bounds for modeled output variables. The model accurately reproduces the key seasonal planktonic patterns in Lake Washington and provides realistic estimates of predictive uncertainty for water quality variables of environmental management interest. A principal component analysis of the annual estimates of the underlying ecological processes highlighted the significant role of the phosphorus recycling stemming from the zooplankton excretion on the planktonic food web variability. We also identified a moderately significant signature of the local climatic conditions (air temperature) on phytoplankton growth (r = 0.41), herbivorous grazing (r = 0.38), and detritus mineralization (r = 0.39). Our study seeks linkages with the conceptual food web model proposed by Hampton et al. [Hampton, S.E., Scheuerell, M.D., Schindler, D.E., 2006b. Coalescence in the Lake Washington story: interaction strengths in a planktonic food web. Limnol. Oceanogr. 51, 2042–2051.] to emphasize the “bottom-up” control of the Lake Washington plankton phenology. The posterior predictive distributions of the plankton model are also used to assess the exceedance frequency and confidence of compliance with total phosphorus (15 μg L⁻¹) and chlorophyll a (4 μg L⁻¹) threshold levels during the summer-stratified period in Lake Washington. Finally, we conclude by underscoring the importance of explicitly acknowledging the uncertainty in ecological forecasts to the management of freshwater ecosystems under a changing global environment.     

Modelling the components of the vertical attenuation of ultraviolet radiation in a wetland lake ecosystem

The extinction of solar UV (290–400 nm) radiation in aquatic ecosystems is a complex phenomena. In this paper, we examine and model the attenuation of UV radiation in a shallow lake ecosystem. In particular we focus our analysis on the specific role of the fractions of dissolved and particulate matter in the water column in the attenuation of radiation. This analysis is aimed by representing the spatial distribution of each fraction making it possible to evidence the spatial variation in habitat quality. In situ and laboratory measurement are used to elaborate a UV attenuation model. The attenuation model distinguishes between the contribution of particulate and dissolved matter in the attenuation of the solar UV flux. In the studied wetland lake (Laguna Iberá, sub-tropical latitude, Argentina) the importance of dissolved matter is dominated in the UVB solar spectrum (290–320 nm) but the effects of the particulate fraction are not negligible, in particular in UVA (320–400 nm). The spatial representation of model results demonstrate the non homogeneous nature between attenuation of the two fractions. Local and global environmental change can have important impacts on dissolved and particulate matter concentrations, which cam have ecological consequences in relation to the high flux of incoming UV radiation. The model developed to examine the relative attenuation of the dissolved and particulate fractions and is a useful instrument to identify the role that these fractions have on the optical characteristics of aquatic ecosystems.     

Phytoplankton patchiness and their role in the modelled productivity of a large, seasonally stratified lake

Phytoplankton concentration in Lake Kinneret (Israel) has varied up to 10-fold in space and time, with horizontal patches ranging from a couple of kilometres to a basin scale. Previous studies have used a 1D model to reproduce the temporal evolution of physical and biogeochemical variables in this lake. The question that arises then is how appropriate is a 1D approach to represent the dynamic of a spatially heterogeneous system, where there are non-linear dependencies between variables. Field data, a N-P-Z model coupled to both a 1D and a 3D hydrodynamic model, a 1D diffusion-reaction equation and scaling analysis are used to understand the role of spatial variability, expressed as phytoplankton patchiness, in the modelling of primary production. The analysis and results are used to investigate the effect of horizontal variability in the forcing and in the free mechanisms that affect the growth of patterns. The study shows that the use of averaged properties in a 1D approach may produce misleading results in the presence of localised patches, in terms of both concentration and composition of phytoplankton. The reason lies in the fact that the calibration process of ecological parameters in the 1D model appears to be site and process specific. That is, it depends on the pattern's characteristics and the underlying physical processes causing them. And this is a critical point for the success of numerical simulations under spatial variability. In this study, it is also shown that a length scale based on diffusion and growth rate of phytoplankton could be used as a criterion to assess the appropriateness of the 1D assumption.     

Seasonal rhythm in a Mediterranean coastal fish community as monitored by a cabled observatory

To adequately describe seasonal rhythms in habitat utilization by fish communities, observations are needed that occur at a relatively high frequency and over large temporal windows. For the first time, images collected from a cabled video-observatory (the western Mediterranean expandable SEAfloor OBservatory; OBSEA) were utilized to examine habitat utilization by coastal fishes. Over the course of a year, hourly digital images of the water column and an artificial reef were obtained during daylight hours. For each image, the total number of species, an estimate of the abundance of individuals, and the Shannon Diversity Index were quantified. A total of 22 fish species commonly associated with nearshore western Mediterranean habitats were identified, and significant spatial (water column vs. artificial reef), seasonal (spring, summer, autumn, and winter), and daily (morning, midday, and sunset) differences in habitat utilization occurred. Four species (Diplodus vulgaris, Diplodus cervinus, Diplodus sargus, and Scorpaena porcus) were associated with the artificial reef, while one species (Pagrus pagrus) chiefly occurred in the water column. Chromis chromis and Diplodus annularis occurred at the site more frequently in the winter and autumn, respectively, while 14 other species utilized the site more frequently in either the spring (8 spp) or summer (6 spp). In addition, Dentex dentex and Spicara maena occurred more frequently at sunset, while D. cervinus displayed a crepuscular rhythm (occurring more frequently in the morning and at sunset). Species diversity was highest in the summer and lowest in the winter. If not taken into account when planning in situ sampling, such seasonal and/or diel differences may lead to spurious estimates of population sizes and biodiversity. We suggest that cabled video-observatories offer a non-invasive and reliable technology for faunistic sampling and population assessment in coastal water of the Mediterranean and likely elsewhere.     

Appreciating Ecological Complexity: Habitat Contours as a Conceptual Landscape Model

Abstract:  Organisms respond to their surroundings at multiple spatial scales, and different organisms respond differently to the same environment. Existing landscape models, such as the "fragmentation model" (or patch-matrix-corridor model) and the "variegation model," can be limited in their ability to explain complex patterns for different species and across multiple scales. An alternative approach is to conceptualize landscapes as overlaid species-specific habitat contour maps. Key characteristics of this approach are that different species may respond differently to the same environmental conditions and at different spatial scales. Although similar approaches are being used in ecological modeling, there is much room for habitat contours as a useful conceptual tool. By providing an alternative view of landscapes, a contour model may stimulate more field investigations stratified on the basis of ecological variables other than human-defined patches and patch boundaries. A conceptual model of habitat contours may also help to communicate ecological complexity to land managers. Finally, by incorporating additional ecological complexity, a conceptual model based on habitat contours may help to bridge the perceived gap between pattern and process in landscape ecology. Habitat contours do not preclude the use of existing landscape models and should be seen as a complementary approach most suited to heterogeneous human-modified landscapes.     

Monitoring and estimating the flow conditions and fish presence probability under various flow conditions at reach scale using genetic algorithms and kriging methods

The combination of current velocity and water depth influences stream flow conditions, and fish activities prefer particular flow conditions. This study develops a novel optimal flow classification method for identifying types of stream flow based on the current velocity and the water depth using a genetic algorithm. It is applied to the Datuan stream in northern Taiwan. Fish were sampled and their habitat investigated at the study site during the spring, summer, fall and winter of 2008-2009. The current velocity, water depth and maps of the presence probability of fish were estimated by ordinary and indicator kriging. The optimal classification results were compared with the classification results obtained using the Froude number and empirical methods. The flow classification results demonstrate that the proposed optimal flow classification method that considers depth-velocity and optimally identified criteria for classifying flow types, yields a current velocity and water depth of 0.32 (m/s) and 0.29 (m), respectively, and classifies the flow conditions in the study area as pool, run, riffle and slack. The variography results of the current velocity and the water depth data reveal that seasonal flows are not spatially stationary among seasons in the study area. Kriging methods and a two-dimensional hydrodynamic model (River 2D) with empirical and optimal flow classification methods are more effective than the Froude number method in classifying flow conditions in the study area. The flow condition classifications and probability maps were generated by River 2D, ordinary kriging and indicator kriging, to quantify the flow conditions preferred by Sicyopterus japonicus in the study area. However, the proposed optimal classification method with kriging and River 2D is an effective alternative method for mapping flow conditions and determining the relationship between flow and the presence probability of target fish in support of stream restoration.     

Two-dimensional habitat modeling of Chinese sturgeon spawning sites

Since the construction of the Gezhouba Dam in the 1980s, the number of Chinese sturgeon in the Yangtze River has been rapidly declining. The Gezhouba Dam has cutoff the migration path of these sturgeon, resulting in an overall reduction of suitable sturgeon habitat. This paper describes a habitat suitability index model that is used to evaluate the impacts of the Gezhouba Dam and Three Gorges Project on Chinese sturgeon spawning sites. Based on research concerning the reproduction characteristics of Chinese sturgeon, ten major ecological factors that influence reproduction were analyzed, including: water temperature, velocity, water depth, substrate, suspended sediment concentration, and the amount of egg predatory fish. The suitability index (SI) curves based on these ten ecological factors were obtained, and a habitat suitability function was developed. A two-dimensional mathematical model was also created to simulate and predict physical habitat situation (such as hydraulic, sediment, and substrate) of the Chinese sturgeon. By coupling the habitat suitability function and a two-dimensional mathematical model, a habitat suitability index model for Chinese sturgeon was established. The habitat suitability index model was validated by comparing measured data with predictions from the model. These comparisons showed that the computed results agreed well with the measured results, and the high calculated habitat suitability index (HSI) corresponded to high measured quantity of eggs per unit (1000 m³) discharge (CPUE_d). The calculated habitat suitability index for Chinese sturgeon also showed that the habitat suitability index was better in 1999, before the impoundment of the Three Gorges Project, compared with the habitat suitability in 2003. Simulation results of different discharges from Gezhouba Dam predicted that flow discharges between 10,000 and 30,000 m³/s were most suitable for sturgeon spawning.     

滇池湖泊生态系统水动力学模拟

针对滇池富营养化状况,依据滇池生态环境建设规划,在以前研究工作的基础上,综合考虑湿地植物、底泥、水域条件变化,根据水动力学原理建立了垂向平均的二维生态系统水动力学模型,模拟了滇池现有水体、拆除防浪堤水域扩展后水体和滇池湿地建成后总水体3种情况的流场和浓度场,以及入滇河道达标排放后的水质浓度场,分析了湿地植物和水域扩展对滇池流场的影响,以及水域扩展、湿地植物、入湖河流达标排放对滇池水质的影响。结果表明,水域拓展对滇池流场的影响是局部的,湿地植物对滇池流场的影响主要发生在湿地区域。水域拓展对滇池水质具有一定影响,TN、TP平均质量浓度由2.08和0.19 mg.L-1(现有水体)降至1.69和0.16 mg.L-1;人工湿地建设对滇池水体TN、TP浓度影响较大,湿地建成后平均质量浓度降至0.76和0.05 mg.L-1,基本达到GB 3838—2002《地表水环境质量标准》Ⅲ类水质标准;人工湿地建设完成并且入滇池河流达标排放后,TN、TP平均质量浓度分别降至0.17和0.01 mg.L-1,水质明显改善。因此,人工湿地建设及入滇池河流的达标排放对滇池生态环境的改善具有重要作用。     

Estimating Plecoglossus altivelis altivelis migration using a mass balance model expressed by hydrological distribution parameters in a major limpid river basin in Japan

We examined the influence of several hydrological and meteorological parameters on the migratory movements of ayu Plecoglossus altivelis altivelis in central Japan. When comprehensively evaluating rivers and ayu behaviour on a catchment scale, the subjects of analysis typically include human activities and hydrological and meteorological phenomena. However, limiting analyses to such factors may be too restrictive when human activities are being conducted. Accordingly, we incorporated a biological viewpoint into the evaluation method, analysing hydrological data (river discharge, river water temperature, sea water temperature) to determine watershed characteristics and examining the relationship between these characteristics and the habitat conditions of ayu. Then we constructed a numerical model for ayu migratory runs that incorporated ayu ecology and watershed characteristics. Analyses of ayu movements from a lower estuarine dam demonstrated that downstream displacements were associated with high water flows of more than 200 m³ s⁻¹ at the beginning of summer. We conclude that it is important to consider the effects of environmental parameters on the movements of different fish species to understand the causes of spatial variation in fish distribution in lowland rivers.     

Crustacean habitat potential mapping in a tidal flat using remote sensing and GIS

To make a macrofaunal (crustacean) habitat potential map, the spatial distribution of ecological variables in the Hwangdo tidal flat, Korea, was explored. Spatial variables were mapped using remote sensing and a geographic information system (GIS) combined with field observations. A frequency ratio (FR) and logistic regression (LR) model were employed to map the macrofauna potential area for the Ilyoplax dentimerosa, a crustacean species. Spatial variables affecting the tidal macrofauna distribution were selected based on abundance and biomass and used within a spatial database derived from remotely sensed data of various types of sensors. The spatial variables included the intertidal digital elevation model (DEM), slope, distance from a tidal channel, tidal channel density, surface sediment facies, spectral reflectance of the near infrared (NIR) bands and the tidal exposure duration. The relation between the I. dentimerosa and each spatial variable was calculated using the FR and LR. The species was randomly divided into a training set (70%) to analyse habitat potential using FR and LR and a test set (30%) to validate the predicted habitat potential map. The relations were overlaid to produce a habitat potential map with the species potential index (SPI) value for each pixel. The potential habitat maps were compared with the surveyed habitat locations such as validation data set. The comparison results showed that the LR model (accuracy is 85.28%) is better in prediction than the FR (accuracy is 78.96%) model. The performance of models gave satisfactory accuracies. The LR provides the quantitative influence of variables on a potential habitat of species; otherwise, the FR shows the quantitative influence of a class in each variable. The combination of a GIS-based frequency ratio and logistic regression models and remote sensing with field observations is an effective method to determine locations favorable for macrofaunal species occurrences in a tidal flat.     

3D model for a secondary facultative pond

This paper describes a comprehensive model of wastewater treatment in secondary facultative ponds, which combines 3D hydrodynamics with a mechanistic water quality model. The hydrodynamics are based on the Navier-Stokes equation for incompressible fluids under shallow water and Boussinesq assumptions capturing the flow dynamics along length, breadth and depth of the pond. The water quality sub model is based on the Activated Sludge Model (ASM) concept, describing COD and nutrient removal as function of bacterial growth following Monod kinetics, except for Escherichia coli removal, which was modelled as first order decay. The model was implemented in the Delft3D software and was used to evaluate the effect of wind and the addition of baffles on the water flow pattern, temperature profiles in the pond and treatment efficiency. In contrast to earlier models reported in the literature, our simulation results did not show any significant improvement in COD removal (based on the ASM concept) with addition of baffles or under intermittent wind-induced mixing. However, E. coli removal efficiency, based on a first order decay approach, showed a fair improvement in the presence of baffles or intermittent wind-induced mixing. Furthermore, simulations with continuous wind effect showed a decrease in removal efficiency for COD but a further increase in E. coli removal efficiency. Such contrasting results for two different approaches in modelling could indicate that the first order decay concept might not be appropriate to describe all the interactions between biochemical processes in a pond. However, these interpretations remain theoretical, as the model needs validation with field data.     

Consideration of fuzziness: Is it necessary in modelling fish habitat preference of Japanese medaka (Oryzias latipes)?

The present study aims to clarify the necessity and effectiveness of considering fuzziness in modelling fish habitat preference, and the advantages which would be achieved by considering it. For this purpose, genetic algorithm (GA) optimized habitat preference models under three different levels of fuzzification were compared with regard to prediction ability of the habitat use of Japanese medaka (Oryzias latipes) dwelling in agricultural canals in Japan. Field surveys were conducted in agricultural canals in Japan to establish a relationship between fish habitat preference and physical environments of water depth, current velocity, lateral cover ratio and percent vegetation coverage. The habitat preference models employed for testing the fuzzy-based approach were category model, fuzzy habitat preference model, and fuzzy habitat preference model with fuzzy inputs. All the models were developed at 50 different initial conditions. The effectiveness of the fuzzification in fish habitat modelling was assessed by comparing mean square error and standard deviation of the models, and fluctuation in habitat preference curves evaluated by each model. As a result, the effect of fuzzification appeared as smoother curves and was found to reduce fluctuation in habitat preference curves in proportion to the level of fuzzification. The smooth curves would be appropriate for expressing uncertainty in habitat preference of the fish, by which fuzzy habitat preference model with fuzzy input achieve the best prediction ability among the models. In conclusion, the present study revealed that there are two advantages of fuzzification: reducing fluctuations in habitat preference evaluation and improving prediction ability of the model. Therefore, the consideration of fuzziness would be appropriate for representing fish habitat preference under natural conditions.     

A spatial model of white sturgeon rearing habitat in the lower Columbia River,USA

Concerns over the potential effects of in-water placement of dredged materials prompted us to develop a GIS-based model that characterizes in a spatially explicit manner white sturgeon Acipenser transmontanus rearing habitat in the lower Columbia River, USA. The spatial model was developed using water depth, riverbed slope and roughness, fish positions collected in 2002, and Mahalanobis distance (D²). We created a habitat suitability map by identifying a Mahalanobis distance under which >50% of white sturgeon locations occurred in 2002 (i.e., high-probability habitat). White sturgeon preferred relatively moderate to high water depths, and low to moderate riverbed slope and roughness values. The eigenvectors indicated that riverbed slope and roughness were slightly more important than water depth, but all three variables were important. We estimated the impacts that fill might have on sturgeon habitat by simulating the addition of fill to the thalweg, in 3-m increments, and recomputing Mahalanobis distances. Channel filling simulations revealed that up to 9 m of fill would have little impact on high-probability habitat, but 12 and 15 m of fill resulted in habitat declines of ∼12% and ∼45%, respectively. This is the first spatially explicit predictive model of white sturgeon rearing habitat in the lower Columbia River, and the first to quantitatively predict the impacts of dredging operations on sturgeon habitat. Future research should consider whether water velocity improves the accuracy and specificity of the model, and to assess its applicability to other areas in the Columbia River.     

Effects of wind-driven circulation on river intrusion in Lake Tegel: modeling study with projection on transport of pollutants

Prediction of mixing intensity of water masses in riverine Lake Tegel (Berlin, Germany) can be used to trace the fate of pollutants that enter the lake through several inflows. Because the contributions of each inflow have not yet been quantified and because the lake features complex bathymetry and numerous islands, a hydrodynamic circulation model with high spatial resolution and dynamic wind forcing is useful. We applied the two-dimensional version of the Princeton Ocean Model to separate the influences of wind and river discharge on the currents and mass transport in Lake Tegel. For model validation, we compared the simulation results with 1?year of electrical conductivity data, which was used as a conservative tracer to distinguish between water from the River Havel and water supplied by a smaller second inflow. Calculation of currents alone is insufficient to investigate water exchanges between rivers and lakes, especially when several islands create multiple pathways for river intrusion. Therefore, mass transport simulations are applied. Our calculations based on archetypical scenarios indicate that the proportion of (polluted) water from the River Havel in the main basin of Lake Tegel fluctuates with river discharge and wind, which either amplify or neutralize each other.     

Meso-scale hydrodynamic and reproductive asynchrony affects the source–sink metapopulation structure of the coastal polychaete Pectinaria koreni

The polychaete Pectinaria koreni exhibits a complex life cycle characterized by non-overlapping generations and widespread larval dispersal. To explore how “local” metapopulation genetic structure varies spatially and temporally during population turnover, we combined observations on demography, larval dispersal through hydrodynamic modelling and population genetics of successive age cohorts in the Baie de Seine (eastern English Channel, France). Mature adults (March), newly settled (July) and later-stage juveniles (September) were sampled in 2003 on the edge and in the main demes of the metapopulation. Demes displayed an asynchronous dynamics due to variations in habitat quality affecting reproductive timing (e.g. three distinct spawning events observed) and in local larval supply linked to temporal fluctuations of hydrodynamism. Two-source populations were identified among dense areas with the greatest larval retention and self-recruitment rates: one with a single recruitment event, stable temporal genetic variation and a strong spatial genetic re-homogenization during turnover, and the other with two recruitment events and significant allele frequency changes over time. Sink habitats displayed single recruitment event and experienced strong inter-generational (juveniles vs adults) genetic changes due to genetic drift associated with strong winter mortalities. Altogether, results suggested that adult spawning asynchrony and sweepstakes reproductive success, together with genetic drift, played a greater role than hydrodynamics itself in determining effective recruitment rates at some sites and generating genetic patchiness.