Understanding the performance of the FLake model over two African Great Lakes

The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu and Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding automatic weather stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration made it possible to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and water temperatures. At Lake Kivu, mixolimnion temperatures predicted by FLake were found to be sensitive both to minimal variations in the external parameters and to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may lead to a regime switch, from the correctly represented seasonal mixed layer deepening to either completely mixed or permanently stratified conditions from similar to 10 m downwards. In contrast, model temperatures were found to be robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterise tropical lake water surface temperatures within atmospheric prediction models. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60m during the main dry season is primarily due to enhanced lake evaporation and secondarily to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.

Morphological Diversity and the Roles of Contingency, Chance and Determinism in African Cichlid Radiations

Background: Deterministic evolution, phylogenetic contingency and evolutionary chance each can influence patterns of morphological diversification during adaptive radiation. In comparative studies of replicate radiations, convergence in a common morphospace implicates determinism, whereas non-convergence suggests the importance of contingency or chance. Methodology/Principal Findings: The endemic cichlid fish assemblages of the three African great lakes have evolved similar sets of ecomorphs but show evidence of non-convergence when compared in a common morphospace, suggesting the importance of contingency and/or chance. We then analyzed the morphological diversity of each assemblage independently and compared their axes of diversification in the unconstrained global morphospace. We find that despite differences in phylogenetic composition, invasion history, and ecological setting, the three assemblages are diversifying along parallel axes through morphospace and have nearly identical variance-covariance structures among morphological elements. Conclusions/Significance: By demonstrating that replicate adaptive radiations are diverging along parallel axes, we have shown that non-convergence in the common morphospace is associated with convergence in the global morphospace. Applying these complimentary analyses to future comparative studies will improve our understanding of the relationship between morphological convergence and non-convergence, and the roles of contingency, chance and determinism in driving morphological diversification.

Holocene climatic change and the development of the lake-effect snowbelt in Michigan, USA

Lake-effect snow is an important constraint on ecological and socio-economic systems near the North American Great Lakes. Little is known about the Holocene history of lake-effect snowbelts, and it is difficult to decipher how lake-effect snowfall abundance affected ecosystem development. We conducted oxygen-isotope analysis of calcite in lake-sediment cores from northern Lower Michigan to infer Holocene climatic variation and assess snowbelt development. The two lakes experience the same synoptic-scale climatic systems, but only one of them (Huffman Lake) receives a significant amount of lake-effect snow. A 177-cm difference in annual snowfall causes groundwater inflow at Huffman Lake to be 18O-depleted by 2.3‰ relative to O'Brien Lake. To assess when the lake-effect snowbelt became established, we compared calcite-δ18O profiles of the last 11,500 years from these two sites. The chronologies are based on accelerator-mass-spectrometry 14C ages of 11 and 17 terrestrial-plant samples from Huffman and O'Brien lakes, respectively. The values of δ18O are low at both sites from 11,500 to 9500 cal yr BP when the Laurentide Ice Sheet (LIS) exerted a dominant control over the regional climate and provided periodic pulses of meltwater to the Great Lakes basin. Carbonate δ18O increases by 2.6‰ at O'Brien Lake and by 1.4‰ at Huffman Lake between 9500 and 7000 cal yr BP, suggesting a regional decline in the proportion of runoff derived from winter precipitation. The Great Lakes snowbelt probably developed between 9500 and 5500 cal yr BP as inferred from the progressive 18O-depletion at Huffman Lake relative to O'Brien Lake, with the largest increase of lake-effect snow around 7000 cal yr BP. Lake-effect snow became possible at this time because of increasing contact between the Great Lakes and frigid arctic air. These changes resulted from enhanced westerly flow over the Great Lakes as the LIS collapsed, and from rapidly rising Great Lakes levels during the Nipissing Transgression. The δ18O difference between Huffman and O'Brien lakes declines after 5500 cal yr BP, probably because of a northward shift of the polar vortex that brought increasing winter precipitation to the entire region. However, δ18O remains depleted at Huffman Lake relative to O'Brien Lake because of the continued production of lake-effect snow.

Footprint of recycled water subsidies downwind of Lake Michigan

Continental evaporation is a significant and dynamic flux within the atmospheric water budget, but few methods provide robust observational constraints on the large-scale hydroclimatological and hydroecological impacts of this ‘recycled-water' flux. We demonstrate a geospatial analysis that provides such information, using stable isotope data to map the distribution of recycled water in shallow aquifers downwind from Lake Michigan. The δ2H and δ18O values of groundwater in the study region decrease from south to north, as expected based on meridional gradients in climate and precipitation isotope ratios. In contrast, deuterium excess (d = δ2H − 8 × δ18O) values exhibit a significant zonal gradient and finer-scale spatially patterned variation. Local d maxima occur in the northwest and southwest corners of the Lower Peninsula of Michigan, where ‘lake-effect' precipitation events are abundant. We apply a published model that describes the effect of recycling from lakes on atmospheric vapor d values to estimate that up to 32% of recharge into individual aquifers may be derived from recycled Lake Michigan water. Applying the model to geostatistical surfaces representing mean d values, we estimate that between 10% and 18% of the vapor evaporated from Lake Michigan is re-precipitated within downwind areas of the Lake Michigan drainage basin. Our approach provides previously unavailable observational constraints on regional land-atmosphere water fluxes in the Great Lakes Basin and elucidates patterns in recycled-water fluxes that may influence the biogeography of the region. As new instruments and networks facilitate enhanced spatial monitoring of environmental water isotopes, similar analyses can be widely applied to calibrate and validate water cycle models and improve projections of regional hydroecological change involving the coupled lake-atmosphere-land system. Read More: http://www.esajournals.org/doi/abs/10.1890/ES12-00062.1

LakeMIP Kivu: evaluating the representation of a large, deep tropical lake by a set of one-dimensional lake models

The African great lakes are of utmost importance for the local economy (fishing), as well as being essential to the survival of the local people. During the past decades, these lakes experienced fast changes in ecosystem structure and functioning, and their future evolution is a major concern. In this study, for the first time a set of one-dimensional lake models are evaluated for Lake Kivu (2.28°S; 28.98°E), East Africa. The unique limnology of this meromictic lake, with the importance of salinity and subsurface springs in a tropical high-altitude climate, presents a worthy challenge to the seven models involved in the Lake Model Intercomparison Project (LakeMIP). Meteorological observations from two automatic weather stations are used to drive the models, whereas a unique dataset, containing over 150 temperature profiles recorded since 2002, is used to assess the model’s performance. Simulations are performed over the freshwater layer only (60 m) and over the average lake depth (240 m), since salinity increases with depth below 60 m in Lake Kivu and some lake models do not account for the influence of salinity upon lake stratification. All models are able to reproduce the mixing seasonality in Lake Kivu, as well as the magnitude and seasonal cycle of the lake enthalpy change. Differences between the models can be ascribed to variations in the treatment of the radiative forcing and the computation of the turbulent heat fluxes. Fluctuations in wind velocity and solar radiation explain inter-annual variability of observed water column temperatures. The good agreement between the deep simulations and the observed meromictic stratification also shows that a subset of models is able to account for the salinity- and geothermal-induced effects upon deep-water stratification. Finally, based on the strengths and weaknesses discerned in this study, an informed choice of a one-dimensional lake model for a given research purpose becomes possible.

Origins of Shared Genetic Variation in African Cichlids

Cichlid fishes have evolved tremendous morphological and behavioral diversity in the waters of East Africa. Within each of the Great Lakes Tanganyika, Malawi, and Victoria, the phenomena of hybridization and retention of ancestral polymorphism explain allele sharing across species. Here, we explore the sharing of single nucleotide polymorphisms (SNPs) between the major East African cichlid assemblages. A set of approximately 200 genic and nongenic SNPs was ascertained in five Lake Malawi species and genotyped in a diverse collection of 160 species from across Africa. We observed segregating polymorphism outside of the Malawi lineage for more than 50% of these loci; this holds similarly for genic versus nongenic SNPs, as well as for SNPs at putative CpG versus non-CpG sites. Bayesian and principal component analyses of genetic structure in the data demonstrate that the Lake Malawi endemic flock is not monophyletic and that river species have likely contributed significantly to Malawi genomes. Coalescent simulations support the hypothesis that river cichlids have transported polymorphism between lake assemblages. We observed strong genetic differentiation between Malawi lineages for approximately 8% of loci, with contributions from both genic and nongenic SNPs. Notably, more than half of these outlier loci between Malawi groups are polymorphic outside of the lake. Cichlid fishes have evolved diversity in Lake Malawi as new mutations combined with standing genetic variation shared across East Africa.

Inheritance of female mating preference in a sympatric sibling species pair of Lake Victoria cichlids: implications for speciation

Female mate choice has often been proposed to play an important role in cases of rapid speciation, in particular in the explosively evolved haplochromine cichlid species flocks of the Great Lakes of East Africa. Little, if anything, is known in cichlid radiations about the heritability of female mating preferences. Entirely sympatric distribution, large ecological overlap and conspicuous differences in male nuptial coloration, and female preferences for these, make the sister species Pundamilia pundamilia and P. nyererei from Lake Victoria an ideally suited species pair to test assumptions on the genetics of mating preferences made in models of sympatric speciation. Female mate choice is necessary and sufficient to maintain reproductive isolation between these species, and it is perhaps not unlikely therefore, that female mate choice has been important during speciation. A prerequisite for this, which had remained untested in African cichlid fish, is that variation in female mating preferences is heritable. We investigated mating preferences of females of these sister species and their hybrids to test this assumption of most sympatric speciation models, and to further test the assumption of some models of sympatric speciation by sexual selection that female preference is a single-gene trait. We find that the differences in female mating preferences between the sister species are heritable, possibly with quite high heritabilities, and that few but probably more than one genetic loci contribute to this behavioural speciation trait with no apparent dominance. We discuss these results in the light of speciation models and the debate about the explosive radiation of cichlid fishes in Lake Victoria.

The genomic substrate for adaptive radiation in African cichlid fish

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Process and pattern in cichlid radiations - inferences for understanding unusually high rates of evolutionary diversification

The cichlid fish radiations in the African Great Lakes differ from all other known cases of rapid speciation in vertebrates by their spectacular trophic diversity and richness of sympatric species, comparable to the most rapid angiosperm radiations. I review factors that may have facilitated these radiations and compare these with insights from recent work on plant radiations. Work to date suggests that it was a coincidence of ecological opportunity, intrinsic ecological versatility and genomic flexibility, rapidly evolving behavioral mate choice and large amounts of standing genetic variation that permitted these spectacular fish radiations. I propose that spatially orthogonal gradients in the fit of phenotypes to the environment facilitate speciation because they allow colonization of alternative fitness peaks during clinal speciation despite local disruptive selection. Such gradients are manifold in lakes because of the interaction of water depth as an omnipresent third spatial dimension with other fitness-relevant variables. I introduce a conceptual model of adaptive radiation that integrates these elements and discuss its applicability to, and predictions for, plant radiations.

Explosive speciation rates and unusual species richness in haplochromine cichlid fishes: Effects of sexual selection

Ancient lakes are often unusually species rich, mostly as a result of radiation and species-flock formation having taken place in only one or a few of many taxa present. Understanding why some taxa radiate and others do not is at the heart of understanding biodiversity. In this chapter I discuss possible explanations for disproportionally large species numbers in some cichlid fish lineages in East African Great Lakes: the halochromine cichlid fishes in Lakes Victoria and Malawi. I show that speciation rates in this group are higher than in any other lacustrine fish radiation. Against this background, I review hypotheses put forward to explain diversity in cichlid species flocks. The evolution of species diversity requires three processes: speciation, ecological radiation and anatomical diversification, and it is wrong to consider hypotheses that are relevant to different processes as alternatives to each other. The African cichlid species flocks show unusually high ecological species packing in several phylogenetic groups and unusually high speciation rates in haplochromines. Therefore, it maybe concluded that at least two evolutionary models are required to explain the difference between cichlid diversity and other fish diversity in East African Lakes: one for speciation in haplochromines and one for coexistence. Subsequently I review work on speciation in haplochromines, and in particular studies aimed at testing the hypothesis of speciation by sexual selection. Haplochromines have a polygynous mating system, conducive to sexual selection, but other polygynous cichlids are not particularly species rich. This suggests that more than just strong sexual selection is required to explain haplochromine species richness. Recent palaeoecological evidence undermines the previously popular hypotheses that explained the species richness of Lake Victoria in terms of speciation under varying natural or sexual selection regimes in satellite lakes or in isolated lake basins. I summarize experimental and comparative studies, which provide evidence for two mechanisms of sympatric speciation by disruptive sexual selection on polymorphic coloration. Such modes of speciation may explain (i) the high speciation rates in colour polymorphic lineages of haplochromine cichlids under conditions where colour variation is visible in clear water, and (ii) in combination with factors that affect population survival, the unusual species richness in haplochromine species flocks. I argue that sexual selection, if disruptive, can accelerate the pace of adaptive radiation because the resultant genetic population fragmentation allows a much increased rate of differential response to disruptive natural selection. Hence, the ecological pattern of diversity resembles that produced by disruptive natural selection, with the difference that disruptive sexual selection continues to cause (gross) speciation even after niche space is saturated. This may explain the unusually high numbers of very closely related and ecologically similar species in haplochromine species flocks. The role of disruptive sexual selection is twofold: it not only causes speciation, but also maintains reproductive isolation in sympatry between species that have evolved in sympatry or allopatry. Therefore, the maintenance of diversity in species flocks that originated through sexual selection depends on the persistence of the selection regime within the environmental signal space under which that diversity evolved.

Local adaptations in populations of rock-dwelling haplochromines (Pisces:Cichlidae) from southern Lake Victoria

The radiation of cichlid fishes in the African great lakes is often described as adaptive, because, at a superficial level, cichlid fishes seem adapted to the ecological niches they occupy. However, adaptiveness has rarely been studied. We investigated to what extent island populations of three species of the rock-dwelling genus Neochromis, endemic to Lake Victoria, are adapted anatomically to exploit locally abundant resources. Specifically, we asked whether different resource environments were reflected in differences in the feeding apparatus, both within species and between species. In populations of two specialized biters, the algae scrapers N. rufocaudalis and N. omnicaeruleus, the biting force of the lower jaw increased with increasing amount of items that require biting in the diet. N. greenwoodi is a less specialized biter; we found differences between two populations in the hyoid position and in the premaxilla that enhance suction feeding. These adaptations were related to the amount of items requiring suction. Comparing across three sympatric pairs of species, in each case different diets were reflected in differences in anatomy.

Can ancient colour polymorphisms explain why some cichlid lineages speciate rapidly under disruptive sexual selection?

It is not sufficiently understood why some lineages of cichlid fishes have proliferated in the Great Lakes of East Africa much more than anywhere else in the world, and much faster than other cichlid lineages or any other group of freshwater fish. Recent field and experimental work on Lake Victoria haplochromines suggests that mate choice-mediated disruptive sexual selection on coloration, that can cause speciation even in the absence of geographical isolation, may explain it. We summarize the evidence and propose a hypothesis for the genetics of coloration that may help understand the phenomenon. By detl ning colour patterns by hue and arrangement of hues on the body, we could assign almost all observed phenotypes of Lake Victoria cichlids to one of three female («plain», «orange blotched», «black and white») and three male («blue», «red-ventrum», «reddorsum») colour patterns. These patterns diagnose species but frequently eo-occur also as morphs within the same population, where they are associated with variation in mate preferences, and appear to be transient stages in speciation. Particularly the male patterns occur in almost every genus of the species flock. We propose that the patterns and their association into polymorphisms express an ancestral trait that is retained across speciation. Our model for male colour pattern assumes two structural loci. When both are switched off, the body is blue. When switched on by a cascade of polymorphic regulatory genes, one expresses a yellow to red ventrum, the other one a yellow to red dorsum. The expression of colour variation initiates speciation. The blue daughter species will inherit the variation at the regulatory genes that can, without new mutational events, purely by recombination, again expose the colour polymorphism, starting the process anew. Very similar colour patterns also dominate among the Mbuna of Lake Malawi. In contrast, similar colour polymorphisms do not exist in the lineages that have not proliferated in the Great Lakes. The colour pattern polymorphism may be an ancient trait in the lineage (or lineages) that gave rise to the two large haplochromine radiations. We propose two tests of our hypothesis.

Lake-effect snow as the dominant control of mesic-forest distribution in Michigan, USA

1Recent studies demonstrated the sensitivity of northern forest ecosystems to changes in the amount and duration of snow cover at annual to decadal time scales. However, the consequences of snowfall variability remain uncertain for ecological variables operating at longer time scales, especially the distributions of forest communities. 2The Great Lakes region of North America offers a unique setting to examine the long-term effects of variable snowfall on forest communities. Lake-effect snow produces a three-fold gradient in annual snowfall over tens of kilometres, and dramatic edaphic variations occur among landform types resulting from Quaternary glaciations. We tested the hypothesis that these factors interact to control the distributions of mesic (dominated by Acer saccharum, Tsuga canadensis and Fagus grandifolia) and xeric forests (dominated by Pinus and Quercus spp.) in northern Lower Michigan. 3We compiled pre-European-settlement vegetation data and overlaid these data with records of climate, water balance and soil, onto Landtype Association polygons in a geographical information system. We then used multivariate adaptive regression splines to model the abundance of mesic vegetation in relation to environmental controls. 4Snowfall is the most predictive among five variables retained by our model, and it affects model performance 29% more than soil texture, the second most important variable. The abundance of mesic trees is high on fine-textured soils regardless of snowfall, but it increases with snowfall on coarse-textured substrates. Lake-effect snowfall also determines the species composition within mesic forests. The weighted importance of A. saccharum is significantly greater than of T. canadensis or F. grandifolia within the lake-effect snowbelt, whereas T. canadensis is more plentiful outside the snowbelt. These patterns are probably driven by the influence of snowfall on soil moisture, nutrient availability and fire return intervals. 5Our results imply that a key factor dictating the spatio-temporal patterns of forest communities in the vast region around the Great Lakes is how the lake-effect snowfall regime responds to global change. Snowfall reductions will probably cause a major decrease in the abundance of ecologically and economically important species, such as A. saccharum.

Are great Cascadia earthquakes recorded in the sedimentary records from small forearc lakes?

Here we investigate sedimentary records from four small inland lakes located in the southern Cascadia forearc region for evidence of earthquakes. Three of these lakes are in the Klamath Mountains near the Oregon–California border, and one is in the central Oregon Coast range. The sedimentary sequences recovered from these lakes are composed of normal lake sediment interbedded with disturbance event layers. The thickest of these layers are graded, and appear to be turbidites or linked debrites (turbidites with a basal debris-flow deposit), suggesting rapid deposition. Variations in particle size and organic content of these layers are reflected in the density and magnetic susceptibility data. The frequency and timing of these events, based on radiocarbon ages from detrital organics, is similar to the offshore seismogenic turbidite record from trench and slope basin cores along the Cascadia margin. Stratigraphic correlation of these anomalous deposits based on radiocarbon ages, down-core density, and magnetic susceptibility data between lake and offshore records suggest synchronous triggering. The areal extent and multiple depositional environments over which these events appear to correlate suggest that these deposits were most likely caused by shaking during great Cascadia earthquakes.

Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on 210Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for 226Ra, 210Pb and 137Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of 210Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the 137Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the 137Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments