Friday, October 14, 2016

Diatoms on slides: They're everywhere and can tell you so much!

Diatoms are one of the most abundant biological microfossils in sediment cores and on smear slides. Diatoms are one group of microscopic algae and characterized by their golden-brown pigmentation and having a cell wall made of biologically produced glass or biogenic opaline silica.  Because they are made of opaline silica, diatom cell walls are fairly resistant to breakage and dissolution and can form up to 60% of the dry weight in some lake cores (Edlund et al. 2000). The cell wall in a living diatom is put together like a Petri dish; the two halves are called valves with one valve slightly larger and overlapping the other valve. There are a few other parts in the cell wall, primarily hoop or band-like structures called girdle bands that form the overlap between the two valves. Collectively, the siliceous parts surrounding a single cell are termed the frustule. Size, shape, and ornamentation on the valves distinguish the more than 40,000 diatom species in the world. For detailed diatom analysis, sediments are specially prepared using peroxide or strong acids to digest the cell contents; the "cleaned" cell walls are mounted on microscope slides using specialized mountants that optimize resolution of valve ornamentation for species-level identification.  More detailed information on diatom preparation and species identification can be found here.


But back to smear slides…you'll find diatoms in your smear slides as single valves, single cells, or many species form characteristic colonies of filaments, bands, zig-zags, or star shapes. Diatoms may be empty, filled with an air bubble, or can even have cell contents remaining (many species of diatoms get buried alive in the sediments as resting cells and can be rejuvenated to grow again even after decades of burial; Sicko-Goad et al. 1989. Scan around your slide at 200-400x (20x or 40x objective) to see how abundant diatoms are in your material. Don't be overwhelmed by their many shapes and forms. Are there lots of different forms or just a few? What kind of shapes do you see: circles, barrel shapes lined up in filaments, banded colonies that look like fences, big ones and little ones, or just lots of funky shapes. There are expert diatomists out there that can help if you want to do a detailed analysis of the diatoms, but learning just a little bit about the diatom flora in your core using the below tutorials can quickly and easily provide some of your first and best paleoenvironmental information. Let's take a look at some simple qualitative and quantitative measures that can be made from smear slides related to abundance, preservation, and diatom ecology.

Diatom abundance
The quantity of diatoms in sediments varies considerably due to the lake's paleoecology, depositional environment, and preservation.  Noting the quantity of diatoms in a smear slide can help focus future analyses, show general shifts in lake productivity (e.g., eutrophication events), and identify changes in the depositional environment (e.g., turbidites). If you have carefully made your smear slides using similar amounts of material for each slide, the relative abundance of diatoms can be determined by simple comparison. Look at a few of your slides under 200-400x magnification. The abundance of diatoms present in a smear slide can be reported using a simple relative scale as abundant, common, few, or absent.



Diatom preservation
The condition of diatom microfossils in a smear slide is a reflection of the depositional environment and will determine their suitability and usefulness as environmental indicators. Several environmental conditions enhance the dissolution of diatoms: breakage, salinity, alkalinity, and silica-poor environments.  Diatom valves are composed of biogenic opaline silica and each valve has areas where silica is deposited more thickly (e.g., central areas, raphe ribs) as well as areas with patterns of ornamentation comprising holes (areolae, pore fields) and lines (striae).  The dissolution of diatoms generally follows a pattern where highly ornamented or thin areas are lost first, followed by sequential loss depending on silica thickness, until only the most heavily silicified areas of the diatom remain. Diatom species are also differentially susceptible to dissolution.  A simple metric (% pristine) based on the condition of diatom microfossils in smear slides can be used to summarize the level of dissolution in your samples.
Examples:
Cyclotella species

Epithemia species
Fragilaria species

Ryves et al. (2001, 2009) discussed issues surrounding diatom dissolution in sediment cores, proposed several indices to describe dissolution artefacts, and offered insights into how dissolution might affect interpretation of paleorecords.

Diatom ecology – plankton:benthic ratio
Diatoms can be ecologically separated in two broad categories: planktonic species that live in the water column and benthic species that live in the sediments or attached to rock, plant, or other substrates.  A ratio called P:B or the planktonic:benthic ratio provides a simple metric of possible changes in the paleoenviroment based on lake level change or nutrients.  The ratio is often expressed as a percentage, i.e. a P:B ratio of 1.0 would equal 50% percent planktonic species. Depending on lake morphometry lake level changes can result in dramatic shifts in the proportion of open water and littoral zone habitat (Wolin and Stone 2010).  Sediments preserve a temporal and spatial integration of a lake's diatom communities so periods with greater littoral zone extent will generally have great benthic diatom numbers. The other major factor that drives changes in P:B is lake eutrophication. Nutrient addition to lakes often results in proportionally greater planktonic diatom growth and higher P:B ratios (Edlund et al. 2009).
Calculating the P:B ratio for your smear slides requires two steps. First let's learn how to recognize planktonic vs benthic diatom, and try out our skills. Then we'll discuss how many fields might be appropriate to count to get a representative value for your smear slide.
Diatoms come in all shapes and sizes. Planktonic species generally have five common shapes: circular, barrel-shaped in filaments, band shaped colonies, star or zig-zag colonies, or long thin species. Here's a montage of many planktonic forms:

Benthic diatoms also vary in shape, but few will be round, rather most have a shape that is symmetric along a line and more similar to the outline of a boat (oval or elliptic), or shaped like the moon (lunate) or a club, and all may have variously shaped ends. Here's a montage of many different benthic species:

Let's look at some pictures from smear slides and see how easy it is to calculate a P:B ratio. This is an image from a eutrophic lake in Minnesota. Can you refer to the montage of planktonic and benthic diatoms and then count how many diatom valves you see from each category?  Remember that each cell of a diatom has two valves—sometimes they remain together in a cell, other times they get dissociated. In real life, you may need to focus up and down a bit to figure out if both valves are there (Hint: an air bubble in the diatom usually means there are two valves). Here's an easy one:



Did you get 1 planktonic and 4 benthic for a 1:4 or 0.25 P:B or 20% planktonic species? Note one of the diatoms has both valves and it's laying on its side, and we usually don't count a broken valve with less than ½ a diatom.  



Heres another image, a little tougher:

Did you find all 16 diatom valves on that last one. Note that when there is an air bubble in the diatom, usually both valves are present.  
Now a really complicated one. There are at least 37 valves in the image.
Did you find 37 valves, about 31 planktonic and 6 benthic to get a P:B ratio of 5.2 with 83.8% planktonic? That band-shaped colony in the middle is the planktonic Fragilaria capucina and each dark line in the colony is where two valves are connected together. Note that you also have to make a decision when counting of what to do with diatoms at the edge of your field of view. Perhaps count them if >1/2 of the valve shows.


You can see that the density of diatoms and the ease of counting varies among slides. To generate good P:B data you will need to count multiple fields on a slide until your P:B number converge on a value. It is easy to just enter your data into a spreadsheet as you count each random field on your smear slide. Keep a running estimate of your P:B scores and consider your count complete once you converge on a stable value. Typically diatomists count between 400 and 600 diatoms to get these values; you may be able to count fewer.

Diatom ecology – diatom:cyst ratio
Another important and easy set of data to recover from smear slides is the ratio of diatom valves to chrysophycean cysts, or the cyst:frustule ratio, which is usally reported as the percentage of cysts to the total of cysts plus frustules. The chrysophytes and synurophytes are two groups of golden-brown microalgae that have incorporated siliceous structures into their vegetative cells (as a coat of siliceous scales and spines) and life histories. As part of their life history, many species produce siliceous resting cysts that are variously referred to as chrysophycean cysts, statocysts, statospores, or hypnospores.  Both the cell coat components (scales and spines; often analyzed in detail to determine pH changes; Zeeb and Smol 2001) and cysts can be found abundantly in sediment cores and smear slides.
The rato of cysts:diatom valves was first recognized by Smol (1985) as a way to distinguish the relative importance of each group in a lake's sediment record. The ratio was later shown to be responsive to several environmental drivers and gradients. For example, in temperate lakes lower numbers of chrysophyte cysts relative to diatoms is a classic sign of eutrophication (Smol 1985).  In arctic lakes, the ratio is responsive to past extent of ice cover, with higher values of cysts:frustule during periods of greater ice cover (reviewed in Smol 1988). The ratio also typically gets lower as lake salinity increases (Cumming et al. 1993). Other applications of the cyst:frustule ratio are discussed in Smol (1995) and Zeeb and Smol (2001).
It's easy to differentiate chysophycean cysts from diatoms. Cysts are generally heavily silicifed, spherical or ovoid in shape, with a single circular pore at one end that may or may not have a collared extension. The cyst exterior is variously ornamented and can be smooth, and with or without spines, bumps, ridges, or reticulate patterns. A working group has created a standardized terminology and method for describing cyst morphotypes (see Cronberg and Sandgren 1986, Duff et al. 1995, Wilkinson et al. 2001).
Let's look at a few images from smear slides and see how to estimate the cyst:frustule ratio in cores. Remember that diatoms have two valves in each frustule, so if you count valves you will need to divide by 2 to convert to frustules. Here's an easy one:

Did you find 4 cysts and 8 diatom valves making a total of 4 diatom frustules? Did you catch the cyst and diatom valve that were hiding behind the big diatom in the upper right? Remember to focus up and down to see those hidden ones. So a cyst:frustule ratio of 1:1 or 50% cysts.  


Let's try another one:
Lots of cysts…7 cysts and only 7 diatom valves for a total of 3.5 diatom frustules? That makes a cyst:frustule ratio of 2.0 and 66.7% cysts in this sediment assemblage.  

One more test of your skills, some tough diatoms on this one:
Only 3 cysts but 10 diatom frustules, mostly tied up in those two band-shaped chains of diatoms. They are laying on their sides (called "girdle view") and you have to focus carefully to determine how many valves are present in each chain.  So a cyst:frustule ratio of 0.33 or 23% cysts in the assemblage


As noted above for estimating the diatom P:B ratio, it is important to count enough fields under the scope to have your estimate of cysts:frustules converge on a value for one smear slide. This may take from 400 to 600 counts, just keep a running total of the ratio as you count your random fields.

Diatom Communities –
Diatom abundance and distribution are strongly controlled by a number of abiotic or environmental gradients (and biogeography-e.g., diatoms in South America are different than diatoms in North America). Among the most important gradients are pH, salinity, nutrient levels, and habitat. We earlier looked at how diatoms could be divided into their habitat preferences of planktonic vs benthic to understand lake level or nutrient changes. Here are some montages of diatoms that provide a representation of diatom communities from various lake types:
Low pH – Diatom communities from lakes with pH <7 are characterized by species that typically also prefer low conductivity and tolerate some level of DOC staining in the water. Typical acidophilic genera include Tabellaria, small-celled Aulacoseira, Asterionella and the benthic genera Eunotia, Neidium, and Brachysira.


Saline lakes – Lakes are saline based on elevated concentrations of various ion combinations of notably Na+, Ca++, Mg++, SO4--, Cl-, CO3--. The transition of diatom communities from fresh water to saline is one of the most dramatic with many fresh water diatoms intolerant of high conductivity and salinity, many saline genera more represented in true marine conditions, and a set of species that are able to tolerate a wide range of salinity. These include the genera Chaetoceros and its resting spores, saline or tolerant Cyclotella, Amphora, Halamphora, Anomoeoneis, Cocconeis, Epithemia, Surirella, Campylodiscus, Ctenophora, and Rhoicosphenia.



Oligotrophic or low nutrient lakes – Low nutrients, high light, deep water, and strong seasonality characterize most oligotrophic lakes.  The diatoms in cores typically reflect a set of planktonic species that represent communities associated with spring and fall turnover, as well as a community that forms the "deep chlorophyll layer" and is adapted to living along the thermocline in the summer months.  These species often include Cyclotella, Handmannia/Puncticulata, Stephanodiscus, Aulacoseira, Asterionella, and long-celled Tabellaria and Ulnaria species.

Eutrophic or high nutrient lakes – Many lakes suffer from excess nutrients and diatom communities shift in both composition and to greater abundance as nutrient levels increase. Diatom taxa that thrive under high nutrients include small Stephanodiscus, Cyclostephanos, some Aulacoseira, Fragilaria, Pseudostaurosira, and Staurosira species.



References cited:

Cronberg, G. & C. D. Sandgren. 1986. A proposal for the develop- ment of standardized nomenclature and terminology for chrys- ophycean statospores. In Kristiansen, J. & R. A. Andersen (eds), Chrysophytes: aspects and problems, Cambridge Univer- sity Press, Cambridge: 317–328.

Cumming, B. F., Wilson, S. E., & Smol, J. P. (1993). Paleolimnological potential of chrysophyte cysts and scales and of sponge spicules as indicators of lake salinity. International Journal of Salt Lake Research, 2(1), 87-92.

Duff, K. E., B. A. Zeeb & J. P. Smol, 1995. Atlas of Chrysophycean cysts. Kluwer Academic Press, Dordrecht, The Netherlands, 189 pp.

Edlund, M. B. and Stoermer, E. F. 2000. A 200,000-year, high-resolution record of diatom productivity and community makeup from Lake Baikal shows high correspondence to the marine oxygen-isotope record of climate change. Limnology and Oceanography 45: 948-962.

Edlund, M. B., Engstrom, D. R., Triplett, L., Lafrancois, B. M. and Leavitt, P. R. 2009. Twentieth-century eutrophication of the St. Croix River (Minnesota-Wisconsin, USA) reconstructed from the sediments of its natural impoundment. Journal of Paleolimnology 41: 641-657. DOI:10.1007/s10933-008-9296-1

Ryves, D. B., Juggins, S., Fritz, S. C., & Battarbee, R. W. (2001). Experimental diatom dissolution and the quantification of microfossil preservation in sediments. Palaeogeography, Palaeoclimatology, Palaeoecology, 172(1), 99-113

Ryves, D. B., Battarbee, R. W., & Fritz, S. C. (2009). The dilemma of disappearing diatoms: Incorporating diatom dissolution data into palaeoenvironmental modelling and reconstruction. Quaternary Science Reviews, 28(1), 120-136

Sicko-Goad, L., Stoermer, E. F., & Kociolek, J. P. (1989). Diatom resting cell rejuvenation and formation: time course, species records and distribution. Journal of Plankton Research, 11(2), 375-389

Smol, J. P. (1985). The ratio of diatom frustules to chrysophycean statospores: a useful paleolimnological index. Hydrobiologia, 123(3), 199-208.

Smol, J. P. (1988). Chrysophycean microfossils in paleolimnological studies. Palaeogeography, Palaeoclimatology, Palaeoecology, 62(1), 287-297.

Smol, J. P. (1995). Application of chrysophytes to problems in paleoecology In: Chrysophyte algae: ecology, phylogeny and development, p. 303.

Wilkinson, A.N., Zeeb, B.A., and Smol, J.P. 2001. Atlas of chrysophycean cysts. Volume II. Dev. Hydrobiol. 157: 1–169.

Wolin, J. A., & Stone, J. R. (2010). Diatoms as indicators of water-level change in freshwater lakes. The diatoms: applications for the environmental and earth sciences, 174-185.

Zeeb, B. A., & Smol, J. P. (2001). Chrysophyte scales and cysts. In Tracking environmental change using lake sediments (pp. 203-223). Springer Netherlands.

Zeeb, B. A., & Smol, J. P. (2001). Chrysophyte scales and cysts. In Tracking environmental change using lake sediments (pp. 203-223). Springer Netherlands.



Saturday, October 18, 2014

Opaque materials

 Opaque materials in smear slides

by Mark Shapley

Most smear slides will include somewhere between a trace and a few percent grains that do not transmit light and appear opaque in plane-polarized light under a petrographic microscope.  These are likely to be of diverse origins, perhaps including charcoal, chitin and other organic debris, sulfide and oxide minerals, and pseudo-opaque slide manufacturing artifacts such as bubbles.

Opaque or near-opaque organic debris are often distinguishable by form and surface texture. Charcoal may show wood grain or the form of plant structure, depending on degree of preservation, and will appear black even along thin edges of a fragment. Unburned but heavily oxidized organic debris may show similar form, but will usually appear brownish along fragment margins.  When probed under a dissecting microscope, charcoal often shows a brittle fracture, while oxidized plant debris tends to crumble; material preserved in a smear slide may also exhibit signs of this mechanical contrast between combusted and un-combusted material.  Chitin is a common constituent of lacustrine sediment and in some cases may be opaque or nearly so.  Chitin preserving the body structures of insects (limbs, eyes, wings, mouth-parts, etc.) is readily recognized.  Dark chitin also occurs as irregular fragments that may superficially resemble charcoal, plant debris, or even strongly colored or opaque minerals.  Most will show some translucence, and (under a dissecting microscope) will flex readily.
Charcoal
Chitin
Thick plant material

Pyrite





Magnetite


Bubbles


Truly opaque mineral grains that commonly occur as components of unlithified sediment are most often iron-bearing sulfide or oxide minerals.  Distinguishing between them can be challenging, but is valuable for the information to be gained about depositional and diagenetic environments.  The presence of sulfide minerals (including pyrite, FeS2, and monosulfides such as greigite, mackinawite, and pyrrhotite) signifies active microbial sulfur cycling, the presence of free sulfide, and a sulfide-based sediment sink for iron and other available cations.  In contrast, opaque Fe-oxides, usually magnetite, are typically of detrital origin, often associated with other grains of high specific gravity in clastic lag deposits comprising heavy mineral assemblages or the high-density mineral-rich fraction of some tephras. Other commonly-seen Fe-oxides (e.g. hematite) may appear strongly colored in smear slides but are typically translucent to transparent, rather than opaque. Biomineralized magnetite produced by magnetotactic bacteria is widespread in sedimentary environments, but in a sub-micron size class well below grains discriminated in smear slides.


While Fe-sulfides and magnetite share opacity as a defining characteristic, their optical characteristics under reflected light can help distinguish the two mineral classes. (Oblique reflected light can be cast on the stage of a petrographic microscope using an overhead fiber-optic light source.) Pyrite and greigite typically show a bronze hue in reflected light, while magnetite usually appears gray. Pyrite and greigite very commonly occur in spheroidal aggregates (framboids) of fine sub-micron crystallites, sometimes as infilling or replacement of diatom frustules and other biological structures and very often associated with organic matter, the oxidation of which promotes low redox potential and liberation of sulfide ions necessary for sulfide formation.  Pyrite also occurs as monocrystalline, euhedral cubes and octahedra, however, and the latter may be confused with magnetite, which commonly shows this form. Here the assemblage of sediment constituents may provide some of the more useful clues to opaque mineral identification, at least as far as mineral class. An opaque mineral grain with cubic or faceted grain boundaries that is associated with mafic silicate minerals, zircons, or other high-density grains is more likely magnetite or a related oxide, while a similar grain morphometry associated with abundant biogenic material is more likely pyrite. Fine opaque grains associated with biotite or other Fe-bearing silicate minerals may indicate alteration of detrital components and Fe mobility within the sedimentary environment.


Macroscopic (core face) observations also give valuable clues to the mineralogy of opaque phases seen in smear slides. Disseminated Fe-monosulfides (greigite, mackinawite) oxidize within minutes or hours of initial exposure to air, dramatically changing the appearance of freshly split cores if the sulfides are abundant.  Such reactive macro-scale behavior from exposed sediments alerts a smear-slide analyst to the presence of Fe-sulfides and/or their oxidized reaction products.        
Freshly split core with weakly crystalline sulfides (black) oxidizing at the edges of the core

Solute evolution and smear slide analysis

Mineral precipitation, solute evolution and balance, and smear slide analysis

by Mark Shapley

Among the many valuable paleolimnological inferences that can be supported by smear slide analysis are those related to the evolving chemical composition of lake waters as chemical sediments are precipitated and as water is evaporated.  Developing a conceptual model of lake solute evolution provides the framework for interpreting changes in the occurrence or dominance of chemically precipitated sediments both spatially (i.e., between cores) and temporally (down-core).

The underlying ideas are relatively simple. 1) Weathering of the landscape delivers a mix of solutes to streams, groundwater, and ultimately to lakes; the composition of this flux of solutes is determined by the relative proportions and solubility of minerals present in or delivered to the watershed.  2) Minerals may achieve a state of saturation in lakes determined by thermodynamics (solubility product) and then can precipitate, removing ions from solution in proportion to the mineral composition. 3) The ratio of solutes (ions) removed to form a precipitating mineral phase need not, and generally will not, be the same as the ratio of those ions in the parent water.  4) The ratio of ions remaining in solution must therefore shift as mineral precipitation proceeds; the compositional trajectory followed by a water mass evolving in this way, and the series of minerals that will be precipitated with increasing evaporative concentration, can be predicted on the basis of its initial ionic composition.    

These concepts were developed in influential review publications by Hardie and Eugster (1970), Eugster and Jones (1979), and others, following upon case studies from regions of evaporative lake occurrence around the world.  Figure 1 encapsulates common alternate pathways of lake solute and mineral-saturation evolution, based on the initial solute ratios of evaporating lake-source waters. Very commonly, initial precipitation of low-Mg calcite (CaCO3) shifts lake water in the direction of more Mg-rich compositions and equilibrium with magnesium-bearing carbonate minerals.  Thermodynamic saturation may not result in the formation of magnesian carbonate minerals, however, as kinetic barriers come into play.  A commonly observed progression in mineral formation occurs from low-Mg calcite, to high-Mg calcite (with several mol % Mg substituted for Ca), to the CaCO3 polymorph aragonite (as higher Mg2+ activity above an uncertainly-defined Mg:Ca threshold interferes with calcite formation).  The path of lake water evolution in the direction of high or low HCO3-/(Ca2+ + Mg2+) compositions as carbonate mineral precipitation proceeds is largely determined by ion ratios of the initial lake inflow.    



Eugster and Jones 1979_no cap.png

Figure 1, modified from Eugster and Hardie, 1978, illustrating the concept of lake solute pathways, determined by initial inflow solute composition and sequential saturation of concentrating water with different mineral species as ion ratios evolve in response to mineral formation and sedimentation.  

Although shown here as progressive evolution from low to very high salinity, hydrologically balanced inflows and outflows may result in a quasi-equilibrium composition precipitating large masses of a relatively simple suite of minerals (Sanford and Wood, 1991).

Eugster and co-authors mostly considered closed-system (terminal lake) evolution of compositions from some starting point, evolving along compositionally determined pathways toward eventual saturation with the most soluble of salts.  Subsequent extension of these concepts considered the more typical situation of ‘leaky’ lakes where outflow of water evolved by in-lake processes to an intermediate composition occurs, allowing the development of steady-state solute compositions. Under conditions like these, persistent precipitation of a given mineral or restricted suite of minerals can occur over a long period of time, resulting in thick accumulations of a chemical sediment type (Sanford and Wood 1991; Teller and Last 1990; Donovan 1994).

The implications for smear slide analysis of lacustrine sediments may point in a couple of useful directions.  Associations of chemically precipitated minerals may be expected to fit a model of chemical evolution along a particular pathway; if they don’t, for example if associations of minerals representing different chemical divide pathways co-occur, then one might infer changing lake source compositions (or perhaps double-check one’s mineral identifications.)  If a given chemically precipitated mineral or mineral assemblage persists over long periods of time, especially in smaller lakes, one implication is that you are looking at a hydrologically ‘open’ system where outflow permits a state of dynamic equilibrium in solute balance.  The integration of down-core smear slide observations with spatial distribution of cores may indicate patterns of solute evolution and mineral endogenesis over time that can guide more quantitative evaluations of lake inflow composition and changing lake water balance.

Evidence of more concentrated lake solute status and composition will not necessarily be expressed in sediment mineral composition observable in smear slides, since brine outflow and/or later dissolution of highly soluble minerals may remove materials accumulated during saline lake phases.  Viewed in company with other compositional indicators (such as clastic components) as well as stratigraphic clues, however, valuable inferences of  lake solute evolution may be drawn  from smear slide observations.  For example, oscillatory transitions in lacustrine CaCO3 mineralogy between calcite and aragonite may imply changing Mg:Ca balance driven by CaCO3 precipitation and solute balance.  Carbonate or chloride salts, if preserved, each imply extensive evolution of continental waters but by progress along different branches of the ‘geochemical divide’, governed by initial solute balance and by preceding mineral precipitation.  Sulfate salts may occur in different associations, influenced both by initial inflow composition and by redox chemistry and its influence on microbial sulfate reduction.  Complex mineral sequences recording repetitive solute evolution can imply cyclic lake volume and solute changes, while persistent accumulations of restricted mineral suites suggest more narrowly varying lake volume and salinity, balanced by solute-bearing outflow.     


Donovan, J. J., 1994, On the measurement of reactive mass fluxes in evaporative groundwater-source lakes, Sedimentology and Geochemistry of Modern and Ancient Saline Lakes: SEPM Special Publication No. 50, SEPM (Society for Sedimentary Geology), p. 33 - 50.
Eugster, H. P., and Jones, B. F., 1979, Behavior of major solutes during closed-basin brine evolution: American Journal of Science, v. 279, p. 609-631.
Hardie, L. A., and Eugster, H. P., 1970, The evolution of closed-basin brines: Mineralogical Society of America Special Publications, v. 3, p. 273-290.
Sanford, W. E., and Wood, W., 1991, Brine evolution and mineral deposition in hydrologically open evaporite basins: American Journal of Science, v. 291, p. 687 - 710.
Teller, J. T., and Last, W. M., 1990, Paleohydrologic indicators in playas and salt lakes, with examples from Canada, Australia and Africa: Paleogeography, Paleoclimatology, Paleoecology, v. 76, p. 215-240.