Search USGSSearch

SEARCH · Search USGS

Results for “Elements”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28Linked to original sources

Elemental composition and potential toxicity of the riverine macrophyte Podostemum ceratophyllum Michx. reflects land use in eastern North America

Land use influences surface water quality, often alleviating stoichiometric constraints on primary production and altering biogeochemical cycling. However, land use effects on nutrient content and potential trace metal accumulation in aquatic plants remain unclear, and high concentrations of metals and altered nutrient ratios could impact the health of herbivores and detritivores. We tested for land use effects on nutrient and trace metal accumulation in a widespread riverine macrophyte, Podostemum ceratophyllum , collected from 91 locations from Georgia to Maine, USA in 2014–2016. We quantified carbon (C), nitrogen (N), phosphorus (P), their molar and mass ratios, N and C stable isotopes, and 17 additional elements in dried plants collected from each location to estimate relationships between plant tissue content and watershed land use, which we quantified as agriculture, forest, and development. Decreasing forest cover was correlated with increasing δ 15 N, Mg, Mn, and P in Podostemum tissue. Increasing urban development was correlated with increasing δ 15 N, Mg and P, while increasing agriculture was correlated with a decrease in C: P and the concentrations of multiple metals, along with increases in P, Mg and δ 15 N. Decreases in ratios of N: P and C:P with increasing agriculture and urban development in the watershed indicate more rapid P storage relative to C and N in plant tissue, and increased resource quality of the plant to consumers in these watersheds. We also observed potentially toxic dietary concentrations of some trace metals (B, Cd, Tl, Zn) in plant tissue which could be related to the plant's natural herbivory defense system or to land use. We conclude that land use influences the elemental composition of P. ceratophyllum , and potentially the quality and toxicity of the plant to herbivores and detritivores in eastern North American rivers.

Science of the Total Environment

Distribution of potentially hazardous trace elements in coals from Shanxi province, China

Shanxi province, located in the center of China, is the biggest coal base of China. There are five coal-forming periods in Shanxi province: Late Carboniferous (Taiyuan Formation), Early Permian (Shanxi Formation), Middle Jurassic (Datong Formation), Tertiary (Taxigou Formation), and Quaternary. Hundred and ten coal samples and a peat sample from Shanxi province were collected and the contents of 20 potentially hazardous trace elements (PHTEs) (As, B, Ba, Cd, Cl, Co, Cr, Cu, F, Hg, Mn, Mo, Ni, Pb, Sb, Se, Th, U, V and Zn) in these samples were determined by instrumental neutron activation analysis, atomic absorption spectrometry, cold-vapor atomic absorption spectrometry, ion chromatography spectrometry, and wet chemical analysis. The result shows that the brown coals are enriched in As, Ba, Cd, Cr, Cu, F and Zn compared with the bituminous coals and anthracite, whereas the bituminous coals are enriched in B, Cl, Hg, and the anthracite is enriched in Cl, Hg, U and V. A comparison with world averages and crustal abundances (Clarke values) shows that the Quaternary peat is highly enriched in As and Mo, Tertiary brown coals are highly enriched in Cd, Middle Jurassic coals, Early Permian coals and Late Carboniferous coals are enriched in Hg. According to the coal ranks, the bituminous coals are highly enriched in Hg, whereas Cd, F and Th show low enrichments, and the anthracite is also highly enriched in Hg and low enrichment in Th. The concentrations of Cd, F, Hg and Th in Shanxi coals are more than world arithmetic means of concentrations for the corresponding elements. Comparing with the United States coals, Shanxi coals show higher concentrations of Cd, Hg, Pb, Se and Th. Most of Shanxi coals contain lower concentrations of PHTEs. ?? 2004 Elsevier Ltd. All rights reserved.

Fuel

Coal quality trends and distribution of potentially hazardous trace elements in eastern Kentucky coals

Coal in the Eastern Kentucky coalfield has been, and continues to be, a valuable energy resource, especially for the electric utility industry. However, Federal mandates in Titles III and IV of the Clean Air Act Amendments of 1990 have placed increasingly stringent demands on the type and grade of coal that can be burnt in an environmentally acceptable manner. Therefore, a greater understanding of the spatial and temporal distribution of thickness and quality parameters, and the geologic factors that control their distribution, is critical if the Eastern Kentucky coalfield is to continue to be a major producer of high-quality coal. Information from the Kentucky Geological Survey's Coal Resource Information System database is used in this paper to document the geographic and stratigraphic distribution of important factors such as bed thickness, calorific value, ash yield and total sulfur content. The distribution of 15 elements that naturally occur in trace amounts in Kentucky coal is also discussed, as these elements may require monitoring with passage of Title III of the Clean Air Act Amendments of 1990.

Kentucky

Environmentally critical elements in channel and cleaned samples of Illinois coals

Sixteen trace and minor elements that occur in coal are listed among 189 substances identified as ‘hazardous air pollutants’ (HAPs) in the US Clean Air Act Amendments of 1990. We investigated the occurrence and cleanability of the 16 HAPs in Illinois coals, as a contribution to the discussion about the potential effect of pending environmental regulations on the future use of these coals in power generation. The average ash content of the samples of conventionally cleaned as-shipped coals is about 20% lower than that of standard channel samples. Conventional cleaning reduces the average concentrations of As, Cd, Co, Hg, Mn, Ni, Pb, Sb and Th in the as-shipped coals by more than 20% relative to channel samples. Thus, basing assessments of health risks from emissions of HAPs during coal combustion on channel samples without appropriate adjustment would overestimate the risk. Additional cleaning by froth-flotation reduces the ash content of finely-ground as-shipped coals by as much as 76% at an 80% combustibles recovery. Although the average froth-flotation cleanability for the majority of HAPs is less than that for ash, the cleanabilities in some individual cases approaches, or even exceeds, the cleanability for ash, depending on the modes of occurrences of the elements.

Illinois

Standard reference water samples for rare earth element determinations

Standard reference water samples (SRWS) were collected from two mine sites, one near Ophir, CO, USA and the other near Redding, CA, USA. The samples were filtered, preserved, and analyzed for rare earth element (REE) concentrations (La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu) by inductively coupled plasma-mass spectrometry (ICP-MS). These two samples were acid mine waters with elevated concentrations of REEs (0.45–161 μg/l). Seventeen international laboratories participated in a ‘round-robin’ chemical analysis program, which made it possible to evaluate the data by robust statistical procedures that are insensitive to outliers. The resulting most probable values are reported. Ten to 15 of the participants also reported values for Ba, Y, and Sc. Field parameters, major ion, and other trace element concentrations, not subject to statistical evaluation, are provided.

Applied Geochemistry

Geochemical and mineralogical controls on trace element release from the Penn Mine base-metal slag dump, California

Base-metal slag deposits at the Penn Mine in Calaveras County, California, are a source of environmental contamination through leaching of potentially toxic elements. Historical Cu smelting at Penn Mine (1865-1919) generated approximately 200,000 m3 of slag. The slag deposits, which are flooded annually by a reservoir used for drinking water and irrigation, also may be in contact with acidic ground waters (pH < 4) from the adjacent mine area. Slags vary from grey to black, are glassy to crystalline, and range in size from coarse sand to large (0.6 ?? 0.7 ?? 1.5 m), tub-shaped casts. Metals are hosted by a variety of minerals and two glass phases. On the basis of mineralogy, slags are characterized by 4 main types: fayalite-rich, glassy, willemite-rich, and sulfide-rich. The ranges in metal and metalloid concentrations of 17 slag samples are: As, 0.0004-0.92; Ba, 0.13-2.9; Cd, 0.0014-1.4; Cu, 0.18-6.4; Pb, 0.02-11; and Zn, 3.2-28 wt.%. Leachates from Toxicity Characteristic Leaching Procedure tests (acetic acid buffered at pH 4.93) on two wiltemite-rich slags contained Cd and Pb concentrations (up to 2.5 and 30 mg/l, respectively) in excess of US Environmental Protection Agency (USEPA) regulatory limits. Analyses of filtered (0.45 ??m) water, collected within the flooded slag dump during reservoir drawdown, reveal concentrations of Cd (1.7 ??g/l), Cu (35 ??g/l), and Zn (250 ??g/l) that exceed USEPA chronic toxicity guidelines for the protection of aquatic life. Data from field and laboratory studies were used to develop geochemical models with the program EQ3/6 that simulate irreversible mass-transfer between slag deposits and reservoir waters. These models include kinetic rate laws for abiotic sulfide oxidation and surface-controlled dissolution of silicates, oxides, and glass. Calculations demonstrate that the main processes controlling dissolved metal concentrations are (1) dissolution of fayalite, willemite, and glass; (2) sulfide oxidation; and (3) secondary phase precipitation. Close agreement between model results and measured concentrations of Al, Ba, Cu, Fe, SiO2, and SO4 in the slag dump pore waters suggests that the dissolved concentrations of these elements are controlled by solubility equilibrium with secondary phases. Differences between predicted and measured Cd and Pb concentrations imply that field weathering rates of glass and sulfides are approximately two orders of magnitude lower than laboratory rates. Overprediction of Pb release may also reflect other attenuation processes in the natural system, such as sorption or coprecipitation. ?? 2001 Elsevier Science Ltd. All rights reserved.

Applied Geochemistry

Inductively coupled plasma-mass spectrometry as an element-specific detector for field-flow fractionation particle separation

An inductively coupled plasma-mass spectrometer was used for the quantitative measurement of trace elements In specific,submicrometer size-fraction particulates, separated by sedimentation field-flow fractionation. Fractions were collected from the eluent of the field-flow fractionation centrifuge and nebulized, with a Babington-type pneumatic nebulizer, into an argon inductively coupled plasma-mass spectrometer. Measured Ion currents were used to quantify the major, minor, and trace element composition of the size-separated colloidal (< 1-microm diameter) particulates. The composition of surface-water suspended matter collected from the Yarra and Darling rivers in Australia is presented to illustrate the usefulness of this tool for characterizing environmental materials. An adsorption experiment was performed using cadmium lon to demonstrate the utility for studying the processes of trace metal-suspended sediment interactions and contaminant transport in natural aquatic systems.

Analytical Chemistry

Modeling precipitation and sorption of elements during mixing of river water and porewater in the Coeur d'Alene River basin

Reddish brown flocs form along the edge of the Coeur d'Alene River when porewater drains into river water during the annual lowering of water level in the basin. The precipitates are efficient scavengers of dissolved elements and have characteristics that may make metals associated with them bioavailable. This work characterizes the geochemistry of the porewater and models the formation and composition of the flocs. Porewater is slightly acidic, has suboxic to anoxic characteristics, tends to have higher alkalinity, and contains elevated concentrations of many constituents relative to river water. Laboratory mixing experiments involving porewater and river water were done to produce the precipitates. Thermodynamic predictions using PHREEQC indicate that predicted amounts of ferrihydrite and gibbsite agree with removal of Fe and Al. Predictions of element removal by adsorption onto ferrihydrite are consistent with observed removal using a combination of surface complexation constants for the generalized two-layer model (As and Se), alternative surface constants derived from experiments at high sorbate-to-sorbent ratios (Cd, Co, Cu, Ni, Pb, and Zn), and adjusted surface constants to fit experimental data (Cr, Mo, and Sb). This new set of surface complexation constants needs further testing in other contaminated systems.

Idaho

Isotopic and trace element variations in the Ruby Batholith, Alaska, and the nature of the deep crust beneath the Ruby and Angayucham Terranes

Thirty-six samples from plutons of the Ruby batholith of central Alaska were collected and analyzed for 22 trace elements, and many were analyzed for the isotopic compositions of Sr, Nd, O, and Pb in order to delimit the processes that produced the diversity of granodioritic to granitic compositions, to deduce the nature of the source of magmas at about 110 Ma, and to characterize the deep crust beneath the Ruby and Angayucham terranes. Plutons of the batholith show a substantial range in initial 87 Sr/ 86 Sr (SIR) of 0.7055&ndash;0.7235 and a general decrease from southwest to northeast. Initial 143 Nd/ 144 Nd (NIR) have a range of 0.51150&ndash;0.51232 and generally increase from southwest to northeast. The &delta; 18 O values for most whole rocks have a range of +8.4 to +11.8 and an average of +10.3&permil;. Rb, Cs, U, and Th show large ranges of concentration, generally increase as SiO 2 increases, and are higher in southwest than in northeast plutons. Sr, Ba, Zr, Hf, Ta, Sc, Cr, Co, and Zr show large ranges of concentration and generally decrease as SiO 2 increases. Rare earth elements (REE) show fractionated patterns and negative Eu anomalies. REE concentrations and anomalies are larger in the southwest than in the northeast plutons. Uniformity of SIR and NIR in Sithylemenkat and Jim River plutons suggests a strong role for fractional crystallization or melting of uniform magma sources at depth. Isotopic variability in Melozitna, Ray Mountains, Hot Springs, and Kanuti plutons suggests complex magmatic processes such as magma mixing and assimilation, probably combined with fractional crystallization, or melting of a complex source at depth. The large variations in SIR and NIR in the batholith require a variation in source materials at depth. The southwestern plutons probably had dominantly siliceous sources composed of metamorphosed Proterozoic and Paleozoic upper crustal rocks. The northeastern plutons probably had Paleozoic sources that were mixtures of siliceous and intermediate to mafic crustal rocks. The inferred sources could well have been the higher-metamorphic-grade lithologic equivalents of the exposed Proterozoic(?) to Paleozoic schists, orthogneisses, and metavolcanic rocks of Ruby terrane, the silicic portions of which are quite radiogenic. The deeper crustal sources that gave rise to most of the batholithic magmas are inferred to be similar under both the Ruby metamorphic terrane and the Angayucham ophiolitic terrane.

Alaska

Remarkable isotopic and trace element trends in potassic through sodic Cretaceous plutons of the Yukon-Koyukuk Basin, Alaska, and the nature of the lithosphere beneath the Koyukuk terrane

During the period from 110 to 80 m.y. ago, a 450-km-long magmatic belt was active along the northern margin of Yukon-Koyukuk basin and on eastern Seward Peninsula. The plutons intruded Upper Jurassic(?) and Lower Cretaceous volcanic arc rocks and Cretaceous sedimentary rocks in Yukon-Koyukuk basin and Proterozoic and lower Paleozoic continental rocks in Seward Peninsula. Within Yukon-Koyukuk basin, the plutons vary in composition from calc-alkalic plutons on the east to potassic and ultrapotassic alkalic plutons on the west. Plutons within Yukon-Koyukuk basin were analyzed for trace element and isotopic compositions in order to discern their origin and the nature of the underling lithosphere. Farthest to the east, the calc-alkalic rocks of Indian Mountain pluton are largely tonalite and sodic granodiorite, and have low Rb (average 82 ppm), high Sr (>600 ppm), high chondrite-normalized (cn) Ce/Yb (16&ndash;37), low &delta; 18 O (+6.5 to +7.1), low initial 87 Sr/ 86 Sr (SIR) (0.704), and high initial 143 Nd/ 144 Nd (NIR) (0.5126). These rocks resemble those modelled elsewhere as partial melts and subsequent fractionates of basaltic or gabbroic metaigneous rocks, and may be products of melting in the deeper parts of the Late Jurassic(?) and Early Cretaceous volcanic arc. Farthest to the west, the two ultrapotassic bodies of Selawik and Inland Lake are high in Cs (up to 93 ppm), Rb (up to 997 ppm), Sr, Ba, Th, and light rare earth elements, have high (Ce/Yb)cn (30, 27), moderate to low &delta; 18 O (+8.4, +6.9), high SIR (0.712, 0.710), and moderate NIR (0.5121&ndash;0.5122). These rocks resemble rocks of Australia and elsewhere that were modelled as melts of continental mantle that had been previously enriched in large cations. This mantle may be Paleozoic or older. The farthest west alkalic pluton of Selawik Hills is largely monzonite, quartz monzonite, and granite; has moderate Rb (average 284 ppm), high Sr (>600 ppm), high (Ce/Yb)cn (15&ndash;25), moderate &delta; 18 O (+8.3 to +8.6), high SIR (0.708&ndash;0.712), and moderate NIR (0.5121&ndash;0.5122). These rocks may be the product of interaction of magma derived from old continental mantle and magma derived from old continental crust. Plutons between eastern and western extremes show completely gradational variations in the concentration of K and Rb and in the isotopic compositions of Sr, Nd, and O. These plutons probably originated either by melting in a mixed source composed of a Paleozoic or older continental section (mantle + crust) overlain by Mesozoic mafic arc rocks, or by mixing of ultrapotassic to potassic magmas from continental sources (mantle + crust), and tonalitic magmas from arc sources. We infer from these results that the northwest portion of Yukon-Koyukuk basin is underlain by a substantial continental basement of Paleozoic or greater age. This basement probably thins out to the east. There is no geochemical evidence for continental basement east of about longitude 157&deg;, or along a belt of at least 50 km width flanking Ruby Geanticline as far to the southwest as about longitude 161&deg;. These areas are probably underlain by oceanic and Mesozoic arc rocks.

Alaska

Stratigraphic, lithologic, and major element geochemical constraints on magmatic evolution at Lassen volcanic center, California

The Lassen volcanic center is the most recent of several long-lived volcanic centers in the southernmost Cascade Range. These centers have erupted products ranging from basaltic andesite to rhyolite and are superimposed on a background of regional basaltic to andesitic volcanism. The evolution of the Lassen volcanic center is described in three stages. Stages I and II comprise the Brokeoff volcano, and 80 km 3 andesitic stratocone, active from 600 to 400 ka. Brokeoff volcano is compositionally equivalent to the regional basaltic andesite to andesite volcanism in the Lassen region and is the result of structurally controlled focusing of the diffuse regional mafic magmatism. Stage III comprises a silicic dome field and adjacent area of hybrid andesites and has a total volume of about 100 km 3 . Volcanism during stage III was episodic and is subdivided into four sequences of lithologically and temporally distinct lavas. Stage III began at 400 ka with a rhyolitic, caldera-forming pyroclastic eruption and chemically related lavas. Additional sequences of dacite erupted between 250–200 ka and 100–0 ka. Hybrid andesites erupted adjacent to the silicic dome field between 300 and 0 ka. Porphyritic andesite and dacite with high Al 2 O 3 , low TiO 2 , medium K 2 O and FeO/MgO ratios of 1.5–2.0 are the most abundant rock types in the Lassen volcanic center. However, the single most voluminous unit is sparsely phyric rhyolite pumice. In general, the lavas of Lassen volcanic center form a single coherent trend on major element variation diagrams and in pseudo-quaternary phase space, consistant with an origin either by fractional crystallization or magma mixing. In detail, however, the lack of systematic temporal change in silica and subtly crossing trends indicate a complex origin. A variety of statistically successful fractional crystallization models can be constructed that derive Brokeoff andesites from regional magmas. An important conclusion of the modeling is that if fractional crystallization is the process responsible for generation of Brokeoff andesite, then the parent magma must be low to medium K in geochemical affinity in order to explain the variation in K 2 O. However, although major element variation can be modeled by fractional crystallization, petrographic and stratigraphic evidence indicates that magma mixing is an important but subtle process in Brokeoff lavas and suggests that lavas evolved in small independent batches. Lavas erupted during stage III, while predominantly silicic, range from 53 to 75% SiO 2 . Disequilibrium mineral assemblages in the stage III lavas indicate that they are not directly derived from Brokeoff andesite by fractional crystallization. Mixing of silicic magma with regional mafic magma and disaggregation of andesitic quenched magmatic inclusions play dominant roles in the compositional diversity of stage III lavas.

Journal of Geophysical Research Solid Earth

Comparison of digestion methods for total elemental analysis of peat and separation of its organic and inorganic components

In order to find the most efficient digestion method for the total elemental recovery in peat, ten samples were subjected to different techniques and analysed for Ca, Mg, Fe, AI, Na, K, Mn, P, Zn, Cu, Li, Cd, Co, Ni, Pb and Si using atomic-absorption spectrophotometry. The most satisfactory procedures were dry ashing followed by hydrofluoric acid treatment and wet digestion using a mixture of hot nitric, perchloric and hydrofluoric acids. The wet digestion offers the advantage of a single decomposition method for the determination of Ca, Mg, Fe, AI, K, Na, Mn, Cu, Li, Zn and P. An alkaline fusion technique was required for the determination of Si. Hydrogen peroxide was used to separate the peat into its organic and inorganic components, leading to the total recovery of the elements for both fractions.

The Analyst

Atomic weights of the elements 1999

The biennial review of atomic-weight, Ar(E), determinations and other cognate data have resulted in changes for the standard atomic weights of the following elements: from to nitrogen 14.006 74??0.000 07 14.0067??0.0002 sulfur 32.066??0.006 32.065??0.005 chlorine 35.4527??0.0009 35.453??0.002 germanium 72.61??0.02 72.64??0.01 xenon 131.29??0.02 131.293??0.006 erbium 167.26??0.03 167.259??0.003 uranium 238.0289??0.0001 238.028 91??0.000 03 Presented are updated tables of the standard atomic weights and their uncertainties estimated by combining experimental uncertainties and terrestrial variabilities. In addition, this report again contains an updated table of relative atomic mass values and half-lives of selected radioisotopes. Changes in the evaluated isotopic abundance values from those published in 1997 are so minor that an updated list will not be published for the year 1999. Many elements have a different isotopic composition in some nonterrestrial materials. Some recent data on parent nuclides that might affect isotopic abundances or atomic-weight values are included in this report for the information of the interested scientific community. ?? 2001 American Institute of Physics.

Journal of Physical and Chemical Reference Data

Characterization of fly ash from low-sulfur and high-sulfur coal sources: Partitioning of carbon and trace elements with particle size

Fly ash samples were collected in November and December of 1994, from generating units at a Kentucky power station using high- and low-sulfur feed coals. The samples are part of a two-year study of the coal and coal combustion byproducts from the power station. The ashes were wet screened at 100, 200, 325, and 500 mesh (150, 75, 42, and 25 mu m, respectively). The size fractions were then dried, weighed, split for petrographic and chemical analysis, and analyzed for ash yield and carbon content. The low-sulfur ''heavy side'' and ''light side'' ashes each have a similar size distribution in the November samples. In contrast, the December fly ashes showed the trend observed in later months, the light-side ash being finer (over 20 % more ash in the 500 mesh [25 mu m] fraction) than the heavy-side ash. Carbon tended to be concentrated in the coarse fractions in the December samples. The dominance of the 325 mesh (42 mum) fractions in the overall size analysis implies, though, that carbon in the fine sizes may be an important consideration in the utilization of the fly ash. Element partitioning follows several patterns. Volatile elements, such as Zn and As, are enriched in the finer sizes, particularly in fly ashes collected at cooler, light-side electrostatic precipitator (ESP) temperatures. The latter trend is a function of precipitation at the cooler-ESP temperatures and of increasing concentration with the increased surface area of the finest fraction. Mercury concentrations are higher in high-carbon fly ashes, suggesting Hg adsorption on the fly ash carbon. Ni and Cr are associated, in part, with the spinel minerals in the fly ash.

Energy Sources

The application of an analytic element model to investigate groundwater-lake interactions at Pretty Lake, Wisconsin

Pretty Lake is a 64 acre, sandy-bottomed groundwater flow-through lake that has a history of hydrologic disturbance. Residents and regulators require a better understanding of lake-groundwater interaction to develop measures to protect the lake's hydrologic system and water quality. A groundwater flow model was constructed as a tool to synthesize field data collected at the site, delineate recharge areas that supply groundwater to the lake, and predict die effect of dredging an adjacent drainage ditch. The one layer, two-dimensional steady-state areal model used analytic element (AE) methods because they are quick to apply and include sophisticated simulation of groundwater-surface water interaction. The model calibrated well to groundwater heads (mean absolute difference = 0.05 m), lake stage (within 0.05 m) and ditch fluxes (mean absolute difference = 0.0023 m 3 &middot;s &minus;1 ). Model results showed that a single 1000 m wide recharge area supplies all the groundwater inflow to the lake. In addition, the model predicted that dredging an adjacent ditch by 3.0 m would lower the lake level by 0.31 m. The analytic element model was verified using a widely accepted finite-difference (FD) code; differences were less than &plusmn;0.015 m near die lake area and reached a maximum of 0.08 m at far corners of the FD grid. These differences are likely a result of die nodal interpolation inherent to FD techniques and error associated with applying a discrete boundary to die AE infinite aquifer. Although developed recently, AE methods have great potential to aid characterizations of groundwater-lake systems.

Wisconsin

Structured populations of Sulfolobus acidocaldarius with susceptibility to mobile genetic elements

The impact of a structured environment on genome evolution can be determined through comparative population genomics of species that live in the same habitat. Recent work comparing three genome sequences of Sulfolobus acidocaldarius suggested that highly structured, extreme, hot spring environments do not limit dispersal of this thermoacidophile, in contrast to other co-occurring Sulfolobus species. Instead, a high level of conservation among these three S. acidocaldarius genomes was hypothesized to result from rapid, global-scale dispersal promoted by low susceptibility to viruses that sets S. acidocaldarius apart from its sister Sulfolobus species. To test this hypothesis, we conducted a comparative analysis of 47 genomes of S. acidocaldarius from spatial and temporal sampling of two hot springs in Yellowstone National Park. While we confirm the low diversity in the core genome, we observe differentiation among S. acidocaldarius populations, likely resulting from low migration among hot spring “islands” in Yellowstone National Park. Patterns of genomic variation indicate that differing geological contexts result in the elimination or preservation of diversity among differentiated populations. We observe multiple deletions associated with a large genomic island rich in glycosyltransferases, differential integrations of the Sulfolobus turreted icosahedral virus, as well as two different plasmid elements. These data demonstrate that neither rapid dispersal nor lack of mobile genetic elements result in low diversity in the S. acidocaldarius genomes. We suggest instead that significant differences in the recent evolutionary history, or the intrinsic evolutionary rates, of sister Sulfolobus species result in the relatively low diversity of the S. acidocaldarius genome.

Genome Biology and Evolution

Post-earthquake relaxation using a spectral element method: 2.5-D case

The computation of quasi-static deformation for axisymmetric viscoelastic structures on a gravitating spherical earth is addressed using the spectral element method (SEM). A 2-D spectral element domain is defined with respect to spherical coordinates of radius and angular distance from a pole of symmetry, and 3-D viscoelastic structure is assumed to be azimuthally symmetric with respect to this pole. A point dislocation source that is periodic in azimuth is implemented with a truncated sequence of azimuthal order numbers. Viscoelasticity is limited to linear rheologies and is implemented with the correspondence principle in the Laplace transform domain. This leads to a series of decoupled 2-D problems which are solved with the SEM. Inverse Laplace transform of the independent 2-D solutions leads to the time-domain solution of the 3-D equations of quasi-static equilibrium imposed on a 2-D structure. The numerical procedure is verified through comparison with analytic solutions for finite faults embedded in a laterally homogeneous viscoelastic structure. This methodology is applicable to situations where the predominant structure varies in one horizontal direction, such as a structural contrast across (or parallel to) a long strike-slip fault.

Geophysical Journal International

A grid-doubling finite-element technique for calculating dynamic three-dimensional spontaneous rupture on an earthquake fault

We present a new finite-element technique for calculating dynamic 3-D spontaneous rupture on an earthquake fault, which can reduce the required computational resources by a factor of six or more, without loss of accuracy. The grid-doubling technique employs small cells in a thin layer surrounding the fault. The remainder of the modelling volume is filled with larger cells, typically two or four times as large as the small cells. In the resulting non-conforming mesh, an interpolation method is used to join the thin layer of smaller cells to the volume of larger cells. Grid-doubling is effective because spontaneous rupture calculations typically require higher spatial resolution on and near the fault than elsewhere in the model volume. The technique can be applied to non-planar faults by morphing, or smoothly distorting, the entire mesh to produce the desired 3-D fault geometry. Using our FaultMod finite-element software, we have tested grid-doubling with both slip-weakening and rate-and-state friction laws, by running the SCEC/USGS 3-D dynamic rupture benchmark problems. We have also applied it to a model of the Hayward fault, Northern California, which uses realistic fault geometry and rock properties. FaultMod implements fault slip using common nodes, which represent motion common to both sides of the fault, and differential nodes, which represent motion of one side of the fault relative to the other side. We describe how to modify the traction-at-split-nodes method to work with common and differential nodes, using an implicit time stepping algorithm.

Geophysical Journal International