Modeling contaminant behavior in Lake Superior : a comparison of PCBs, PBDEs, and mercury

A mass‐balance model for Lake Superior was applied to polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), and mercury to determine the major routes of entry and the major mechanisms of loss from this ecosystem as well as the time required for each contaminant class to approach steady state. A two‐box model (water column, surface sediments) incorporating seasonally adjusted environmental parameters was used. Both numerical (forward Euler) and analytical solutions were employed and compared. For validation, the model was compared with current and historical concentrations and fluxes in the lake and sediments. Results for PCBs were similar to prior work showing that air‐water exchange is the most rapid input and loss process. The model indicates that mercury behaves similarly to a moderately‐chlorinated PCB, with air‐water exchange being a relatively rapid input and loss process. Modeled accumulation fluxes of PBDEs in sediments agreed with measured values reported in the literature. Wet deposition rates were about three times greater than dry particulate deposition rates for PBDEs. Gas deposition was an important process for tri‐ and tetra‐BDEs (BDEs 28 and 47), but not for higher‐brominated BDEs. Sediment burial was the dominant loss mechanism for most of the PBDE congeners while volatilization was still significant for tri‐ and tetra‐BDEs. Because volatilization is a relatively rapid loss process for both mercury and the most abundant PCBs (tri‐ through penta‐), the model predicts that similar times (from 2 ‐ 10 yr) are required for the compounds to approach steady state in the lake. The model predicts that if inputs of Hg(II) to the lake decrease in the future then concentrations of mercury in the lake will decrease at a rate similar to the historical decline in PCB concentrations following the ban on production and most uses in the U.S. In contrast, PBDEs are likely to respond more slowly if atmospheric concentrations are reduced in the future because loss by volatilization is a much slower process for PBDEs, leading to lesser overall loss rates for PBDEs in comparison to PCBs and mercury. Uncertainties in the chemical degradation rates and partitioning constants of PBDEs are the largest source of uncertainty in the modeled times to steady‐state for this class of chemicals. The modeled organic PBT loading rates are sensitive to uncertainties in scavenging efficiencies by rain and snow, dry deposition velocity, watershed runoff concentrations, and uncertainties in air‐water exchange such as the effect of atmospheric stability.

An investigation of phenolic glycoside and condensed tannin homeostasis in Populus by salicyl alcohol feeding to cell cultures and by transgenic manipulation of the sucrose transporter, PTSUT4, in planta

Secondary metabolites play an important role in plant protection against biotic and abiotic stress. In Populus, phenolic glycosides (PGs) and condensed tannins (CTs) are two such groups of compounds derived from the common phenylpropanoid pathway. The basal levels and the inducibility of PGs and CTs depend on genetic as well as environmental factors, such as soil nitrogen (N) level. Carbohydrate allocation, transport and sink strength also affect PG and CT levels. A negative correlation between the levels of PGs and CTs was observed in several studies. However, the molecular mechanism underlying such relation is not known. We used a cell culture system to understand negative correlation of PGs and CTs. Under normal culture conditions, neither salicin nor higher-order PGs accumulated in cell cultures. Several factors, such as hormones, light, organelles and precursors were discussed in the context of aspen suspension cells’ inability to synthesize PGs. Salicin and its isomer, isosalicin, were detected in cell cultures fed with salicyl alcohol, salicylaldehyde and helicin. At higher levels (5 mM) of salicyl alcohol feeding, accumulation of salicins led to reduced CT production in the cells. Based on metabolic and gene expression data, the CT reduction in salicin-accumulating cells is partly a result of regulatory changes at the transcriptional level affecting carbon partitioning between growth processes, and phenylpropanoid CT biosynthesis. Based on molecular studies, the glycosyltransferases, GT1-2 and GT1-246, may function in glycosylation of simple phenolics, such as salicyl alcohol in cell cultures. The uptake of such glycosides into vacuole may be mediated to some extent by tonoplast localized multidrug-resistance associated protein transporters, PtMRP1 and PtMRP6. In Populus, sucrose is the common transported carbohydrate and its transport is possibly regulated by sucrose transporters (SUTs). SUTs are also capable of transporting simple PGs, such as salicin. Therefore, we characterized the SUT gene family in Populus and investigated, by transgenic analysis, the possible role of the most abundantly expressed member, PtSUT4, in PG-CT homeostasis using plants grown under varying nitrogen regimes. PtSUT4 transgenic plants were phenotypically similar to the wildtype plants except that the leaf area-to-stem volume ratio was higher for transgenic plants. In SUT4 transgenics, levels of non-structural carbohydrates, such as sucrose and starch, were altered in mature leaves. The levels of PGs and CTs were lower in green tissues of transgenic plants under N-replete, but were higher under N-depleted conditions, compared to the levels in wildtype plants. Based on our results, SUT4 partly regulates N-level dependent PG-CT homeostasis by differential carbohydrate allocation.

DEVELOPMENT OF A CONCEPTUAL MODEL AND PARAMETERIZATION OF A MATHEMATICAL MODEL FOR SEDIMENT PHOSPHORUS DIAGENESIS, SORPTION AND RELEASE FROM ONONDAGA LAKE SYRACUSE, NEW YORK

Eutrophication is a persistent problem in many fresh water lakes. Delay in lake recovery following reductions in external loading of phosphorus, the limiting nutrient in fresh water ecosystems, is often observed. Models have been created to assist with lake remediation efforts, however, the application of management tools to sediment diagenesis is often neglected due to conceptual and mathematical complexity. SED2K (Chapra et al. 2012) is proposed as a "middle way", offering engineering rigor while being accessible to users. An objective of this research is to further support the development and application SED2K for sediment phosphorus diagenesis and release to the water column of Onondaga Lake. Application of SED2K has been made to eutrophic Lake Alice in Minnesota. The more homogenous sediment characteristics of Lake Alice, compared with the industrially polluted sediment layers of Onondaga Lake, allowed for an invariant rate coefficient to be applied to describe first order decay kinetics of phosphorus. When a similar approach was attempted on Onondaga Lake an invariant rate coefficient failed to simulate the sediment phosphorus profile. Therefore, labile P was accounted for by progressive preservation after burial and a rate coefficient which gradual decreased with depth was applied. In this study, profile sediment samples were chemically extracted into five operationally-defined fractions: CaCO3-P, Fe/Al-P, Biogenic-P, Ca Mineral-P and Residual-P. Chemical fractionation data, from this study, showed that preservation is not the only mechanism by which phosphorus may be maintained in a non-reactive state in the profile. Sorption has been shown to contribute substantially to P burial within the profile. A new kinetic approach involving partitioning of P into process based fractions is applied here. Results from this approach indicate that labile P (Ca Mineral and Organic P) is contributing to internal P loading to Onondaga Lake, through diagenesis and diffusion to the water column, while the sorbed P fraction (Fe/Al-P and CaCO3-P) is remaining consistent. Sediment profile concentrations of labile and total phosphorus at time of deposition were also modeled and compared with current labile and total phosphorus, to quantify the extent to which remaining phosphorus which will continue to contribute to internal P loading and influence the trophic status of Onondaga Lake. Results presented here also allowed for estimation of the depth of the active sediment layer and the attendant response time as well as the sediment burden of labile P and associated efflux.

Spatial isolation and genetic differentiation in naturally fragmented plant populations of the Swiss Alps

Aims The effect Of anthropogenic landscape fragmentation on the genetic diversity and adaptive potential of plant populations is a major issue in conservation biology. However, little is known about the partitioning of genetic diversity in alpine species, which occur in naturally fragmented habitats. Here, we, investigate molecular patterns of three alpine plants (Epilobium fleischeri, Geum reptans and Campanula thyrsoides) across Switzerland and ask whether Spatial isolation has led to high levels of populations differentiation, increasing over distance, and a decrease of within-population variability. We further hypothesize that file contrasting potential for long-distance dispersal (LDD) of Seed in these Species will considerably influence and explain diversity partitioning. Methods For each study species, we Sampled 20-23 individuals from each of 20-32 populations across entire Switzerland. We applied Random Amplified Polymorphic Dimorphism markers to assess genetic diversity within (Nei's expected heterozygosity, H-e; percentage of polymorphic hands, P-P) and among (analysis of molecular variance, Phi(st)) populations and correlated population size and altitude with within-populalion diversity. Spatial patterns of genetic relatedness were investigated using Mantel tests and standardized major axis regression as well as unweighted pair group method with arithmetic mean cluster analyses and Monmonier's algorithm. To avoid known biases, We standardized the numbers of populations, individuals and markers using multiple random reductions. We modelled LDD with a high alpine wind data set using the terminal velocity and height of seed release as key parameters. Additionally, we assessed a number of important life-history traits and factors that potentially influence genetic diversity partitioning (e.g. breeding system, longevity and population size). Important findings For all three species, We found a significant isolation-by-distance relationship but only a moderately high differentiation among populations (Phi(st): 22.7, 48 and 16.8%, for E. fleischeri, G. reptans and C. thyrsoides, respectively). Within-population diversity (H-c: 0.19-0.21, P-p: 62-75%) was not reduced in comparison to known results from lowland species and even small populations with < 50 reproductive individuals contained high levels of genetic diversity. We further found no indication that a high long-distance seed dispersal potential enhances genetic connectivity among populations. Gene flow seems to have a strong stochastic component causing large dissimilarity between population pairs irrespective of the spatial distance. Our results suggest that other life-history traits, especially the breeding System, may play an important role in genetic diversity partitioning. We conclude that spatial isolation in the alpine environment has a strong influence on population relatedness but that a number of factors can considerably influence the strength of this relationship.

The role of landuse and natural determinants for grassland vegetation composition in the Swiss Alps

The Alps provide a high habitat diversity for plant species, structured by broad- and fine-scale abiotic site conditions. In man-made grasslands, vegetation composition is additionally affected by the type of landuse. We recorded vegetation composition in 216 parcels of grassland in 12 municipalities representing an area of 170 x 70 km in the south-eastern part of the Swiss Alps. Each parcel was characterized by a combination of altitudinal level (valley, intermediate, alp). traditional landuse (mown. grazed), current management (mown, grazed, abandoned). and Fertilization (unfertilized, fertilized). For each parcel we also assessed the abiotic factors aspect, slope, pH value, and geographic coordinates, and for each municipality annual precipitation and its cultural tradition. We analysed vegetation composition using (i) variation partitioning in RDA. (ii) cover of graminoids. non-legume forbs, and legumes, and (iii) dominance and frequency of species. Species composition was determined by, in decreasing order of variation explained. landuse, broad-scale abiotic factors, fine-scale abiotic factors. and cultural tradition. Current socio-economically motivated landuse changes, such as grazing of unfertilized former meadows or their abandonment, strongly affect vegetation composition. In our study, the frequency of characteristic meadow species was significantly smaller in grazed and even smaller in abandoned parcels than in still mown ones, suggesting less severe consequences of grazing for vegetation composition than of abandonment. Therefore. low-intensity grazing and mowing every few years should be considered valuable conservation alternatives to abandonment. Furthermore. because each landuse type was characterized by different species. a high variety of landuse types should be promoted to preserve plant species diversity in Alpine grasslands. (C) 2007 Gesellschaft fur Okologie. Published by Elsevier GmbH. All rights reserved.

Using P-GENESIS for Parallel Simulation of GENESIS Models

P-GENESIS is an extension to the GENESIS neural simulator that allows users to take advantage of parallel machines to speed up the simulation of their network models or concurrently simulate multiple models. P-GENESIS adds several commands to the GENESIS script language that let a script running on one processor execute remote procedure calls on other processors, and that let a script synchronize its execution with the scripts running on other processors. We present here some brief comments on the mechanisms underlying parallel script execution. We also offer advice on parallelizing parameter searches, partitioning network models, and selecting suitable parallel hardware on which to run P-GENESIS.

Evidence for past variations in methane availability in a Siberian thermokarst lake based on δ13C of chitinous invertebrate remains

Understanding past methane dynamics in arctic wetlands and lakes is crucial for estimating future methane release. Methane fluxes from lake ecosystems have increasingly been studied, yet only few reconstructions of past methane emissions from lakes are available. In this study, we develop an approach to assess changes in methane availability in lakes based on δ13C of chitinous invertebrate remains and apply this to a sediment record from a Siberian thermokarst lake. Diffusive methane fluxes from the surface of ten newly sampled Siberian lakes and seven previously studied Swedish lakes were compared to taxon-specific δ13C values of invertebrate remains from lake surface sediments to investigate whether these invertebrates assimilated 13C-depleted carbon typical for methane. Remains of chironomid larvae of the tribe Orthocladiinae that, in the study lakes, mainly assimilate plant-derived carbon had higher δ13C than other invertebrate groups. δ13C of other invertebrates such as several chironomid groups (Chironomus, Chironomini, Tanytarsini, and Tanypodinae), cladocerans (Daphnia), and ostracods were generally lower. δ13C of Chironomini and Daphnia, and to a lesser extent Tanytarsini was variable in the lakes and lower at sites with higher diffusive methane fluxes. δ13C of Chironomini, Tanytarsini, and Daphnia were correlated significantly with diffusive methane flux in the combined Siberian and Swedish dataset (r = −0.72, p = 0.001, r = −0.53, p = 0.03, and r = −0.81, p < 0.001, respectively), suggesting that δ13C in these invertebrates was affected by methane availability. In a second step, we measured δ13C of invertebrate remains from a sediment record of Lake S1, a shallow thermokarst lake in northeast Siberia. In this record, covering the past ca 1000 years, δ13C of taxa most sensitive to methane availability (Chironomini, Tanytarsini, and Daphnia) was lowest in sediments deposited from ca AD 1250 to ca AD 1500, and after AD 1970, coinciding with warmer climate as indicated by an independent local temperature record. As a consequence the offset in δ13C between methane-sensitive taxa and bulk organic matter was higher in these sections than in other parts of the core. In contrast, δ13C of other invertebrate taxa did not show this trend. Our results suggest higher methane availability in the study lake during warmer periods and that thermokarst lakes can respond dynamically in their methane output to changing environmental conditions.

Isotopes Trace Biogeochemistry and Sources of Cu and Zn in an intertidal soil

River floodplain soils are sinks and potential sources for toxic trace metals like Cu and Zn. We hypothesize that stable Cu and Zn isotope ratios reflect both the mobilization and the sources of metals. We determined the soil properties, the concentrations and partitioning of Cu and Zn, and variations in δ65Cu and δ66Zn values in a core obtained from an Aquic Udifluvent developed on a freshwater intertidal mudflat of the River Elbe, Germany. The core was sampled at 2 cm intervals to a depth of 34 cm, which corresponds to approximately 9 yr of sedimentation. Elevated concentrations of Cu (up to 320 μg g−1) and Zn (up to 2080 μg g−1) indicated anthropogenic pollution. At the time of sampling the redox conditions changed from oxic (Eh 200 to 400 mV, above 22 cm deep) to strongly anoxic conditions (-100 to -200 mV, below 22 cm deep). The δ65Cu values varied systematically with depth (from -0.02 to 0.16‰) and were correlated with the Fe, C, and N concentrations. Although pre-depositional variations cannot be ruled out, the systematic variation with depth suggests post-sedimentation fractionation of δ65Cu in response to seasonally variable organic matter deposition and redox conditions. In contrast, the δ66ZnIRMM values were uniform (from -0.07 to 0.01‰) throughout the core, indicating that the Zn isotopes did not significantly fractionate after deposition and that the Zn sources were homogeneous throughout the sedimentation.

Habitat selection of three cryptic Plecotus bat species in the European Alps reveals contrasting implications for conservation

Assessing the ecological requirements of species coexisting within a community is an essential requisite for developing sound conservation action. A particularly interesting question is what mechanisms govern the stable coexistence of cryptic species within a community, i.e. species that are almost impossible to distinguish. Resource partitioning theory predicts that cryptic species, like other sympatric taxa, will occupy distinct ecological niches. This prediction is widely inferred from eco-morphological studies. A new cryptic long-eared bat species, Plecotus macrobullaris, has been recently discovered in the complex of two other species present in the European Alps, with even evidence for a few mixed colonies. This discovery poses challenges to bat ecologists concerned with planning conservation measures beyond roost protection. We therefore tested whether foraging habitat segregation occurred among the three cryptic Plecotus bat species in Switzerland by radiotracking 24 breeding female bats (8 of each species). We compared habitat features at locations visited by a bat versus random locations within individual home ranges, applying mixed effects logistic regression. Distinct, species-specific habitat preferences were revealed. P. auritus foraged mostly within traditional orchards in roost vicinity, with a marked preference for habitat heterogeneity. P. austriacus foraged up to 4.7 km from the roost, selecting mostly fruit tree plantations, hedges and tree lines. P. macrobullaris preferred patchy deciduous and mixed forests with high vertical heterogeneity in a grassland dominated-matrix. These species-specific habitat preferences should inform future conservation programmes. They highlight the possible need of distinct conservation measures for species that look very much alike.

Designing forest biodiversity experiments: general considerations illustrated by a new large experiment in subtropical China

1. Biodiversity-ecosystem functioning (BEF) experiments address ecosystem-level consequences of species loss by comparing communities of high species richness with communities from which species have been gradually eliminated. BEF experiments originally started with microcosms in the laboratory and with grassland ecosystems. A new frontier in experimental BEF research is manipulating tree diversity in forest ecosystems, compelling researchers to think big and comprehensively. 2. We present and discuss some of the major issues to be considered in the design of BEF experiments with trees and illustrate these with a new forest biodiversity experiment established in subtropical China (Xingangshan, Jiangxi Province) in 2009/2010. Using a pool of 40 tree species, extinction scenarios were simulated with tree richness levels of 1, 2, 4, 8 and 16 species on a total of 566 plots of 25.8x25.8m each. 3. The goal of this experiment is to estimate effects of tree and shrub species richness on carbon storage and soil erosion; therefore, the experiment was established on sloped terrain. The following important design choices were made: (i) establishing many small rather than fewer larger plots, (ii) using high planting density and random mixing of species rather than lower planting density and patchwise mixing of species, (iii) establishing a map of the initial ecoscape' to characterize site heterogeneity before the onset of biodiversity effects and (iv) manipulating tree species richness not only in random but also in trait-oriented extinction scenarios. 4. Data management and analysis are particularly challenging in BEF experiments with their hierarchical designs nesting individuals within-species populations within plots within-species compositions. Statistical analysis best proceeds by partitioning these random terms into fixed-term contrasts, for example, species composition into contrasts for species richness and the presence of particular functional groups, which can then be tested against the remaining random variation among compositions. 5. We conclude that forest BEF experiments provide exciting and timely research options. They especially require careful thinking to allow multiple disciplines to measure and analyse data jointly and effectively. Achieving specific research goals and synergy with previous experiments involves trade-offs between different designs and requires manifold design decisions.

Comparison of inorganic nitrogen uptake dynamics following snowmelt and at peak biomass in subalpine grasslands

Subalpine grasslands are highly seasonal environments and likely subject to strong variability in nitrogen (N) dynamics. Plants and microbes typically compete for N acquisition during the growing season and particularly at plant peak biomass. During snowmelt, plants could potentially benefit from a decrease in competition by microbes, leading to greater plant N uptake associated with active growth and freeze-thaw cycles restricting microbial growth. In managed subalpine grasslands, we expect these interactions to be influenced by recent changes in agricultural land use, and associated modifications in plant and microbial communities. At several subalpine grasslands in the French Alps, we added pulses of 15N to the soil at the end of snowmelt, allowing us to compare the dynamics of inorganic N uptake in plants and microbes during this period with that previously reported at the peak biomass in July. In all grasslands, while specific shoot N translocation (per g of biomass) of dissolved inorganic nitrogen (DIN) was two to five times greater at snowmelt than at peak biomass, specific microbial DIN uptakes were similar between the two sampling dates. On an area basis, plant communities took more DIN than microbial communities at the end of snowmelt when aboveground plant biomasses were at least two times lower than at peak biomass. Consequently, inorganic N partitioning after snowmelt switches in favor of plant communities, allowing them to support their growing capacities at this period of the year. Seasonal differences in microbial and plant inorganic N-related dynamics were also affected by past (terraced vs. unterraced) rather than current (mown vs. unmown) land use. In terraced grasslands, microbial biomass N remained similar across seasons, whereas in unterraced grasslands, microbial biomass N was higher and microbial C : N lower at the end of snowmelt as compared to peak biomass. Further investigations on microbial community composition and their organic N uptake dynamics are required to better understand the decrease in microbial DIN uptake.

The Wabar impact craters, Saudi Arabia, revisited

The very young Wabar craters formed by impact of an iron meteorite and are known to the scientific community since 1933. We describe field observations made during a visit to the Wabar impact site, provide analytical data on the material collected, and combine these data with poorly known information discovered during the recovery of the largest meteorites. During our visit in March 2008, only two craters (Philby-B and 11 m) were visible; Philby-A was completely covered by sand. Mapping of the ejecta field showed that the outcrops are strongly changing over time. Combining information from different visitors with our own and satellite images, we estimate that the large seif dunes over the impact site migrate by approximately 1.0–2.0 m yr␣1 southward. Shock lithification took place even at the smallest, 11 m crater, but planar fractures (PFs) and undecorated planar deformation features (PDFs), as well as coesite and stishovite, have only been found in shock-lithified material from the two larger craters. Shock-lithified dune sand material shows perfectly preserved sedimentary structures including cross-bedding and animal burrows as well as postimpact structures such as open fractures perpendicular to the bedding, slickensides, and radiating striation resembling shatter cones. The composition of all impact melt glasses can be explained as mixtures of aeolian sand and iron meteorite. We observed a partial decoupling of Fe and Ni in the black impact glass, probably due to partitioning of Ni into unoxidized metal droplets. The absence of a Ca-enriched component demonstrates that the craters did not penetrate the bedrock below the sand sheet, which has an estimated thickness of 20–30 m.

Electron density building block approach for metal organic frameworks

A general introduction to the state of the art in modeling metal organic materials using transferable atomic multipoles is provided. The method is based on the building block partitioning of the electron density, which is illustrated with some examples of potential applications and with detailed discussions of the advantages and pitfalls. The interactions taking place between building blocks are summarized and are used to discuss the properties that can be calculated.

Diversity versus disparity and the role of ecological opportunity in a continental bird radiation

Aim The Neotropical parrots (Arini) are an unusually diverse group which colonized South America in the Oligocene. The newly invaded Neotropics may have functioned as an underused adaptive zone and provided novel ecological opportunities that facilitated diversification. Alternatively, diversification may have been driven by ecological changes caused by Andean uplift and/or climate change from the Miocene onwards. Our aim was to find out whether Arini diversified in a classical adaptive radiation after their colonization of South America, or whether their diversification occurred later and was influenced by more recent environmental change. Location Neotropics. Methods We generated a time-calibrated phylogeny of more than 80% of all Arini species in order to analyse lineage diversification. This chronogram was also used as the basis for the reconstruction of morphological evolution within Arini using a multivariate ratio analysis of three size measurements. Results We found a concentration of size evolution and partitioning of size niches in the early history of Arini consistent with the process of adaptive radia- tion, but there were no signs of an early burst of speciation or a decrease in speci- ation rates through time. Although we detected no overall temporal shifts in diversification rates, we discovered two young, unexpectedly species-rich clades. Main conclusions Arini show signs of an early adaptive radiation, but we found no evidence of the slowdown in speciation rate generally considered a feature of island or lake radiations. Historical processes and environmental change from the Miocene onwards may have kept diversification rates roughly constant ever since the colonization of the Neotropics. Thus, Arini may not yet have reached equilibrium diversity. The lack of diversity-dependent speciation might be a general feature of adaptive radiations on a continental scale, and diversification processes on continents might therefore not be as ecologically limited as in isolated lakes or on oceanic islands.

From transient to steady state deformation and grain size: A thermodynamic approach using elasto-visco-plastic numerical modeling

Numerical simulation experiments give insight into the evolving energy partitioning during high-strain torsion experiments of calcite. Our numerical experiments are designed to derive a generic macroscopic grain size sensitive flow law capable of describing the full evolution from the transient regime to steady state. The transient regime is crucial for understanding the importance of micro structural processes that may lead to strain localization phenomena in deforming materials. This is particularly important in geological and geodynamic applications where the phenomenon of strain localization happens outside the time frame that can be observed under controlled laboratory conditions. Ourmethod is based on an extension of the paleowattmeter approach to the transient regime. We add an empirical hardening law using the Ramberg-Osgood approximation and assess the experiments by an evolution test function of stored over dissipated energy (lambda factor). Parameter studies of, strain hardening, dislocation creep parameter, strain rates, temperature, and lambda factor as well asmesh sensitivity are presented to explore the sensitivity of the newly derived transient/steady state flow law. Our analysis can be seen as one of the first steps in a hybrid computational-laboratory-field modeling workflow. The analysis could be improved through independent verifications by thermographic analysis in physical laboratory experiments to independently assess lambda factor evolution under laboratory conditions.